! $Id: strat_nml.h,v 20.0 2013/12/13 23:22:19 fms Exp $ ! $Name: tikal $ !------------------------------------------------------------------------ !---namelist------ !------------------------------------------------------------------------ ! the namelist variables are stored in a strat_nml_type derived type ! variable so that they may be conveniently passed to other modules. !------------------------------------------------------------------------ ! ! ! a netcdf restart file should be written ? ! ! ! threshold rel hum for cloud formation by large-scale condensation. ! ! ! the low-level u00 ECMWF profile should be applied? ! ! ! liquid cloud drop radius threshold for autoconversion. ! ! ! the Khairoutdinov and Kogan (2000) autoconversion should be used? ! ! ! critical relative humidity above which rain does not evaporate. ! ! ! normal erosion rate constant for cloud destruction (default = 1.E-06) ! ! ! enhanced erosion should occur in turbulent or convective conditions ? ! ! ! erosion rate constant for cloud destruction for convective conditions. ! ! ! erosion rate constant for cloud destruction for turbulent conditions. ! ! ! convective mass-flux threshold for enhanced erosion to turn on. ! ! ! diffusion coefficient threshold for enhanced erosion to turn on. ! ! ! excess vapor in supersaturated conditions should be put into ! cloud water (true) or precipitation fluxes (false) ? ! ! ! cloud liquid, ice and fraction are advected by the grid resolved motion? ! ! ! minimum permissible value of cloud liquid, cloud ice, saturated volume ! fraction, or rain and snow areas. ! NOTE: qmin should be chosen such that the range of {qmin, max(qa,ql,qi)} ! is resolved by the precision of the numbers used. ! ! ! minimum permissible dissipation in analytic integration of qa, ql, qi ! equations. This constant only affects the method by which the ! prognostic equations are integrated. ! NOTE: Dmin will be MACHINE DEPENDENT and occur when ! a) 1. -exp(-Dmin) = 0. instead of Dmin in the limit of very small Dmin. ! AND ! b) 1. - exp(-D) < D for all D > Dmin ! ! ! factor used in formula to enhance cloud erosion as a function of height ! ! ! factor used to enhance ice fall velocity. ! ! ! factor in bergeron process calculation when the effect of dust ! is not considered; values > 1 enhance the conversion to ice. ! ! ! the old version of snow melting, which has a bug, should be run? ! ! ! critical ice-water content below which to apply alternate fall ! speed formula ! ! ! factor for alternate fall speed formula ! ! ! exponent for alternate fall speed formula ! ! ! ! ! ! ! ! ! the microphysics scheme being used (currently either ! "morrison_gettelman" or "rotstayn_klein") ! ! ! the macrophysics scheme being used (currently either ! "tiedtke" or " ") ! ! ! the aerosol activation scheme being used (currently either ! "dqa" or "total") ! ! ! should we use ensure water mass conservation by adjusting precip to ! balance column water mass change ? ! ! ! should all potentially available ice particles be activated under all ! conditions ? (or only when dqa is increasing) ! ! ! the hallet-mossop process should be included in the NCAR microphysics? ! ! ! the minor bug present in CM3, in which several small terms in qv and ! temp equations were retained while corresponding terms in ql and qi ! were not, is retained? ! ! ! flag to indicate how to treat supersaturation; 0 => don't allow. ! ! ! should we use activate ice nuclei by assuming a pdf ? ! ! ! the online aerosol fields should be used for nucleation source ! (rather than climo values) ? ! ! ! only sub-micron seasalt particles should be used for nucleation ? ! ! ! scaling factor used when offline seasalt aerosol is used for ! nucleation ! ! ! scaling factor used to convert offline organic matter to organic ! carbon for nucleation ! ! ! assumed number of cloud drops per unit volume in liquid clouds ! over land when droplet number is not prdicted. ! ! ! assumed number of cloud drops per unit volume in liquid clouds ! over ocean when droplet number is not predicted. ! ! ! minimum value of the variance in the vertical velocity pdf ! ! ! the prognostic droplet number option is activated ? ! ! ! sub-micron dust particles are used as ice nuclei for bergeron ! process evaluation? ! ! ! minimum number of droplets allowed in a grid box when predicted ! droplet number code is activated ! ! ! the statistical cloud scheme should be run? ! ! ! p-parameter to the beta distribution - used when do_pdf_clouds is true ! ! ! half-width to the qt PDF - used when do_pdf_clouds is true and ! diagnostic variance ! The fraction of qtbar (mean total water in the grid box) that the ! maximum and minimum of the distribution differ from qtbar. That is, ! total water at the sub-grid scale may take on values anywhere between ! (1.-qthalfwidth)*qtbar and (1.+qthalfwidth)*qtbar ! ! ! number of sublevels to vertical sub-grid cloud structure - used when ! do_pdf_cloud is true ! ! ! PPM partial remap integer - used when do_pdf_cloud is true and if ! vertical subgrid structure is used ! ! ! PPM method number - used when do_pdf_cloud is true and if vertical ! subgrid structure is used ! ! ! define pdf based on input water fields (rather than updated ones) ? ! ! ! number of grid points where instantaneous output will be saved to ! file strat.data. ! num_strat_pts must be <= max_strat_pts. ! NOTE: this option is currently not supported. ! ! ! "num_strat_pts" pairs of grid indices, i.e., the global indices for i,j. ! NOTE: this option is currently not supported. ! ! ! should we activated debug diagnostics ? ! ! ! i coordinate of grid point for which debug data is to be written ! ! ! j coordinate of grid point for which debug data is to be written ! ! ! k coordinate of grid point for which debug data is to be written ! ! INTEGER, PARAMETER :: max_strat_pts = 5 logical :: do_netcdf_restart = .true. real :: U00 = 0.80 logical :: u00_profile = .false. real :: rthresh = 10. logical :: use_kk_auto = .false. real :: U_evap = 1.0 real :: eros_scale = 1.E-06 logical :: eros_choice = .false. real :: eros_scale_c = 8.E-06 real :: eros_scale_t = 5.E-05 real :: mc_thresh = 0.001 real :: diff_thresh = 1.0 logical :: super_choice = .false. logical :: tracer_advec = .false. real :: qmin = 1.E-10 real :: Dmin = 1.E-08 real :: efact = 0.0 real :: vfact = 1.0 real :: cfact = 1.0 logical :: do_old_snowmelt= .false. real :: iwc_crit = 0. real :: vfall_const2 = 3.29 real :: vfall_exp2 = 0.16 real :: num_mass_ratio1= 1. real :: num_mass_ratio2= 1. character(len=64) :: microphys_scheme = 'rotstayn_klein' character(len=64) :: macrophys_scheme = 'tiedtke' character(len=64) :: aerosol_activation_scheme = 'dqa' logical :: mass_cons = .true. logical :: activate_all_ice_always= .true. logical :: do_hallet_mossop = .false. logical :: retain_cm3_bug = .false. integer :: super_ice_opt = 0 logical :: do_ice_nucl_wpdf = .false. logical :: use_online_aerosol = .false. logical :: use_sub_seasalt = .true. real :: sea_salt_scale = 0.1 real :: om_to_oc = 1.67 real :: N_land = 250.E+06 real :: N_ocean = 100.E+06 real :: var_limit = 0.0 logical :: do_liq_num = .false. logical :: do_dust_berg = .false. real :: N_min = 1.E6 logical :: do_pdf_clouds = .false. integer :: betaP = 5 real :: qthalfwidth = 0.1 integer :: nsublevels = 1 integer :: kmap = 1 integer :: kord = 7 logical :: pdf_org = .true. integer :: num_strat_pts = 0 integer,dimension(2,max_strat_pts) :: strat_pts = 0 logical :: debugo = .false. integer :: isamp = 1 integer :: jsamp = 1 integer :: ksamp = 1 ! 1 / relative variance of sub-grid cloud water distribution ! see morrison and gettelman, 2007, J. Climate for details real :: qcvar = 2. namelist / strat_cloud_nml / & do_netcdf_restart, U00, u00_profile, rthresh, use_kk_auto, & U_evap, eros_scale, eros_choice, eros_scale_c, eros_scale_t, & mc_thresh, diff_thresh, super_choice, tracer_advec, qmin, Dmin, & efact, vfact, cfact, do_old_snowmelt, iwc_crit, vfall_const2, & vfall_exp2, num_mass_ratio1, num_mass_ratio2, & microphys_scheme, macrophys_scheme, aerosol_activation_scheme, & mass_cons, activate_all_ice_always, do_hallet_mossop, & retain_cm3_bug, super_ice_opt, do_ice_nucl_wpdf, & use_online_aerosol, use_sub_seasalt, sea_salt_scale, om_to_oc, & N_land, N_ocean, var_limit, do_liq_num, do_dust_berg, N_min, & do_pdf_clouds, betaP, qthalfwidth, nsublevels, kmap, kord, pdf_org, & num_strat_pts, strat_pts, debugo, isamp, jsamp, ksamp, & qcvar