pkg/aim_v23/aim_readparms.F

b7411f1a84 Jean*0001 #undef ALLOW_AIM_TAVE
d676f916b2 Jean*0002 #include "AIM_OPTIONS.h"
                0003
                0004 CBOP
                0005 C     !ROUTINE: AIM_READPARMS
                0006 C     !INTERFACE:
                0007       SUBROUTINE AIM_READPARMS( myThid )
                0008
                0009 C     !DESCRIPTION: \bv
                0010 C     *==========================================================*
                0011 C     | S/R AIM_READPARMS
                0012 C     | o Read AIM physics package parameters
                0013 C     *==========================================================*
d24aab1d8d Jean*0014 C     | Initialized parameter in common blocks:
d676f916b2 Jean*0015 C     |  FORCON, SFLCON, CNVCON, LSCCON, RADCON, VDICON
                0016 C     *==========================================================*
                0017 C     \ev
d24aab1d8d Jean*0018
d676f916b2 Jean*0019 C     !USES:
                0020       IMPLICIT NONE
                0021
                0022 C     == Global variables ===
                0023 #include "AIM_SIZE.h"
                0024 #include "EEPARAMS.h"
a27dd2281d Jean*0025 #include "PARAMS.h"
65007c221b Jean*0026 #include "AIM_PARAMS.h"
2a4f7f3ad0 Jean*0027 #include "AIM_CO2.h"
d0a9461855 Jean*0028 #ifdef ALLOW_MNC
                0029 #include "MNC_PARAMS.h"
                0030 #endif
d676f916b2 Jean*0031
                0032 C-   Physical constants + functions of sigma and latitude
                0033 c #include "com_physcon.h"
                0034 C-   Constants for sub-grid-scale physics
                0035 #include "com_forcon.h"
                0036 #include "com_sflcon.h"
                0037 #include "com_cnvcon.h"
                0038 #include "com_lsccon.h"
d24aab1d8d Jean*0039 #include "com_radcon.h"
                0040 #include "com_vdicon.h"
d676f916b2 Jean*0041
                0042 C     !INPUT/OUTPUT PARAMETERS:
d24aab1d8d Jean*0043 C     == Routine Arguments ==
                0044 C     myThid :: my Thread Id number
d676f916b2 Jean*0045       INTEGER myThid
                0046 CEOP
                0047
                0048 #ifdef ALLOW_AIM
                0049
d24aab1d8d Jean*0050 C     == Local Variables ==
2a4f7f3ad0 Jean*0051 C     msgBuf         :: Informational/error message buffer
                0052 C     iUnit          :: Work variable for IO unit number
b7411f1a84 Jean*0053 C     errCount       :: error counter
2a4f7f3ad0 Jean*0054 C     aim_fixed_pCO2 :: fixed and uniform pCO2 value for AIM radiation
                0055 C     atmpCO2init    :: Atmospheric pCO2 initial value of well-mixed atm. box
d676f916b2 Jean*0056       CHARACTER*(MAX_LEN_MBUF) msgBuf
2a4f7f3ad0 Jean*0057       CHARACTER*10 endList
                0058       INTEGER buffI(1)
d676f916b2 Jean*0059       INTEGER iUnit
2a4f7f3ad0 Jean*0060       INTEGER errCount
                0061       _RL aim_fixed_pCO2
                0062       _RL atmpCO2init
b7411f1a84 Jean*0063 C-    retired parameters:
                0064 C     aim_taveFreq :: retired; Frequency^-1 for time-average output (s)
                0065       _RL aim_taveFreq
                0066       LOGICAL aim_timeave_mnc
d676f916b2 Jean*0067
                0068 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
                0069
65007c221b Jean*0070 C--   AIM interface parameter (common AIM_PARM_?):
                0071 C     aim_useFMsurfBC :: select surface B.C. from Franco Molteni
                0072 C     aim_useMMsurfFc :: select Monthly Mean surface forcing (e.g., NCEP)
d0a9461855 Jean*0073 C     aim_surfForc_TimePeriod :: Length of forcing time period (e.g. 1 month)
                0074 C     aim_surfForc_NppCycle   :: Number of time period per Cycle (e.g. 12)
                0075 C     aim_surfForc_TransRatio :: transition ratio from one month to the next
65007c221b Jean*0076 C     aim_surfPotTemp :: surf.Temp input file is in Pot.Temp (aim_useMMsurfFc)
b3097ed02d Jean*0077 C     aim_energPrecip :: account for energy of precipitation (snow & rain temp)
                0078 C     aim_splitSIOsFx :: compute separately Sea-Ice & Ocean surf. Flux
                0079 C                 (also land SW & LW) ; default=F as in original version
65007c221b Jean*0080 C     aim_MMsufx   :: sufix for all Monthly Mean surface forcing files
                0081 C     aim_MMsufxLength :: Length of sufix (Monthly Mean surf. forc. files)
                0082 C     aim_LandFile :: file name for Land fraction       (F.M. surfBC)
                0083 C     aim_albFile  :: file name for Albedo input file   (F.M. surfBC)
                0084 C     aim_vegFile  :: file name for vegetation fraction (F.M. surfBC)
                0085 C     aim_sstFile  :: file name for  Sea.Surf.Temp      (F.M. surfBC)
                0086 C     aim_lstFile  :: file name for Land.Surf.Temp      (F.M. surfBC)
                0087 C     aim_oiceFile :: file name for Sea Ice fraction    (F.M. surfBC)
                0088 C     aim_snowFile :: file name for Snow depth          (F.M. surfBC)
                0089 C     aim_swcFile  :: file name for Soil Water content  (F.M. surfBC)
2a4f7f3ad0 Jean*0090 C     aim_dragStrato  :: stratospheric-drag damping time scale (s)
                0091 C     aim_select_pCO2 :: select AIM CO2 formulation (= 0,1,2 or 3)
                0092 C     aim_selectOceAlbedo :: Ocean albedo(latitude) scheme (integer)
                0093 C                     :: 0 = constant (default); 1 = Briegleb et al. 1986
e749d70ece Jean*0094 C     aim_clrSkyDiag :: compute clear-sky radiation for diagnostics
65007c221b Jean*0095 C     aim_diagFreq :: Frequency^-1 for diagnostic output (s)
                0096 C     aim_tendFreq :: Frequency^-1 for tendencies output (s)
                0097       NAMELIST /AIM_PARAMS/
d24aab1d8d Jean*0098      &    aim_useFMsurfBC, aim_useMMsurfFc,
d0a9461855 Jean*0099      &    aim_surfForc_TimePeriod, aim_surfForc_NppCycle,
                0100      &    aim_surfForc_TransRatio, aim_surfPotTemp,
                0101      &    aim_energPrecip, aim_splitSIOsFx,
65007c221b Jean*0102      &    aim_MMsufx, aim_MMsufxLength,
                0103      &    aim_LandFile, aim_albFile, aim_vegFile,
                0104      &    aim_sstFile, aim_lstFile, aim_oiceFile, aim_snowFile,
5e328a6c4a Davi*0105      &    aim_swcFile, aim_qfxFile,
2a4f7f3ad0 Jean*0106      &    aim_dragStrato, aim_selectOceAlbedo,
                0107      &    aim_select_pCO2, aim_abs_pCO2, aim_fixed_pCO2, atmpCO2init,
ad0ca66649 Ed H*0108      &    aim_clrSkyDiag, aim_taveFreq, aim_diagFreq, aim_tendFreq,
d0a9461855 Jean*0109      &    aim_timeave_mnc, aim_snapshot_mnc,
2a4f7f3ad0 Jean*0110      &    aim_pickup_write_mnc, aim_pickup_read_mnc
65007c221b Jean*0111
d676f916b2 Jean*0112 C--   Physical constants (common PHYCON) :
b3097ed02d Jean*0113 C       P0    = reference pressure                 [Pa=N/m2]
                0114 C       GG    = gravity accel.                     [m/s2]
                0115 C       RD    = gas constant for dry air           [J/kg/K]
                0116 C       CP    = specific heat at constant pressure [J/kg/K]
                0117 C       ALHC  = latent heat of condensation        [J/g]
                0118 C       ALHF  = latent heat of freezing            [J/g]
                0119 C       SBC   = Stefan-Boltzmann constant
                0120 C     rainCP  = heat capacity of liquid water      [J/g/K]
                0121 C     tFreeze = freezing temperature of pure water [K]
d676f916b2 Jean*0122
                0123 C--   Constants for forcing fields (common FORCON) :
                0124 C      SOLC   = Solar constant (area averaged) in W/m^2
7f98c35e47 Davi*0125 C      OBLIQ  = Obliquity (in degree) only used if #define ALLOW_INSOLATION
d676f916b2 Jean*0126 C      ALBSEA = Albedo over sea
                0127 C      ALBICE = Albedo over sea ice (for ice fraction = 1)
                0128 C      ALBSN  = Albedo over snow (for snow depth > SDALB)
                0129 C      SDALB  = Snow depth (mm water) corresponding to maximum albedo
                0130 C      SWCAP  = Soil wetness at field capacity (volume fraction)
                0131 C      SWWIL  = Soil wetness at wilting point  (volume fraction)
bea6259b65 Jean*0132 C      hSnowWetness :: snow depth (m) corresponding to maximum wetness
d676f916b2 Jean*0133       NAMELIST /AIM_PAR_FOR/
                0134      &    SOLC, ALBSEA, ALBICE, ALBSN,
7f98c35e47 Davi*0135      &    SDALB, SWCAP, SWWIL, hSnowWetness, OBLIQ
d676f916b2 Jean*0136
d24aab1d8d Jean*0137 C--   Constants for surface fluxes (common SFLCON) :
d676f916b2 Jean*0138 C      FWIND0 = ratio of near-sfc wind to lowest-level wind
                0139 C      FTEMP0 = weight for near-sfc temperature extrapolation (0-1) :
                0140 C               1 : linear extrapolation from two lowest levels
                0141 C               0 : constant potential temperature ( = lowest level)
                0142 C      FHUM0  = weight for near-sfc specific humidity extrapolation (0-1) :
                0143 C               1 : extrap. with constant relative hum. ( = lowest level)
                0144 C               0 : constant specific hum. ( = lowest level)
                0145 C      CDL    = drag coefficient for momentum over land
                0146 C      CDS    = drag coefficient for momentum over sea
                0147 C      CHL    = heat exchange coefficient over land
                0148 C      CHS    = heat exchange coefficient over sea
                0149 C      VGUST  = wind speed for sub-grid-scale gusts
                0150 C      CTDAY  = daily-cycle correction (dTskin/dSSRad)
                0151 C      DTHETA = Potential temp. gradient for stability correction
e749d70ece Jean*0152 C      dTstab = potential temp. increment for stability function derivative
d676f916b2 Jean*0153 C      FSTAB  = Amplitude of stability correction (fraction)
                0154 C      HDRAG  = Height scale for orographic correction
                0155 C      FHDRAG = Amplitude of orographic correction (fraction)
                0156       NAMELIST /AIM_PAR_SFL/
                0157      &    FWIND0, FTEMP0, FHUM0,
                0158      &    CDL, CDS, CHL, CHS, VGUST, CTDAY,
e749d70ece Jean*0159      &    DTHETA, dTstab, FSTAB, HDRAG, FHDRAG
d676f916b2 Jean*0160
                0161 C--   Convection constants (common CNVCON) :
                0162 C      PSMIN  = minimum (norm.) sfc. pressure for the occurrence of convection
                0163 C      TRCNV  = time of relaxation (in hours) towards reference state
                0164 C      QBL    = specific hum. threshold in the boundary layer
                0165 C      RHBL   = relative hum. threshold in the boundary layer
                0166 C      RHIL   = rel. hum. threshold in intermed. layers for secondary mass flux
                0167 C      ENTMAX = max. entrainment as a fraction of cloud-base mass flux
                0168 C      SMF    = ratio between secondary and primary mass flux at cloud-base
                0169       NAMELIST /AIM_PAR_CNV/
                0170      &    PSMIN, TRCNV, QBL, RHBL, RHIL, ENTMAX, SMF
                0171
d24aab1d8d Jean*0172 C--   Constants for large-scale condendation (common LSCCON) :
d676f916b2 Jean*0173 C      TRLSC  = Relaxation time (in hours) for specific humidity
                0174 C      RHLSC  = Maximum relative humidity threshold (at sigma=1)
                0175 C      DRHLSC = Vertical range of relative humidity threshold
                0176 C      QSMAX  = used to define the maximum latent heat release
                0177       NAMELIST /AIM_PAR_LSC/
                0178      &    TRLSC, RHLSC, DRHLSC, QSMAX
                0179
                0180 C--    Radiation constants (common RADCON) :
                0181 C      RHCL1  = relative hum. corresponding to cloud cover = 0
                0182 C      RHCL2  = relative hum. corresponding to cloud cover = 1
                0183 C      QACL1  = specific hum. threshold for cloud cover in the upper troposphere
                0184 C      QACL2  = specific hum. threshold for cloud cover in the upper troposphere
                0185 C      ALBCL  = cloud albedo (for cloud cover = 1)
                0186 C      EPSSW  = fraction of incoming solar radiation absorbed by ozone
                0187 C      EPSLW  = fraction of surface LW radiation emitted directly to space
                0188 C      EMISFC = longwave surface emissivity
                0189 C--:            shortwave absorptivities (for dp = 10^5 Pa) :
                0190 C      ABSDRY = abs. of dry air      (visible band)
                0191 C      ABSAER = abs. of aerosols     (visible band)
                0192 C      ABSWV1 = abs. of water vapour (visible band, for dq = 1 g/kg)
                0193 C      ABSWV2 = abs. of water vapour (near IR band, for dq = 1 g/kg)
                0194 C      ABSCL1 = abs. of clouds       (visible band, constant term)
                0195 C      ABSCL2 = abs. of clouds       (visible band, for dw = 1 g/kg)
                0196 C--:            longwave absorptivities (per dp = 10^5 Pa) :
                0197 C      ABLWIN = abs. of air in "window" band
                0198 C      ABLCO2 = abs. of air in CO2 band
                0199 C      ABLWV1 = abs. of water vapour in H2O band 1 (weak),   for dq = 1 g/kg
                0200 C      ABLWV2 = abs. of water vapour in H2O band 2 (strong), for dq = 1 g/kg
                0201 C      ABLCL1 = abs. of clouds       in "window" band,       constant term
                0202 C      ABLCL2 = abs. of clouds       in "window" band,       for dw = 1 g/kg
                0203       NAMELIST /AIM_PAR_RAD/
                0204      &    RHCL1,  RHCL2,  QACL1,  QACL2,  ALBCL,
                0205      &    EPSSW,  EPSLW,  EMISFC,
                0206      &    ABSDRY, ABSAER, ABSWV1, ABSWV2, ABSCL1, ABSCL2,
d24aab1d8d Jean*0207      &    ABLWIN, ABLCO2, ABLWV1, ABLWV2, ABLCL1, ABLCL2
d676f916b2 Jean*0208
d24aab1d8d Jean*0209 C--   Constants for vertical dif. and sh. conv. (common VDICON) :
d676f916b2 Jean*0210 C      TRSHC  = relaxation time (in hours) for shallow convection
                0211 C      TRVDI  = relaxation time (in hours) for moisture diffusion
                0212 C      TRVDS  = relaxation time (in hours) for super-adiab. conditions
                0213 C      RHGRAD = maximum gradient of relative humidity (d_RH/d_sigma)
                0214 C      SEGRAD = minimum gradient of dry static energy (d_DSE/d_phi)
                0215       NAMELIST /AIM_PAR_VDI/
                0216      &    TRSHC, TRVDI, TRVDS, RHGRAD, SEGRAD
                0217
                0218 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
                0219
                0220       _BEGIN_MASTER(myThid)
d24aab1d8d Jean*0221
d0a9461855 Jean*0222 C--  Set default value for AIM interface code (AIM_PARAMS.h):
                0223       aim_useFMsurfBC = .TRUE.
                0224       aim_useMMsurfFc = .FALSE.
                0225       aim_surfPotTemp = .FALSE.
                0226       aim_energPrecip = .FALSE.
                0227       aim_splitSIOsFx = .FALSE.
                0228       aim_clrSkyDiag  = .FALSE.
                0229 #ifdef ALLOW_MNC
                0230       aim_snapshot_mnc      = snapshot_mnc .AND. useMNC
                0231       aim_pickup_write_mnc  = pickup_write_mnc .AND. useMNC
                0232       aim_pickup_read_mnc   = pickup_read_mnc .AND. useMNC
                0233 #else
                0234       aim_snapshot_mnc      = .FALSE.
                0235       aim_pickup_write_mnc  = .FALSE.
                0236       aim_pickup_read_mnc   = .FALSE.
                0237 #endif
                0238       aim_MMsufx = '.bin'
                0239       aim_MMsufxLength = 4
                0240 C-    Length (s) of one time period:
                0241       aim_surfForc_TimePeriod = 30.*86400.
                0242 C-    Number of time period per Cycle:
                0243       aim_surfForc_NppCycle   = 12
                0244 C-    define how fast the (linear) transition is from one month to the next
                0245 C      = 1                 -> linear between 2 midle month
                0246 C      > TimePeriod/deltaT -> jump from one month to the next one
                0247       aim_surfForc_TransRatio = 1.
                0248       aim_LandFile = ' '
                0249       aim_albFile  = ' '
                0250       aim_vegFile  = ' '
                0251       aim_sstFile  = ' '
                0252       aim_lstFile  = ' '
                0253       aim_oiceFile = ' '
                0254       aim_snowFile = ' '
                0255       aim_swcFile  = ' '
5e328a6c4a Davi*0256       aim_qfxFile  = ' '
d0a9461855 Jean*0257       aim_dragStrato = 0.
2a4f7f3ad0 Jean*0258       aim_selectOceAlbedo = 0
                0259       aim_select_pCO2 = 0
                0260       aim_ref_pCO2   = 320. _d -6
455d18a2e4 Jean*0261 C- Note: aim_abs_pCO2 default value has been calibrated (David Ferreira)
                0262 C        for pCO2 values within [ 100.ppm, 1000.ppm ]
                0263       aim_abs_pCO2   = 1.2616 _d 0
2a4f7f3ad0 Jean*0264       aim_fixed_pCO2 = UNSET_RL
                0265       atmpCO2init    = UNSET_RL
                0266 C-
d0a9461855 Jean*0267       aim_diagFreq = dumpFreq
                0268       aim_tendFreq = 0.
7f98c35e47 Davi*0269       OBLIQ = 23.45 _d 0
d0a9461855 Jean*0270
b7411f1a84 Jean*0271 C-    retired params default:
                0272       aim_taveFreq    = UNSET_RL
                0273       aim_timeave_mnc = .FALSE.
                0274
d0a9461855 Jean*0275 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
                0276
2a4f7f3ad0 Jean*0277       WRITE(msgBuf,'(A)') ' AIM_READPARMS: opening data.aimphys'
                0278       CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
                0279      &                    SQUEEZE_RIGHT, myThid )
d676f916b2 Jean*0280
2a4f7f3ad0 Jean*0281       CALL OPEN_COPY_DATA_FILE( 'data.aimphys', 'AIM_READPARMS',
d676f916b2 Jean*0282      O                          iUnit, myThid )
                0283
                0284 C--   Read parameters from open data file:
                0285
65007c221b Jean*0286 C-    Parameters for AIM interface code:
                0287       READ(UNIT=iUnit,NML=AIM_PARAMS)
                0288
d676f916b2 Jean*0289 C-    Constants for boundary forcing
                0290       READ(UNIT=iUnit,NML=AIM_PAR_FOR)
                0291
                0292 C-    Constants for surface fluxes
                0293       READ(UNIT=iUnit,NML=AIM_PAR_SFL)
d24aab1d8d Jean*0294
d676f916b2 Jean*0295 C-    Constants for convection
                0296       READ(UNIT=iUnit,NML=AIM_PAR_CNV)
                0297
                0298 C-    Constants for large-scale condensation
                0299       READ(UNIT=iUnit,NML=AIM_PAR_LSC)
                0300
                0301 C-    Constants for radiation
                0302       READ(UNIT=iUnit,NML=AIM_PAR_RAD)
                0303
                0304 C-    Constants for vertical diffusion and sh. conv.
                0305       READ(UNIT=iUnit,NML=AIM_PAR_VDI)
                0306
d24aab1d8d Jean*0307       WRITE(msgBuf,'(A)')
2a4f7f3ad0 Jean*0308      &   ' AIM_READPARMS: finished reading data.aimphys'
                0309       CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
                0310      &                    SQUEEZE_RIGHT, myThid )
d24aab1d8d Jean*0311
d676f916b2 Jean*0312 C--   Close the open data file
7a77863887 Mart*0313 #ifdef SINGLE_DISK_IO
d676f916b2 Jean*0314       CLOSE(iUnit)
7a77863887 Mart*0315 #else
                0316       CLOSE(iUnit,STATUS='DELETE')
                0317 #endif /* SINGLE_DISK_IO */
d676f916b2 Jean*0318
                0319 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
3a474aa258 Jean*0320 C--   Check parameters and model configuration
2a4f7f3ad0 Jean*0321       errCount = 0
                0322
b7411f1a84 Jean*0323 C-    Retired params:
                0324       IF ( aim_taveFreq .NE. UNSET_RL ) THEN
                0325         WRITE(msgBuf,'(2A)') 'AIM_READPARMS: ',
                0326      &    '"aim_taveFreq" is no longer allowed in file "data.aimphys"'
                0327         CALL PRINT_ERROR( msgBuf, myThid )
                0328         WRITE(msgBuf,'(2A)') 'AIM_READPARMS: ',
                0329      &    ' since "pkg/timeave" has been removed.'
                0330         CALL PRINT_ERROR( msgBuf, myThid )
                0331         errCount = errCount + 1
                0332       ENDIF
                0333       IF ( aim_timeave_mnc ) THEN
                0334         WRITE(msgBuf,'(2A)') 'AIM_READPARMS: "aim_timeave_mnc"',
                0335      &    ' is no longer allowed in file "data.aimphys"'
                0336         CALL PRINT_ERROR( msgBuf, myThid )
                0337         WRITE(msgBuf,'(2A)') 'AIM_READPARMS: ',
                0338      &    ' since "pkg/timeave" has been removed.'
                0339         CALL PRINT_ERROR( msgBuf, myThid )
                0340         errCount = errCount + 1
                0341       ENDIF
                0342 #ifdef ALLOW_AIM_TAVE
                0343         WRITE(msgBuf,'(2A)') 'AIM_READPARMS: :'
                0344      &    '#define ALLOW_AIM_TAVE not allowed in AIM_OPTIONS.h'
                0345         CALL PRINT_ERROR( msgBuf, myThid)
                0346         WRITE(msgBuf,'(2A)') 'AIM_READPARMS: :'
                0347      &    ' since "pkg/timeave" has been removed.'
                0348         CALL PRINT_ERROR( msgBuf, myThid)
                0349         errCount = errCount + 1
                0350 #endif /* ALLOW_AIM_TAVE */
                0351
                0352 C-    Consistency check:
2a4f7f3ad0 Jean*0353       IF ( aim_select_pCO2.LT.0 .OR. aim_select_pCO2.GT.3 ) THEN
                0354         WRITE(msgBuf,'(2A)') 'AIM_READPARMS: ',
                0355      &       'Invalid value for "aim_select_pCO2"'
                0356         CALL PRINT_ERROR( msgBuf, myThid )
                0357         errCount = errCount + 1
                0358       ENDIF
                0359 #ifndef ALLOW_AIM_CO2
                0360       IF ( aim_select_pCO2.GE.2 .AND. aim_select_pCO2.LE.3 ) THEN
                0361         WRITE(msgBuf,'(2A)') 'AIM_READPARMS: ',
                0362      &       'aim_select_pCO2 > 1 but ALLOW_AIM_CO2 undef'
                0363         CALL PRINT_ERROR( msgBuf, myThid )
                0364         WRITE(msgBuf,'(2A)')
                0365      &       'Re-compile with: #define ALLOW_AIM_CO2 (AIM_OPTIONS.h)'
                0366         CALL PRINT_ERROR( msgBuf, myThid )
                0367         errCount = errCount + 1
                0368       ENDIF
                0369 #endif /* ndef ALLOW_AIM_CO2 */
                0370
                0371       IF ( aim_fixed_pCO2.EQ.UNSET_RL ) THEN
                0372         aim_fixed_pCO2 = aim_ref_pCO2
                0373       ELSEIF ( aim_select_pCO2.NE.1 .AND.
                0374      &         aim_select_pCO2.GE.0 .AND. aim_select_pCO2.LE.3 ) THEN
                0375         WRITE(msgBuf,'(2A)') 'AIM_READPARMS: ',
                0376      &       'aim_fixed_pCO2 only used if aim_select_pCO2=1'
                0377         CALL PRINT_ERROR( msgBuf, myThid )
                0378         errCount = errCount + 1
                0379       ENDIF
                0380       IF ( atmpCO2init.EQ.UNSET_RL ) THEN
                0381         atmpCO2init = aim_ref_pCO2
                0382       ELSEIF ( aim_select_pCO2.GE.0 .AND. aim_select_pCO2.LT.2 ) THEN
                0383         WRITE(msgBuf,'(2A)') '** WARNING ** AIM_READPARMS: ',
                0384      &       'atmpCO2init only used if aim_select_pCO2 > 1'
                0385         CALL PRINT_MESSAGE( msgBuf, errorMessageUnit,
                0386      &                      SQUEEZE_RIGHT, myThid )
                0387       ENDIF
3a474aa258 Jean*0388
e749d70ece Jean*0389 #ifndef ALLOW_CLR_SKY_DIAG
                0390 C-    If aim_clrSkyDiag is set, then must compile the Clear-Sky Diagnostic code
                0391       IF ( aim_clrSkyDiag ) THEN
2a4f7f3ad0 Jean*0392         WRITE(msgBuf,'(2A)') 'AIM_READPARMS: ',
                0393      &       'aim_clrSkyDiag=T but ALLOW_CLR_SKY_DIAG undef'
                0394         CALL PRINT_ERROR( msgBuf, myThid )
e749d70ece Jean*0395         WRITE(msgBuf,'(2A)')
                0396      &  'Re-compile with: #define ALLOW_CLR_SKY_DIAG (AIM_OPTIONS.h)'
2a4f7f3ad0 Jean*0397         CALL PRINT_ERROR( msgBuf, myThid )
                0398         errCount = errCount + 1
e749d70ece Jean*0399       ENDIF
                0400 #endif
                0401
ad0ca66649 Ed H*0402 #ifdef ALLOW_MNC
d24aab1d8d Jean*0403       aim_snapshot_mnc       = useMNC .AND. aim_snapshot_mnc
                0404       aim_pickup_write_mnc   = useMNC .AND. aim_pickup_write_mnc
                0405       aim_pickup_read_mnc    = useMNC .AND. aim_pickup_read_mnc
                0406 #else /* ALLOW_MNC */
78ad6fc820 Ed H*0407       aim_snapshot_mnc       = .FALSE.
                0408       aim_pickup_write_mnc   = .FALSE.
                0409       aim_pickup_read_mnc    = .FALSE.
ad0ca66649 Ed H*0410 #endif /* ALLOW_MNC */
d24aab1d8d Jean*0411       aim_snapshot_mdsio     = (.NOT. aim_snapshot_mnc)
                0412      &                     .OR. outputTypesInclusive
                0413       aim_pickup_write_mdsio = (.NOT. aim_pickup_write_mnc)
                0414      &                     .OR. outputTypesInclusive
                0415       aim_pickup_read_mdsio  = (.NOT. aim_pickup_read_mnc)
                0416      &                     .OR. outputTypesInclusive
ad0ca66649 Ed H*0417
3a474aa258 Jean*0418 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
d676f916b2 Jean*0419 C--   Print out parameter values :
2a4f7f3ad0 Jean*0420       iUnit = standardMessageUnit
                0421       endList = '    ;     '
d676f916b2 Jean*0422
                0423       WRITE(msgBuf,'(A)') ' '
2a4f7f3ad0 Jean*0424       CALL PRINT_MESSAGE( msgBuf, iUnit, SQUEEZE_RIGHT, myThid )
d676f916b2 Jean*0425       WRITE(msgBuf,'(A)') '// ==================================='
2a4f7f3ad0 Jean*0426       CALL PRINT_MESSAGE( msgBuf, iUnit, SQUEEZE_RIGHT, myThid )
d676f916b2 Jean*0427       WRITE(msgBuf,'(A)') '// AIM physics parameters :'
2a4f7f3ad0 Jean*0428       CALL PRINT_MESSAGE( msgBuf, iUnit, SQUEEZE_RIGHT, myThid )
d676f916b2 Jean*0429       WRITE(msgBuf,'(A)') '// ==================================='
2a4f7f3ad0 Jean*0430       CALL PRINT_MESSAGE( msgBuf, iUnit, SQUEEZE_RIGHT, myThid )
d676f916b2 Jean*0431
a27dd2281d Jean*0432        CALL WRITE_0D_L( useLand, INDEX_NONE,
                0433      &                 'useLand =',
d24aab1d8d Jean*0434      &'   /* use Land package for surf.BC over land  */')
a27dd2281d Jean*0435
65007c221b Jean*0436 C- namelist AIM_PARAMS:
                0437        CALL WRITE_0D_L( aim_useFMsurfBC, INDEX_NONE,
                0438      &                 'aim_useFMsurfBC =',
d24aab1d8d Jean*0439      &'   /* F.Molteni surf. BC scheme on/off flag */')
65007c221b Jean*0440        CALL WRITE_0D_L( aim_useMMsurfFc, INDEX_NONE,
                0441      &                 'aim_useMMsurfFc =',
d24aab1d8d Jean*0442      &'   /* Monthly Mean surf. Forcing on/off flag */')
4da4b49499 Jean*0443        CALL WRITE_0D_RL(aim_surfForc_TimePeriod, INDEX_NONE,
d0a9461855 Jean*0444      &                 'aim_surfForc_TimePeriod =',
                0445      &  ' /* Length of forcing time period (s) */')
                0446        CALL WRITE_0D_I( aim_surfForc_NppCycle, INDEX_NONE,
                0447      &                 'aim_surfForc_NppCycle   =',
                0448      &  ' /* Number of time period per Cycle */')
4da4b49499 Jean*0449        CALL WRITE_0D_RL(aim_surfForc_TransRatio, INDEX_NONE,
d0a9461855 Jean*0450      &                 'aim_surfForc_TransRatio =',
                0451      &  ' /* transition ratio from one month to the next */')
65007c221b Jean*0452        CALL WRITE_0D_L( aim_surfPotTemp, INDEX_NONE,
                0453      &                 'aim_surfPotTemp =',
d24aab1d8d Jean*0454      &'   /* Surf Temp file in Pot.Temp on/off flag */')
b3097ed02d Jean*0455        CALL WRITE_0D_L( aim_energPrecip, INDEX_NONE,
                0456      &                 'aim_energPrecip =',
d24aab1d8d Jean*0457      &  ' /* account for energy of precip. on/off flag */')
b3097ed02d Jean*0458        CALL WRITE_0D_L( aim_splitSIOsFx, INDEX_NONE,
                0459      &                 'aim_splitSIOsFx =',
d24aab1d8d Jean*0460      &  ' /* separate Sea-Ice & Ocean Flux on/off flag */')
4da4b49499 Jean*0461        CALL WRITE_0D_RL(aim_dragStrato,INDEX_NONE,'aim_dragStrato=',
299f32bec2 Jean*0462      &    ' /* stratospheric-drag damping time scale (s) */')
2a4f7f3ad0 Jean*0463        CALL WRITE_0D_I( aim_selectOceAlbedo, INDEX_NONE,
                0464      &  'aim_selectOceAlbedo =', ' /* Ocean albedo(latitude) flag */')
                0465 C-
                0466        WRITE(msgBuf,'(2A)') 'aim_select_pCO2=',
                0467      &              ' /* select AIM CO2 option */'
                0468        CALL PRINT_MESSAGE( msgBuf, iUnit, SQUEEZE_RIGHT, myThid )
                0469        buffI(1) = aim_select_pCO2
                0470        CALL PRINT_LIST_I( buffI, 1, 1, INDEX_NONE,
                0471      &                    .FALSE., .TRUE., iUnit )
                0472        WRITE(msgBuf,'(2A)')
                0473      &      '   =0: prescribed LW absorption in CO2 band (=ABLCO2)'
                0474        CALL PRINT_MESSAGE( msgBuf, iUnit, SQUEEZE_RIGHT, myThid )
                0475        WRITE(msgBuf,'(2A)')
                0476      &      '   =1: absorption fct of prescribed pCO2 (=aim_fixed_pCO2)'
                0477        CALL PRINT_MESSAGE( msgBuf, iUnit, SQUEEZE_RIGHT, myThid )
                0478        WRITE(msgBuf,'(2A)')
                0479      &      '   =2: same as 0 with passive,  well-mixed pCO2 atm box'
                0480        CALL PRINT_MESSAGE( msgBuf, iUnit, SQUEEZE_RIGHT, myThid )
                0481        WRITE(msgBuf,'(2A)')
                0482      &      '   =3: same as 1 with evolving, well-mixed pCO2 atm-box'
                0483        CALL PRINT_MESSAGE( msgBuf, iUnit, SQUEEZE_RIGHT, myThid )
                0484        CALL PRINT_MESSAGE(endList, iUnit, SQUEEZE_RIGHT, myThid )
                0485        IF ( aim_select_pCO2.EQ.1 .OR. aim_select_pCO2.EQ.3 )
                0486      & CALL WRITE_0D_RL( aim_abs_pCO2, INDEX_NONE,
                0487      &  'aim_abs_pCO2 =', ' /* pCO2 dep. coeff. of CO2 absortion */')
                0488        IF ( aim_select_pCO2.EQ.1 )
                0489      & CALL WRITE_0D_RL( aim_fixed_pCO2, INDEX_NONE,
                0490      &  'aim_fixed_pCO2 =', ' /* fixed pCO2 for AIM radiation */')
                0491 #ifdef ALLOW_AIM_CO2
                0492        CALL WRITE_0D_RL( atmpCO2init, INDEX_NONE,
                0493      &  'atmpCO2init =', ' /* Initial atm-box pCO2 */')
                0494 #endif
e749d70ece Jean*0495        CALL WRITE_0D_L( aim_clrSkyDiag, INDEX_NONE,
                0496      &                 'aim_clrSkyDiag =',
d24aab1d8d Jean*0497      &  ' /* do clear-sky radiation diagnostics */')
4da4b49499 Jean*0498        CALL WRITE_0D_RL( aim_diagFreq, INDEX_NONE,'aim_diagFreq =',
65007c221b Jean*0499      &  '   /* Frequency^-1 for diagnostic output (s) */')
4da4b49499 Jean*0500        CALL WRITE_0D_RL( aim_tendFreq, INDEX_NONE,'aim_tendFreq =',
65007c221b Jean*0501      &  '   /* Frequency^-1 for tendencies output (s) */')
                0502
d676f916b2 Jean*0503 C- namelist AIM_PAR_FOR:
4da4b49499 Jean*0504        CALL WRITE_0D_RL( SOLC,  INDEX_NONE,'AIM_FOR:   SOLC =',
d676f916b2 Jean*0505      &  '   /* Solar constant (area averaged) in W/m2 */')
7f98c35e47 Davi*0506        CALL WRITE_0D_RL( OBLIQ, INDEX_NONE,'AIM_FOR:  OBLIQ =',
                0507      &  '   /* Obliquity (degree) used with ALLOW_INSOLATION */')
4da4b49499 Jean*0508        CALL WRITE_0D_RL( ALBSEA,INDEX_NONE,'AIM_FOR: ALBSEA =',
d676f916b2 Jean*0509      &  '   /* Albedo over sea [0-1]      */')
4da4b49499 Jean*0510        CALL WRITE_0D_RL( ALBICE,INDEX_NONE,'AIM_FOR: ALBICE =',
d676f916b2 Jean*0511      &  '   /* Albedo over sea ice (for ice fraction =1) */')
4da4b49499 Jean*0512        CALL WRITE_0D_RL( ALBSN, INDEX_NONE,'AIM_FOR:  ALBSN =',
d676f916b2 Jean*0513      &  '   /* Albedo over snow (for snow depth > SDALB) */')
4da4b49499 Jean*0514        CALL WRITE_0D_RL( SDALB, INDEX_NONE,'AIM_FOR:  SDALB =',
d676f916b2 Jean*0515      &  '   /* Snow depth (mm H2O) corresp. maximum albedo */')
4da4b49499 Jean*0516        CALL WRITE_0D_RL( SWCAP, INDEX_NONE,'AIM_FOR:  SWCAP =',
d676f916b2 Jean*0517      &  '   /* Soil wetness at field capacity (Vol.fract.) */')
4da4b49499 Jean*0518        CALL WRITE_0D_RL( SWWIL, INDEX_NONE,'AIM_FOR:  SWWIL =',
d676f916b2 Jean*0519      &  '   /* Soil wetness at wilting point  (Vol.fract.) */')
4da4b49499 Jean*0520        CALL WRITE_0D_RL( hSnowWetness, INDEX_NONE,
bea6259b65 Jean*0521      &  'AIM_FOR: hSnowWetness=',
                0522      &  '   /* snow depth corresp. maximum wetness (m) */')
d676f916b2 Jean*0523
                0524 C- namelist AIM_PAR_SFL:
4da4b49499 Jean*0525        CALL WRITE_0D_RL( FWIND0,INDEX_NONE,'AIM_SFL: FWIND0 =',
d676f916b2 Jean*0526      &  '   /* ratio of near-sfc wind to lowest-level wind */')
4da4b49499 Jean*0527        CALL WRITE_0D_RL( FTEMP0,INDEX_NONE,'AIM_SFL: FTEMP0 =',
d676f916b2 Jean*0528      &  '  /* weight for near-sfc temp. extrapolation (0-1)*/')
4da4b49499 Jean*0529        CALL WRITE_0D_RL( FHUM0, INDEX_NONE,'AIM_SFL:  FHUM0 =',
d676f916b2 Jean*0530      &  '  /* weight for near-sfc spec.humid. extrap. (0-1)*/')
4da4b49499 Jean*0531        CALL WRITE_0D_RL( CDL,   INDEX_NONE,'AIM_SFL:    CDL =',
d676f916b2 Jean*0532      &  '   /* drag coefficient for momentum over land */')
4da4b49499 Jean*0533        CALL WRITE_0D_RL( CDS,   INDEX_NONE,'AIM_SFL:    CDS =',
d676f916b2 Jean*0534      &  '   /* drag coefficient for momentum over sea  */')
4da4b49499 Jean*0535        CALL WRITE_0D_RL( CHL,   INDEX_NONE,'AIM_SFL:    CHL =',
d676f916b2 Jean*0536      &  '   /* heat exchange coefficient over land */')
4da4b49499 Jean*0537        CALL WRITE_0D_RL( CHS,   INDEX_NONE,'AIM_SFL:    CHS =',
d676f916b2 Jean*0538      &  '   /* heat exchange coefficient over sea  */')
4da4b49499 Jean*0539        CALL WRITE_0D_RL( VGUST, INDEX_NONE,'AIM_SFL:  VGUST =',
d676f916b2 Jean*0540      &  '   /* wind speed [m/s] for sub-grid-scale gusts */')
4da4b49499 Jean*0541        CALL WRITE_0D_RL( CTDAY, INDEX_NONE,'AIM_SFL:  CTDAY =',
d676f916b2 Jean*0542      &  '   /* daily-cycle correction (dTskin/dSSRad) */')
4da4b49499 Jean*0543        CALL WRITE_0D_RL( DTHETA,INDEX_NONE,'AIM_SFL: DTHETA =',
d676f916b2 Jean*0544      &  '   /* Pot.Temp. gradient for stability correction */')
4da4b49499 Jean*0545        CALL WRITE_0D_RL( dTstab,INDEX_NONE,'AIM_SFL: dTstab =',
e749d70ece Jean*0546      &  '   /* Pot.Temp. increment for stab.funct. derivative */')
4da4b49499 Jean*0547        CALL WRITE_0D_RL( FSTAB, INDEX_NONE,'AIM_SFL:  FSTAB =',
d676f916b2 Jean*0548      &  '   /* Amplitude of stability correction (fract.) */')
4da4b49499 Jean*0549        CALL WRITE_0D_RL( HDRAG, INDEX_NONE,'AIM_SFL:  HDRAG =',
d676f916b2 Jean*0550      &  '   /* Height scale for orographic correction */')
4da4b49499 Jean*0551        CALL WRITE_0D_RL( FHDRAG,INDEX_NONE,'AIM_SFL: FHDRAG =',
d676f916b2 Jean*0552      &  '   /* Amplitude of orographic correction (fract.)*/')
                0553
                0554 C- namelist AIM_PAR_CNV:
4da4b49499 Jean*0555        CALL WRITE_0D_RL( PSMIN, INDEX_NONE,'AIM_CNV:  PSMIN =',
d676f916b2 Jean*0556      &  '   /* min norm.sfc.P. for occurrence of convect. */')
4da4b49499 Jean*0557        CALL WRITE_0D_RL( TRCNV, INDEX_NONE,'AIM_CNV:  TRCNV =',
d676f916b2 Jean*0558      &  '   /* time of relaxation [h] towards ref state */')
4da4b49499 Jean*0559        CALL WRITE_0D_RL( QBL,   INDEX_NONE,'AIM_CNV:    QBL =',
d676f916b2 Jean*0560      &  '   /* specific hum. threshold in the Bound.layer */')
4da4b49499 Jean*0561        CALL WRITE_0D_RL( RHBL,  INDEX_NONE,'AIM_CNV:   RHBL =',
d676f916b2 Jean*0562      &  '   /* relative hum. threshold in the Bound.layer */')
4da4b49499 Jean*0563        CALL WRITE_0D_RL( RHIL,  INDEX_NONE,'AIM_CNV:   RHIL =',
d676f916b2 Jean*0564      &  '   /* rel.hum. threshold (intern) for 2nd mass.Flx*/')
4da4b49499 Jean*0565        CALL WRITE_0D_RL( ENTMAX,INDEX_NONE,'AIM_CNV: ENTMAX =',
d676f916b2 Jean*0566      &  '   /* time of relaxation [h] towards neutral eq. */')
4da4b49499 Jean*0567        CALL WRITE_0D_RL( SMF,   INDEX_NONE,'AIM_CNV:    SMF =',
d676f916b2 Jean*0568      &  '   /* ratio 2ndary/primary mass.Flx at cloud-base*/')
                0569
                0570 C- namelist AIM_PAR_LSC:
4da4b49499 Jean*0571        CALL WRITE_0D_RL( TRLSC, INDEX_NONE,'AIM_LSC:  TRLSC =',
d676f916b2 Jean*0572      &  '   /* relaxation time [h] for supersat. spec.hum. */')
4da4b49499 Jean*0573        CALL WRITE_0D_RL( RHLSC, INDEX_NONE,'AIM_LSC:  RHLSC =',
d676f916b2 Jean*0574      &  '   /* Max rel.humidity threshold (at sigma=1) */')
4da4b49499 Jean*0575        CALL WRITE_0D_RL( DRHLSC,INDEX_NONE,'AIM_LSC: DRHLSC =',
d676f916b2 Jean*0576      &  '   /* Vertical range of rel.humidity threshold */')
4da4b49499 Jean*0577        CALL WRITE_0D_RL( QSMAX, INDEX_NONE,'AIM_LSC:  QSMAX =',
d676f916b2 Jean*0578      &  '   /* Define the maximum latent heat release */')
                0579
                0580 C- namelist AIM_PAR_RAD:
4da4b49499 Jean*0581        CALL WRITE_0D_RL( RHCL1, INDEX_NONE,'AIM_RAD:  RHCL1 =',
d676f916b2 Jean*0582      &  '   /*  rel.hum. corresponding to cloud cover = 0 */')
4da4b49499 Jean*0583        CALL WRITE_0D_RL( RHCL2, INDEX_NONE,'AIM_RAD:  RHCL2 =',
d676f916b2 Jean*0584      &  '   /*  rel.hum. corresponding to cloud cover = 1 */')
4da4b49499 Jean*0585        CALL WRITE_0D_RL( QACL1, INDEX_NONE,'AIM_RAD:  QACL1 =',
d676f916b2 Jean*0586      &  '   /* spec.hum. threshold for cloud, upper trop. */')
4da4b49499 Jean*0587        CALL WRITE_0D_RL( QACL2, INDEX_NONE,'AIM_RAD:  QACL2 =',
d676f916b2 Jean*0588      &  '   /* spec.hum. threshold for cloud, upper trop. */')
4da4b49499 Jean*0589        CALL WRITE_0D_RL( ALBCL, INDEX_NONE,'AIM_RAD:  ALBCL =',
d676f916b2 Jean*0590      &  '   /*  cloud albedo (for cloud cover = 1) */')
4da4b49499 Jean*0591        CALL WRITE_0D_RL( EPSSW, INDEX_NONE,'AIM_RAD:  EPSSW =',
d24aab1d8d Jean*0592      &  '   /*  fract. of inc.solar rad. absorbed by ozone */')
4da4b49499 Jean*0593        CALL WRITE_0D_RL( EPSLW, INDEX_NONE,'AIM_RAD:  EPSLW =',
d24aab1d8d Jean*0594      &  '   /*  fract. of sfc LW emitted directly to space */')
4da4b49499 Jean*0595        CALL WRITE_0D_RL( EMISFC,INDEX_NONE,'AIM_RAD: EMISFC =',
d24aab1d8d Jean*0596      &  '   /*  longwave surface emissivity */')
d676f916b2 Jean*0597
                0598        WRITE(msgBuf,'(A)')
                0599      &  ' AIM_RAD : ShortWave absorptivities (for dp = 10^5 Pa) :'
2a4f7f3ad0 Jean*0600       CALL PRINT_MESSAGE( msgBuf, iUnit, SQUEEZE_RIGHT, myThid )
4da4b49499 Jean*0601        CALL WRITE_0D_RL( ABSDRY,INDEX_NONE,'AIM_RAD: ABSDRY =',
d676f916b2 Jean*0602      &  '   /* abs. of dry air (visible band) */')
4da4b49499 Jean*0603        CALL WRITE_0D_RL( ABSAER,INDEX_NONE,'AIM_RAD: ABSAER =',
d676f916b2 Jean*0604      &  '   /* abs. of aerosols (visible band) */')
4da4b49499 Jean*0605        CALL WRITE_0D_RL( ABSWV1,INDEX_NONE,'AIM_RAD: ABSWV1 =',
d676f916b2 Jean*0606      &  '   /* abs. of water vap. (vis. band)(/dq, 1g/kg) */')
4da4b49499 Jean*0607        CALL WRITE_0D_RL( ABSWV2,INDEX_NONE,'AIM_RAD: ABSWV2 =',
d676f916b2 Jean*0608      &  '   /* abs. of water vap.(near IR bd)(/dq, 1g/kg) */')
4da4b49499 Jean*0609        CALL WRITE_0D_RL( ABSCL1,INDEX_NONE,'AIM_RAD: ABSCL1 =',
d676f916b2 Jean*0610      &  '   /* abs. of clouds (vis.band)(constant term) */')
4da4b49499 Jean*0611        CALL WRITE_0D_RL( ABSCL2,INDEX_NONE,'AIM_RAD: ABSCL2 =',
d676f916b2 Jean*0612      &  '   /* abs. of clouds (vis.band) (/dw, 1g/kg) */')
                0613
                0614        WRITE(msgBuf,'(A)')
                0615      &  ' AIM_RAD : LongWave absorptivities (per dp = 10^5 Pa) :'
4da4b49499 Jean*0616        CALL WRITE_0D_RL( ABLWIN,INDEX_NONE,'AIM_RAD: ABLWIN =',
d676f916b2 Jean*0617      &  '   /* abs. of air in "window" band */')
4da4b49499 Jean*0618        CALL WRITE_0D_RL( ABLCO2,INDEX_NONE,'AIM_RAD: ABLCO2 =',
d676f916b2 Jean*0619      &  '   /* abs. of air in CO2 band */')
4da4b49499 Jean*0620        CALL WRITE_0D_RL( ABLWV1,INDEX_NONE,'AIM_RAD: ABLWV1 =',
d676f916b2 Jean*0621      &  '   /* abs. of Water vap. H2O bd-1(weak)(/dq,1g/kg)*/')
4da4b49499 Jean*0622        CALL WRITE_0D_RL( ABLWV2,INDEX_NONE,'AIM_RAD: ABLWV2 =',
d676f916b2 Jean*0623      &  '   /* abs. of W. vap., H2O bd-2(strong)(/dq,1g/kg)*/')
4da4b49499 Jean*0624        CALL WRITE_0D_RL( ABLCL1,INDEX_NONE,'AIM_RAD: ABLCL1 =',
d676f916b2 Jean*0625      &  '   /* abs. of clouds in window band (const term) */')
4da4b49499 Jean*0626        CALL WRITE_0D_RL( ABLCL2,INDEX_NONE,'AIM_RAD: ABLCL2 =',
d676f916b2 Jean*0627      &  '   /* abs. of clouds in window band (/dw, 1g/kg) */')
                0628
                0629 C- namelist AIM_PAR_VDI:
4da4b49499 Jean*0630        CALL WRITE_0D_RL( TRSHC, INDEX_NONE,'AIM_VDI:  TRSHC =',
d676f916b2 Jean*0631      &  '   /* relaxation time [h] for shallow convection */')
4da4b49499 Jean*0632        CALL WRITE_0D_RL( TRVDI, INDEX_NONE,'AIM_VDI:  TRVDI =',
d676f916b2 Jean*0633      &  '   /* relaxation time [h] for moisture diffusion */')
4da4b49499 Jean*0634        CALL WRITE_0D_RL( TRVDS, INDEX_NONE,'AIM_VDI:  TRVDS =',
d676f916b2 Jean*0635      &  '   /* relaxation time [h] for super-adiab. cond. */')
4da4b49499 Jean*0636        CALL WRITE_0D_RL( RHGRAD,INDEX_NONE,'AIM_VDI: RHGRAD =',
d676f916b2 Jean*0637      &  '   /* max gradient of rel.humidity (d_RH/d_sigma)*/')
4da4b49499 Jean*0638        CALL WRITE_0D_RL( SEGRAD,INDEX_NONE,'AIM_VDI: SEGRAD =',
d676f916b2 Jean*0639      &  '   /* max grad. of dry static Energy(d_DSE/d_phi)*/')
                0640
2a4f7f3ad0 Jean*0641       IF ( errCount.GE.1 ) THEN
                0642         WRITE(msgBuf,'(A,I3,A)')
                0643      &       'AIM_READPARMS: detected', errCount,' fatal error(s)'
                0644         CALL PRINT_ERROR( msgBuf, myThid )
                0645         CALL ALL_PROC_DIE( 0 )
                0646         STOP 'ABNORMAL END: S/R AIM_READPARMS'
                0647       ENDIF
4c105e462e Davi*0648
d676f916b2 Jean*0649 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
                0650
2a4f7f3ad0 Jean*0651 C- Initialise CO2 concentration (a better place would be from AIM_INIT_VARIA):
                0652       aim_pCO2 = aim_fixed_pCO2
                0653 #ifdef ALLOW_AIM_CO2
                0654       atm_pCO2 = atmpCO2init
                0655 #endif /* ALLOW_AIM_CO2 */
                0656
d676f916b2 Jean*0657       _END_MASTER(myThid)
d24aab1d8d Jean*0658
d676f916b2 Jean*0659 C--   Everyone else must wait for the parameters to be loaded
                0660       _BARRIER
                0661
                0662 #endif /* ALLOW_AIM */
                0663
                0664       RETURN
                0665       END