CCPP SciDoc v7.0.0  v7.0.0
Common Community Physics Package Developed at DTC
 
Loading...
Searching...
No Matches
micro_mg3_0 Module Reference

MG microphysics version 3.0 - Update of MG microphysics with prognostic hail OR graupel. More...

Functions/Subroutines

subroutine, public micro_mg_init (kind, gravit, rair, rh2o, cpair, eps, tmelt_in, latvap, latice, rhmini_in, micro_mg_dcs, ts_auto, mg_qcvar, micro_mg_do_hail_in, micro_mg_do_graupel_in, microp_uniform_in, do_cldice_in, use_hetfrz_classnuc_in, micro_mg_precip_frac_method_in, micro_mg_berg_eff_factor_in, allow_sed_supersat_in, do_sb_physics_in, do_ice_gmao_in, do_liq_liu_in, nccons_in, nicons_in, ncnst_in, ninst_in, ngcons_in, ngnst_in)
 This subroutine initializes the microphysics and needs to be called once at start of simulation.
 
subroutine calc_rercld (lamr, n0r, lamc, pgam, qric, qcic, ncic, rercld, mgncol, nlev)
 This subroutine calculates effective radii for rain and cloud.
 
subroutine, public micro_mg_tend (mgncol, nlev, deltatin, t, q, qcn, qin, ncn, nin, qrn, qsn, nrn, nsn, qgr, ngr, relvar, accre_enhan_i, p, pdel, cldn, liqcldf, icecldf, qsatfac, qcsinksum_rate1ord, naai, npccnin, rndst, nacon, tlat, qvlat, qctend, qitend, nctend, nitend, qrtend, qstend, nrtend, nstend, qgtend, ngtend, effc, effc_fn, effi, sadice, sadsnow, prect, preci, nevapr, evapsnow, am_evp_st, prain, prodsnow, cmeout, deffi, pgamrad, lamcrad, qsout, dsout, qgout, ngout, dgout, lflx, iflx, gflx, rflx, sflx, qrout, reff_rain, reff_snow, reff_grau, qcsevap, qisevap, qvres, cmeitot, vtrmc, vtrmi, umr, ums, umg, qgsedten, qcsedten, qisedten, qrsedten, qssedten, pratot, prctot, mnuccctot, mnuccttot, msacwitot, psacwstot, bergstot, bergtot, melttot, homotot, qcrestot, prcitot, praitot, qirestot, mnuccrtot, mnuccritot, pracstot, meltsdttot, frzrdttot, mnuccdtot, pracgtot, psacwgtot, pgsacwtot, pgracstot, prdgtot, qmultgtot, qmultrgtot, psacrtot, npracgtot, nscngtot, ngracstot, nmultgtot, nmultrgtot, npsacwgtot, nrout, nsout, refl, arefl, areflz, frefl, csrfl, acsrfl, fcsrfl, rercld, ncai, ncal, qrout2, qsout2, nrout2, nsout2, drout2, dsout2, qgout2, ngout2, dgout2, freqg, freqs, freqr, nfice, qcrat, prer_evap, xlat, xlon, lprnt, iccn, nlball)
 This subroutine calculates the MG3 microphysical processes.
 

Variables

logical nccons
 nccons = .true. to specify constant cloud droplet number
 
logical nicons
 nicons = .true. to specify constant cloud ice number
 
logical ngcons
 ngcons = .true. to specify constant graupel number
 
real(r8) ncnst
 droplet num concentration when nccons=.true. (m-3)
 
real(r8) ninst
 ice num concentration when nicons=.true. (m-3)
 
real(r8) ngnst
 graupel num concentration when ngcons=.true. (m-3)
 
real(r8), parameter csmin = -30._r8
 Range of cloudsat reflectivities (dBz) for analytic simulator.
 
real(r8), parameter csmax = 26._r8
 
real(r8), parameter mindbz = -99._r8
 
real(r8), parameter minrefl = 1.26e-10_r8
 
real(r8) dcs
 
real(r8) ts_au
 
real(r8) ts_au_min
 
real(r8), public qcvar
 
real(r8), parameter mi0l_min = 4._r8/3._r8*pi*rhow*(4.e-6_r8)**3
 
real(r8), parameter sublim_factor = 0.0_r8
 
real(r8), parameter zero =0.0_r8
 
real(r8), parameter one =1.0_r8
 
real(r8), parameter two =2.0_r8
 
real(r8), parameter three =3.0_r8
 
real(r8), parameter four =4.0_r8
 
real(r8), parameter five =5.0_r8
 
real(r8), parameter six =6._r8
 
real(r8), parameter half =0.5_r8
 
real(r8), parameter ten =10.0_r8
 
real(r8), parameter forty =40.0_r8
 
real(r8), parameter oneo6 =one/six
 
real(r8) g
 gravity
 
real(r8) r
 dry air gas constant
 
real(r8) rv
 water vapor gas constant
 
real(r8) cpp
 specific heat of dry air
 
real(r8) tmelt
 freezing point of water (K)
 
real(r8) xxlv
 vaporization
 
real(r8) xlf
 freezing
 
real(r8) xxls
 sublimation
 
real(r8) rhmini
 Minimum rh for ice cloud fraction > 0.
 
logical microp_uniform
 
logical do_cldice
 
logical use_hetfrz_classnuc
 
logical do_hail
 
logical do_graupel
 
real(r8) rhosu
 typical 850mn air density
 
real(r8) icenuct
 ice nucleation temperature: currently -5 degrees C
 
real(r8) snowmelt
 what temp to melt all snow: currently 2 degrees C
 
real(r8) rainfrze
 what temp to freeze all rain: currently -5 degrees C
 
real(r8) rhogtmp
 hail or graupel density (kg m-3)
 
real(r8) agtmp
 tmp ag/ah parameter
 
real(r8) bgtmp
 tmp fall speed parameter
 
real(r8) gamma_br_plus1
 
real(r8) gamma_bs_plus1
 
real(r8) gamma_bi_plus1
 
real(r8) gamma_bj_plus1
 
real(r8) gamma_bg_plus1
 
real(r8) gamma_br_plus4
 
real(r8) gamma_bs_plus4
 
real(r8) gamma_bi_plus4
 
real(r8) gamma_bj_plus4
 
real(r8) gamma_bg_plus4
 
real(r8) xxlv_squared
 
real(r8) xxls_squared
 
real(r8) omeps
 
real(r8) epsqs
 
character(len=16) micro_mg_precip_frac_method
 type of precipitation fraction method
 
real(r8) micro_mg_berg_eff_factor
 berg efficiency factor
 
logical allow_sed_supersat
 Allow supersaturated conditions after sedimentation loop.
 
logical do_sb_physics
 do SB 2001 autoconversion or accretion physics
 
logical do_ice_gmao
 
logical do_liq_liu
 

Detailed Description

Authors
Andrew Gettelman, Hugh Morrison
Version
3 history: Sep 2016: development begun for hail, graupel This version:https://svn-ccsm-models.cgd.ucar.edu/cam1/branch_tags/mg3_tags/mg3_33_cam5_4_153/
2 history: Sep 2011: Development begun.
Feb 2013: Added of prognostic precipitation.
Aug 2015: Published and released version

Contributions from: Sean Santos, Peter Caldwell, Xiaohong Liu and Steve Ghan

  • Anning Cheng adopted mg2 for FV3GFS 9/29/2017
    add GMAO ice conversion and Liu et. al liquid water
    conversion in 10/12/2017
  • Anning showed promising results for FV3GFS on 10/15/2017
  • S. Moorthi - Oct/Nov 2017 - optimized the MG2 code
  • S. Moorthi - Nov 2017 - made the sedimentation quasi-implicit
  • S. Moorthi - Feb 2018 - updated to MG3 - modified graupel sedimentation other modifications to eliminate blowup.
  • S. Moorthi - Mar 2018 - fixed a few bugs and added option to run as MG2
  • S. Moorthi - Oct,29,2018 - change nlb from nlev/3 to levels with p/ps < 0.05 (nlball)

invoked in CAM by specifying -microphys=mg3

References:

      Gettelman, A. and H. Morrison, Advanced Two-Moment Microphysics for Global Models.
      Part I: Off line tests and comparisons with other schemes.
      J. Climate, 28, 1268-1287. doi: 10.1175/JCLI-D-14-00102.1, 2015.

      Gettelman, A., H. Morrison, S. Santos, P. Bogenschutz and P. H. Caldwell
      Advanced Two-Moment Microphysics for Global Models.
      Part II: Global model solutions and Aerosol-Cloud Interactions.
      J. Climate, 28, 1288-1307. doi:10.1175/JCLI-D-14-00103.1 , 2015.

for questions contact Hugh Morrison, Andrew Gettelman

e-mail: morri.nosp@m.son@.nosp@m.ucar..nosp@m.edu, andre.nosp@m.w@uc.nosp@m.ar.ed.nosp@m.u

NOTE: Modified to allow other microphysics packages (e.g. CARMA) to do ice microphysics in cooperation with the MG liquid microphysics. This is controlled by the do_cldice variable.

If do_cldice is false, then MG microphysics should not update CLDICE or NUMICE; it is assumed that the other microphysics scheme will have updated CLDICE and NUMICE. The other microphysics should handle the following processes that would have been done by MG:

  • Detrainment (liquid and ice)
  • Homogeneous ice nucleation
  • Heterogeneous ice nucleation
  • Bergeron process
  • Melting of ice
  • Freezing of cloud drops
  • Autoconversion (ice -> snow)
  • Growth/Sublimation of ice
  • Sedimentation of ice

This option has not been updated since the introduction of prognostic precipitation, and probably should be adjusted to cover snow as well.