search

NCAR / amwg_dev

Repo to store model sandboxes and cases used for CAM development
9 stars 2 forks source link

b.e23_alpha16b.BLT1850.ne30_t232.044 #399

Open cecilehannay opened 1 year ago

cecilehannay commented 1 year ago

Description:

Same as https://github.com/NCAR/amwg_dev/issues/398, but using the 37 tuning for the non-clubb tunings and also for clubb_gamma_coef and clubb_gamma_coefb 

clubb_gamma_coef       = 0.25->0.3
clubb_gamma_coefb      = 0.32->0.3

micro_mg_dcs           = 650->500.D-6
micro_mg_vtrmi_factor  = 1.25->1.0
microp_aero_wsub_scale = 0.01->0.25

cldfrc_dp2             =  750->500.0D0 
zmconv_dmpdz           =  -0.0005->-1.0E-3

This run includes namelist changes and SourceMods. All of these described below.

Case directory: Locally (if still available): /glade/p/cesmdata/cseg/runs/cesm2_0/b.e23_alpha16b.BLT1850.ne30_t232.044

On github: https://github.com/NCAR/amwg_dev/tree/b.e23_alpha16b.BLT1850.ne30_t232.044

Sandbox: Locally (if still available): /glade/work/hannay/cesm_tags/cesm2_3_alpha16b_taus

The details of how to create this sandbox are available below.

Diagnostics: Diags (if available)

https://webext.cgd.ucar.edu/BLT1850/b.e23_alpha16b.BLT1850.ne30_t232.044/

Output: Output (if still available): /glade/derecho/scratch/hannay/archive/b.e23_alpha16b.BLT1850.ne30_t232.044/

Contacts: @bstephens82, @vlarson, @JulioTBacmeister, @PeterHjortLauritzen, @adamrher, @cecilehannay

Extra details:

Sandbox

# Clone the repository
cd /glade/work/hannay/cesm_tags 

# get the cesm tag
git clone https://github.com/ESCOMP/cesm cesm2_3_alpha16b_taus

# check out the tag
cd  /glade/work/hannay/cesm_tags/cesm2_3_alpha16b_taus
git checkout  cesm2_3_alpha16b
./manage_externals/checkout_externals -o

# Modify the clubb external 
cd /glade/work/hannay/cesm_tags/cesm2_3_alpha16b_taus/components/cam

In Externals_CAM.cfg

##replace:

[clubb]
local_path = src/physics/clubb
protocol = svn
repo_url = https://github.com/larson-group/clubb_release/tags/
tag = clubb_4ncar_20221129_59cb19f_20230330_branchtag/src/CLUBB_core
required = True

## with
[clubb]
local_path = src/physics/clubb
protocol = svn
repo_url = https://github.com/bstephens82/clubb_src/tags/
tag = bstephens82.clubb_src.aa1d359
required = True

# checkout the updated clubb extrenal
cd  /glade/work/hannay/cesm_tags/cesm2_3_alpha16b_taus
rm -rf components/cam/src/physics/clubb
cd components/cam/
../../manage_externals/checkout_externals -e Externals_CAM.cfg

# modify the namelist files 
cd /glade/work/hannay/cesm_tags/cesm2_3_alpha16b_taus/components/cam
## you need to modify:
bld/build-namelist   
bld/namelist_files/namelist_defaults_cam.xml  
bld/namelist_files/namelist_definition.xml

Namelist changes

user_nl_cam

mfilt           =       0,       5,     20,      40,      12,       120,      1,   1
nhtfrq              =       0,     -24,    -24,      -3,       0,       -2,      0,  -8760
ndens               =       2,       2,      2,       2,       2,       1,      2,   1
interpolate_output  =  .true.,  .true., .true., .false., .false., .true.,  .true.
interpolate_nlat    =     192,     192,    192,     192,     192,     192,   192
interpolate_nlon    =     288,     288,    288,     288,     288,     288,   288 

empty_htapes = .true.

fincl1 = 'ACTNI', 'ACTNL', 'ACTREI', 'ACTREL', 'AODDUST', 'AODVIS', 'AODVISdn','BURDENBC', 'BURDENDUST', 'BURDENPOM', 'BURDENSEASALT', 
'BURDENSO4', 'BURDENSOA', 'CAPE', 'CCN3', 'CDNUMC', 'CH4', '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', 'RCMTEND_CLUBB', 'RELHUM', 'RVMTEND_CLUBB', 'SHFLX', 'SOLIN', 'SST', 'STEND_CLUBB', 'SWCF', 
'T', 'TAUX', 'TAUY', 'TFIX', 'TGCLDIWP', 'TGCLDLWP', 'TMQ', 'TREFHT', 'TS', 'TTGW', 'U', 'U10', 'UBOT', 'UTGWORO', 'UTGW_TOTAL', 
'V', 'VBOT', 'VTGWORO', 'VTGW_TOTAL', 'WPRTP_CLUBB', 'WPTHLP_CLUBB', 'Z3', 'ZMDQ', 'ZMDT', 'N2O', 'CO2','CFC11','CFC12', 
'CLD_MISR','FISCCP1_COSP','CLD_CAL','CLD_MISR','CLDTOT_CAL','CLDHGH_CAL', 'CLDMED_CAL','CLDLOW_CAL','CLMODIS', 'AODVISdn','AODDUSTdn',
'CCN3', 'CDNUMC', 'H2O', 'NUMICE', 'NUMLIQ','OMEGA500'

fincl3 = 'PRECT', 'PRECC', 'FLUT', 'U850', 'U200', 'V850', 'V200', 'OMEGA', 'PSL'

fincl4 =  'PRECC','PRECL'

fincl5 = 'Uzm','Vzm','Wzm','THzm', 'VTHzm','WTHzm','UVzm','UWzm'
phys_grid_ctem_nfreq=-6
phys_grid_ctem_zm_nbas=120
phys_grid_ctem_za_nlat=90

fincl7= 'AQSO4_H2O2','AQSO4_O3', 'bc_a1', 'bc_a4', 'dst_a1', 'dst_a2', 'dst_a3', 'ncl_a1',
'ncl_a1', 'ncl_a2', 'ncl_a3', 'pom_a1', 'pom_a4', 'so4_a1', 'so4_a2', 'so4_a3',
'soa_a1', 'num_a1', 'num_a2', 'num_a3', 'num_a4',
'bc_a1SFWET', 'bc_a4SFWET', 'dst_a1SFWET', 'dst_a2SFWET', 'dst_a3SFWET', 'ncl_a1SFWET',
'ncl_a2SFWET', 'ncl_a3SFWET', 'pom_a1SFWET', 'pom_a4SFWET', 'so4_a1SFWET', 'so4_a2SFWET', 'so4_a3SFWET', 'soa_a1SFWET',
'soa_a2SFWET', 'bc_c1SFWET', 'bc_c4SFWET', 'dst_c1SFWET', 'dst_c2SFWET', 'dst_c3SFWET', 'ncl_c1SFWET', 'ncl_c2SFWET',
'ncl_c3SFWET', 'pom_c1SFWET', 'pom_c4SFWET', 'so4_c1SFWET', 'so4_c2SFWET', 'so4_c3SFWET', 'soa_c1SFWET', 'soa_c2SFWET',
'bc_a1DDF', 'bc_a4DDF', 'dst_a1DDF', 'dst_a2DDF', 'dst_a3DDF', 'ncl_a1DDF', 'ncl_a2DDF', 'ncl_a3DDF',
'pom_a1DDF', 'pom_a4DDF', 'so4_a1DDF', 'so4_a2DDF', 'so4_a3DDF', 'soa_a1DDF', 'soa_a2DDF',
'so4_a1_CLXF', 'so4_a2_CLXF', 'SFbc_a4', 'SFpom_a4', 'SFso4_a1', 'SFso4_a2',
'so4_a1_sfgaex1', 'so4_a2_sfgaex1', 'so4_a3_sfgaex1', 'soa_a1_sfgaex1', 'soa_a2_sfgaex1',
'SFdst_a1','SFdst_a2', 'SFdst_a3', 'SFncl_a1', 'SFncl_a2', 'SFncl_a3',
'num_a2_sfnnuc1', 'SFSO2', 'OCN_FLUX_DMS', 'SAD_SULFC', 'SAD_TROP', 'SAD_AERO'

srf_emis_specifier      = 'bc_a4 -> /glade/p/cesmdata/inputdata/atm/cam/chem/emis/historical_ne30pg3/emissions-cmip6_bc_a4_anthro_surface_175001-201412_ne30pg3_c20191118.nc',
         'bc_a4 -> /glade/p/cesmdata/inputdata/atm/cam/chem/emis/historical_ne30pg3/emissions-cmip6_bc_a4_bb_surface_175001-201512_ne30pg3_c20191118.nc',
         'DMS -> /glade/p/cesmdata/inputdata/atm/cam/chem/emis/historical_ne30pg3/emissions-cmip6_DMS_bb_surface_175001-201512_ne30pg3_c20191118.nc',
         'num_a1 -> /glade/p/cesmdata/inputdata/atm/cam/chem/emis/historical_ne30pg3/emissions-cmip6_num_so4_a1_bb_surface_175001-201512_ne30pg3_c20191118.nc',
         'num_a1 -> /glade/p/cesmdata/inputdata/atm/cam/chem/emis/historical_ne30pg3/emissions-cmip6_num_so4_a1_anthro-ag-ship_surface_mol_175001-201412_ne30pg3_c20200103.nc',
         'num_a2 -> /glade/p/cesmdata/inputdata/atm/cam/chem/emis/historical_ne30pg3/emissions-cmip6_num_so4_a2_anthro-res_surface_mol_175001-201412_ne30pg3_c20200103.nc',
         'num_a4 -> /glade/p/cesmdata/inputdata/atm/cam/chem/emis/historical_ne30pg3/emissions-cmip6_num_bc_a4_bb_surface_175001-201512_ne30pg3_c20191118.nc',
         'num_a4 -> /glade/p/cesmdata/inputdata/atm/cam/chem/emis/historical_ne30pg3/emissions-cmip6_num_bc_a4_anthro_surface_175001-201412_ne30pg3_c20191118.nc',
         'num_a4 -> /glade/p/cesmdata/inputdata/atm/cam/chem/emis/historical_ne30pg3/emissions-cmip6_num_pom_a4_anthro_surface_175001-201412_ne30pg3_c20191118.nc',
         'num_a4 -> /glade/p/cesmdata/inputdata/atm/cam/chem/emis/historical_ne30pg3/emissions-cmip6_num_pom_a4_bb_surface_175001-201512_ne30pg3_c20191118.nc',
         'pom_a4 -> /glade/p/cesmdata/inputdata/atm/cam/chem/emis/historical_ne30pg3/emissions-cmip6_pom_a4_anthro_surface_175001-201412_ne30pg3_c20191118.nc',
         'pom_a4 -> /glade/p/cesmdata/inputdata/atm/cam/chem/emis/historical_ne30pg3/emissions-cmip6_pom_a4_bb_surface_175001-201512_ne30pg3_c20191118.nc',
         'SO2 -> /glade/p/cesmdata/inputdata/atm/cam/chem/emis/historical_ne30pg3/emissions-cmip6_SO2_anthro-ag-ship-res_surface_mol_175001-201412_ne30pg3_c20200103.nc',
         'SO2 -> /glade/p/cesmdata/inputdata/atm/cam/chem/emis/historical_ne30pg3/emissions-cmip6_SO2_anthro-ene_surface_mol_175001-201412_ne30pg3_c20200103.nc',
         'SO2 -> /glade/p/cesmdata/inputdata/atm/cam/chem/emis/historical_ne30pg3/emissions-cmip6_SO2_bb_surface_175001-201512_ne30pg3_c20191118.nc',
         'so4_a1 -> /glade/p/cesmdata/inputdata/atm/cam/chem/emis/historical_ne30pg3/emissions-cmip6_so4_a1_anthro-ag-ship_surface_mol_175001-201412_ne30pg3_c20200103.nc',
         'so4_a1 -> /glade/p/cesmdata/inputdata/atm/cam/chem/emis/historical_ne30pg3/emissions-cmip6_so4_a1_bb_surface_175001-201512_ne30pg3_c20191118.nc',
         'so4_a2 -> /glade/p/cesmdata/inputdata/atm/cam/chem/emis/historical_ne30pg3/emissions-cmip6_so4_a2_anthro-res_surface_mol_175001-201412_ne30pg3_c20200103.nc',
         'SOAE -> 2.5592*/glade/p/cesmdata/cseg/inputdata/atm/cam/chem/emis/CMIP6_emissions_1750_2015/emissions-cmip6_BENZENE_anthro_surface_1750-2015_0.9x1.25_c20170608.nc',
         'SOAE -> 2.5592*/glade/p/cesmdata/cseg/inputdata/atm/cam/chem/emis/CMIP6_emissions_1750_2015/emissions-cmip6_BENZENE_bb_surface_1750-2015_0.9x1.25_c20170322.nc',
         'SOAE -> 0.5954*/glade/p/cesmdata/cseg/inputdata/atm/cam/chem/emis/CMIP6_emissions_1750_2015/emissions-cmip6_ISOP_bb_surface_1750-2015_0.9x1.25_c20170322.nc',
         'SOAE -> 5.1004*/glade/p/cesmdata/cseg/inputdata/atm/cam/chem/emis/CMIP6_emissions_1750_2015/emissions-cmip6_MTERP_bb_surface_1750-2015_0.9x1.25_c20170322.nc',
         'SOAE -> 8.2367*/glade/p/cesmdata/cseg/inputdata/atm/cam/chem/emis/CMIP6_emissions_1750_2015/emissions-cmip6_TOLUENE_anthro_surface_1750-2015_0.9x1.25_c20170608.nc',
         'SOAE -> 8.2367*/glade/p/cesmdata/cseg/inputdata/atm/cam/chem/emis/CMIP6_emissions_1750_2015/emissions-cmip6_TOLUENE_bb_surface_1750-2015_0.9x1.25_c20170322.nc',
         'SOAE -> 6.5013*/glade/p/cesmdata/cseg/inputdata/atm/cam/chem/emis/CMIP6_emissions_1750_2015/emissions-cmip6_XYLENES_anthro_surface_1750-2015_0.9x1.25_c20170608.nc',
         'SOAE -> 6.5013*/glade/p/cesmdata/cseg/inputdata/atm/cam/chem/emis/CMIP6_emissions_1750_2015/emissions-cmip6_XYLENES_bb_surface_1750-2015_0.9x1.25_c20170322.nc',
         'SOAE -> 8.5371*/glade/p/cesmdata/cseg/inputdata/atm/cam/chem/emis/CMIP6_emissions_1750_2015/emissions-cmip6_IVOC_anthro_surface_1750-2015_0.9x1.25_c20170608.nc',
         'SOAE -> 8.5371*/glade/p/cesmdata/cseg/inputdata/atm/cam/chem/emis/CMIP6_emissions_1750_2015/emissions-cmip6_IVOC_bb_surface_1750-2015_0.9x1.25_c20170322.nc',
         'SOAE -> 16.6500*/glade/p/cesmdata/cseg/inputdata/atm/cam/chem/emis/CMIP6_emissions_1750_2015/emissions-cmip6_SVOC_anthro_surface_1750-2015_0.9x1.25_c20170608.nc',
         'SOAE -> 16.6500*/glade/p/cesmdata/cseg/inputdata/atm/cam/chem/emis/CMIP6_emissions_1750_2015/emissions-cmip6_SVOC_bb_surface_1750-2015_0.9x1.25_c20170322.nc'

  megan_specifier = 'SOAE = 0.5954*isoprene + 5.1004*(carene_3 + pinene_a + thujene_a + bornene +',
      ' terpineol_4 + terpineol_a + terpinyl_ACT_a + myrtenal + sabinene + pinene_b + camphene +',
      ' fenchene_a + limonene + phellandrene_a + terpinene_a + terpinene_g + terpinolene +',
      ' phellandrene_b + linalool + ionone_b + geranyl_acetone + neryl_acetone + jasmone +',
      ' verbenene + ipsenol + myrcene + ocimene_t_b + ocimene_al + ocimene_c_b + 2met_nonatriene) + ',
      ' 12.3942*(farnescene_a + caryophyllene_b + acoradiene + aromadendrene + bergamotene_a +',
      ' bergamotene_b + bisabolene_a + bisabolene_b + bourbonene_b + cadinene_d + cadinene_g +',
      ' cedrene_a + copaene_a + cubebene_a + cubebene_b + elemene_b + farnescene_b +',
      ' germacrene_B + germacrene_D + gurjunene_b + humulene_a + humulene_g + isolongifolene +',
      ' longifolene + longipinene + muurolene_a + muurolene_g + selinene_b + selinene_d +',
      ' nerolidol_c + nerolidol_t)'

!----non-CLUBB parameters
dust_emis_fact = 0.80D0
seasalt_emis_scale=1.5

micro_mg_dcs=500.D-6
micro_mg_vtrmi_factor = 1.0
microp_aero_wsubi_scale  = 1.0
microp_aero_wsub_scale  = 0.25
cldfrc_dp2             =  500.0D0 
zmconv_dmpdz           =  -1.0E-3 

!----CLUBB parameters
clubb_gamma_coef = 0.3
clubb_gamma_coefb = 0.3

clubb_C1 = 0.80000000000000
clubb_C1b = 0.80000000000000
clubb_C1c = 0.750000000000000
clubb_C2rt = 1.00000000000000
clubb_C2thl = 1.00000000000000
clubb_C2rtthl = 1.00000000000000
clubb_C4 = 1.00000000000000
clubb_c_uu_shr = 0.1
clubb_C_uu_buoy = 0.000000000000000E+000
clubb_C6rt = 1.00000000000000
clubb_C6rtb = 1.00000000000000
clubb_C6rtc = 0.500000000000000
clubb_C6thl = 1.00000000000000
clubb_C6thlb = 1.00000000000000
clubb_C6thlc = 0.500000000000000
clubb_c7 = 0.1
clubb_C7b = 0.700000000000000
clubb_C7c = 0.500000000000000
clubb_C8 = 0.80000000000
clubb_C8b = 0.000000000000000E+000
clubb_C11 = 0.500000000000000
clubb_C11b = 0.800000000000000
clubb_C11c = 0.850000000000000
clubb_C12 = 1.00000000000000
clubb_C14 = 0.80000000000000
clubb_C_wp3_pr_turb = 1.00000000000000
clubb_C_wp2_splat = 0.000000000000000E+000
clubb_C6rt_Lscale0 = 14.0000000000000
clubb_C6thl_Lscale0 = 14.0000000000000
clubb_C7_Lscale0 = 0.850000000000000
clubb_wpxp_L_thresh = 100.000000000000
clubb_c_K = 0.200000000000000
clubb_c_K1 = 1.00000000000000
clubb_nu1 = 10.0000000000000
clubb_c_K2 = 0.100000000000000
clubb_nu2 = 1.00000000000000
clubb_c_K6 = 0.375000000000000
clubb_nu6 = 5.00000000000000
clubb_c_K8 = 10.0000000000000
clubb_nu8 = 20.0000000000000
clubb_c_K9 = 0.250000000000000
clubb_nu9 = 20.0000000000000
clubb_nu10 = 0.000000000000000E+000
clubb_c_K_hm = 0.750000000000000
clubb_c_K_hmb = 0.750000000000000
clubb_K_hm_min_coef = 0.100000000000000
clubb_nu_hm = 1.50000000000000
clubb_slope_coef_spread_DG_means_w = 21.0000000000000
clubb_pdf_component_stdev_factor_w = 1.00000000000000
clubb_coef_spread_DG_means_rt = 0.800000000000000
clubb_coef_spread_DG_means_thl = 0.800000000000000
clubb_gamma_coefc = 1.20000000000000
clubb_mu = 5.000000000000000E-004
clubb_beta = 2.00000000000000
clubb_lmin_coef = 0.100000000000000
clubb_omicron = 0.500000000000000
clubb_zeta_vrnce_rat = 0.000000000000000E+000
clubb_upsilon_precip_frac_rat = 0.550000000000000
clubb_lambda0_stability_coef = 3.000000000000000E-002
clubb_mult_coef = 1.00000000000000
clubb_taumin = 90.0000000000000
clubb_taumax = 3600.00000000000
clubb_Lscale_mu_coef = 2.00000000000000
clubb_Lscale_pert_coef = 0.100000000000000
clubb_alpha_corr = 0.150000000000000
clubb_Skw_denom_coef = 0.000000000000000E+000
clubb_c_K10 = 0.300000000000000
clubb_c_K10h = 0.350000000000000
clubb_thlp2_rad_coef = 1.00000000000000
clubb_thlp2_rad_cloud_frac_thresh = 0.100000000000000
clubb_up2_sfc_coef = 2.00000000000000
clubb_Skw_max_mag = 4.50000000000000
clubb_C_invrs_tau_bkgnd = 1.10000000000000
clubb_C_invrs_tau_sfc = 5.000000000000000E-002
clubb_C_invrs_tau_shear = 0.220000000000000
clubb_C_invrs_tau_N2 = 0.300000000000000
clubb_C_invrs_tau_N2_wp2 = 0.100000000000000
clubb_C_invrs_tau_N2_xp2 = 0.000000000000000E+000
clubb_C_invrs_tau_N2_wpxp = 2.00000000000000
clubb_C_invrs_tau_N2_clear_wp3 = 2.00000000000000
clubb_C_invrs_tau_wpxp_Ri = 3.00000000000000
clubb_C_invrs_tau_wpxp_N2_thresh = 3.500000000000000E-004
clubb_xp3_coef_base = 0.250000000000000
clubb_xp3_coef_slope = 1.000000000000000E-002
clubb_altitude_threshold = 150.000000000000
clubb_rtp2_clip_coef = 0.500000000000000
clubb_Cx_min = 0.330000000000000
clubb_Cx_max = 0.950000000000000
clubb_Richardson_num_min = 0.250000000000000
clubb_Richardson_num_max = 400.000000000000
clubb_a3_coef_min = 1.60000000000000
clubb_a_const = 1.6000000000000
clubb_c_wp2_pr_dfsn = 0.
clubb_c_wp3_pr_dfsn = 0.
clubb_c_wp3_pr_tp = 0.

!-----CLUBB flags
clubb_l_use_precip_frac = .true.
clubb_l_predict_upwp_vpwp = .true.
clubb_l_min_wp2_from_corr_wx = .false.
clubb_l_min_xp2_from_corr_wx = .false.
clubb_l_C2_cloud_frac = .false.
clubb_l_diffuse_rtm_and_thlm = .false.
clubb_l_stability_correct_Kh_N2_zm = .false.
clubb_l_calc_thlp2_rad = .true.
clubb_l_upwind_xpyp_ta = .true.
clubb_l_upwind_xm_ma = .true.
clubb_l_uv_nudge = .false.
clubb_l_rtm_nudge = .false.
clubb_l_tke_aniso = .true.
clubb_l_vert_avg_closure = .false.
clubb_l_trapezoidal_rule_zt = .true.
clubb_l_trapezoidal_rule_zm = .true.
clubb_l_call_pdf_closure_twice = .true.
clubb_l_standard_term_ta = .false.
clubb_l_partial_upwind_wp3 = .false.
clubb_l_godunov_upwind_wpxp_ta = .false.
clubb_l_godunov_upwind_xpyp_ta = .false.
clubb_l_use_cloud_cover = .false.
clubb_l_diagnose_correlations = .false.
clubb_l_calc_w_corr = .false.
clubb_l_const_Nc_in_cloud = .false.
clubb_l_fix_w_chi_eta_correlations = .true.
clubb_l_stability_correct_tau_zm = .false.
clubb_l_damp_wp2_using_em = .true.
clubb_l_do_expldiff_rtm_thlm = .false.
clubb_l_Lscale_plume_centered = .false.
clubb_l_diag_Lscale_from_tau = .true.
clubb_l_use_C7_Richardson = .true.
clubb_l_use_C11_Richardson = .false.
clubb_l_use_shear_Richardson = .false.
clubb_l_brunt_vaisala_freq_moist = .false.
clubb_l_use_thvm_in_bv_freq = .false.
clubb_l_rcm_supersat_adj = .true.
clubb_l_damp_wp3_Skw_squared = .true.
clubb_l_prescribed_avg_deltaz = .false.
clubb_l_lmm_stepping = .false.
clubb_l_e3sm_config = .false.
clubb_l_vary_convect_depth = .false.
clubb_l_use_tke_in_wp3_pr_turb_term = .false.
clubb_l_use_tke_in_wp2_wp3_K_dfsn = .false.
clubb_l_smooth_Heaviside_tau_wpxp = .false.
clubb_l_enable_relaxed_clipping = .true.
clubb_l_linearize_pbl_winds = .false.
clubb_l_mono_flux_lim_thlm = .true.
clubb_l_mono_flux_lim_rtm = .true.
clubb_l_mono_flux_lim_um = .true.
clubb_l_mono_flux_lim_vm = .true.
clubb_l_mono_flux_lim_spikefix = .true.

user_nl_clm

use_init_interp = .true.
fsurdat = '/glade/work/slevis/git/mksurfdata_toolchain/tools/mksurfdata_esmf/surfdata_ne30np4.pg3_SSP5-8.5_78pfts_CMIP6_1850-2100_c230227.nc'
finidat='/glade/p/cgd/tss/people/oleson/CLM5_restarts/ctsm51_cesm23a14a_ne30pg3ne30pg3mg17_CPLHIST_1850pAD.clm2.r.0561-01-01-00000.nc'

user_nl_cice

f_fresh = 'mxxxx'
f_fresh_ai = 'mxxxx'
ice_ic = '/glade/p/cesmdata/cseg/inputdata/cesm2_init/gmom.e23.GJRAv4.TL319_t061_zstar_N65.tx2_3v2.001/0061-01-01/gmom.e23.GJRAv4.TL319_t061_zstar_N65.tx2_3v2.001.cice.r.0061-01-01-00000.nc'

user_nl_mom

MEKE_GEOMETRIC_ALPHA = 0.10

SourceMods

Looks at the SourceMods in: /glade/p/cesmdata/cseg/runs/cesm2_0/b.e23_alpha16b.BLT1850.ne30_t232.043/SourceMods

src.cam

clubb_intr.F90
hb_diff.F90
mixing_length.F90
vertical_diffusion.F90

src.clm

clm_varpar.F90 
surfrdMod.F90

src.mom

diag_table  
input.nml
dabail10 commented 1 year ago

I have added sea ice diagnostics for this run.

https://webext.cgd.ucar.edu/BLT1850/b.e23_alpha16b.BLT1850.ne30_t232.044/ice/b.e23_alpha16b.BLT1850.ne30_t232.044-obs/yrs1-22/

olyson commented 1 year ago

Land diagnostics:

https://webext.cgd.ucar.edu/BLT1850/b.e23_alpha16b.BLT1850.ne30_t232.044/lnd/b.e23_alpha16b.BLT1850.ne30_t232.044.13_22-b.cesm3_cam058_mom_e.B1850WscMOM.ne30_L58_t061.camdev_cice5.026g.13_22/setsIndex.html

bstephens82 commented 1 year ago

I added atm diagnostics here:

https://webext.cgd.ucar.edu/BLT1850/b.e23_alpha16b.BLT1850.ne30_t232.044/atm/b.e23_alpha16b.BLT1850.ne30_t232.044_1_22_vs_Obs/website/

On the plots it just says years 1-2, which I don't understand since there are 22 years available and I processed all of them, as you can see in the name of the folder above (ie "1_22_vs_Obs"). Maybe that's a bug? Is there a way to know if it's really 2 or 22 years?

vlarson commented 1 year ago

I added atm diagnostics here:

webext.cgd.ucar.edu/BLT1850/b.e23_alpha16b.BLT1850.ne30_t232.044/atm/b.e23_alpha16b.BLT1850.ne30_t232.044_1_22_vs_Obs/website

On the plots it just says years 1-2, which I don't understand since there are 22 years available and I processed all of them, as you can see in the name of the folder above (ie "1_22_vs_Obs"). Maybe that's a bug? Is there a way to know if it's really 2 or 22 years?

Thanks, @bstephens82. Are results posted anywhere from a run that has a similar configuration but uses CLUBB's default Lscale code? Also, I didn't find output of PRECT.

From the comparison with obs, the taus run appears to have near-coastal Sc that are too dim and excessively bright sub-equatorial regions. The bright sub-equatorial regions may be manifesting themselves as a strong double-ITCZ.

bstephens82 commented 1 year ago

@vlarson I'm sure we could run some diagnostics comparing this case to a somewhat similar Lscale case. I will also be uploading diags for the taus run with my non-CLUBB parameters (still running) and we can see if these two taus runs are significantly different as well.

bstephens82 commented 1 year ago

@vlarson here are some additional diagnostics comparing this case (test 044) to test 037, which should be identical apart from the CLUBB setup; i.e. 37 is using the Lscale code, but 044 is using taus along with the other flags/paremeters I got from Zhun, for the most part. I tweaked a few things like gamma and C8. https://webext.cgd.ucar.edu/BLT1850/b.e23_alpha16b.BLT1850.ne30_t232.044/atm/b.e23_alpha16b.BLT1850.ne30_t232.044_1_40_vs_b.e23_alpha16b.BLT1850.ne30_t232.037_1_40/website/

bstephens82 commented 1 year ago

Some diagnostics comparing the first 25 years of this experiment and 043: https://webext.cgd.ucar.edu/BLT1850/b.e23_alpha16b.BLT1850.ne30_t232.044/atm/b.e23_alpha16b.BLT1850.ne30_t232.044_1_25_vs_b.e23_alpha16b.BLT1850.ne30_t232.043_1_25/website/

vlarson commented 12 months ago

Some diagnostics comparing the first 25 years of this experiment and 043: webext.cgd.ucar.edu/BLT1850/b.e23_alpha16b.BLT1850.ne30_t232.044/atm/b.e23_alpha16b.BLT1850.ne30_t232.044_1_25_vs_b.e23_alpha16b.BLT1850.ne30_t232.043_1_25/website

This experiment (044) has brighter (i.e., improved) near-coastal Sc than does (043). It is unclear to me which of the tuning differences between 044 and 043 leads to the improvement, although the tuner suggests that maybe it's the changes in vtrmi or wsub_scale.

olyson commented 11 months ago

Updated land diagnostics here:

https://webext.cgd.ucar.edu/BLT1850/b.e23_alpha16b.BLT1850.ne30_t232.044/lnd/b.e23_alpha16b.BLT1850.ne30_t232.044.25_44-b.cesm3_cam058_mom_e.B1850WscMOM.ne30_L58_t061.camdev_cice5.026g.25_44/setsIndex.html