f.cam6_3_106.FLTHIST_v0a.ne30.dcs_effgw_rdg.001

cecilehannay commented 1 year ago

Description:

Same as https://github.com/NCAR/amwg_dev/issues/235:

start the run in 1995

dcs = 500->250 microns

micro_mg_dcs            =           250.D-6

gw tuning:

use_gw_convect_dp    = .false.
use_gw_front         = .false.
effgw_rdg_beta       = 2.0D0
effgw_rdg_beta_max   = 2.0D0

Case directory: Locally (if still available): /glade/p/cesmdata/cseg/runs/cesm2_0/f.cam6_3_106.FLTHIST_v0a.ne30.dcs_effgw_rdg.001

On github: https://github.com/NCAR/amwg_dev/tree/$CASE

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

On github: hash: 0470324 (https://github.com/ESCOMP/CAM/tree/cam6_3_106)

Diagnostics: ADF diags (if available) https://webext.cgd.ucar.edu/FLTHIST/f.cam6_3_106.FLTHIST_v0a.ne30.dcs_effgw_rdg.001/atm/

Contacts: @cacraigucar, @cecilehannay, @julio, @adamrher, @brianpm, @dan800, @tilmes, @swrneale, @bstephens82

cecilehannay commented 1 year ago

Here is the complete user_nl_cam

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

fexcl1 = ' '  

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'

fincl2 = 'OMEGA', 'PMID', 'PS', 'Q', 'QRL', 'QRS', 'T', 'TROP_P', 'TROP_T', 'U', 'V', 'Z3'

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

fincl4 =  'PRECC','PRECL'

fincl5 = 'Uzm','Vzm','Wzm','THzm', 'VTHzm','WTHzm','UVzm','UWzm'

fincl6 =  'WV_phBF','WL_phBF','WI_phBF','SE_phBF','KE_phBF',
'WV_phBP','WL_phBP','WI_phBP','SE_phBP','KE_phBP',
'WV_phAP','WL_phAP','WI_phAP','SE_phAP','KE_phAP',
'WV_phAM','WL_phAM','WI_phAM','SE_phAM','KE_phAM',
'WV_dyBF','WL_dyBF','WI_dyBF','SE_dyBF','KE_dyBF',
'WV_dyBP','WL_dyBP','WI_dyBP','SE_dyBP','KE_dyBP',
'WV_dyAP','WL_dyAP','WI_dyAP','SE_dyAP','KE_dyAP',
'WV_dyAM','WL_dyAM','WI_dyAM','SE_dyAM','KE_dyAM',
'WV_dBF','WL_dBF','WI_dBF','SE_dBF','KE_dBF',
'WV_dED','WL_dED','WI_dED','SE_dED','KE_dED',
'WV_dAD','WL_dAD','WI_dAD','SE_dAD','KE_dAD',
'WV_dAF','WL_dAF','WI_dAF','SE_dAF','KE_dAF',
'WV_dBD','WL_dBD','WI_dBD','SE_dBD','KE_dBD',
'WV_dAR','WL_dAR','WI_dAR','SE_dAR','KE_dAR',
'WV_dAH','WL_dAH','WI_dAH','SE_dAH','KE_dAH',
'WV_dBH','WL_dBH','WI_dBH','SE_dBH','KE_dBH',
'WV_dCH','WL_dCH','WI_dCH','SE_dCH','KE_dCH',
'WV_dAS','WL_dAS','WI_dAS','SE_dAS','KE_dAS',
'WV_dBS','WL_dBS','WI_dBS','SE_dBS','KE_dBS',
'EFIX'

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_a2', 
'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'

phys_grid_ctem_nfreq=-6
phys_grid_ctem_zm_nbas=120
phys_grid_ctem_za_nlat=90

ncdata = '/glade/p/cesm/amwg_dev/juliob/FWsc_ne30pg3_58L_GRID_48_taperstart10km_lowtop_BL10_v3_beta1p75_Top_43km.nc'

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            =           250.D-6
use_gw_convect_dp        = .false.
use_gw_front        = .false.
effgw_rdg_beta        = 2.0D0
effgw_rdg_beta_max        = 2.0D0

cecilehannay commented 1 year ago

LWCF is very low with Dcs = 250 microns Screen Shot 2023-04-11 at 1 57 13 PM

brianpm commented 1 year ago

The total cloud cover is fairly similar to CAM6 (the old, CMIP6 version of CAM6). At least based on CALIPSO observations and simulator: Global averages:

OBS: 0.6705907397435672,
CAM6: 0.573204437986628,
DEV: 0.5615052033358574 Near-global (65S-65N) averages:
OBS: 0.6699508555382104,
CAM6: 0.5635547349994824,
DEV: 0.5629364973276819 Bias (model minus obs):
REF: 0.10639612053872784,
DEV: 0.10701435821052874 The normalized mean square error components w.r.t. observed climo: CAM6: (0.525, 0.329, 0.218, 2.28) DEV: (0.532, 0.350, 0.195, 2.38) Where the values are the unconditional bias (basically mean bias), the conditional bias, the phase error (basically pattern correlation), and the scaled variance ratio. They are quite close together, but this simulation has a slightly better spatial pattern than the default model.

That said, the actual errors in cloud cover remain very large.

brianpm commented 1 year ago

Comparison with CALIPSO total, high, middle, low cloud amounts.

And here's the GCSS cross section, this shows how the stratus deck breaks up too soon and the shallow cumulus tend to be too low. Based on annual average, so take with a grain of salt, but I think is indicative of general bias in low clouds.

adamrher commented 1 year ago

thanks @brianpm! I had previously shown that the L32->L58 improves the cloud fraction fields in the GCSS transect -- at least Marine Sc side of things. Are you able to plot L32 CAM6 transect as well?

brianpm commented 1 year ago

Yes. (edit: improved by making colors the same / edit 2: one colorbar, better labeling)

swrneale commented 1 year ago

Wow! Stop everything now.

adamrher commented 1 year ago

When I've looked at this before, the breakup into Cu clouds as you go south is primarily instigated by ZM. I'd argue it's an uphill battle to tame ZM in this regime, and that maybe an increase in dp1 and dp2 could recover wider Cu clouds?

I'd be curious to see Trude's PPE plots showing the sensitivity to dp1 and dp2.

justin-richling commented 1 year ago

Updated ADF diags vs f.cam6_3_105.FLTHIST_v0a.ne30.001 are here

Zonal log-p plots are under Log-P category FYI

NCAR / amwg_dev

f.cam6_3_106.FLTHIST_v0a.ne30.dcs_effgw_rdg.001 #239