search

NCAR / wawg_dev

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

f.e30.FWscHIST.ne30pg3_mg17_L135_cam6_4_038_beres0.70 #63

Open mijeong135 opened 1 month ago

mijeong135 commented 1 month ago

Case Name f.e30.FWscHIST.ne30pg3_mg17_L135_cam6_4_038_beres0.70 Keywords 1deg, ne30, WACCM, 135L, cam6_4_038, derecho, cam7

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

cd /glade/u/home/mijeong/cesm_tags git clone https://github.com/ESCOMP/CAM.git cam6_4_038 cd cam6_4_038 git checkout cam6_4_038 ./bin/git-fleximod update

Set up a WACCM run (135L), Tag: cam6_4_038 use specified chemistry (FWscHIST)

ssh derecho

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

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

cd /glade/u/home/mijeong/cesm3/cases/f.e30.FWscHIST.ne30pg3_mg17_L135_cam6_4_038_beres0.70

./case.setup

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

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

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.70D0 se_nsplit=4

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=6 ./xmlchange CONTINUE_RUN=FALSE ./xmlchange RESUBMIT=2 ./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

PeterHjortLauritzen commented 1 month ago

If you want to match Cecile's setup (see https://github.com/NCAR/amwg_dev/issues/583) a couple of things need to be added:

./xmlchange --append CAM_CONFIG_OPTS="-rad rrtmgp"

and

 seasalt_emis_scale             = 0.75D0
clubb_c8        =  4.95

to user_nl_cam.

I am not sure if this is needed in WACCM (in user_nl_cam):

 rad_climate        = 'A:Q:H2O', 'N:O2:O2',
         'A:CO2:CO2', 'N:ozone:O3',
         'A:N2O:N2O', 'A:CH4:CH4',
         'N:CFC11STAR:CFC11', 'A:CFC12:CFC12',
         'M:mam4_mode1:/glade/campaign/cesm/cesmdata/inputdata/atm/cam/physprops/mam4_mode1_rrtmg_aeronetdust_sig1.6_dgnh.48_c140304.nc',
     'M:mam4_mode2:/glade/campaign/cesm/cesmdata/inputdata/atm/cam/physprops/mam4_mode2_rrtmg_aitkendust_c141106.nc',
         'M:mam4_mode3:/glade/campaign/cesm/cesmdata/inputdata/atm/cam/physprops/mam4_mode3_rrtmg_aeronetdust_c141106.nc',
     'M:mam4_mode4:/glade/campaign/cesm/cesmdata/inputdata/atm/cam/physprops/mam4_mode4_rrtmg_c130628.nc',
         'N:VOLC_MMR1:/glade/campaign/cesm/cesmdata/inputdata/atm/cam/physprops/volc_camRRTMG_byradius_sigma1.6_mode1_c210211.nc',
     'N:VOLC_MMR2:/glade/campaign/cesm/cesmdata/inputdata/atm/cam/physprops/volc_camRRTMG_byradius_sigma1.6_mode2_c210211.nc',
         'N:VOLC_MMR3:/glade/campaign/cesm/cesmdata/inputdata/atm/cam/physprops/volc_camRRTMG_byradius_sigma1.2_mode3_c210211.nc'

Also, Cecile does not have:

dust_emis_fact = 1.3

Maybe @tilmes can comment on that?

tilmes commented 1 month ago

Regarding the dust factor, I would do what is used in Cecile's run for now.

PeterHjortLauritzen commented 1 month ago

Here are code mods for an increased del4 sponge:

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

Added code to increase del4 in top 5 layers:

Screenshot 2024-10-15 at 1 41 57 PM

If you are using se_nsplit=4 then you need to set

se_hypervis_subcycle=17

If you are trying se_nsplit=3 then

se_hypervis_subcycle=23

(=17*4/3)

mijeong135 commented 1 month ago

Adding this to my run:

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

!In user_nl_cam:

se_nsplit=4 se_hypervis_subcycle=17

mijeong135 commented 1 month ago

Hi Peter,

Francis said remove the 'rad_climate =...' part in user_nl_cam. In WACCM, CO2 is not advected and without including this, the model will use the default setting. We will see if this helps.

-Mijeong

On Tue, Oct 15, 2024 at 1:44 PM Peter Hjort Lauritzen < @.***> wrote:

Here are code mods for an increased del4 sponge:

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

Added code to increase del4 in top 5 layers:

Screenshot.2024-10-15.at.1.41.57.PM.png (view on web) https://github.com/user-attachments/assets/b062926b-e0b7-47be-9800-efe3b26c9064

If you are using se_nsplit=4 then you need to set

se_hypervis_subcycle=17

If you are trying se_nsplit=3 then

se_hypervis_subcycle=23

(=17*4/3)

— Reply to this email directly, view it on GitHub https://github.com/NCAR/wawg_dev/issues/63#issuecomment-2414866323, or unsubscribe https://github.com/notifications/unsubscribe-auth/BCCRGOLMA6LEVOSDVDPVHCDZ3VWCTAVCNFSM6AAAAABPTIVQI6VHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDIMJUHA3DMMZSGM . You are receiving this because you authored the thread.Message ID: @.***>

mijeong135 commented 2 weeks ago

Adding a bug fix on the vertical sponge layer diffusion Adam found (https://github.com/ESCOMP/CAM/issues/1182)

download vertical_diffusion.F90 (here)

cp vertical_diffusion.F90 ./SourceMods/src.cam change the line explained below

Below is an email from Adam ;---------- 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)

This pertains to line 1106 in the most recent cam tag, here. You can either change your source tree vertical_diffusion.F90, or you can copy the original into SourceMods/src.cam/ of an individual case, and implement the bug fix for an individual case.