f.e30_cam6_4_036.FMTHIST.ne120_L93.001

cecilehannay commented 1 month ago

Purpose:

cam7 high resolution run: ne120 L93

Description:

pe_layout
```
./xmlchange NTASKS=6144
```

rrtmgp

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

start in 2000
```
./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', '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', 'AODVISdn','CCN3', 'CDNUMC', 'H2O', 'NUMICE', 'NUMLIQ','OMEGA500', 'TROP_P','TROP_T','TROP_Z', 'num_a1','num_a2','num_a3','num_a4', 'dgnd_a01','dgnd_a02','dgnd_a03','dgnd_a04', 'ncl_a1','ncl_a2','ncl_a3', 'bc_a1','bc_a4', 'dst_a1','dst_a3', 'pom_a1','pom_a4', 'so4_a1','so4_a2','so4_a3', 'soa_a1','soa_a2', 'CCN3','CCN4', 'BUTGWSPEC', 'BVTGWSPEC','BVTGWSPEC','TTGWORO','BTTGWSPEC','TAUGWX', 'TAUGWY','TAUBLJX','TAUBLJY', 'UTEND_PHYSTOT', 'UTEND_GWDTOT', 'UTEND_DCONV', 'UTEND_CLUBB', 'UTEND_CORE', 'VTEND_PHYSTOT', 'VTEND_GWDTOT', 'VTEND_DCONV', 'VTEND_CLUBB', 'VTEND_CORE'

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

!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'

!6 hourly - instantaneous with input from Ben (h8) fincl9 = 'um:I','upwp:I','up2_tp:I','UPWP_UG_CLUBB:I','UPWP_UG_CT_CLUBB:I','DUDZ_CLUBB:I','UPWP_UG_MEAN_LEV:I','UPWP_UG_MIN_LEV:I','UPWP_UG_MAX_LEV:I','UPWP_UG_MEAN_HGT:I' clubb_vars_zt='um' clubb_vars_zm='upwp','up2_tp' clubb_history=.true.

!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'

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' seasalt_emis_scale = 0.75D0

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=576

ncdata = '/glade/campaign/cgd/amp/aherring/restart/inic/ne120/ERA5_x_ne120np4_L93_rgC2_WO.2000-01-01-00000.nc' bnd_topo = '/glade/campaign/cgd/amp/aherring/restart/TopoCESM2/cube_to_target/output/ne120pg3_gmted2010_modis_bedmachine_nc3000_Laplace0025_noleak_20240913.nc'

use_gw_front = .true. taubgnd = 1.5D-6

- **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

flanduse_timeseries = '/glade/derecho/scratch/aherring/restart/ctsm5.2.028/tools/mksurfdata_esmf/landuse.timeseries_ne120pg3_SSP3-7.0_1979-2026_78pfts_c240925.nc' fsurdat = '/glade/derecho/scratch/aherring/restart/ctsm5.2.028/tools/mksurfdata_esmf/surfdata_ne120pg3_SSP3-7.0_1979_78pfts_c240925.nc' finidat = '/glade/campaign/cesm/cesmdata/inputdata/lnd/clm2/initdata_map/clmi.F2000.2000-01-01.ne120pg3_mt13_simyr2000_c200728.nc' use_init_interp = .true. use_excess_ice = .false.


- **SourceMods** 
add Z850, Ben's output

/glade/campaign/cesm/cesmdata/cseg/runs/cesm2_0/f.e30_cam6_4_036.FMTHIST.ne120_L93.001/SourceMods/src.cam/



___
**Case directory:**
- Locally (if still available):
`/glade/campaign/cesm/cesmdata/cseg/runs/cesm2_0/f.e30_cam6_4_036.FMTHIST.ne120_L93.001`

- On github:
https://github.com/NCAR/amwg_dev/tree/f.e30_cam6_4_036.FMTHIST.ne120_L93.001
___
**Sandbox:**
- Locally (if still available):
/glade/work/hannay/cesm_tags/cam6_4_036

- On github:
https://github.com/ESCOMP/CAM/tree/cam6_4_036
___
**Diagnostics:**
- ADF diags (if available)
[https://webext.cgd.ucar.edu/FMTHIST/f.e30_cam6_4_036.FMTHIST.ne120_L93.001/atm/](https://webext.cgd.ucar.edu/FMTHIST/f.e30_cam6_4_036.FMTHIST.ne120_L93.001/atm/)
___
**Output:**
- Raw output on derecho (if still available)
`/glade/derecho/scratch/hannay/archive/f.e30_cam6_4_036.FMTHIST.ne120_L93.001/`
- CAM Climos: `/glade/campaign/cgd/amp/amwg/climo/f.e30_cam6_4_036.FMTHIST.ne120_L93.001`
___
**Contacts:**
@JulioTBacmeister, @adamrher, @islasimpson, @swrneale, @PeterHjortLauritzen, @cmcclus, @bstephens82, @megandevlan, @jedwards4b, @billsack, @wwieder, @ekluzek, @shieldsca, @slevis-lmwg
___
**Extra details:**

cecilehannay commented 1 month ago

From @JulioTBacmeister

First JFM and JJA zonal mean U and T

Memory of the initial conditions here seems remarkable, i.e., our zonal means for the 3-month JFM 2000 look a lot like the ERA JFM 2000 mean.

JJA U shows clearly that summertime stratospheric jet is not closing off. This is also evident in JFM. This is almost certainly caused by the fact that we have shut off the frontal gravity waves. Not too worried about it. More worried about the warming that seems to be taking place at the top of the model. Maybe it's just seasonal, but are we sure have turned on radiation above 1Pa?

@brianpm confirmed the radiaton change went in with cam6_4_016 (https://github.com/ESCOMP/CAM/pull/1100).

@JulioTBacmeister thinks we might be just seeing seasonal variation. He will keep an eye on the model top T.

cecilehannay commented 1 month ago

Note we restarted the simulation clean => this is not relevant anymore.

change in gw parameters during the run

from: 2000-01-01 to 2003-05
```
use_gw_front        = .true.
```

from: 2003-05-01

use_gw_front        = .true.
taubgnd        = 1.5D-6

@JulioTBacmeister

cecilehannay commented 1 month ago

@PeterHjortLauritzen: The run is crashing with the error:

dec2012.hsn.de.hpc.ucar.edu 5720:  ERROR: negative moist layer thickness.  timestep or remap time too large

Should I adjust the sponge?

cecilehannay commented 1 month ago

From: @PeterHjortLauritzen:
We are getting this warning from the estimates Screenshot 2024-10-21 at 11 26 41 AM (1)

can you try se_hypervis_subcycle_sponge=9

cecilehannay commented 1 month ago

I successfully run with se_hypervis_subcycle_sponge=9 for a couple of years but the run crashed again with:

 ERROR: negative moist layer thickness.  timestep or remap time too large

@PeterHjortLauritzen @adamrher : please advice if I can change the sponge further.

adamrher commented 1 month ago

It's crashing at layer 8 with 250+ m/s winds. @PeterHjortLauritzen should we increase se_nsplit or se_rsplit?

k,dp_star_moist,u,v,T    8  -0.015685742200    253.63     91.88    185.27

NCAR / amwg_dev

f.e30_cam6_4_036.FMTHIST.ne120_L93.001 #590

Note we restarted the simulation clean => this is not relevant anymore.

change in gw parameters during the run