!-----------------------------------------------------------------------
!
! This README details namelist parameters for RegCM V4
!
!-----------------------------------------------------------------------
!
! DOMAIN dimension Stanza
!
 &dimparam 
 iy     = 34,   ! This is number of points in the N/S direction
 jx     = 48,   ! This is number of points in the E/W direction
 kz     = 18,   ! Number of vertical levels
 dsmin  = 0.01, ! Minimum sigma spacing (only used if kz is not 14, 18, 23, 41)
 dsmax  = 0.05, ! Maximum sigma spacing (only used if kz is not 14, 18, 23, 41)
 nsg    = 1,    ! For subgridding, number of points to decompose. If nsg=1,
                ! no subgridding is performed. CLM does NOT work as of now with
                ! subgridding enabled.
 njxcpus = -1,  ! Number of CPUS to be used in the jx (lon) dimension.
                ! If <=0 , the executable will try to figure out a suitable
                ! decomposition.
 niycpus = -1,  ! Number of CPUS to be used in the iy (lat) dimension.
                ! If <=0 , the executable will try to figure out a suitable
                ! decomposition.
 /
!
! Choice of the dynamical core
!
 &coreparam
 idynamic = 1,  ! Choice of dynamical core
                ! 1 = MM4 hydrostatical core
                ! 2 = MM5 NON hydrostatical core
 /
!
! DOMAIN geolocation Stanza
!
 &geoparam
 iproj = 'LAMCON', ! Domain cartographic projection. Supported values are:
                   ! 'LAMCON', Lambert conformal.
                   ! 'POLSTR', Polar stereographic.
                   ! 'NORMER', Normal  Mercator.
                   ! 'ROTMER', Rotated Mercator.
 ds = 60.0,        ! Grid point horizontal resolution in km
 ptop = 5.0,       ! Pressure of model top in cbar
 clat = 45.39,     ! Central latitude  of model domain in degrees
                   ! North hemisphere is positive
 clon = 13.48,     ! Central longitude of model domain in degrees
                   ! West is negative.
 plat = 45.39,     ! Pole latitude (only for rotated Mercator Proj)
 plon = 13.48,     ! Pole longitude (only for rotated Mercator Proj)
 truelatl = 30.0,  ! Lambert true latitude (low latitude side)
 truelath = 60,    ! Lambert true latitude (high latitude side)
 i_band = 0,       ! Use this to enable a tropical band. In this case the ds,
                   ! iproj, clat, clon parameters are not considered.
 /
!
! DOMAIN terrain generation parameters
!
 &terrainparam
 domname  = 'AQWA',           ! Name of the domain/experiment.
                              ! Controls naming of input files
 smthbdy = .true.,            ! Smoothing Control flag
                              !  true -> Perform extra smoothing in boundaries
 lakedpth    = .false.,       ! If using lakemod (see below), produce from
                              ! terrain program the domain bathymetry
 lsmoist     = .false.,       ! Use Satellite Soil Moisture Dataset for
                              ! initialization of soil moisture.
 fudge_lnd   = .false.,       ! Fudging Control flag, for landuse of grid 
 fudge_lnd_s = .false.,       ! Fudging Control flag, for landuse of subgrid
 fudge_tex   = .false.,       ! Fudging Control flag, for texture of grid
 fudge_tex_s = .false.,       ! Fudging Control flag, for texture of subgrid
 fudge_lak   = .false.,       ! Fudging Control flag, for lake of grid
 fudge_lak_s = .false.,       ! Fudging Control flag, for lake of subgrid
 h2opct = 50.,                ! Surface min H2O percent to be considered water
 h2ohgt = .true.,             ! Allow water points to have hgt greater than 0
 ismthlev = 1,                ! How many times apply the 121 smoothing
 dirter = 'input/',           ! Output directory for terrain files
 inpter = 'globdata/',        ! Input directory for SURFACE dataset
 moist_filename = 'moist.nc', ! Read initial moisture and snow from this file
 /
!
! Debug I/O printout level.
!
 &debugparam
 debug_level = 0, ! Currently value of 2 and 3 control previous DIAG flag
 dbgfrq = 3,      ! Interval for printout if debug_level >= 3
 /
!
! Buffer Zone Control relaxation + diffusion term
!
 &boundaryparam
 nspgx  = 12, ! nspgx-1 represent the number of cross point slices on
              ! the boundary sponge or relaxation boundary conditions.
 nspgd  = 12, ! nspgd-1 represent the number of dot point slices on 
              ! the boundary sponge or relaxation boundary conditions.
 high_nudge =    3.0D0, ! Nudge value high range
 medium_nudge =  2.0D0, ! Nudge value medium range
 low_nudge =     1.0D0  ! Nudge value low range
 /
!
! ICBC Global data input control
!
 &globdatparam
 ibdyfrq =     6,            ! boundary condition interval (hours)
 ssttyp = 'OI_WK',           ! Type of Sea Surface Temperature used
                             !  One in: GISST, OISST, OI2ST, OI_WK, OI2WK,
                             !          FV_A2, FV_B2, EH5A2, EH5B1, EHA1B,
                             !          EIN75, EIN15, CCSST, CA_XX, HA_XX,
                             !          EC_XX, IP_XX, GF_XX, CN_XX, MP_XX
 dattyp = 'EIN15',           ! Type of global analysis datasets used
                             !  One in: ECMWF, ERA40, EIN75, EIN15, EIN25,
                             !          ERAHI, NNRP1, NNRP2, NRP2W, GFS11,
                             !          FVGCM, FNEST, EH5A2, EH5B1, EHA1B,
                             !          CCSMN, ECEXY, CA_XX, HA_XX, EC_XX,
                             !          IP_XX, GF_XX, CN_XX, MP_XX
                             ! with XX for CMIP5 datasets in 26, 45, 85
 chemtyp = 'MZ6HR',          ! Type of Global Chemistry boundary conditions
                             ! One in : MZ6HR, 6 hours MOZART output
                             !        : MZCLM, MOZART climatology 1999-2009
 gdate1 = 1990060100,        ! Start date for ICBC data generation
 gdate2 = 1990070100,        ! End data for ICBC data generation
 calendar = 'gregorian',     ! Calendar type : gregorian, noleap, 360_day
 dirglob = 'input/',         ! Path for ICBC produced input files
 inpglob = 'globdata/',      ! Path for ICBC global input datasets.
 ensemble_run = .false.,     ! If this is a member of a perturbed ensemble
                             ! run. Activate random noise added to input
                             ! ICBC controlled by the perturbparam stanza
                   ! Look http://users.ictp.it/~pubregcm/RegCM4/globedat.htm
                   ! on how to download them.
 /
!
! Nesting control
!
 &fnestparam
 coarse_outdir = 'globdata/RegCM', ! Coarse domain output dir if FNEST
 coarse_domname = 'EUROPE',        ! Coarse domain domname
 /
!
! Perturbation control for ensembles
!
 &perturbparam
 lperturb_topo = .false.,     ! Add perturbation to surface elevation
 perturb_frac_topo = 0.001D0, ! Fractional value of the perturbation on topo
 lperturb_ts = .false.,       ! Add perturbation to surface temeprature
 perturb_frac_ts = 0.001D0,   ! Fractional value of the perturbation on ts
 lperturb_ps = .false.,       ! Add perturbation to surface pressure
 perturb_frac_ps = 0.001D0,   ! Fractional value of the perturbation on ps
 lperturb_t  = .false.,       ! Add perturbation to temperature
 perturb_frac_t  = 0.001D0,   ! Fractional value of the perturbation on t
 lperturb_q  = .false.,       ! Add perturbation to humidity mixing ratio
 perturb_frac_q  = 0.001D0,   ! Fractional value of the perturbation on q
 lperturb_u  = .false.,       ! Add perturbation to zonal velocity
 perturb_frac_u  = 0.001D0,   ! Fractional value of the perturbation on u
 lperturb_v  = .false.,       ! Add perturbation to meridional velocity
 perturb_frac_v  = 0.001D0,   ! Fractional value of the perturbation on v
 /
!
! Model start/restart control
!
 &restartparam
 ifrest  = .false. ,   ! If a restart
 mdate0  = 1990060100, ! Global start (is globidate1)
 mdate1  = 1990060100, ! Start date of this run
 mdate2  = 1990060200, ! End date for this run
 /
!
! Model timing parameters
! All values except dt can be zero. Convenient values are chosen by the model.
!
 &timeparam
 dt     =   150.,   ! time step in seconds
 dtrad  =     0.,   ! time interval solar radiation calculated (minutes)
 dtabem =     0.,   ! time interval absorption-emission calculated (hours)
 dtsrf  =     0.,   ! time interval at which land model is called (seconds)
 dtcum  =   -1.0,   ! time interval at which cumuls is called (seconds)
                    ! By default use same timestep as dynamical core.
 dtche  =   900.,   ! time interval at which chem model is called (seconds)
 /
!
! Model Output control
!
 &outparam
 ifsave  = .true. ,           ! Create SAV files for restart
 savfrq  =     0.,            ! Frequency in hours to create them (0 = monthly)
 ifatm   = .true. ,           ! Output ATM ?
 atmfrq  =     6.,            ! Frequency in hours to write to ATM
 ifrad   = .true. ,           ! Output RAD ?
 radfrq  =     6.,            ! Frequency in hours to write to RAD
 ifsts   = .true. ,           ! Output STS (frequence is daily) ?
 ifsrf   = .true. ,           ! Output SRF ?
 srffrq  =     3.,            ! Frequency in hours to write to SRF
 ifsub   = .true. ,           ! Output SUB ?
 subfrq  =     6.,            ! Frequency in hours to write to SUB
 iflak   = .true.,            ! Output LAK ?
 lakfrq  =     6.,            ! Frequency in hours to write to LAK
 ifchem  = .true.,            ! Output CHE ?
 ifopt   = .false.,           ! Output OPT ?
 chemfrq =     6.,            ! Frequency in hours to write to CHE
 enable_atm_vars = 67*.true., ! Mask to eventually disable variables ATM
 enable_srf_vars = 35*.true., ! Mask to eventually disable variables SRF
 enable_rad_vars = 25*.true., ! Mask to eventually disable variables RAD
 enable_sub_vars = 18*.true., ! Mask to eventually disable variables SUB
 enable_sts_vars = 18*.true., ! Mask to eventually disable variables STS
 enable_lak_vars = 18*.true., ! Mask to eventually disable variables LAK
 enable_opt_vars = 19*.true., ! Mask to eventually disable variables OPT
 enable_che_vars = 26*.true., ! Mask to eventually disable variables CHE
 dirout  = './output',        ! Path where all output will be placed
 lsync   = .false.,           ! If sync of output files at every timestep is
                              ! requested. Note, it has a performance impact.
                              ! Enabled by default if debug_level > 2
 icosp = 0,                   ! Enable COSP needed fields in output in the ATM
 idiag = 0,                   ! Enable tendency diagnostic output in the ATM
                              ! file. NOTE: output file gets HUGE.
 do_parallel_netcdf_in  = .false., ! This enables paralell input
                                   ! Each processors reads its slice in the
                                   ! input file. Enable ONLY in case of
                                   ! HUGE input bandwidth,
 do_parallel_netcdf_out = .false., ! This enables paralell output if the 
                                   ! hdf5/netcdf libraries support it and
                                   ! the model is compiled with :
                                   !    --enable-nc4-parallel
 /
!
! Model Physics
!
 &physicsparam
 iboudy  =          5,  ! Lateral Boundary conditions scheme
                        !   0 => Fixed
                        !   1 => Relaxation, linear technique.
                        !   2 => Time-dependent
                        !   3 => Time and inflow/outflow dependent.
                        !   4 => Sponge (Perkey & Kreitzberg, MWR 1976)
                        !   5 => Relaxation, exponential technique.
 isladvec =         0,  ! Semilagrangian advection scheme for tracers and
                        ! humidity
                        !   0 => Disabled
                        !   1 => Enable Semi Lagrangian Scheme
   iqmsl =          1,  ! Quasi-monotonic Semi Lagrangian
                        !   0 => Standard Semi-Lagrangian
                        !   1 => Bermejo and Staniforth 1992 QMSL scheme
 ibltyp  =          1,  ! Boundary layer scheme
                        !   0 => Frictionless
                        !   1 => Holtslag PBL (Holtslag, 1990)
                        !   2 => UW PBL (Bretherton and McCaa, 2004)
 icup_lnd =          4, ! Cumulus convection scheme Over Land
 icup_ocn =          4, ! Cumulus convection scheme Over Icean
                        !   1 => Kuo
                        !   2 => Grell
                        !   3 => Betts-Miller (1986) DOES NOT WORK !!!
                        !   4 => Emanuel (1991)
                        !   5 => Tiedtke (1996)
                        !   6 => Kain-Fritsch (1990), Kain (2004)
 ipptls  =          1,  ! Moisture scheme
                        !   1 => Explicit moisture (SUBEX; Pal et al 2000)
                        !   2 => Explicit moisture Nogherotto/Tompkins
 iocncpl =          0,  ! Ocean SST from coupled Ocean Model through RegESM
                        !   1 => Coupling activated
 iwavcpl =          0,  ! Ocean roughness from coupled Wave Model through RegESM
                        !   1 => Coupling activated
 iocnflx =          2,  ! Ocean Flux scheme
                        !   1 => Use BATS1e Monin-Obukhov
                        !   2 => Zeng et al (1998)
                        !   3 => Coare bulk flux algorithm
   iocnrough =      1,  ! Zeng Ocean model roughness formula to use.
                        !   1 => (0.0065*ustar*ustar)/egrav
                        !   2 => (0.013*ustar*ustar)/egrav + 0.11*visa/ustar
                        !   3 => (0.017*ustar*ustar)/egrav
                        !   4 => Huang 2012 free convection and swell effects
                        !   5 => four regime formulation
   iocnzoq =        1,  ! Zeng Ocean model factors for t,q roughness
                        !   1 => 2.67*(re**d_rfour) - 2.57
                        !   2 => min(4.0e-4, 2.0e-4*re**(-3.3))
                        !   3 => COARE formulation as in bulk flux above
 ipgf    =          0,  ! Pressure gradient force scheme
                        !   0 => Use full fields
                        !   1 => Hydrostatic deduction with pert. temperature
 iemiss  =          0,  ! Use computed long wave emissivity
 lakemod =          0,  ! Use lake model
 ichem   =          0,  ! Use active aerosol chemical model
 scenario =      'A1B', ! IPCC Scenario to use in A1B,A2,B1,B2
                        ! RCP Scenarios in RPC2.6,RCP4.5,RCP6,RCP8.5
                        ! CONST scenario at year ghg_year_const
 ghg_year_const = 1950, ! Year to use for a constant GHG concentration values
 idcsst   =          0, ! Use diurnal cycle sst scheme
 iseaice  =          0, ! Model seaice effects
 idesseas =          0, ! Model desert seasonal albedo variability
 iconvlwp =          1, ! Use convective algo for lwp in the large-scale
                        ! This is reset to zero if using ipptls = 2
 icldfrac =          0, ! Cloud fraction algorithm
                        !   0 : Original SUBEX
                        !   1 : Xu-Randall empirical
 icldmstrat =        1, ! Simulate stratocumulus clouds
 icumcloud =         1, ! Formulas to use for cumulus clouds (cf and lwc)
                        ! Cloud fractions, only if mass fluxes are not
                        ! available (Kuo and BM):
                        !   0,1 => cf = 1-(1-clfrcv)**(1/kdepth)
                        !   2   => cf = cloud profile
                        ! Liquid water content:
                        !   0   => constant in cloud
                        !   1,2 => function of temperature
 irrtm    =          0, ! Use RRTM radiation scheme instead of CCSM
 iclimao3 =          0, ! Use O3 climatic dataset from SPARC CMIP5
 isolconst =         0, ! Use a constant 1367 W/m^2 instead of the prescribed
                        ! TSI recommended CMIP5 solar forcing data.
 islab_ocean =       0, ! Activate the SLAB ocean model
 itweak =            0, ! Enable tweak scenario
 /
!
! Non-hydrostatic core option
!
 &nonhydroparam
 base_state_pressure = 101325.0,  ! Base state reference pressure
 base_state_temperature = 288.15, ! Base state reference temperature
 logp_lrate = 50.0,               ! Logp lapse rate d(T)/d(ln P) [K/ln(Pa)]
 ifupr = 0,         ! Upper radiative boundary condition (Klemp and Durran,
                    ! Bougeault, 1983)
 ckh = 1.0,         ! Background diffusion multiplication factor
 diffu_hgtf = 1,    ! Add topographic effect to diffusion
 nhbet = 0.4,       ! Ikawa beta parameter (0.=centered, 1.=backward)
                    ! determines the time-weighting, where zero gives a
                    ! time-centered average and positive values give a bias
                    ! towards the future time step that can be used for
                    ! acoustic damping. In practice, values of
                    ! nhbet = 0.2 - 0.4 are used (MM5 manual, Sec. 2.5.1)
 nhxkd = 0.1,       ! Time weighting for weighting old/new pp
 /
!
! Moisture scheme comtrol (ipptls = 1)
!
 &subexparam
 ncld      = 1,      ! # of bottom model levels with no clouds (rad only)
 qck1land  = 0.0005, ! Autoconversion Rate for Land
 qck1oce   = 0.0005, ! Autoconversion Rate for Ocean
 gulland   = 0.65,   ! Fract of Gultepe eqn (qcth) when prcp occurs (land)
 guloce    = 0.30,   ! Fract of Gultepe eqn (qcth) for ocean
 rhmax     = 1.01,   ! RH at whicn FCC = 1.0
 rhmin     = 0.01,   ! RH min value
 rh0land   = 0.80,   ! Relative humidity threshold for land
 rh0oce    = 0.90,   ! Relative humidity threshold for ocean
 tc0       = 238.0,  ! Below this temp, rh0 begins to approach unity
 cevaplnd  = 1.0e-5, ! Raindrop evap rate coef land [[(kg m-2 s-1)-1/2]/s]
 cevapoce  = 1.0e-5, ! Raindrop evap rate coef ocean [[(kg m-2 s-1)-1/2]/s]
 caccrlnd  = 6.0,    ! Raindrop accretion rate land  [m3/kg/s]
 caccroce  = 6.0,    ! Raindrop accretion rate ocean [m3/kg/s]
 cllwcv    = 0.3e-3, ! Cloud liquid water content for convective precip.
 clfrcvmax = 0.75,   ! Max cloud fractional cover for convective precip.
 cftotmax  = 0.75,   ! Max total cover cloud fraction for radiation
 conf      = 1.00,   ! Condensation efficiency
 rcrit     = 13.5,   ! Mean critical radius
 coef_ccn  = 2.0,    ! Geometric mean Diameter and standard deviation
 abulk     = 0.9,    ! Bulk activation ratio
 lsrfhack  = .false. ! Surface radiation hack
 /
!
! New Microphysics Option (ipptls = 2)
!
 &microparam
 stats = .false.,           ! Produce debug variables in output files
 budget_compute = .false.,  ! Verify enthalpy and moisture conservation
 nssopt = 1,                ! Supersaturation Computation
                            ! 0 => No scheme
                            ! 1 => Tompkins
                            ! 2 => Lohmann and Karcher
                            ! 3 => Gierens
 iautoconv = 4,             !  Choose the autoconversion paramaterization
                            ! => 1 Klein & Pincus (2000)
                            ! => 2 Khairoutdinov and Kogan (2000)
                            ! => 3 Kessler (1969)
                            ! => 4 Sundqvist
 rsemi = 1.0,               ! Implicit/Explicit control 
                            ! NOT ACTIVATED YET - IT DOES NOT WORK!
                            ! rsemi == 0 => scheme is fully explicit
                            ! rsemi == 1 => scheme is fully implicit
                            ! 0<rsemi<1  => scheme is semi-implicit
 vfqr = 4.0,                ! Rain fall speed (default is 4 m/s)
 vfqi = 0.15,               ! Ice fall speed (default is 0.15 m/s)
 vfqs = 1.0,                ! Snow fall speed (default is 1 m/s)
 auto_rate_khair = 0.355,   ! Autoconversion coefficient for kautoconv=2
 auto_rate_kessl = 1.e-3,   ! Autoconversion coefficient for kautoconv=3
 auto_rate_klepi = 0.5e-3,  ! Autoconversion coefficient for kautoconv=1
 rkconv = 1.666e-4,         ! Autoconversion coefficient for kautoconv=4
 rcovpmin = 0.1,            ! Minimum precipitation coverage
 rpecons = 5.547e-5,        ! Evaporation constant Kessler
 /
!
! Grell Cumulus scheme control
!
 &grellparam
 igcc  = 2,          ! Cumulus closure scheme
                     !   1 => Arakawa & Schubert (1974)
                     !   2 => Fritsch & Chappell (1980)
 gcr0 = 0.0020,      ! Conversion rate from cloud to rain
 edtmin      = 0.20, ! Minimum Precipitation Efficiency land
 edtmin_ocn  = 0.20, ! Minimum Precipitation Efficiency ocean
 edtmax      = 0.80, ! Maximum Precipitation Efficiency land
 edtmax_ocn  = 0.80, ! Maximum Precipitation Efficiency ocean
 edtmino     = 0.20, ! Minimum Tendency Efficiency (o var) land
 edtmino_ocn = 0.20, ! Minimum Tendency Efficiency (o var) ocean
 edtmaxo     = 0.80, ! Maximum Tendency Efficiency (o var) land
 edtmaxo_ocn = 0.80, ! Maximum Tendency Efficiency (o var) ocean
 edtminx     = 0.20, ! Minimum Tendency Efficiency (x var) land
 edtminx_ocn = 0.20, ! Minimum Tendency Efficiency (x var) ocean
 edtmaxx     = 0.80, ! Maximum Tendency Efficiency (x var) land
 edtmaxx_ocn = 0.80, ! Maximum Tendency Efficiency (x var) ocean
 shrmin      = 0.30, ! Minimum Shear effect on precip eff. land
 shrmin_ocn  = 0.30, ! Minimum Shear effect on precip eff. ocean
 shrmax      = 0.90, ! Maximum Shear effect on precip eff. land
 shrmax_ocn  = 0.90, ! Maximum Shear effect on precip eff. ocean
 pbcmax = 50.0,      ! Max depth (mb) of stable layer b/twn LCL & LFC
 mincld = 150.0,     ! Min cloud depth (mb).
 htmin = -250.0,     ! Min convective heating
 htmax = 500.0,      ! Max convective heating
 skbmax = 0.4,       ! Max cloud base height in sigma
 dtauc = 30.0D0      ! Fritsch & Chappell (1980) ABE Removal Timescale (min)
 /
!
! Emanuel Cumuls scheme control
!
 &emanparam
 minsig = 0.95,       ! Lowest sigma level from which convection can originate
 elcrit_ocn = 0.0011, ! Autoconversion threshold water content (g/g) (ocean)
 elcrit_lnd = 0.0011, ! Autoconversion threshold water content (g/g) (land)
 tlcrit = -55.0,      ! Below tlcrit auto-conversion threshold is zero
 entp = 1.5,          ! Coefficient of mixing in the entrainment formulation
 sigd = 0.05,         ! Fractional area covered by unsaturated dndraft
 sigs = 0.12,         ! Fraction of precipitation falling outside of cloud
 omtrain = 50.0,      ! Fall speed of rain (Pa/s)
 omtsnow = 5.5,       ! Fall speed of snow (Pa/s)
 coeffr = 1.0,        ! Coefficient governing the rate of rain evaporation
 coeffs = 0.8,        ! Coefficient governing the rate of snow evaporation
 cu = 0.7,            ! Coefficient governing convective momentum transport
 betae = 10.0,        ! Controls downdraft velocity scale
 dtmax = 0.9,         ! Max negative parcel temperature perturbation below LFC
 alphae = 0.2,        ! Controls the approach rate to quasi-equilibrium
 damp = 0.1,          ! Controls the approach rate to quasi-equilibrium
 epmax_ocn = 0.999,   ! Maximum precipitation efficiency (ocean)
 epmax_lnd = 0.999,   ! Maximum precipitation efficiency (land)
 /
!
! Tiedtke Cumulus scheme control
!
 &tiedtkeparam
 iconv = 4,           ! Actual used scheme.
 entrmax = 1.75e-3,   ! Max entrainment iconv=[1,2,3]
 entrdd  = 3.0e-4,    ! Entrainment rate for cumulus downdrafts
 entrpen = 1.75e-3,   ! Entrainment rate for penetrative convection
 entrscv = 3.0e-4,    ! Entrainment rate for shallow convection iconv=[1,2,3]
 entrmid = 1.0e-4,    ! Entrainment rate for midlevel convection iconv=[1,2,3]
 cprcon = 1.0e-4,     ! Coefficient for determining conversion iconv=[1,2,3]
 detrpen = 0.75e-4,   ! Detrainment rate for penetrative convection iconv=4
 entshalp = 2.0,      ! shallow entrainment factor for entrorg iconv=4
 rcuc_lnd = 0.05,     ! Convective cloud cover for rain evporation iconv=4
 rcuc_ocn = 0.05,     ! Convective cloud cover for rain evporation iconv=4
 rcpec_lnd = 5.55e-5, ! Coefficient for rain evaporation below cloud iconv=4
 rcpec_ocn = 5.55e-5, ! Coefficient for rain evaporation below cloud iconv=4
 rhebc_lnd = 0.7,     ! Critical rh below cloud for evaporation iconv=4
 rhebc_ocn = 0.9,     ! Critical rh below cloud for evaporation iconv=4
 rprc_lnd = 1.4e-3,   ! conversion coefficient from cloud water iconv=4
 rprc_ocn = 1.4e-3,   ! conversion coefficient from cloud water iconv=4
 cmtcape = 3600.0,    ! CAPE adjustment timescale iconv=4
 lmfpen    = .true.,  ! penetrative conv is switched on
 lmfmid    = .true.,  ! midlevel conv is switched on
 lmfdd     = .true.,  ! cumulus downdraft is switched on
 lepcld    = .true.,  ! prognostic cloud scheme is on
 lmfdudv   = .true.,  ! cumulus friction is switched on
 lmfscv    = .true.,  ! shallow convection is switched on
 lmfuvdis  = .true.,  ! use kinetic energy dissipation
 lmftrac   = .true.,  ! chemical tracer transport is on
 lmfsmooth = .false., ! smoot of mass fluxes for tracers
 lmfwstar  = .false., ! Grant w* closure for shallow conv
 /
!
! Kain Fritsch cumulus parameters
! From WRF KFeta module
!
 &kfparam
 kf_entrate = 0.03,      ! Entrainment rate
 kf_min_pef = 0.2,       ! Minimum precipitation efficiency
 kf_max_pef = 0.9,       ! Maximum precipitation efficiency
 kf_dpp     = 150.0,     ! Start elevation for downdraft above cloud base (mb)
 kf_tkemax = 5.0,        ! Maximum turbolent kinetic energy in sub cloud layer
 kf_min_dtcape = 1800.0, ! Consumption time of CAPE low limit
 kf_max_dtcape = 3600.0, ! Consumption time of CAPE high limit
 /
!
! Holtslag PBL scheme control
!
 &holtslagparam
 ricr_ocn = 0.25,  ! Critical Richardson Number over Ocean
 ricr_lnd = 0.25,  ! Critical Richardson Number over Land
 zhnew_fac = 0.25, ! Multiplicative factor for zzhnew in holtpbl
 ifaholtth10 = 1,  ! First approximation for obhukov length, th10 formula:
                   !      1 => 0.5 * (t+tg) * (1+0.61*q)
                   !      2 => (0.25*t + 0.75*tg) * (1+0.61*q)
                   !      3 => theta + hf/(k*us)*log(z/10)
                   ! t = air temp., tg = ground temp., q = wv mix. ratio
                   ! hf = total heat flux, z = elevation
                   ! theta = virt. pot. t
 ifaholt = 1,      ! th10 final adjustment:
                   !      0 => no adjustment
                   !      1 => max(th10,tg)
                   !      2 => min(th10,tg)
 /
!
! UW PBL scheme control
!
 &uwparam
 iuwvadv = 0,   ! 0=standard T/QV/QC advection, 1=GB01-style advection
                ! 1 is ideal for MSc simulation, but may have stability issues
 atwo = 15.0,   ! Efficiency of enhancement of entrainment by cloud evap.
                ! see Grenier and Bretherton (2001) Mon. Wea. Rev.
                ! and Bretherton and Park (2009) J. Clim.
 rstbl = 1.5,   ! Scaling parameter for stable boundary layer eddy length
                ! scale. Higher values mean stronger mixing in stable
                ! conditions
 czero = 5.869, ! Czero constant in UW PBL (eqn 44a and pgs 856-857)
 nuk = 5.0D0,   ! Multiplication factor for diffusion coefficients
 /
!
! SLAB Ocean model, or q-flux adjusted mixed layer model
!
 &slabocparam
 do_qflux_adj  = .false.,  ! Activate SLAB Ocean model surface fluxes adjust
                           ! from an already created climatology
 do_restore_sst = .true.,  ! Create during the run the SLAB Ocean model surface
                           ! fluxes climatology to be used in a subsequent run
 sst_restore_timescale = 5.0D0,  ! Time interval in days in building the
                                 ! q-flux adjustment
 mixed_layer_depth     = 50.0D0, ! Depth in meters of the Ocean mixed layer.
 /
!
! RRTM Radiation scheme parameters
!
 &rrtmparam
 inflgsw  = 2, ! 0 = use the optical properties calculated in prep_dat_rrtm
               !     (same as standard radiation)
               ! 2 = use RRTM option to calculate cloud optical properties
               !     from water path and cloud drop radius
 iceflgsw = 3, ! Flag for ice particle specification
               !   0 => ice effective radius, r_ec, (Ebert and Curry, 1992),
               !        r_ec must be >= 10.0 microns
               !   1 => ice effective radius, r_ec, (Ebert and Curry, 1992),
               !        r_ec range is limited to 13.0 to 130.0 microns
               !   2 => ice effective radius, r_k, (Key, Streamer Ref. Manual,
               !        1996), r_k range is limited to 5.0 to 131.0 microns
               !   3 => generalized effective size, dge, (Fu, 1996),
               !        dge range is limited to 5.0 to 140.0 microns
               !        [dge = 1.0315 * r_ec]
 liqflgsw = 1, ! Flag for liquid droplet specification
               !   0 => Compute the optical depths due to water clouds in ATM
               !        (currently not supported)
               !   1 => The water droplet effective radius (microns) is input
               !        and the optical depths due to water clouds are computed
               !        as in Hu and Stamnes, J., Clim., 6, 728-742, (1993).
 inflglw  = 2, ! Flag for cloud optical properties as above but for LW
 iceflglw = 3, ! Flag for ice particle specification as above but for LW
 liqflglw = 1, ! Flag for liquid droplet specification as above but for LW
 icld  = 1,    ! Cloud Overlap hypothesis
 irng = 1,     ! mersenne twister random generator for McICA COH
 /
!
! Chemical control parameter
!
 &chemparam
 chemsimtype = 'CBMZ    ', ! Which chemical tracers to be activated.
                           ! One in :
                           !   DUST   : Activate 4 dust bins scheme
                           !   SSLT   : Activate 2 bins Sea salt scheme
                           !   DUSS   : Activate DUST + SSLT
                           !   DU12   : Activate 12 dust bins scheme
                           !   CARB   : Activate 4 species black/anthropic
                           !            carbon simulations
                           !   SULF   : Activate SO2 and SO4 tracers
                           !   SUCA   : Activate both SUKF and CARB
                           !   AERO   : Activate all DUST, SSLT, CARB and SULF
                           !   CBMZ   : Activate gas phase and sulfate
                           !   DCCB   : Activate CBMZ +DUST +CARB
                           !   POLLEN : Activate POLLEN transport scheme
 ichsolver = 1,  ! Activate the chemical solver
 ichsursrc = 1,  ! Enable the emissions
 ichdrdepo = 1,  ! 1 = enable tracer surface dry deposition. For dust,
                 !     it is calculated by a size settling and dry
                 !     deposition scheme. For other aerosol,a dry
                 !     deposition velocity is simply prescribed further.
 ichebdy = 1,    ! Enable boundary conditions read
 ichcumtra = 1,  ! 1 = enable tracer convective transport and mixing.
 ichremlsc = 1,  ! 1 = wet removal of chemical species (washout and rainout
                 !     by total rain) is enabled
 ichremcvc = 1,  ! 1 = wet removal of chemical species (washout and rainout
                 !     by convective rain) is enabled
 ichdustemd = 1, ! Choice for parametrisation of dust emission size distribution
                 ! 1 = use the standard scheme (Alfaro et al., Zakey et al.)
                 ! 2 = use the the revised soil granulometry + Kok et al., 2011
 ichdiag = 0,    ! 1 = enable writing of additional diagnostics in the output
 idirect   = 1,  ! Choice to enable or not aerosol feedbacks on radiation and
                 ! dynamics (aerosol direct and semi direct effcts):
                 ! 1 = no coupling to dynamic and thermodynamic. However
                 !     the clear sky surface and top of atmosphere 
                 !     aerosol radiative forcings are diagnosed.
                 ! 2 = allows aerosol feedbacks on radiative,
                 !     thermodynamic and dynamic fields.
 iindirect = 0,  ! Enable sulfate indirect effect in radiation scheme
 rdstemfac = 1.0,! Aerosol correction factor (Laurent et al, 2008)
 ichjphcld = 1,  ! Impact of cloud aod on photolysis coef
 ichbion = 0,    ! ???????????????????????????????????///
 /
!
! CLM parameters
!
 &clmparam
 dirclm = 'input/', ! CLM path to Input data produced by clm2rcm. If 
                    ! relative, It should be how to reach the Input dir
                    ! from the Run dir.
 clmfrq =  12.,     ! Frequency for CLM own output write
 imask  =   1,      ! For CLM, Type of land surface parameterization
                    !   1 => using DOMAIN.INFO for landmask (same as BATS)
                    !   2 => using mksrf_navyoro file landfraction for
                    !        landmask and perform a weighted average over
                    !        ocean/land gridcells; for example:
                    ! tgb = tgb_ocean*(1-landfraction)+tgb_land*landfraction
 /
!
! Coupled model parameters
!
 &cplparam
 cpldt        = 10800., ! Coupling time step in sec. 
 zomax        = 0.005,  ! Maximum surface roughness from wave model coupling
 ustarmax     = 2.0,    ! Maximum friction velocity from wave model coupling
 cpldbglevel  = 1,      ! Debug level
                        ! 1 => basic output
                        ! 2 => grid information of each component is written
                        !      to VTK file
                        ! 3 => exchange fields are written to netcdf file at
                        !      each coupling time step 
                        ! 4 => exchange fields are written to a ASCII file at
                        !      each coupling time step
 /
!
! Tweaking parameters
!
 &tweakparam
 itweak_sst = 0,              ! Enable adding sst_tweak to input TS
 itweak_temperature = 0,      ! Enable adding temperature_tweak to input T
 itweak_solar_irradiance = 0, ! Add solar_tweak to solar constant
 itweak_greenhouse_gases = 0, ! Multiply gas_tweak_factors to GG concentrations
 sst_tweak = 0.0D0,           ! In K
 temperature_tweak = 0.0D0,   ! In K
 solar_tweak = 0.0D0,         ! In W m-2 (1367.0 is default solar)
 gas_tweak_factors = 1.0D0, 1.0D0 , 1.0D0 , 1.0D0 , 1.0D0,
 !                   CO2    CH4     N2O     CFC11   CFC12
 /
!
! CLM 4.5 parameters
!
 &clm_inparm
 fpftcon = 'pft-physiology.c130503.nc',
 fsnowoptics = 'snicar_optics_5bnd_c090915.nc',
 fsnowaging = 'snicar_drdt_bst_fit_60_c070416.nc',
 /
 &clm_soilhydrology_inparm
 h2osfcflag = 1,
 origflag = 0,
 /
 &clm_hydrology1_inparm
 oldfflag = 0,
 /
! ####################################################################
