SHDOM: Surface File Format

There are currently four types of surface property files: Lambertian, Fresnel, RPV, and Ocean. The file format allows the temperature and reflection properties to vary arbitrarily across the surface. The surface properties are specified on a regular, evenly spaced grid. There should be Nx*Ny points in the surface file. The properties are bilinearly interpolated to the domain bottom grid point locations. The surface domain is also periodic in X and Y, so that the domain is delX*Nx by delY*Ny in size. Thus there are Nx cells across, and the routine that reads in the surface properties assigns the properties of the left side of the domain (IX=1) to the right side of the last cell (IX=NX+1). A uniform surface can be specified by having Nx=1 and Ny=1 and making delX and delY large enough to cover the whole domain.

The Lambertian surface file format is as follows:

    L
    Nx  Ny   delX  delY
    IX IY  SfcTemp SfcAlbedo  
      . . .
 

The Fresnel surface is specified with the complex index of refraction of the dielectric surface (Mre, Mim):

    F
    Nx  Ny   delX  delY
    IX IY  SfcTemp  Mre  Mim
      . . .
 

The RPV surface is specified with three parameters (rho0, k, Theta):

    R
    Nx  Ny   delX  delY
    IX IY  SfcTemp  rho0  k  Theta
      . . .
 

The Ocean surface is specified with the near surface wind speed (m/s) and the chlorophyl-alpha pigment concentration (mg/m^3):

    O
    Nx  Ny   delX  delY
    IX IY  SfcTemp  WindSpeed  PigmentConc
      . . .
 
The IX and IY are the surface grid numbers (1 to Nx, 1 to Ny). The surface temperature is specified in Kelvin.