--- MITgcm/verification/aim.5l_cs/input/data.aimphys	2002/12/10 04:19:37	1.3
+++ MITgcm/verification/aim.5l_cs/input/data.aimphys	2004/11/14 20:09:46	1.10
@@ -4,9 +4,11 @@
 #     aim_useFMsurfBC :: select surface B.C. from Franco Molteni
 #     aim_useMMsurfFc :: select Monthly Mean surface forcing (e.g., NCEP)
 #     aim_surfPotTemp :: surf.Temp input file is in Pot.Temp (aim_useMMsurfFc)
+#     aim_energPrecip :: account for energy of precipitation (snow & rain temp)
+#     aim_splitSIOsFx :: compute separately Sea-Ice & Ocean surf. Flux
+#     aim_LandFile :: file name for Land fraction [0-1]
 #     aim_MMsufx   :: sufix for all Monthly Mean surface forcing files
 #     aim_MMsufxLength :: Length of sufix (Monthly Mean surf. forc. files)
-#     aim_LandFile :: file name for Land fraction       (F.M. surfBC)
 #     aim_albFile  :: file name for Albedo input file   (F.M. surfBC)
 #     aim_vegFile  :: file name for vegetation fraction (F.M. surfBC)
 #     aim_sstFile  :: file name for  Sea.Surf.Temp      (F.M. surfBC)
@@ -14,23 +16,40 @@
 #     aim_oiceFile :: file name for Sea Ice fraction    (F.M. surfBC)
 #     aim_snowFile :: file name for Snow depth          (F.M. surfBC)
 #     aim_swcFile  :: file name for Soil Water content  (F.M. surfBC)
+#     aim_dragStrato :: stratospheric-drag damping time scale (s)
+#     aim_clrSkyDiag :: compute clear-sky radiation for diagnostics
+#     aim_taveFreq :: Frequency^-1 for time-average output (s)
 #     aim_diagFreq :: Frequency^-1 for diagnostic output (s)
 #     aim_tendFreq :: Frequency^-1 for tendencies output (s)
  &AIM_PARAMS
-   aim_useFMsurfBC=.FALSE.,
-   aim_useMMsurfFc=.TRUE.,
+#  aim_useFMsurfBC=.FALSE.,
+#  aim_useMMsurfFc=.TRUE.,
 #  aim_surfPotTemp=.TRUE.,
-   aim_MMsufx = '.cs.bin',
-   aim_MMsufxLength = 7,
+#  aim_energPrecip=.TRUE.,
+#  aim_splitSIOsFx=.TRUE.,
+   aim_LandFile= 'landFrc.2f2.bin',
+   aim_albFile = 'albedo.FM.bin',
+   aim_vegFile = 'vegetFrc.FM.bin',
+   aim_sstFile = 'seaSurfT.FM.bin',
+   aim_lstFile = 'lndSurfT.2f2.bin',
+   aim_oiceFile= 'seaIce.2f2.bin',
+   aim_snowFile= 'snowDepth.FM.bin',
+   aim_swcFile = 'soilMoist.FM.bin',
+#  aim_dragStrato= 2592000.,
+#  aim_clrSkyDiag=.TRUE.,
+#  aim_diagFreq=3600.,
  &
 
 #--   Physical constants
-#      P0   = reference pressure
-#      GG   = gravity accel.
-#      RD   = gas constant for dry air
-#      CP   = specific heat at constant pressure
-#      ALHC = latent heat of condensation
-#      SBC  = Stefan-Boltzmann constant            
+#       P0    = reference pressure                 [Pa=N/m2]
+#       GG    = gravity accel.                     [m/s2]
+#       RD    = gas constant for dry air           [J/kg/K]
+#       CP    = specific heat at constant pressure [J/kg/K]
+#       ALHC  = latent heat of condensation        [J/g]
+#       ALHF  = latent heat of freezing            [J/g]
+#       SBC   = Stefan-Boltzmann constant
+#     rainCP  = heat capacity of liquid water      [J/g/K]
+#     tFreeze = freezing temperature of pure water [K]
 
 #--   Constants for forcing fields (namelist AIM_PAR_FOR):
 #      SOLC   = Solar constant (area averaged) in W/m^2
@@ -40,10 +59,12 @@
 #      SDALB  = Snow depth (mm water) corresponding to maximum albedo
 #      SWCAP  = Soil wetness at field capacity (volume fraction)
 #      SWWIL  = Soil wetness at wilting point  (volume fraction)
+#      hSnowWetness :: snow depth (m) corresponding to maximum wetness
  &AIM_PAR_FOR
+   hSnowWetness = 0.,
  &
 
-#--   Constants for surface fluxes (namelist AIM_PAR_SFL): 
+#--   Constants for surface fluxes (namelist AIM_PAR_SFL):
 #      FWIND0 = ratio of near-sfc wind to lowest-level wind
 #      FTEMP0 = weight for near-sfc temperature extrapolation (0-1) :
 #               1 : linear extrapolation from two lowest levels
@@ -58,10 +79,13 @@
 #      VGUST  = wind speed for sub-grid-scale gusts
 #      CTDAY  = daily-cycle correction (dTskin/dSSRad)
 #      DTHETA = Potential temp. gradient for stability correction
+#      dTstab = potential temp. increment for stability function derivative
 #      FSTAB  = Amplitude of stability correction (fraction)
 #      HDRAG  = Height scale for orographic correction
 #      FHDRAG = Amplitude of orographic correction (fraction)
  &AIM_PAR_SFL
+#  FTEMP0 = 0.,
+#  FHUM0  = 0.,
  &
 
 #--   Convection constants (namelist AIM_PAR_CNV):
@@ -75,7 +99,7 @@
  &AIM_PAR_CNV
  &
 
-#--   Constants for large-scale condendation (namelist AIM_PAR_LSC): 
+#--   Constants for large-scale condensation (namelist AIM_PAR_LSC):
 #      TRLSC  = Relaxation time (in hours) for specific humidity
 #      RHLSC  = Maximum relative humidity threshold (at sigma=1)
 #      DRHLSC = Vertical range of relative humidity threshold
@@ -109,7 +133,7 @@
  &AIM_PAR_RAD
  &
 
-#--   Constants for vertical dif. and sh. conv. (namelist AIM_PAR_VDI): 
+#--   Constants for vertical dif. and sh. conv. (namelist AIM_PAR_VDI):
 #      TRSHC  = relaxation time (in hours) for shallow convection
 #      TRVDI  = relaxation time (in hours) for moisture diffusion
 #      TRVDS  = relaxation time (in hours) for super-adiab. conditions