Skip to content

Commit

Permalink
Add 3-moment option to NSSL microphysics (plus other minor scheme upd…
Browse files Browse the repository at this point in the history 
…ates) (#1876)

TYPE: enhancement

KEYWORDS: microphysics, NSSL scheme

SOURCE: Ted Mansell (NOAA/NSSL) (MicroTed)

DESCRIPTION OF CHANGES:

3-moment bulk microphysics provides more accurate size-sorting of fast-falling hydrometeors (rain, graupel, hail) by effectively predicting the shape parameter of the gamma function particle size distribution. Also various updates to increase run-time flexibility, for example, it is now possible to predict CCN concentration but have hail turned off. Some bug fixes and changes in default settings, as well. The previous mp_physics options (17,19,21,22) still work as expected, but can be replicated with mp_physics=18 plus new namelist options. A new README.NSSLmp gives an overview of the scheme and tips on usage and some details on an internal namelist for setting additional options within the module. The 5th-order WENO advection (module_advect_em.F) has a decreased 'epsilon' value (1.d-40) to accommodate the order of magnitude of reflectivity moments. This should not noticeably affect scalar advection but will cause round-off differences. Another change fixes an NSSL-MP bug where sedimentation did not work for k > 128 (i.e., if number of model levels exceeds 128, sedimentation only worked for k <= 128).

LIST OF MODIFIED FILES:

Registry/Registry.EM_COMMON
Registry/registry.var
Registry/wrfplus
chem/chemics_init.F
dyn_em/module_advect_em.F
dyn_em/solve_em.F
dyn_em/start_em.F
phys/module_microphysics_driver.F
phys/module_mp_nssl_2mom.F
phys/module_physics_init.F
phys/module_diagnostics_driver.F
phys/module_diag_nwp.F
wrftladj/module_microphysics_driver_ad.F
wrftladj/module_microphysics_driver_tl.F
run/README.namelist
share/module_check_a_mundo.F
doc/README.NSSLmp (new file)

TESTS CONDUCTED:
It passed the regression tests.

RELEASE NOTE: Adds 3-moment option to NSSL microphysics scheme (for rain, graupel, and hail) (mp_physics=18). Changes to snow aggregation parameters reduce stratiform reflectivity bias. Also adds namelist options to enable/disable bulk CCN prediction and graupel/hail density prediction, and deprecates old mp_physics options (17, 19, 21, 22). Main default parameter changes are for graupel/hail fall speed options (icdx, icdxhl; changed from 3 to 6), and default maximum gr/hail droplet collection efficiencies (ehw0/ehlw0 changed from 0.5/0.75 to 0.9/0.9) See README.NSSLmp for details. Also fixes a problem for high resolution vertical domains with nz > 128, where hydrometeor sedimentation did not work at k > 128 (this bug affected most or all previous releases).
  • Loading branch information
@MicroTed
MicroTed authored Jan 23, 2024
1 parent 5da1d24 commit e18082c
Show file tree
Hide file tree
Showing 17 changed files with 5,821 additions and 1,214 deletions.
70 changes: 51 additions & 19 deletions Registry/Registry.EM_COMMON
View file
Original file line number Diff line number Diff line change
Expand Up @@ -543,6 +543,12 @@ state real qvolg ikjftb scalar 1 - \
i0rhusdf=(bdy_interp:dt) "QVGRAUPEL" "Graupel Particle Volume" "m(3) kg(-1)"
state real qvolh ikjftb scalar 1 - \
i0rhusdf=(bdy_interp:dt) "QVHAIL" "Hail Particle Volume" "m(3) kg(-1)"
state real qzr ikjftb scalar 1 - \
i0rhusdf=(bdy_interp:dt) "QZRAIN" "Sixth moment rain" "m(6) kg(-1)"
state real qzg ikjftb scalar 1 - \
i0rhusdf=(bdy_interp:dt) "QZGRAUPEL" "Sixth moment graupel" "m(6) kg(-1)"
state real qzh ikjftb scalar 1 - \
i0rhusdf=(bdy_interp:dt) "QZHAIL" "Sixth moment hail" "m(6) kg(-1)"
state real qrimef ikjftb scalar 1 - \
i0rhusdf=(bdy_interp:dt) "QRIMEF" "rime factor * qi" "kg kg-1"
state real qir ikjftb scalar 1 - \
Expand Down Expand Up @@ -591,6 +597,12 @@ state real dfi_qnn ikjftb dfi_scalar 1 - \
rusdf=(bdy_interp:dt) "DFI_QNCC" "DFI CNN Number concentration" "# kg(-1)"
state real dfi_qnc ikjftb dfi_scalar 1 - \
rusdf=(bdy_interp:dt) "DFI_QNCLOUD" "DFI Cloud Number concentration" "# kg(-1)"
state real dfi_qzr ikjftb dfi_scalar 1 - \
rhusdf=(bdy_interp:dt) "DFI_QZRAIN" "DFI Rain Reflectivity" "m(6) kg(-1)"
state real dfi_qzg ikjftb dfi_scalar 1 - \
rhusdf=(bdy_interp:dt) "DFI_QZGRAUPEL" "DFI Graupel Reflectivity" "m(6) kg(-1)"
state real dfi_qzh ikjftb dfi_scalar 1 - \
rhusdf=(bdy_interp:dt) "DFI_QZHAIL" "DFI Hail Reflectivity" "m(6) kg(-1)"
state real dfi_qvolg ikjftb dfi_scalar 1 - \
rhusdf=(bdy_interp:dt) "DFI_QVGRAUPEL" "DFI Graupel Particle Volume" "m(3) kg(-1)"
state real dfi_qvolh ikjftb dfi_scalar 1 - \
Expand Down Expand Up @@ -2385,16 +2397,27 @@ rconfig logical write_thompson_tables namelist,physics 1 .tru
rconfig logical write_thompson_mp38table namelist,physics 1 .false.
rconfig integer mp_physics namelist,physics max_domains -1 irh "mp_physics" "" ""
#rconfig integer milbrandt_ccntype namelist,physics max_domains 0 rh "milbrandt select maritime(1)/continental(2)" "" ""
rconfig real nssl_cccn namelist,physics max_domains 0.5e9 rh "Base CCN concentration for NSSL microphysics" "" ""
rconfig real nssl_alphah namelist,physics max_domains 0 rh "Graupel PSD shape paramter" "" ""
rconfig real nssl_alphahl namelist,physics max_domains 1 rh "Hail PSD shape paramter" "" ""
rconfig real nssl_cnoh namelist,physics max_domains 4.e5 rh "Graupel intercept paramter" "" ""
rconfig real nssl_cnohl namelist,physics max_domains 4.e4 rh "Hail intercept paramter" "" ""
rconfig real nssl_cnor namelist,physics max_domains 8.e5 rh "Rain intercept paramter" "" ""
rconfig real nssl_cnos namelist,physics max_domains 3.e6 rh "Snow intercept paramter" "" ""
rconfig real nssl_rho_qh namelist,physics max_domains 500. rh "Graupel particle density" "" ""
rconfig real nssl_rho_qhl namelist,physics max_domains 900. rh "Hail particle density" "" ""
rconfig real nssl_rho_qs namelist,physics max_domains 100. rh "Snow particle density" "" ""
rconfig real nssl_cccn namelist,physics 1 0.5e9 rh "Base CCN concentration for NSSL microphysics" "" ""
rconfig real nssl_alphah namelist,physics 1 0 rh "Graupel PSD shape paramter" "" ""
rconfig real nssl_alphahl namelist,physics 1 1 rh "Hail PSD shape paramter" "" ""
rconfig real nssl_cnoh namelist,physics 1 4.e5 rh "Graupel intercept paramter" "" ""
rconfig real nssl_cnohl namelist,physics 1 4.e4 rh "Hail intercept paramter" "" ""
rconfig real nssl_cnor namelist,physics 1 8.e5 rh "Rain intercept paramter" "" ""
rconfig real nssl_cnos namelist,physics 1 3.e6 rh "Snow intercept paramter" "" ""
rconfig real nssl_rho_qh namelist,physics 1 500. rh "Graupel particle density" "" ""
rconfig real nssl_rho_qhl namelist,physics 1 900. rh "Hail particle density" "" ""
rconfig real nssl_rho_qs namelist,physics 1 100. rh "Snow particle density" "" ""
rconfig integer nssl_icdx namelist,physics 1 6 rh "NSSL Graupel fall speed option" "" ""
rconfig integer nssl_icdxhl namelist,physics 1 6 rh "NSSL Hail fall speed option" "" ""
rconfig integer nssl_hail_on namelist,physics max_domains -1 rh "NSSL Hail flag" "" ""
rconfig integer nssl_ccn_on namelist,physics 1 -1 rh "NSSL CCN flag" "" ""
rconfig integer nssl_ccn_is_ccna namelist,physics 1 0 rh "NSSL flag that CCN is CCNA" "" ""
rconfig integer nssl_2moment_on namelist,physics 1 -1 rh "NSSL 2-moment flag" "" ""
rconfig integer nssl_3moment namelist,physics 1 0 rh "NSSL 3-moment flag" "" ""
rconfig integer nssl_density_on namelist,physics 1 -1 rh "NSSL graupel/hail density flag" "" ""



rconfig integer CCNTY namelist,physics 1 2 rh "Aerosol background type for NTU microphysics" "" ""

# Lightning Qv Nudging
Expand Down Expand Up @@ -2995,11 +3018,8 @@ package cammgmpscheme mp_physics==11 - moist:qv,qc
package sbu_ylinscheme mp_physics==13 - moist:qv,qc,qr,qi,qs;state:rimi
package wdm5scheme mp_physics==14 - moist:qv,qc,qr,qi,qs;scalar:qnn,qnc,qnr;state:re_cloud,re_ice,re_snow
package wdm6scheme mp_physics==16 - moist:qv,qc,qr,qi,qs,qg;scalar:qnn,qnc,qnr;state:re_cloud,re_ice,re_snow
package nssl_2mom mp_physics==17 - moist:qv,qc,qr,qi,qs,qg,qh;scalar:qndrop,qnr,qni,qns,qng,qnh,qvolg,qvolh;state:re_cloud,re_ice,re_snow
package nssl_2momccn mp_physics==18 - moist:qv,qc,qr,qi,qs,qg,qh;scalar:qnn,qndrop,qnr,qni,qns,qng,qnh,qvolg,qvolh;state:re_cloud,re_ice,re_snow
package nssl_1mom mp_physics==19 - moist:qv,qc,qr,qi,qs,qg,qh;scalar:qvolg
package nssl_1momlfo mp_physics==21 - moist:qv,qc,qr,qi,qs,qg
package nssl_2momg mp_physics==22 - moist:qv,qc,qr,qi,qs,qg;scalar:qndrop,qnr,qni,qns,qng,qvolg;state:re_cloud,re_ice,re_snow
# Note: Options 17, 19, 21, 22 are deprecated but still reserved for compatibility -- for now
package nssl_2mom mp_physics==18 - moist:qv,qc,qr,qi,qs,qg
package wsm7scheme mp_physics==24 - moist:qv,qc,qr,qi,qs,qg,qh;state:re_cloud,re_ice,re_snow
package wdm7scheme mp_physics==26 - moist:qv,qc,qr,qi,qs,qg,qh;scalar:qnn,qnc,qnr;state:re_cloud,re_ice,re_snow
package thompsonaero mp_physics==28 - moist:qv,qc,qr,qi,qs,qg;scalar:qni,qnr,qnc,qnwfa,qnifa,qnbca;state:re_cloud,re_ice,re_snow,qnwfa2d,qnifa2d,taod5503d,taod5502d
Expand All @@ -3015,6 +3035,16 @@ package etampnew mp_physics==95 - moist:qv,qc
package gsfcgcescheme mp_physics==97 - moist:qv,qc,qr,qi,qs,qg
package madwrf_mp mp_physics==96 - moist:qv,qc,qi,qs

package nssl2mconc nssl_2moment_on==1 - scalar:qndrop,qnr,qni,qns,qng;state:re_cloud,re_ice,re_snow
package nssl3mg nssl_3moment==1 - scalar:qzr,qzg
package nssl3m nssl_3moment==2 - scalar:qzr,qzg,qzh
package nssl_hail nssl_hail_on==1 - moist:qh;scalar:qnh
package nssl_hail1m nssl_hail_on==2 - moist:qh;
package nssl_ccn_opt nssl_ccn_on==1 - scalar:qnn
package nssl_graupelvol nssl_density_on==1 - scalar:qvolg
package nssl_hailvol nssl_density_on==2 - scalar:qvolg,qvolh


package radar_refl compute_radar_ref==1 - state:refl_10cm,refd_max
endif

Expand All @@ -3038,10 +3068,12 @@ package morr_two_moment_dfi mp_physics_dfi==10 - dfi_moist:dfi
#package sbu_ylinscheme_dfi mp_physics==13 - dfi_moist:dfi_qv,dfi_qc,dfi_qr,dfi_qi,dfi_qs;state:rimi
package wdm5scheme_dfi mp_physics_dfi==14 - dfi_moist:dfi_qv,dfi_qc,dfi_qr,dfi_qi,dfi_qs;dfi_scalar:dfi_qnn,dfi_qnc,dfi_qnr;state:dfi_re_cloud,dfi_re_ice,dfi_re_snow
package wdm6scheme_dfi mp_physics_dfi==16 - dfi_moist:dfi_qv,dfi_qc,dfi_qr,dfi_qi,dfi_qs,dfi_qg;dfi_scalar:dfi_qnn,dfi_qnc,dfi_qnr;state:dfi_re_cloud,dfi_re_ice,dfi_re_snow
package nssl_2mom_dfi mp_physics_dfi==17 - dfi_moist:dfi_qv,dfi_qc,dfi_qr,dfi_qi,dfi_qs,dfi_qg,dfi_qh;dfi_scalar:dfi_qndrop,dfi_qnr,dfi_qni,dfi_qns,dfi_qng,dfi_qnh,dfi_qvolg,dfi_qvolh;state:dfi_re_cloud,dfi_re_ice,dfi_re_snow
package nssl_2mom_dficcn mp_physics_dfi==18 - dfi_moist:dfi_qv,dfi_qc,dfi_qr,dfi_qi,dfi_qs,dfi_qg,dfi_qh;dfi_scalar:dfi_qndrop,dfi_qnn,dfi_qnr,dfi_qni,dfi_qns,dfi_qng,dfi_qnh,dfi_qvolg;state:dfi_re_cloud,dfi_re_ice,dfi_re_snow
package nssl_1mom_dfi mp_physics_dfi==19 - dfi_moist:dfi_qv,dfi_qc,dfi_qr,dfi_qi,dfi_qs,dfi_qg,dfi_qh;dfi_scalar:dfi_qvolg
package nssl_1momlfo_dfi mp_physics_dfi==21 - dfi_moist:dfi_qv,dfi_qc,dfi_qr,dfi_qi,dfi_qs,dfi_qg
#package nssl_2mom_dfi mp_physics_dfi==17 - dfi_moist:dfi_qv,dfi_qc,dfi_qr,dfi_qi,dfi_qs,dfi_qg,dfi_qh;dfi_scalar:dfi_qndrop,dfi_qnr,dfi_qni,dfi_qns,dfi_qng,dfi_qnh,dfi_qvolg,dfi_qvolh;state:dfi_re_cloud,dfi_re_ice,dfi_re_snow
#package nssl_2mom_dficcn mp_physics_dfi==18 - dfi_moist:dfi_qv,dfi_qc,dfi_qr,dfi_qi,dfi_qs,dfi_qg,dfi_qh;dfi_scalar:dfi_qndrop,dfi_qnn,dfi_qnr,dfi_qni,dfi_qns,dfi_qng,dfi_qnh,dfi_qvolg;state:dfi_re_cloud,dfi_re_ice,dfi_re_snow
package nssl_2mom_dfi mp_physics_dfi==18 - dfi_moist:dfi_qv,dfi_qc,dfi_qr,dfi_qi,dfi_qs,dfi_qg,dfi_qh;dfi_scalar:dfi_qndrop,dfi_qnn,dfi_qnr,dfi_qni,dfi_qns,dfi_qng,dfi_qnh,dfi_qvolg;state:dfi_re_cloud,dfi_re_ice,dfi_re_snow
#package nssl_1mom_dfi mp_physics_dfi==19 - dfi_moist:dfi_qv,dfi_qc,dfi_qr,dfi_qi,dfi_qs,dfi_qg,dfi_qh;dfi_scalar:dfi_qvolg
#package nssl_1momlfo_dfi mp_physics_dfi==21 - dfi_moist:dfi_qv,dfi_qc,dfi_qr,dfi_qi,dfi_qs,dfi_qg
#package nssl_2momg_dfi mp_physics_dfi==22 - dfi_moist:dfi_qv,dfi_qc,dfi_qr,dfi_qi,dfi_qs,dfi_qg;dfi_scalar:dfi_qndrop,dfi_qnr,dfi_qni,dfi_qns,dfi_qng,dfi_qvolg;state:dfi_re_cloud,dfi_re_ice,dfi_re_snow
package wsm7scheme_dfi mp_physics_dfi==24 - dfi_moist:dfi_qv,dfi_qc,dfi_qr,dfi_qi,dfi_qs,dfi_qg,dfi_qh;state:dfi_re_cloud,dfi_re_ice,dfi_re_snow
package wdm7scheme_dfi mp_physics_dfi==26 - dfi_moist:dfi_qv,dfi_qc,dfi_qr,dfi_qi,dfi_qs,dfi_qg,dfi_qh;dfi_scalar:dfi_qnn,dfi_qnc,dfi_qnr;state:dfi_re_cloud,dfi_re_ice,dfi_re_snow
package thompsonaero_dfi mp_physics_dfi==28 - dfi_moist:dfi_qv,dfi_qc,dfi_qr,dfi_qi,dfi_qs,dfi_qg;dfi_scalar:dfi_qni,dfi_qnr,dfi_qnc,dfi_qnwfa,dfi_qnifa,dfi_qnbca;state:dfi_re_cloud,dfi_re_ice,dfi_re_snow
Expand Down
15 changes: 5 additions & 10 deletions Registry/registry.var
View file
Original file line number Diff line number Diff line change
Expand Up @@ -593,11 +593,8 @@ package cammgmpscheme mp_physics==11 - moist:qv,qc
package sbu_ylinscheme mp_physics==13 - moist:qv,qc,qr,qi,qs
package wdm5scheme mp_physics==14 - moist:qv,qc,qr,qi,qs
package wdm6scheme mp_physics==16 - moist:qv,qc,qr,qi,qs,qg
package nssl_2mom mp_physics==17 - moist:qv,qc,qr,qi,qs,qg,qh
package nssl_2momccn mp_physics==18 - moist:qv,qc,qr,qi,qs,qg,qh
package nssl_1mom mp_physics==19 - moist:qv,qc,qr,qi,qs,qg,qh
package nssl_1momlfo mp_physics==21 - moist:qv,qc,qr,qi,qs,qg
package nssl_2momg mp_physics==22 - moist:qv,qc,qr,qi,qs,qg
# Note: Options 17, 19, 21, 22 are deprecated but still reserved for compatibility -- for now
package nssl_2mom mp_physics==18 - moist:qv,qc,qr,qi,qs,qg
package thompsonaero mp_physics==28 - moist:qv,qc,qr,qi,qs,qg
package p3_1category mp_physics==50 - moist:qv,qc,qr,qi
package p3_1category_nc mp_physics==51 - moist:qv,qc,qr,qi
Expand All @@ -607,6 +604,7 @@ package ntu mp_physics==56 - moist:qv,qc
package etampnew mp_physics==95 - moist:qv,qc,qr,qs
package lscondscheme mp_physics==98 - moist:qv
package mkesslerscheme mp_physics==99 - moist:qv,qc,qr

#
package mpnotset_4dvar mp_physics_4dvar==-1 - g_moist:g_qv;a_moist:a_qv
package passiveqv_4dvar mp_physics_4dvar==0 - g_moist:g_qv;a_moist:a_qv
Expand All @@ -626,11 +624,8 @@ package cammgmp_4dvar mp_physics_4dvar==11 - g_moist:g_q
package sbu_ylin_4dvar mp_physics_4dvar==13 - g_moist:g_qv,g_qc,g_qr,g_qi,g_qs;a_moist:a_qv,a_qc,a_qr,a_qi,a_qs
package wdm5_4dvar mp_physics_4dvar==14 - g_moist:g_qv,g_qc,g_qr,g_qi,g_qs;a_moist:a_qv,a_qc,a_qr,a_qi,a_qs
package wdm6_4dvar mp_physics_4dvar==16 - g_moist:g_qv,g_qc,g_qr,g_qi,g_qs,g_qg;a_moist:a_qv,a_qc,a_qr,a_qi,a_qs,a_qg
package nssl_2mom_4dvar mp_physics_4dvar==17 - g_moist:g_qv,g_qc,g_qr,g_qi,g_qs,g_qg,g_qh;a_moist:a_qv,a_qc,a_qr,a_qi,a_qs,a_qg,a_qh
package nssl_2momccn_4dvar mp_physics_4dvar==18 - g_moist:g_qv,g_qc,g_qr,g_qi,g_qs,g_qg,g_qh;a_moist:a_qv,a_qc,a_qr,a_qi,a_qs,a_qg,a_qh
package nssl_1mom_4dvar mp_physics_4dvar==19 - g_moist:g_qv,g_qc,g_qr,g_qi,g_qs,g_qg,g_qh;a_moist:a_qv,a_qc,a_qr,a_qi,a_qs,a_qg,a_qh
package nssl_1momlfo_4dvar mp_physics_4dvar==21 - g_moist:g_qv,g_qc,g_qr,g_qi,g_qs,g_qg;a_moist:a_qv,a_qc,a_qr,a_qi,a_qs,a_qg
package nssl_2momg_4dvar mp_physics_4dvar==22 - g_moist:g_qv,g_qc,g_qr,g_qi,g_qs,g_qg;a_moist:a_qv,a_qc,a_qr,a_qi,a_qs,a_qg
# Note: Options 17, 19, 21, 22 are deprecated but still reserved for compatibility -- for now
package nssl_2mom_4dvar mp_physics_4dvar==18 - g_moist:g_qv,g_qc,g_qr,g_qi,g_qs,g_qg,g_qh;a_moist:a_qv,a_qc,a_qr,a_qi,a_qs,a_qg,a_qh
package thompsonaero_4dvar mp_physics_4dvar==28 - g_moist:g_qv,g_qc,g_qr,g_qi,g_qs,g_qg;a_moist:a_qv,a_qc,a_qr,a_qi,a_qs,a_qg
package p3_1category_4dvar mp_physics_4dvar==50 - g_moist:g_qv,g_qc,g_qr,g_qi;a_moist:a_qv,a_qc,a_qr,a_qi
package p3_1category_nc_4dvar mp_physics_4dvar==51 - g_moist:g_qv,g_qc,g_qr,g_qi;a_moist:a_qv,a_qc,a_qr,a_qi
Expand Down
6 changes: 1 addition & 5 deletions Registry/registry.wrfplus
View file
Original file line number Diff line number Diff line change
Expand Up @@ -872,11 +872,7 @@ package cammgmp_plus mp_physics_plus==11 - g_moist:g_qv,
package sbu_ylin_plus mp_physics_plus==13 - g_moist:g_qv,g_qc,g_qr,g_qi,g_qs;a_moist:a_qv,a_qc,a_qr,a_qi,a_qs
package wdm5_plus mp_physics_plus==14 - g_moist:g_qv,g_qc,g_qr,g_qi,g_qs;a_moist:a_qv,a_qc,a_qr,a_qi,a_qs
package wdm6_plus mp_physics_plus==16 - g_moist:g_qv,g_qc,g_qr,g_qi,g_qs,g_qg;a_moist:a_qv,a_qc,a_qr,a_qi,a_qs,a_qg
package nssl_2mom_plus mp_physics_plus==17 - g_moist:g_qv,g_qc,g_qr,g_qi,g_qs,g_qg,g_qh;a_moist:a_qv,a_qc,a_qr,a_qi,a_qs,a_qg,a_qh
package nssl_2momccn_plus mp_physics_plus==18 - g_moist:g_qv,g_qc,g_qr,g_qi,g_qs,g_qg,g_qh;a_moist:a_qv,a_qc,a_qr,a_qi,a_qs,a_qg,a_qh
package nssl_1mom_plus mp_physics_plus==19 - g_moist:g_qv,g_qc,g_qr,g_qi,g_qs,g_qg,g_qh;a_moist:a_qv,a_qc,a_qr,a_qi,a_qs,a_qg,a_qh
package nssl_1momlfo_plus mp_physics_plus==21 - g_moist:g_qv,g_qc,g_qr,g_qi,g_qs,g_qg;a_moist:a_qv,a_qc,a_qr,a_qi,a_qs,a_qg
package nssl_2momg_plus mp_physics_plus==22 - g_moist:g_qv,g_qc,g_qr,g_qi,g_qs,g_qg;a_moist:a_qv,a_qc,a_qr,a_qi,a_qs,a_qg
package nssl_2mom_plus mp_physics_plus==18 - g_moist:g_qv,g_qc,g_qr,g_qi,g_qs,g_qg,g_qh;a_moist:a_qv,a_qc,a_qr,a_qi,a_qs,a_qg,a_qh
package thompsonaero_plus mp_physics_plus==28 - g_moist:g_qv,g_qc,g_qr,g_qi,g_qs,g_qg;a_moist:a_qv,a_qc,a_qr,a_qi,a_qs,a_qg
package p3_1category_plus mp_physics_plus==50 - g_moist:g_qv,g_qc,g_qr,g_qi;a_moist:a_qv,a_qc,a_qr,a_qi
package p3_1category_nc_plus mp_physics_plus==51 - g_moist:g_qv,g_qc,g_qr,g_qi;a_moist:a_qv,a_qc,a_qr,a_qi
Expand Down
4 changes: 2 additions & 2 deletions chem/chemics_init.F
View file
Original file line number Diff line number Diff line change
Expand Up @@ -337,8 +337,8 @@ subroutine chem_init (id,chem,emis_ant,scalar,dt,bioemdt,photdt,chemdt,stepbioe,
call wrf_error_fatal("ERROR: wet scavenging option requires chem_opt = 8 through 13 or 31 to 36 or 41 to 42 or 109 or 503 or 504 or 601 or 611 to function.")
endif
if ( config_flags%mp_physics /= 2 .and. config_flags%mp_physics /= 10 .and. config_flags%mp_physics /= 11 &
.and. config_flags%mp_physics /= 17 .and. config_flags%mp_physics /= 18 .and. config_flags%mp_physics /= 22) then
call wrf_error_fatal("ERROR: wet scavenging option requires mp_phys = 2 (Lin et al.) or 10 (Morrison) or 11 (CAMMGMP) or 17/18/22 NSSL_2mom to function.")
.and. .not. ( config_flags%mp_physics == 18 .and. config_flags%nssl_2moment_on == 1 ) ) then
call wrf_error_fatal("ERROR: wet scavenging option requires mp_phys = 2 (Lin et al.) or 10 (Morrison) or 11 (CAMMGMP) or 18 NSSL_2mom to function.")
endif
elseif( id == 1 ) then
if ( config_flags%mp_physics /= 6 .and. config_flags%mp_physics /= 8 .and. config_flags%mp_physics /= 10 .and. config_flags%mp_physics /= 17 &
Expand Down
Loading

0 comments on commit e18082c

Please sign in to comment.