&version ver="vx.x.x" / ! ---------------------------------------------------------------------------------- &model level = 4 ! Thermodynamical level nxp = 68 ! Number of points in x direction nyp = 68 ! Number of points in y direction nzp = 120 ! Number of vertical levels deltax = 30. ! Grid spacing in x deltay = 30. ! Grid spacing in y deltaz = 10. ! Grid spacing in the vertical nxpart = .TRUE. ! dzmax = 1500. ! Height above which start stretching vertical grid dzrat = 1.0 ! Factor for vertical grid stretching dtlong = 1. ! Max. timestep distim = 100. ! Timescale for the dissipation in sponge layer timmax = 30600. ! Length of the simulation runtype = "INITIAL" ! INITIAL or HISTORY (restart) run CCN = 100.e6 ! corflg = .false. ! Apply coriolis force prndtl = -0.333333 ! filprf = 'cloud_20201031_forces_L4b088OC_23' ! Output filename profile hfilin = 'UAE_Cu_restart_6.014400s' ! History file name frqhis = 36000. !mcflg = .FALSE. ! Do mass conservation statistics sed_aero%switch = .FALSE. ! Calculate sedimentation of aerosol particles sed_cloud%switch = .TRUE. ! - '' - cloud droplets sed_cloud%delay = 3600. sed_precp%switch = .TRUE. ! precipitation sed_precp%delay = 3600. sed_ice%switch = .TRUE. ! ice particles bulk_autoc%switch = .FALSE. bulk_autoc%delay = 3600. ! Autoconversion switch for level = 1-3 lnudging = .FALSE. ! Master swtich for nudging scheme lemission = .FALSE. ! Master switch for aerosol emissions lpback = .FALSE. iradtyp = 3 ! Radiation/large scale forcing strtim = 303.9410 ! Start time case_name = 'default' ! Case name for large-scale forcing schemes div = 6.0e-6 ! Large-scale divergence cntlat = 62.53 / ! ----------------------------------------------------------------------------------- &output breakUndefOutput = .TRUE. ps_intvl = 60. ts_intvl = 60. main_intvl = 60. varlist_main = 'temp','theta','wwind','rp','rc','rrate','Naba','Dwaba','Ncba','Dwcba','Npba','Dwpba','lwp','rh','Np','Dwpa','Naa','Dwaa','Nca','Dwca','press','rflx','sflx','dn','s_n_activ','srp','aSO4a','aOCa','cSO4a','cOCa','pSO4a','pOCa' !varlist_main = 'temp','theta','wwind','rp','rc','rrate','lwp','rh','rflx','dn', 'press', !varlist_ps = '' !varlist_ts = '' TH_rrate = 30. TH_rc = 1.e-5 TH_rr = 1.e-6 / ! ----------------------------------------------------------------------------------- &initialization init_type = 1 ! 1: random perturbations, 2: warm bubble bubble_center = 1000.,0.,0. bubble_diameter = 1500.,6000.,6000. bubble_temp_ampl = 1.5 iseed = 3 th00 = 275. ! Reference temperature umean = 0. vmean = 0. itsflg = 1 ! Flag for temperature type in input sounding ipsflg = 1 ! - '' - for pressure / ! ------------------------------------------------------------------------------------- &surface isfctyp = 0 ! This is not entirely accurate for momentum vs vanZanten et al. 2011 zrough = 0.0002 ubmin = -0.25 sst = 272.85 dthcon = 0. !0.001094 ! with isfctyp these are the bulk formula coefficients drtcon = 0. !0.00133 !C_heat = 2.e6 !deepSoilTemp = 299.8 !W1 = 0.1 !W2 = 0.1 !W3 = 0.1 lConstSoilWater = .TRUE. lConstSoilHeatCap = .TRUE. / ! ------------------------------------------------------------------------------------ &radiation radsounding = 'datafiles/kmls.lay' sfc_albedo = 0.1 laerorad = .FALSE. useMcICA = .TRUE. / ! ------------------------------------------------------------------------------------- ! With lnudging = .TRUE. &nudge nudge_time = 3600. ! Overall time for nudging from the start of the simulation ndg_theta%nudgetype = 1 ndg_theta%tau_type = 2 ! Type of relaxation time (0:constant, 1-3: increasing) ndg_theta%tau_min = 60. ! Min relaxation time (with tau_type=1-3 and constant tau) ndg_theta%tau_max = 3600. ! Max relaxation time (with tau_type=1-3) ndg_theta%tau_max_continue = .FALSE. / ! --------------------------------------------------------------------------------------- ! With lemission = .TRUE. and level >= 4 &emission nEmissionModes = 2 !emitModes(1)%emitType = 3 !emitModes(1)%regime = 2 !emitModes(1)%start_time = 200. !emitModes(1)%end_time = 7201. !emitModes(1)%species = "SS" !emitModes(1)%emitHeightMin = 600. !emitModes(1)%emitHeightMax = 650. !emitModes(1)%emitSizeDistType = 2 !emitModes(1)%emitDiam = 0.6e-6 !emitModes(1)%emitNum = 10.5e7 !emitModes(1)%emitSigma = 1.1 !emitModes(1)%emitMap = 'emitMap_1' !emitModes(1)%scS = 60. !emitModes(1)%z_expan_up = 3 !emitModes(1)%z_expan_dw = 2 emitModes(1)%emitType = 2 emitModes(1)%regime = 2 emitModes(1)%start_time = 15000. emitModes(1)%end_time = 15001. emitModes(1)%species = "SO4" emitModes(1)%emitHeightMin = 5450. emitModes(1)%emitHeightMax = 5550. emitModes(1)%emitSizeDistType = 2 emitModes(1)%emitDiam = 1.5e-6 emitModes(1)%emitNum = 5.e6 emitModes(1)%emitSigma = 1.3 emitModes(2)%emitType = 2 emitModes(2)%regime = 2 emitModes(2)%start_time = 15000. emitModes(2)%end_time = 15001. emitModes(2)%species = "SO4" emitModes(2)%emitHeightMin = 5450. emitModes(2)%emitHeightMax = 5550. emitModes(2)%emitSizeDistType = 2 emitModes(2)%emitDiam = 8.0e-6 emitModes(2)%emitNum = 5.e6 emitModes(2)%emitSigma = 1.6 / ! ----------------------------------------------------------------------------------------- ! With level >= 4 &salsa ! Master process switches lscoag%switch = .TRUE. lscoag%delay = 3600. lscoag%mode = 2 ! 1: update kernels every timestep, 2: reduced update freq cgintvl = 10. ! for mode = 2, update interval in seconds lscnd%switch = .TRUE. lsauto%switch = .TRUE. lsauto%delay = 3600. lsauto%mode = 1 ! 1: precip formation based on coagulation, 2: parameterized autoconversion lsactiv%switch = .TRUE. lsactiv%mode = 1 ! 1: aerosol growth based activation, 2: parameterized cloud base activation lsicenucl%switch = .TRUE. lsicenucl%delay = 100. lsicemelt%switch = .TRUE. ! Ice melting lssecice%switch = .TRUE. ! Secondary ice production lssecice%delay = 200. ! Sub-process switches: Coagulation lscgcc = .TRUE. ! Self-collection of cloud droplets lscgpp = .TRUE. ! Self-collection of rain drops lscgpc = .TRUE. ! Rain collection of cloud droplets lscgaa = .FALSE. ! Aerosol coagulation lscgca = .TRUE. ! Cloud collection of aerosols lscgpa = .TRUE. ! Rain collection of aerosols lscgia = .TRUE. ! Ice collection of aerosols lscgic = .TRUE. ! Ice collection of cloud droplets lscgii = .TRUE. ! Self-collection of ice lscgip = .TRUE. ! Ice collection of rain drops ! Sub-process switches: Condensation lscndgas = .FALSE. ! --Aerosol precursor gas codensation lscndh2oae = .TRUE. ! --Condensation of water on aerosols (if FALSE, equilibrium assumed) lscndh2ocl = .TRUE. ! --Condensation of water on cloud droplets (and drizzle) lscndh2oic = .TRUE. ! --Condensation of water on ice particles ! Sub-process switches: Ice nucleation ice_hom = .TRUE. ! Homogeneous freezing ice_imm = .TRUE. ! Immersion freezing ice_dep = .TRUE. ! Deposition freezing ice_theta_dist = .TRUE. ! Use contact angle distribution for the IN population ! Sub-process swithces: Secondary ice ice_halmos = .FALSE. ! Rime splintering; Hallet-Mossop ! Misc switches lscheckarrays = .FALSE. lsdistupdate = .TRUE. ! Some useful switches for spinup and initilization lsfreeRH%switch = .TRUE. ! Use freely evolving RH. If FALSE, RH is limited by rhlim through entire simulation lsfreeRH%delay = 3600. ! Time until which RH is limited by rhlim, in case lsfreeRH%swithc == TRUE rhlim = 1.001 ! RH limit for SALSA during initialization and spinup lsfreeTheta%switch = .TRUE. ! Use freely evolving lower limit in contact angle integration for ice nucleation lsfreeTheta%delay = 0. ! Time until which lower limit theta is set to be at least initMinTheta, in case lsfreeTheta%switch == TRUE initMinTheta = 0. ! Lower limit for contact angle integration in degrees for initialization and spinup ! Parameters for ice nucleation contact angle distributions (used with ice_theta_dist == TRUE) mean_theta = 79.!46.! 132. ! Distribution mean sigma_theta = 12.!20.! 0.25 ! 20. ! Standard deviation ! m-D and A-D parameters for non-spherical ice iceShapeAlpha = 15.56999e-3 ! m = ALPHA * D**beta iceShapeBeta = 2.02 ! m = alpha * D**BETA iceShapeGamma = 0.55 ! A = GAMMA * D**sigma iceShapeSigma = 1.97 ! A = gamma * D**SIGMA ! Precipitation bin layout parameters bloPrc%nbins = 20 ! Number of bins bloPrc%dlo = 20.e-6 ! Lower limit of the smallest bin bloPrc%vol_ratio = 2. ! Volume ratio between bins ! Ice bin layout parameters bloIce%nbins = 20 ! Number of bins bloIce%dlo = 2.e-6 ! Lower limit of the smallest bin bloIce%vol_ratio = 3. ! Volume ratio between bins ! Initial aerosol size distributions nspec_dry = 2 listspec = 'SO4','OC','','','','','' isdtyp = 0 ! The number of entries must comply with the max number of _compounds_ (excluding water, currently 7) volDistA = 0.12, 0.88, 0.0, 0.0, 0.0, 0.0, 0.0 ! The number of entries must comply with the max number of _modes_ (currently 7) sigmagA = 1.5249, 1.5826,1.1811, 2.0, 2.0, 2.0, 2.0 ! Geometric stdev for initial aerosol size distribution for isdtyp == 0 (uniform) dpgA = 0.039144, 0.21519, 0.73482, 0.2, 0.2, 0.2, 0.2 ! Mode mean diameters in micrometers nA = 456, 155., 25., 0., 0., 0. ,0. ! Mode number concentrations in #/mg ! The number of entries must comply with the max number of _compounds_ (excluding water, currently 7) volDistB = 0.905, 0.0, 0.095, 0.0, 0.0, 0.0, 0.0 ! The number of entries must comply with the max number of _modes_ (currently 7) sigmagB = 1.5249, 1.5826,1.1811, 2.0, 2.0, 2.0, 2.0 ! Geometric stdev for initial aerosol size distribution for isdtyp == 0 (uniform) dpgB = 0.039144, 0.21519, 0.73482, 0.2, 0.2, 0.2, 0.2 ! Mode mean diameters in micrometers nB = 0., 0., 0., 0., 0., 0., 0. ! Mode number concentrations in #/mg /