&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 = 0.5 ! Max. timestep distim = 100. ! Timescale for the dissipation in sponge layer timmax = 25200. ! Length of the simulation runtype = "INITIAL" ! INITIAL or HISTORY (restart) run CCN = 100.e6 ! corflg = .false. ! Apply coriolis force prndtl = -0.333333 ! filprf = 'cloud_20200924_forces_L4b060OC_2ms5' ! Output filename profile hfilin = 'UAE_Cu_restart_6.014400s' ! History file name frqhis = 3600. !mcflg = .FALSE. ! Do mass conservation statistics sed_aero%switch = .FALSE. ! Calculate sedimentation of aerosol particles sed_cloud%switch =.TRUE.!.FALSE. ! - '' - cloud droplets sed_cloud%delay = 3600. sed_precp%switch = .TRUE. ! precipitation sed_precp%delay = 3600. sed_ice%switch = .FALSE. ! ice particles bulk_autoc%switch = .FALSE. bulk_autoc%delay = 3600. ! Autoconversion switch for level = 1-3 lnudging =.TRUE.!.FALSE. ! Master swtich for nudging scheme lemission = .FALSE. ! Master switch for aerosol emissions lpback = .FALSE. iradtyp = 3 ! Radiation/large scale forcing strtim = 267.25 ! Start time case_name = 'default' ! Case name for large-scale forcing schemes div = 1.e-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','rh','dn','press','srp','rrate','lwp','Np','sflx','rflx','albedo','Naba','Dwaba','Ncba','Dwcba','Npba','Dwpba','Nca','Dwca','Dwaa','Naa','aSO4a','cSO4a','pSO4a','aOCa','cOCa','pOCa', 's_n_activ', 'Ncbb','Dwcbb','Nabb','Dwabb','aSO4b','cSO4b','aOCb','cOCb','Ncb','Dwcb','Nab','Dwab' !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 = 3. itsflg = 1 ! Flag for temperature type in input sounding ipsflg = 1 ! - '' - for pressure / ! ------------------------------------------------------------------------------------- &surface isfctyp =5! change 27,04,2021 0 ! This is not entirely accurate for momentum vs vanZanten et al. 2011 zrough = 0.0002 ubmin = -0.25 sst = 283.15 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 = 7200. ! 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 = .TRUE. ! 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 = 132. ! Distribution mean sigma_theta = 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.35, 0.65, 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.8205, 1.494, 2.0, 2.0, 2.0, 2.0, 2.0 ! Geometric stdev for initial aerosol size distribution for isdtyp == 0 (uniform) dpgA = 0.0755, 0.1560, 0.7352, 0.2, 0.2, 0.2, 0.2 ! Mode mean diameters in micrometers nA =615., 927., 0., 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.03, 0.97, 0.0, 0.0, 0.0, 0.0, 0.0 ! The number of entries must comply with the max number of _modes_ (currently 7) sigmagB = 1.8205, 1.494, 1.6, 2.0, 2.0, 2.0, 2.0 ! Geometric stdev for initial aerosol size distribution for isdtyp == 0 (uniform) dpgB = 0.0755, 0.1560, 1.8, 0.2, 0.2, 0.2, 0.2 ! Mode mean diameters in micrometers nB = 264., 398., 0., 0., 0., 0., 0. ! Mode number concentrations in #/mg /