f.cam6_3_109.FLTHIST_v0b.ne30.tuningF.001

[cecilehannay]

Description: FLTHIST with f.cam6_3_109 + tuning "F"

compset: FLTHIST_v0b add cosp output use clm5.2 dataset

Tuning:

clubb_l_predict_upwp_vpwp=.true. 
clubb_l_mono_flux_lim_um   = .true.  
clubb_l_mono_flux_lim_vm   = .true.  
clubb_c_uu_shr = 0.1
clubb_c7=0.1

dust_emis_fact         = 0.80D0

micro_mg_dcs=0.00025
micro_mg_vtrmi_factor=0.5
microp_aero_wsubi_scale=2.5
seasalt_emis_scale=1.5

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',
'CCN3', 'CDNUMC'

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

fincl6 =  'LGHTNG_CLD2GRND:A', 'FLASHFRQ:A'

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'

clubb_l_predict_upwp_vpwp=.true. 
clubb_l_mono_flux_lim_um   = .true.  
clubb_l_mono_flux_lim_vm   = .true.  
clubb_c_uu_shr = 0.1
clubb_c7=0.1

dust_emis_fact         = 0.80D0

micro_mg_dcs=0.00025
micro_mg_vtrmi_factor=0.5
microp_aero_wsubi_scale=2.5
seasalt_emis_scale=1.5

user_nl_clm

use_init_interp = .true.
check_finidat_year_consistency = .false.
fsurdat = '/glade/work/slevis/git/mksurfdata_toolchain/tools/mksurfdata_esmf/surfdata_ne30np4.pg3_SSP5-8.5_78pfts_CMIP6_1850-2100_c230227.nc'
flanduse_timeseries = '/glade/work/slevis/git/mksurfdata_toolchain/tools/mksurfdata_esmf/landuse.timeseries_ne30np4.pg3_SSP5-8.5_78_CMIP6_1850-2100_c230227.nc'

Case directory: Locally (if still available): /glade/p/cesmdata/cseg/runs/cesm2_0/f.cam6_3_109.FLTHIST_v0b.ne30.tuningF.001

On github: https://github.com/NCAR/amwg_dev/tree/f.cam6_3_109.FLTHIST_v0b.ne30.tuningF.001

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

On github: hash: 962921a (https://github.com/ESCOMP/CAM/tree/cam6_3_109)

Diagnostics: ADF diags (if available) https://webext.cgd.ucar.edu/FLTHIST/f.cam6_3_109.FLTHIST_v0b.ne30.tuningF.001/atm/

Contacts:

[cecilehannay]

I am getting an error:

(GETFIL): using /glade/work/slevis/git/mksurfdata_toolchain/tools/mksurfdata_esmf/surfdata_ne30np4.pg3_SSP5-8.5_78pfts_CMIP6_1850-2100_c230227.nc
 UrbanInput                      
 /glade/work/slevis/git/mksurfdata_toolchain/tools/mksurfdata_esmf/surfdata_ne30
 np4.pg3_SSP5-8.5_78pfts_CMIP6_1850-2100_c230227.nc
 UrbanInput: parameter nlevurb=            5 
 does not equal input dataset nlevurb=           10
 ENDRUN:
 ERROR: ERROR in UrbanParamsType.F90 at line 476

@olyson

[adamrher]

DId you add the clm source mods from https://github.com/NCAR/amwg_dev/issues/250? You can omit the SnowHydrology.F90 one, but the rest you need for the new surface dataset your pointing to.

[cecilehannay]

Yes, the user_nl_clm were included.

[olyson]

@adamrher is correct about needing to include the SourceMods for the surface dataset (clm_varpar.F90 and surfrdMod.F90). However, I see that cam6_3_109 has ctsm5.1.dev120 as its external, which does not have the robust snow hydrology fix (which came in with ctsm5.1.dev122), so you would need to use the SnowHydrologyMod.F90 sourcemod. Unless I'm mistaken about the CTSM external in cam6_3_109.

[adamrher]

you're correct, we're still at dev120 so we need the SnowHydrology.F90 mods. Cecile the three modules you need are here:

/glade/u/home/oleson/run_hist_1850_files/SPINUP/casefiles_cecile/SourceMods/ctsm5.1.dev120/src.clm/

[cecilehannay]

Sorry, I thought this went on the trunk when CLM5.1 so I had removed those.

On Fri, May 5, 2023 at 4:25 PM Adam Herrington @.***> wrote:

you're correct, we're still at dev120 so we need the SnowHydrology.F90 mods. Cecile the three modules you need are here:

/glade/u/home/oleson/run_hist_1850_files/SPINUP/casefiles_cecile/SourceMods/ctsm5.1.dev120/src.clm/

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

[olyson]

@cecilehannay , can you include FLASHFREQ in that 2-hourly stream you are generating for fincl6.

[edit] Sorry, that should be FLASHFRQ.

[adamrher]

@cecilehannay , can you include FLASHFREQ in that 2-hourly stream you are generating for fincl6.

also, add :A to both variables, or set the avgflag_pertape(6) = 'A'

[JulioTBacmeister]

@dan800 This is an FLTHIST case. What can we transfer to BLT1850?

[cecilehannay]

This run completed. @justin-richling: could you run the ADF for

f.cam6_3_109.FLTHIST_v0b.ne30.tuningF.001

• versus f.cam6_3_107.FLTHIST_v0a.ne30.clm5_1_ebudget.001
• versus observations

[justin-richling]

This run completed. @justin-richling: could you run the ADF for

f.cam6_3_109.FLTHIST_v0b.ne30.tuningF.001

• versus f.cam6_3_107.FLTHIST_v0a.ne30.clm5_1_ebudget.001
• versus observations

On it now

[justin-richling]

ADF diags vs f.cam6_3_107.FLTHIST_v0a.ne30.clm5_1_ebudget.001 here

ADF vs Obs here

[dan800]

I looked at the ADF vs obs.

• There's a 10 hPa surface pressure bias over the pole in NH winter. Was that better in the run with non-orographic GWD turned on?
• There is no QBO of any form is this run
• The U, T zonal mean plots on logP were not made, so hard to gauge the stratosphere climatology.

[justin-richling]

• The U, T zonal mean plots on logP were not made, so hard to gauge the stratosphere climatology.

U should be in the log_p plots, but it looks like T isn't because it doesn't have an observation file in the variable defaults yaml. I can re-run it again real quick and add ERAI_all_climo.nc for T if that is a reasonable obs file for it.

[dan800]

• The U, T zonal mean plots on logP were not made, so hard to gauge the stratosphere climatology.

U should be in the log_p plots, but it looks like T isn't because it doesn't have an observation file in the variable defaults yaml. I can re-run it again real quick and add ERAI_all_climo.nc for T if that is a reasonable obs file for it.

OK, found it - for some reason I didn't see them before. I do worry how the extension of the strong strat bias in U extending to surface in SON impacts ocean circulation and sea ice extent. Also, if we want to keep CAM-CHEM consistent, this looks like it could be a late breakup of the winter vortex which means the ozone hole will be way off.

[adamrher]

There's a 10 hPa surface pressure bias over the pole in NH winter. Was that better in the run with non-orographic GWD turned on?

It does look a bit low compared to obs here, if that's what you're referring to. Not much difference compared to the non tuned control.

I do recall you showed a ~10 hPa increase in PS from non-oro settings here here. So it might be useful to turn them on.

[dan800]

It does look a bit low compared to obs here, if that's what you're referring to.

Yes.

I do recall you showed a ~10 hPa increase in PS from non-oro settings here here. So it might be useful to turn them on.

I couldn't recall the sign of the change - it does look like it's in the right direction to help alleviate this PS bias. And it would tame that stratospheric jet bias as well.

[tilmes]

@JulioTBacmeister Do you have ideas for improving this?

[JulioTBacmeister]

@JulioTBacmeister Do you have ideas for improving this?

Yes, a couple. But I don't think we should wait to start the coupled testing.
I'll submit two run requests (FLTHIST) when I have some time: 1) frontal and convective turned on in FLT (I think they're still off by default @cecilehannay ?) 2) A run with effgw_rdg==> 2.0 but with low level drag reduced. I'll put the setttings in request.

We also should start an FMTHIST run soon.

[tilmes]

@JulioTBacmeister @cecilehannay I would be happy to help run some of the runs, including the mid-top version. What variables have to be changed from the LT version for 1. and 2.?

[JulioTBacmeister]

@islasimpson This is the FLTHIST control for the other runs