search

NCAR / wawg_dev

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

f.e30.FWmaHIST.ne30pg3_mg17_L135_cam6_4_038_beres0.25 #78

Open mijeong135 opened 5 days ago

mijeong135 commented 5 days ago

Case Name f.e30.FWmaHIST.ne30pg3_mg17_L135_cam6_4_038_beres0.25 Keywords 1deg, ne30, WACCM, 135L, cam6_4_038, derecho, cam7, effgw_beres_dp=0.25, Yuanpu's bug fix (Beres scheme), Adam's bug fix (vertical diffusion coefficient), Peter's vertical sponge layer mod

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

Tag cam6_4_038 use middle atmosphere chemistry (FWmaHIST) add Yuanpu's SourceMods (gw_convect.F90)

Fix a bug in vertical_diffusion.F90 In src/physics/cam/vertical_diffusion.F90 you just need to change the hard coded 1 in the kvm arrays to k: !!!kvm(:ncol,1) = kvm(:ncol,1)+kvm_sponge(k) !+++bug! kvm(:ncol,k) = kvm(:ncol,k)+kvm_sponge(k)

cp /glade/u/home/pel/src/for-mijeong2/global_norms_mod.F90 ./SourceMods/src.cam

Changing effgw_beres_dp = 0.25

ssh derecho

cd /glade/u/home/mijeong/cesm_tags/cam6_4_038/cime/scripts

./create_newcase --compset FWmaHIST --res ne30pg3_ne30pg3_mg17 --case /glade/u/home/mijeong/cesm3/cases/f.e30.FWmaHIST.ne30pg3_mg17_L135_cam6_4_038_beres0.25 --run-unsupported --project P93300607

cd /glade/u/home/mijeong/cesm3/cases/f.e30.FWmaHIST.ne30pg3_mg17_L135_cam6_4_038_beres0.25

./case.setup

Edit env_build.xml: ./xmlchange CAM_CONFIG_OPTS="-phys cam7 -nlev 135 -age_of_air_trcs -chem waccm_ma_mam5 -cppdefs -Dwaccm_debug -pcols 9"

Turn on RRTMGP ./xmlchange --append CAM_CONFIG_OPTS="-rad rrtmgp"

;sourcemods cp ~/sourcemods/gw_convect.F90 ./SourceMods/src.cam

cp /glade/u/home/mijeong/cesm_tags/cam6_4_038/src/physics/cam/vertical_diffusion.F90 ./SourceMods/src.cam !!!kvm(:ncol,1) = kvm(:ncol,1)+kvm_sponge(k) !+++bug! kvm(:ncol,k) = kvm(:ncol,k)+kvm_sponge(k)

cp /glade/u/home/pel/src/for-mijeong2/global_norms_mod.F90 ./SourceMods/src.cam

Edit user_nl_cam:

ncdata='/glade/campaign/cgd/amp/pel/inic/L135_ne30pg3_ne30pg3_mg17.cam.i.1979-01-11-00000.nc' bnd_topo= '/glade/work/juliob/bndtopo/ne30pg3_gmted2010_modis_bedmachine_nc3000_Laplace0100_noleak_20240720.nc'

interpolate_output = .true.,.true.,.true.,.true. interpolate_nlat = 192,192,192,192 interpolate_nlon = 288,288,288,288

empty_htapes=.true.

fincl1 = 'U','V','Q','T','PSL','OMEGA','PS','PRECT','OMEGA500','OMEGA850', 'FREQZM','CAPE','CMFMC_DP','ZMDT','PRECC','PRECL' fincl2 = 'U','V','Q','T','PSL','OMEGA','PS','PRECT','OMEGA500','OMEGA850', 'UTGWORO','VTGWORO','UTGWSPEC','BUTGWSPEC','O3','N2O','H2O' fincl3 = 'NETDT', 'HDEPTH', 'MAXQ0', 'PRECT', 'FLUT' 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 = -6 phys_grid_ctem_za_nlat = 90 phys_grid_ctem_zm_nbas = 120

inithist = 'YEARLY' se_statefreq = 144 effgw_beres_dp = 0.25D0 se_nsplit=4 se_hypervis_subcycle=17

seasalt_emis_scale = 0.75D0 clubb_c8 = 4.95

XMLchanges: ./xmlchange RUN_REFDATE=1980-01-01 ./xmlchange RUN_STARTDATE=1980-01-01 ./xmlchange STOP_OPTION=nmonths ./xmlchange STOP_N=12 ./xmlchange CONTINUE_RUN=FALSE ./xmlchange RESUBMIT=5 ./xmlchange DOUT_S_SAVE_INTERIM_RESTART_FILES=TRUE ./xmlchange DOUT_S_ROOT=/glade/campaign/acom/acom-climate/mijeong/archive/\$CASE ./xmlchange DOUT_S=TRUE

;env_mach_pes.xml (ntasks=-20 ; 20*128=2160) ./xmlchange NTASKS=-20

./case.build

Check env_workflow.xml for clock settings ./xmlchange JOB_WALLCLOCK_TIME=12:00:00 --subgroup case.run ./xmlchange JOB_WALLCLOCK_TIME=04:00:00 --subgroup case.st_archive

./preview_namelists ./preview_run ./case.submit