Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

f.e30_cam6_4_036.FLTHIST.ne30_L32.cam5.001 #603

Open
cecilehannay opened this issue Oct 30, 2024 · 0 comments
Open

f.e30_cam6_4_036.FLTHIST.ne30_L32.cam5.001 #603

cecilehannay opened this issue Oct 30, 2024 · 0 comments
Labels
FLTHIST FLHIST compset L32 Vertical grid: L32

Comments

@cecilehannay
Copy link
Collaborator

Purpose:
Testing recent tags with cam5 physics with ne30_L32


Description:
Same as #596 but using cam5 physics instead of cam7

  • use cam5 physics and L32
    ./xmlchange CAM_CONFIG_OPTS="-phys cam5 -nlev 32"

  • pe_layout
    ./xmlchange NTASKS=1280

  • start date
    ./xmlchange RUN_STARTDATE=2000-01-01

  • user_nl_cam

! file                      h0      h1     h2       h3          h4        h5       h6       h7        h8      h9 
mfilt		        =       0,      30,    0,        100,       30,       4,       80,      40,        4,     48
nhtfrq              =       0,     -24,    0,        -24,       -24,      -6,      -3,      -3,       -6,     -1
ndens               =       2,       2,     2,       2,          2,        2,       2,       2,        2,      2
interpolate_output  =  .true.,  .true., .false., .false.,   .false.,     .false.,  .false., .false., .false.,  .false.
interpolate_nlat    =     768,      768,    192,     192,      192,       192,     192,      192,     192,   192
interpolate_nlon    =     1152,    1152,    288,     288,      288,       288,      288,     288,     288,   288

empty_htapes = .true.

! interpolated output
! output for ADF  monthly-interpolated (h0)    
fincl1 = 'AODDUST', 'AODVIS', 'AODDUSTdn','AODVISdn','BURDENBC', 'BURDENDUST', 'BURDENPOM', 'BURDENSEASALT', 
'BURDENSO4', 'BURDENSOA', 'CDNUMC', 'CLDHGH', 'CLDICE', 'CLDLIQ', 'CLDLOW', 'CLDMED', 'CLDTOT', 'CLOUD', 'FLDS', 'FLNS', 'FLNSC', 'FLNT', 'FLNTC', 'FLUT', 
'FSNS', 'FSNSC', 'FSNT', 'FSNTC', 'FSNTOA', 'ICEFRAC', 'LHFLX', 'LWCF','OMEGA500', 'PBLH','PRECC', 'PRECL', 'PRECSC', 'PRECSL', 'PRECT', 'PS', 'PSL','Q', 
'QFLX', 'RELHUM', 'SHFLX', 'SOLIN', 'SST', 'SWCF', 'T', 'TAUX', 'TAUY', 'TGCLDIWP', 'TGCLDLWP', 'TMQ', 'TREFHT', 'TS', 'U', 'U10', 'V'

!output for MDTF (h1)
fincl2='PRECT', 'PRECC', 'FLUT', 'U850', 'U200', 'V850', 'V200', 'OMEGA500', 'TS', 'SST', 'PSL'

! native grid
!monthly (h2)
fincl3='ACTNI', 'ACTNL', 'ACTREI', 'ACTREL', 'AODDUST', 'AODVIS', 'AODVISdn','BURDENBC', 'BURDENDUST', 'BURDENPOM', 'BURDENSEASALT', 
'BURDENSO4', 'BURDENSOA', 'CAPE', 'CCN3', 'CDNUMC', 'CLDHGH', 'CLDICE', 'CLDLIQ', 'CLDLOW', 'CLDMED', 'CLDTOT', 'CLOUD', 'CMFMC_DP', 
'CT_H2O', 'DCQ', 'DQCORE', 'DTCOND', 'DTCORE', 'DTV', 'EVAPPREC', 'EVAPSNOW', 'FCTI', 'FCTL', 'FICE', 'FLDS', 'FLNS', 'FLNSC', 'FLNT', 'FLNTC', 'FLUT', 
'FREQZM', 'FSDS', 'FSDSC', 'FSNS', 'FSNSC', 'FSNT', 'FSNTC', 'FSNTOA', 'ICEFRAC', 'LANDFRAC', 'LHFLX', 'LWCF', 'MPDICE', 'MPDLIQ', 'MPDQ', 'MPDT', 
'OCNFRAC', 'OMEGA', 'OMEGA500', 'PBLH', 'PHIS', 'PINT', 'PMID', 'PRECC', 'PRECL', 'PRECSC', 'PRECSL', 'PRECT', 'PS', 'PSL', 'PTEQ', 'PTTEND', 'Q', 
'QFLX', 'QRL', 'QRS', 'QTGW',  'RELHUM', 'SHFLX', 'SOLIN', 'SST', 'SWCF', 
'T', 'TAUX', 'TAUY', 'TFIX', 'TGCLDIWP', 'TGCLDLWP', 'TMQ', 'TREFHT', 'TS', 'TTGW', 'U', 'U10', 'UBOT', 'UTGWORO', 'UTGW_TOTAL', 
'V', 'VBOT', 'VTGWORO', 'VTGW_TOTAL', 'Z3', 'ZMDQ', 'ZMDT', 
'AODVISdn','CCN3', 'CDNUMC', 'H2O', 'NUMICE', 'NUMLIQ','OMEGA500',
'TROP_P','TROP_T','TROP_Z',
'num_a1','num_a2','num_a3',
'dgnd_a01','dgnd_a02','dgnd_a03',
'ncl_a1','ncl_a2','ncl_a3',
'bc_a1',
'dst_a1','dst_a3',
'pom_a1',
'so4_a1','so4_a2','so4_a3',
'soa_a1','soa_a2',
'CCN3','CCN4',
'TTGWORO','TAUGWX', 'TAUGWY', 'UTEND_PHYSTOT', 'UTEND_GWDTOT', 'UTEND_DCONV', 'UTEND_CORE',
'VTEND_PHYSTOT', 'VTEND_GWDTOT', 'VTEND_DCONV', 'VTEND_CORE'


!monthly  (h3)
fincl4 = 'Uzm','Vzm','Wzm','THzm', 'VTHzm','WTHzm','UVzm','UWzm'    
phys_grid_ctem_nfreq=-6
phys_grid_ctem_zm_nbas=120
!ne30 value
phys_grid_ctem_za_nlat=90 

!daily - average (h4)
fincl5 = 'SHFLX', 'LHFLX', 'PBLH',      
'PRECT', 'PRECC', 'FLUT', 'U850', 'U200', 'V850', 'V200', 'OMEGA500', 'TS', 'SST', 'PSL'

!6 hourly - instantaneous with input from Julio (h5)
fincl6 = 'PS:I','U:I','V:I','T:I','OMEGA:I', 'PMID:I', 'PSDRY:I'

!3 hourly - average (h6)
fincl7= 'PRECC:A', 'PRECL:A'	

!3 hourly - instantaneous (h7)
fincl8= 'PRECC:I','PRECL:I',
'U850:I','V850:I','T500:I','T200:I','Z850:I','Z500:I','Z300:I','PSL:I','PS:I','U10:I','UBOT:I','VBOT:I',
'FLUT:I','TMQ:I','OMEGA500:I','uIVT:I', 'vIVT:I' 

!1 hourly - instantaneous for ming (h9)
fincl10 = 'Z500:I','PSL:I','U10:I','Z250:I','PRECT:A','PRECC:A','uIVT:I','vIVT:I','U850:I','U200:I','V850:I','V200:I','Q850:I','T850:I','FLUT:I'



seasalt_emis_scale             = 0.75D0

clubb_c8		=  4.95 

ubc_file_path=  "/glade/campaign/cesm/cesmdata/inputdata/atm/cam/chem/ubc/b.e21.BWHIST.f09_g17.CMIP6-historical-WACCM.ensAvg123.cam.h0zm.H2O.1849-2014_c240604.nc"

se_statefreq=144

ncdata = '/glade/campaign/cgd/amp/juliob/ERA5/ne30np4/L32_cdf5/ERA5_x_ne30np4_L32_rgC2_WO_aero.2000-01-01-00000.nc'

use_gw_front        = .true.
taubgnd        = 1.5D-6

tracer_cnst_datapath    = '/glade/campaign/cesm/cesmdata/cseg/inputdata/atm/cam/chem/trop_mozart_aero/oxid'
tracer_cnst_file		= 'oxid_rcp45_v1_1.9x2.5_L26_1995-2105_c100316.nc'
tracer_cnst_filelist	= 'oxid_1.9x2.5_L26_clim_list.c090805.txt'
tracer_cnst_specifier	= 'O3','OH','NO3','HO2'
tracer_cnst_type		= 'INTERP_MISSING_MONTHS'

 ext_frc_specifier		= 'SO2         -> /glade/campaign/cesm/cesmdata/cseg/inputdata/atm/cam/chem/trop_mozart_aero/emis/rcp45_1999-2100/RCP45_mam3_so2_elev_2000-2100_c20111116.nc',
  'bc_a1       -> /glade/campaign/cesm/cesmdata/cseg/inputdata/atm/cam/chem/trop_mozart_aero/emis/rcp45_1999-2100/RCP45_mam3_bc_elev_2000-2100_c20111116.nc',
  'num_a1      -> /glade/campaign/cesm/cesmdata/cseg/inputdata/atm/cam/chem/trop_mozart_aero/emis/rcp45_1999-2100/RCP45_mam3_num_a1_elev_2000-2100_c20111116.nc',
  'num_a2      -> /glade/campaign/cesm/cesmdata/cseg/inputdata/atm/cam/chem/trop_mozart_aero/emis/rcp45_1999-2100/RCP45_mam3_num_a2_elev_2000-2100_c20111116.nc',
  'pom_a1      -> /glade/campaign/cesm/cesmdata/cseg/inputdata/atm/cam/chem/trop_mozart_aero/emis/rcp45_1999-2100/RCP45_mam3_oc_elev_2000-2100_c20111116.nc',
  'so4_a1      -> /glade/campaign/cesm/cesmdata/cseg/inputdata/atm/cam/chem/trop_mozart_aero/emis/rcp45_1999-2100/RCP45_mam3_so4_a1_elev_2000-2100_c20111116.nc',
  'so4_a2      -> /glade/campaign/cesm/cesmdata/cseg/inputdata/atm/cam/chem/trop_mozart_aero/emis/rcp45_1999-2100/RCP45_mam3_so4_a2_elev_2000-2100_c20111116.nc'
 srf_emis_specifier		= 'DMS       -> /glade/campaign/cesm/cesmdata/cseg/inputdata/atm/cam/chem/trop_mozart_aero/emis/aerocom_mam3_dms_surf_1849-2100_c111017.nc',
  'SO2       -> /glade/campaign/cesm/cesmdata/cseg/inputdata/atm/cam/chem/trop_mozart_aero/emis/rcp45_1999-2100/RCP45_mam3_so2_surf_2000-2100_c20111116.nc',
  'SOAG      -> /glade/campaign/cesm/cesmdata/cseg/inputdata/atm/cam/chem/trop_mozart_aero/emis/rcp45_1999-2100/RCP45_soag_1.5_surf_2000-2100_c20111116.nc',
  'bc_a1     -> /glade/campaign/cesm/cesmdata/cseg/inputdata/atm/cam/chem/trop_mozart_aero/emis/rcp45_1999-2100/RCP45_mam3_bc_surf_2000-2100_c20111116.nc',
  'num_a1    -> /glade/campaign/cesm/cesmdata/cseg/inputdata/atm/cam/chem/trop_mozart_aero/emis/rcp45_1999-2100/RCP45_mam3_num_a1_surf_2000-2100_c20111116.nc',
  'num_a2    -> /glade/campaign/cesm/cesmdata/cseg/inputdata/atm/cam/chem/trop_mozart_aero/emis/rcp45_1999-2100/RCP45_mam3_num_a2_surf_2000-2100_c20111116.nc',
  'pom_a1    -> /glade/campaign/cesm/cesmdata/cseg/inputdata/atm/cam/chem/trop_mozart_aero/emis/rcp45_1999-2100/RCP45_mam3_oc_surf_2000-2100_c20111116.nc',
  'so4_a1    -> /glade/campaign/cesm/cesmdata/cseg/inputdata/atm/cam/chem/trop_mozart_aero/emis/rcp45_1999-2100/RCP45_mam3_so4_a1_surf_2000-2100_c20111116.nc',
  'so4_a2    -> /glade/campaign/cesm/cesmdata/cseg/inputdata/atm/cam/chem/trop_mozart_aero/emis/rcp45_1999-2100/RCP45_mam3_so4_a2_surf_2000-2100_c20111116.nc'

uwshcu_rpen = 10
zmconv_ke = 1e-6
zmconv_ke_lnd = 1e-6
use_gw_convect_dp = .true.
  • user_nl_clm
hist_fincl1 = 'QICE_FORC','TSRF_FORC','TOPO_FORC'                      
hist_fincl2 = 'PCT_LANDUNIT','ATM_TOPO','TSA','LWup','LWdown',         
             'ALBD','ALBI','ALBGRD','ALBGRI','SWdown','SWup',          
             'RAIN_ICE','SNOW_ICE','QSNOMELT_ICE','QICE_MELT',         
             'QSNOFRZ_ICE','QRUNOFF_ICE','QICE','QSOIL_ICE',           
             'SNOW_DEPTH','SNOWDP','H2OSNO','FSNO','QICE_FRZ', 
             'QSNOMELT','SOILWATER_10CM','TWS','QFLX_EVAP_TOT',
             'H2OSOI','TSOI','TSOI_10CM','TSOI_ICE','H2OSNO_ICE' 
             'QRUNOFF','SNOW','RAIN','FSH','TREFMNAV','TREFMXAV'                             
hist_avgflag_pertape(1) = 'A'                                           
hist_avgflag_pertape(2) = 'A'                                           
hist_nhtfrq = 0,-24                                                     
hist_mfilt  = 1,365                                                     
hist_ndens  = 2,2  
fsurdat = '/glade/derecho/scratch/aherring/restart/ctsm5.2.028/tools/mksurfdata_esmf/surfdata_ne30pg3_SSP3-7.0_1979_78pfts_c241014.nc'
flanduse_timeseries = '/glade/derecho/scratch/aherring/restart/ctsm5.2.028/tools/mksurfdata_esmf/landuse.timeseries_ne30pg3_SSP3-7.0_1979-2026_78pfts_c241014.nc' 
finidat = '/glade/campaign/cesm/cesmdata/inputdata/lnd/clm2/initdata_esmf/ctsm5.3/ctsm53n04ctsm52028_ne30pg3t232_BgcCrop_exice_hist.clm60.r.2000-01-01.nc'
use_init_interp = .true.
use_excess_ice = .false.
  • SourceMods
/glade/campaign/cesm/cesmdata/cseg/runs/cesm2_0/f.e30_cam6_4_036.FLTHIST.ne30_L32.cam5.001/SourceMods/src.cam/cam_diagnostics.F90 

Case directory:


Sandbox:


Diagnostics:


Output:

  • Raw output on derecho (if still available)
    /glade/derecho/scratch/hannay/archive/f.e30_cam6_4_036.FLTHIST.ne30_L32.cam5.001/
  • CAM Climos: /glade/campaign/cgd/amp/amwg/climo/f.e30_cam6_4_036.FLTHIST.ne30_L32.cam5.001

Contacts:
@cecilehannay @JulioTBacmeister


Extra details:

To set this run, we identified the differences in tuning between cesm2-cam5 and current version.

cesm2-cam5 -> cesm3-cam5

Differences we don't care about:

  • seasalt_emis_scale=1.35->0.75
  • topo: ok to use current cam7 topo
  • co2: ghg->flbc
  • gw_dc = 0-> 2.5
  • pgwv = 0->32
  • tau_0_ubc = false ->true
  • use_gw_front = false->true
  • UBC_FILE

Differences that might matter (for these parameters, we set the values from cam5)

  • uwshcu_rpen = 10->5
  • zmconv_ke = 1->5e-6
  • zmconv_ke_lnd = 1->3e-6

Julio would like also to set:

  • use_gw_convect_dp = .true.

Because cam5 stoppes in 20025, we use RCP45 files to run after 2005.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
FLTHIST FLHIST compset L32 Vertical grid: L32
Projects
None yet
Development

No branches or pull requests

1 participant