NAM_PARAM_LIMA

It contains the options for the 2 moment mixed phase cloud parameterizations used by the model (CCLOUD = “LIMA” in NAM_PARAMn).

NAM_PARAM_LIMA content

Fortran name

Fortran type

Default value

NMOM_C

INTEGER

2

NMOM_R

INTEGER

2

NMOM_I

INTEGER

2

NMOM_S

INTEGER

1

NMOM_G

INTEGER

1

NMOM_H

INTEGER

0

LNUCL

LOGICAL

.TRUE.

LSEDI

LOGICAL

.TRUE.

LHHONI

LOGICAL

.FALSE.

LMEYERS

LOGICAL

.FALSE.

NMOD_IFN

INTEGER

1

XIFN_CONC

REAL

100.0

LIFN_HOM

LOGICAL

.TRUE.

CIFN_SPECIES

CHARACTER(LEN=8)

‘PHILLIPS’

CINT_MIXING

CHARACTER(LEN=8)

‘DM2 ‘

NMOD_IMM

INTEGER

0

NIND_SPECIE

INTEGER

1

CPRISTINE_ICE_LIMA

CHARACTER(LEN=4)

‘PLAT’

XFACTNUC_DEP

REAL

1.0

XFACTNUC_CON

REAL

1.0

NPHILLIPS

INTEGER

8

LACTI

LOGICAL

.TRUE.

LSEDC

LOGICAL

.TRUE.

LACTIT

LOGICAL

.FALSE.

NMOD_CCN

INTEGER

1

XCCN_CONC

REAL

300.0

LCCN_HOM

LOGICAL

.TRUE.

CCCN_MODES

CHARACTER(LEN=8)

‘COPT ‘

HINI_CCN

CHARACTER(LEN=3)

‘AER’

HTYPE_CCN

CHARACTER(LEN=10)

‘M’

XALPHAC

REAL

3.0

XNUC

REAL

1.0

XALPHAR

REAL

1.0

XNUR

REAL

2.0

XFSOLUB_CCN

REAL

1.0

XACTEMP_CCN

REAL

280.0

LSCAV

LOGICAL

.FALSE.

LAERO_MASS

LOGICAL

.FALSE.

LACTTKE

LOGICAL

.TRUE.

LDEPOC

LOGICAL

.TRUE.

XVDEPOC

REAL

0.02

LPTSPLIT

LOGICAL

.TRUE.

LFEEDBACKT

LOGICAL

.TRUE.

NMAXITER

INTEGER

5

XMRSTEP

REAL

0.005

XTSTEP_TS

REAL

20.0

LADJ

LOGICAL

.TRUE.

LSPRO

LOGICAL

.FALSE.

LKHKO

LOGICAL

.FALSE.

LCIBU

LOGICAL

.FALSE.

LRDSF

LOGICAL

.FALSE.

XNDEBRIS_CIBU

REAL

50.0

LMURAKAMI

LOGICAL

.TRUE.

LSNOW_T

LOGICAL

.FALSE.

LKESSLERAC

LOGICAL

.FALSE.

LICE3

LOGICAL

.FALSE.

LSIGMOIDE_G

LOGICAL

.FALSE.

LSIGMOIDE_NG

LOGICAL

.FALSE.

XSIGMOIDE_G

LOGICAL

1E8

XMVDMIN_G

LOGICAL

125E-6

LCRIAUTI

LOGICAL

.FALSE.

XPSH_MAX_RDSF

REAL

0.2

XT0CRIAUTI

LOGICAL

(LOG10(XCRIAUTI)-XBCRIAUTI)/0.06

XCRIAUTI

REAL

0.2E-4

XCRIAUTC

REAL

0.5E-3

XACRIAUTI

REAL

0.06

XBCRIAUTI

REAL

-3.5

CSUBG_PR_PDF

CHARACTER(LEN=4)

‘SIGM’

CSUBG_AUCV_RC

CHARACTER(LEN=4)

‘NONE’

CSUBG_AUCV_RI

CHARACTER(LEN=4)

‘NONE’

LCRYSTAL_SHAPE

LOGICAL

.FALSE.

NNB_CRYSTAL_SHAPE

INTEGER

1

HTYPE_CRYSTAL

CHARACTER(LEN=4)(:)

NNB_CRYSTAL_SHAPE * ‘’

LICE_ISC

LOGICAL

.FALSE.

LINITORILAM

LOGICAL

.FALSE.

LINTERP_CAMS

LOGICAL

.FALSE.

LFREEZ_RATE

LOGICAL

.TRUE.

  • NMOM_C : number of prognostic moments for cloud droplets.

  • NMOM_R : number of prognostic moments for rain drops.

  • NMOM_I : number of prognostic moments for ice crystals.

  • NMOM_S : number of prognostic moments for snow/aggregates.

  • NMOM_G : number of prognostic moments for graupel.

  • NMOM_H : number of prognostic moments for hail.

Note

Note that the full flexibility is only available with the time-splitted version of LIMA (with LPTSPLIT=T). With the original version, two configurations only are available : the original parameterization with NMOM_C=NMOM_R=NMOM_I=2 and NMOM_S=NMOM_G=1; and 2 moments for all species with NMOM_x=2.

  • LNUCL : Switch to activate pristine ice crystals nucleation (both from IFN and homogeneous freezing)

  • LSEDI : Switch to activate the sedimentation of pristine ice crystals

  • LSNOW_T : Switch to activate the representation of snow proposed by Wurtz et al. 2023 which improves the extension and cloud composition of anvils in convective systems.

  • LMURAKAMI : Switch to activate the snow riming and conversion to graupel process can be computed following Murakami (1990). Only avaiable with LPTSPLIT=T.

  • LCIBU : swith to activate the representation of collisional ice break-up (CIBU, Hoarau et al., 2018), a secondary ice production mechanism (ice is produced by the fragmentation of snow upon impact with graupel). The number of fragments formed per collision can be fixed by XNDEBRIS_CIBU.

  • XNDEBRIS_CIBU : number of fragments formed per collision in the ice collisional break-up mechanism. Negative values result in a random number of fragments formed for each collision.

  • LRDSF : switch to activate raindrop shattering freezing, a secondary ice production mechanism.

  • LKHKO : switch to activate a replicate behaviour of the KHKO scheme (useful for stratocumulus clouds and drizzle formation) with LIMA code.

  • LHHONI : Switch to activate CCN homogeneous freezing

  • LMEYERS : Switch to activate the ice nucleation parameterization by Meyers (1992) instead of using the IFN

  • NMOD_IFN : Number of IFN modes

  • XIFN_CONC : Initial reference number concentration for each IFN mode (verb?#/L?)

  • LIFN_HOM : If set to true, the initial concentration of IFN is homogeneous on the vertical. If set to false, the IFN concentration is equal to the reference value XIFN_CONC below 1000m and exponentially decreasing above.

  • CINT_MIXING : String to select the proportion of each IFN type in each IFN mode. Possible values :

    • ‘DM1’ pure small dust particles

    • ‘DM2’ pure large dust particles

    • ‘BC’ pure black carbon

    • ‘O’ pure organics

    • ‘CAMS’, ‘CAMS_JPP’, ‘CAMS_AIT’, ‘CAMS_ACC’, ‘MOCAGE’, mix for use of CAMS or MOCAGE aerosols

    • ‘DEFAULT’ mix as in Phillips et al 2008 or 2013

  • CIFN_SPECIES : String to select the IFN modes size distribution parameters. Available options are :

    • ‘MOCAGE’, ‘CAMS_JPP’, ‘CAMS_AIT’, ‘CAMS_ACC’ for use with MOCAGE/CAMS aerosols

    • ‘DEFAULT’ to use the same parameters as in Phillips et al. 2008 or 2023

  • NMOD_IMM :Number of “coated IFN” modes

  • NIND_SPECIE : Type of the “coated IFN” mode. 1 for dust, 2 for black carbon or 3 for organics

  • CPRISTINE_ICE_LIMA : Select the shape of pristine ice among:

    • ‘PLAT’ : plates

    • ‘COLU’ : columns

    • ‘BURO’ : bullet rosettes

    • ‘YPLA’ : plates from Yang et al. 2013

    • ‘YCOL’ : column from Yang et al. 2013

    • ‘YBUR’ : solid bullet rosette from Yang et al. 2013

    • ‘YDRO’ : droxtal from Yang et al. 2013

    • ‘YHCO’ : hollow column from Yang et al. 2013

    • ‘YHBU’ : hollow bullet rosette from Yang et al. 2013

  • XFACTNUC_DEP : Amplification factor for IN nucleation by deposition (only used if LMEYERS=T)

  • XFACTNUC_CON : Amplification factor for IN nucleation by contact (only used if LMEYERS=T)

  • NPHILLIPS : Version of the Phillips parameterization : 8 for the 2008 paper ; 13 for the 2013 paper

  • LACTI : Switch to activate the CCN activation

  • LSEDC : Switch to activate the cloud droplets sedimentation

  • LACTIT : Switch to activate the representation of radiative cooling in the diagnostic maximum supersaturation computation

  • NMOD_CCN : Number of CCN modes

  • XCCN_CONC : Reference concentration for each CCN mode (verb?#/cm3?)

  • LCCN_HOM : If set to true, the initial concentration of CCN is homogeneous on the vertical. If set to false, the CCN concentration is equal to the reference value XCCN_CONC below 1000m and exponentially decreasing above.

  • CCCN_MODES : Select the size distribution of CCN modes (‘JUNGFRAU’,’COPT’,’CAMS’, ‘CAMS_JPP’,’CAMS_ACC’,’CAMS_AIT’,’SIRTA’,’CPS00’,’FREETROP’)

  • HINI_CCN : Switch to use aerosols as CCN or to describe directly the CCN activation spectrum : ‘AER’ (use aerosols) or ‘CCN’ (use the CCN activation spectrum directly)

  • HTYPE_CCN : Switch to affect maritime or continental activation properties to each CCN mode : ‘C’ continental or ‘M’ maritime. NH42SO4 (=C), NH4NO3, NaCl (=M), H2SO4, NaNO3, NaHSO4, Na2SO4, NH43HSO42, SOA

  • XALPHAC,XNUC : Droplets size distribution parameter

  • XALPHAR,XNUR : Rain size distribution parameter

  • XFSOLUB_CCN : Fractional solubility of the CCN

  • XACTEMP_CCN : Expected temperature of CCN activation

  • LSCAV : Switch to activate below cloud scavenging of aerosols by rain

  • LAERO_MASS : Switch to track the mass of scavenged aerosols

  • LACTTKE : flag to use TKE in the CCN activation formulation

  • LDEPOC : flag to activate droplet deposition

  • XVDEPOC : Droplet deposition velocity

  • LPTSPLIT : flag to activate the time-splitting scheme

  • LFEEDBACKT : Flag to recompute tendencies if temperature reaches 0 (for the time-splitting scheme)

  • NMAXITER : Maximum number of iterations (for the time-splitting scheme)

  • XMRSTEP : Recompute tendencies if any mixing ratio changes by more than XMRSTEP (0=no limit) (for the time-splitting scheme)

  • XTSTEP_TS : Maximum length of sub-time-steps (for the time-splitting scheme)

  • LSPRO : flag to activate the saturation adjustement from Thouron et al. 2012 (diagnostic supersaturation)

  • LADJ : flag to use a saturation adjustment for cloud droplets (if set to T) or the “diagnostic supersaturation” from Thouron et al. 2012 (if set to F and LSPRO=F).

  • LKESSLERAC : Set to T to use the Kessler autoconversion of cloud droplets into rain drops, based on the droplets mixing ratio. Useful if NMOM_C=1 to have results closer to ICE3 simulations.

Note

Use of CAMS with LIMA. You must have set HCAMSFILE and HCAMSFILETYPE at the PREP_REAL_CASE step).

&NAM_PARAM_LIMA  HTYPE_CCN(1)='NaCl',
                  HTYPE_CCN(2)='NH42SO4',
                  HTYPE_CCN(3)='SOA',
                  NMOD_CCN=3,
                  CCCN_MODES='CAMS_AIT',
                  NMOD_IFN=2,
                  CIFN_SPECIES='CAMS_AIT',
                  CINT_MIXING='CAMS',
                  NMOD_IMM=1 /

  • LICE3: Use to mimic the ICE3 scheme. If set to .TRUE., some parameters are set :

NMOM_C=1
NMOM_R=1
NMOM_I=1
NMOM_S=1
NMOM_G=1
NMOM_H=MIN(NMOM_H,1)
NMOD_CCN=0
NMOD_IFN=0
LMURAKAMI=.TRUE.
LKESSLERAC=.TRUE.
XALPHAR=1.
XNUR=1.

  • LSIGMOIDE_G: true to limit graupel growth by XSIGMOIDE_G

  • XSIGMOIDE_G: sigmoide parameter for graupel growth limitation

  • LSIGMOIDE_NG: true to force lambda to be < lambda(Dmin)

  • XMVDMIN_G: minimum MVD for graupel growth lim or lambda(Dmin) calculation

  • LCRIAUTI: true to compute XACRIAUTI and XBCRIAUTI (from XCRIAUTI and XT0CRIAUTI). If false, XT0CRIAUTI is computed from XCRIAUTI and XBCRIAUTI.

  • XPSH_MAX_RDSF: shattering probability normal distribution maximum

  • XT0CRIAUTI: threshold temperature for the ice->snow autoconversion threshold

  • XCRIAUTI: minimum value for the ice \(\rightarrow\) snow autoconversion threshold

  • XACRIAUTI: A parameter for the ice \(\rightarrow\) snow autoconversion 10**(aT+b) law

  • XBCRIAUTI: B parameter for the ice \(\rightarrow\) snow autoconversion 10**(aT+b) law

  • CSUBG_PR_PDF: PDF for subgrid precipitation. Options are the same as in NAM_PARAM_ICEn.

  • CSUBG_AUCV_RC: type of subgrid rc->rr autoconversion method. Options are the same as in NAM_PARAM_ICEn.

  • CSUBG_AUCV_RI: type of subgrid ri->rs autoconversion method. Options are the same as in NAM_PARAM_ICEn.

  • LCRYSTAL_SHAPE: true to enable several ice crystal shapes. It can only be used if LPTSPLIT=T.

  • NNB_CRYSTAL_SHAPE: number of ice crystal shapes ; taken into account if LCRYSTAL_SHAPE=T. For the moment, only 4 ice crystal shapes are allowed.

  • HTYPE_CRYSTAL_SHAPE: ice crystal shapes if LCRYSTAL_SHAPE=T. Can be set to YPLA, YCOL, YBUR or YDRO.

  • LICE_ISC: true to enable self collection of ice crystals

  • LINITORILAM: true to initialize CCN and IF by ORILAM

  • LINTERP_CAMS: true to interpolate CAMS data at each time step (from Large-Scale fields)

  • LFREEZ_RATE: true to limit riming efficiency (heat budget on icy hydrometeors collecting supercooled liquid water)