NCAR / wawg_dev

Repository for tracking CESM/WAWG development simulations
0 stars 0 forks source link

f.e23.FWmaHIST.ne16pg3_mg17_L110_tag132_beres0.85_HBMixMod_ZM1 #42

Open mijeong135 opened 1 week ago

mijeong135 commented 1 week ago

Case Name f.e23.FWmaHIST.ne16pg3_mg17_L110_tag132_beres0.85_HBMixMod_ZM1 Keywords 2deg, ne16, WACCM, 110L, tag132, derecho, cam_dev, HBMixMod, CTEM, ZM1


Case Dir /glade/u/home/mijeong/cesm2/cases/$CASE Run Dir /glade/derecho/scratch/mijeong/$CASE/run Archive Dir /glade/campaign/acom/acom-climate/mijeong/archive/$CASE


This is a test run with cam_dev with 110L with ZM1

Tag #132 (2 deg) cloned f.e23.FWmaHIST.ne16pg3_mg17_L135_tag132_beres0.65_HBMixMod (https://github.com/NCAR/wawg_dev/issues/40) HB Diffusion on but with Modified Mixing (https://github.com/NCAR/amwg_dev/issues/537)

110L ncdata file (mapped from the 135L to 110L): /glade/derecho/scratch/aherring/restart/inic/acom_ne16pg3_ne16pg3_mg17_1536_long2.cam.i.1980-01-01-00000.map_TO_L110.nc

./xmlchange CAM_CONFIG_OPTS="-phys cam_dev -nlev 110 -age_of_air_trcs -chem waccm_ma_mam5 -cppdefs -Dwaccm_debug -pcols 9"

Turn on ZM1: set zmconv_parcel_pbl = .false. in user_nl_cam

Add the three beres output fields (instantaneous - 6-hourly): NETDT, HDEPTH, MAXQ0


ssh derecho cd /glade/work/pel/cesm_tags/cam6_3_132/cime/scripts ./create_clone --case /glade/u/home/mijeong/cesm2/cases/f.e23.FWmaHIST.ne16pg3_mg17_L110_tag132_beres0.85_HBMixMod_ZM1 --clone /glade/u/home/mijeong/cesm2/cases/f.e23.FWmaHIST.ne16pg3_mg17_L135_tag132_beres0.65 --cime-output-root /glade/derecho/scratch/mijeong cd /glade/u/home/mijeong/cesm2/cases/f.e23.FWmaHIST.ne16pg3_mg17_L110_tag132_beres0.85_HBMixMod_ZM1 ./case.setup

HB diffusion - with modified mixing cp /glade/u/home/pel/src/for-unstable-Ri-hb/pbl_utils.F90 ./SourceMods/src.cam

HB Mixing Modified (https://github.com/NCAR/amwg_dev/issues/537)

edit user_nl_cam:

ncdata = '/glade/derecho/scratch/aherring/restart/inic/acom_ne16pg3_ne16pg3_mg17_1536_long2.cam.i.1980-01-01-00000.map_TO_L110.nc'

interpolate_output = .true.,.true.,.true.,.true. empty_htapes=.true.

fincl1 = 'U','V','Q','T','PSL','OMEGA','PS','PRECT','OMEGA500','OMEGA850', 'UTGWORO','VTGWORO','UTEND_TOT','VTEND_TOT','UTEND_CLUBB','VTEND_CLUBB', 'UTEND_PHYSTOT','VTEND_PHYSTOT','UTGWSPEC', 'BUTGWSPEC', 'UTEND_VDIFF','VTEND_VDIFF','UTEND_GWDTOT','VTEND_GWDTOT', 'UTEND_LUNART','VTEND_LUNART','UTEND_CORE','VTEND_CORE', 'O3','CO2','CO','O','N2O','NO2','H2O','UTGWSPEC','BUTGWSPEC','UU', 'VV','OMGAOMGA','VT','VU','OMEGAU','FRONTGF','KVH','KVM','KVH_CLUBB','ABS_dPSdt'

fincl2 = 'U','V','Q','T','PSL','OMEGA','PS','PRECT','OMEGA500','OMEGA850', 'UTGWORO','VTGWORO','UTGWSPEC','BUTGWSPEC','O3','CO2','N2O','H2O'

fincl5 = 'NETDT', 'HDEPTH', 'MAXQ0'

fincl8 = 'Uzm','Vzm','Wzm','THzm', 'VTHzm','WTHzm','UVzm','UWzm'

mfilt = 0, 5, 20, 30 avgflag_pertape = 'A','A','I','A' nhtfrq = 0,-24,-6,-24

phys_grid_ctem_nfreq = -24 phys_grid_ctem_za_nlat = 90 phys_grid_ctem_zm_nbas = 20

qbo_use_forcing = .false. use_gw_rdg_beta = .true. use_gw_convect_dp = .true. use_gw_front = .true. zmconv_parcel_pbl = .false.

se_statefreq = 144 effgw_beres_dp = 0.85D0 gw_qbo_hdepth_scaling = 0.25D0 effgw_cm = 1.0D0 effgw_rdg_beta = 1.2D0 effgw_rdg_beta_max = 1.2D0 gw_top_taper = .false. frontgfc = 1.75D-15 taubgnd = 1.25D-3 nlte_limit_co2=.true.

se_hypervis_subcycle_sponge = 2 se_nu_top = 1E6 se_hypervis_subcycle = 9 se_rsplit = 2 se_nsplit = 3 se_large_courant_incr = .true. se_nu_div = 6.0E15 se_nu = 6.0E15 se_sponge_del4_lev = 20 se_sponge_del4_nu_div_fac = 8

XMLchanges: ./xmlchange CAM_CONFIG_OPTS="-phys cam_dev -nlev 110 -age_of_air_trcs -chem waccm_ma_mam5 -cppdefs -Dwaccm_debug -pcols 9" ./xmlchange RESUBMIT=10 ./xmlchange DOUT_S_SAVE_INTERIM_RESTART_FILES=TRUE ./xmlchange DOUT_S_ROOT=/glade/campaign/acom/acom-climate/mijeong/archive/\$CASE ./xmlchange DOUT_S=TRUE ./xmlchange STOP_OPTION=nmonths ./xmlchange STOP_N=12

Check env_workflow.xml for clock settings <"JOB_WALLCLOCK_TIME" value = "12:00:00"> <"JOB_WALLCLOCK_TIME" value = "03:00:00">

./case.build --clean-all ./case.build ./preview_namelists ./preview_run ./case.submit