search

NCAR / wawg_dev

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

f.e23.FWmaHIST.ne16pg3_mg17_L135_tag132_beres0.65_HBdiff50off #33

Open mijeong135 opened 2 months ago

mijeong135 commented 2 months ago

Case Name f.e23.FWmaHIST.ne16pg3_mg17_L135_tag132_beres0.65_HBdiff50off Keywords 2deg, ne16, WACCM, 135L, tag132, derecho, cam_dev, HBdiffMods, CTEM, CF20

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

Tag #132 (2 deg) Clone of https://github.com/NCAR/wawg_dev/issues/24 with HBdiff off from the surface to below layer 50 (approx. 5 .455 hPa, 35.17 km)

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_L135_tag132_beres0.65_HBdiff50off --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_L135_tag132_beres0.65_HBdiff50off ./case.setup

copy vertical_diffusion.F90 from https://github.com/ESCOMP/CAM/blob/cam_development/src/physics/cam/vertical_diffusion.F90 cp this to ./SourceMods/src.cam

Set tendencies to zero from the surface to layer 50: Replace the line 1075 "do k=clubbtop(i),pverp" with "do k=50,pverp"

edit user_nl_cam:

interpolate_output = .true.,.true.,.true. interpolate_nlat=96 interpolate_nlon=144 empty_htapes=.true.

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

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'

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

phys_grid_ctem_nfreq = -6 phys_grid_ctem_za_nlat = 90 phys_grid_ctem_zm_nbas = 90

se_statefreq = 144 effgw_beres_dp = 0.65D0 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 = 2 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 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 ./case.build ./preview_namelists ./preview_run ./case.submit