search

NCAR / amwg_dev

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

b.e23_alpha17f.BLTHIST.ne30_t232.092 #543

Open cecilehannay opened 6 months ago

cecilehannay commented 6 months ago

Purpose: First attempt to run a 20th century with current cesm3 development version (alpha17f)

Description:

Use alpha17f with HYCOM1. Starting from b.e23_alpha17f.BLT1850.ne30_t232.092 at year 101.

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', 'AODVISdn','CCN3', 'CDNUMC', 'H2O', 'NUMICE', 'NUMLIQ','OMEGA500', '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', 'SFSO2', 'OCN_FLUX_DMS', 'SAD_SULFC', 'SAD_TROP', 'SAD_AERO'

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

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


- **user_nl_clm**

use_init_interp = .true.


- **user_nl_cice**

dumpfreq = 'x','x','x','x','x' f_fresh = 'mxxxx' f_fresh_ai = 'mxxxx'


- **user_nl_mom**

MIN_SALINITY = 1.e-06 INPUTDIR = "INPUT" MAX_LAYER_THICKNESS_CONFIG = "FILE:dz_max_90th_quantile.nc,dz"

- **`xmlchange` changes**

./xmlchange MOM6_VERTICAL_GRID=hycom1 ./xmlchange RUN_TYPE=hybrid ./xmlchange RUN_REFCASE=b.e23_alpha17f.BLT1850.ne30_t232.092 ./xmlchange RUN_REFDATE=0101-01-01 ./xmlchange GET_REFCASE=TRUE ./xmlchange RUN_REFDIR=cesm2_init


___
**Case directory:**
- Locally (if still available):
`/glade/campaign/cesm/cesmdata/cseg/runs/cesm2_0/b.e23_alpha17f.BLTHIST.ne30_t232.092`

- On github:
https://github.com/NCAR/amwg_dev/tree/b.e23_alpha17f.BLTHIST.ne30_t232.092
___
**Sandbox:**
- Locally (if still available): `/glade/work/hannay/cesm_tags/cesm2_3_alpha17f`

- On github: https://github.com/ESCOMP/cesm/tree/cesm2_3_alpha17f
___
**Diagnostics:**
- Diags (if available):
[https://webext.cgd.ucar.edu/BLTHIST/b.e23_alpha17f.BLTHIST.ne30_t232.092/](https://webext.cgd.ucar.edu/BLTHIST/b.e23_alpha17f.BLTHIST.ne30_t232.092/)
___
**Output:**
- Raw output on derecho (if still available): `/glade/derecho/scratch/hannay/archive/b.e23_alpha17f.BLTHIST.ne30_t232.092/`
- CAM Climos: `/glade/campaign/cgd/amp/amwg/climo/b.e23_alpha17f.BLTHIST.ne30_t232.092`
___
**Contacts:**
@cecilehannay @dlawrenncar 
@PeterHjortLauritzen, @adamrher, @brianpm, @JulioTBacmeister, @swrneale, 
@gustavo-marques, @fobryan3, @klindsay28,
@wwieder, @slevis-lmwg, @olyson
___
**Extra details:**
cecilehannay commented 5 months ago

@fvitt: We are starting to start a 20th century. Right now, we get out of the box:

ubc_file_path       = '/glade/campaign/cesm/cesmdata/inputdata/atm/cam/chem/ubc/f.e21.FWHISTBgcCrop.f09_f09_mg17.CMIP6-AMIP-WACCM.ensAvg123.cam.h0zm.UBC.195001-201412_c220322.nc'

which doesn't include the time period we need.

We need a ubc_file_path starting in 1850. I see over files in the same directory but I have no clue there are ok to use.

Please advice. Thanks.

klindsay28 commented 5 months ago

Could you add USE_CFC_CAP = True to user_nl_mom?

On Sat, Jun 1, 2024 at 10:13 AM Cecile Hannay @.***> wrote:

@fvitt https://github.com/fvitt: We are starting to start a 20th century. Right now, we get out of the box:

ubc_file_path = '/glade/campaign/cesm/cesmdata/inputdata/atm/cam/chem/ubc/f.e21.FWHISTBgcCrop.f09_f09_mg17.CMIP6-AMIP-WACCM.ensAvg123.cam.h0zm.UBC.195001-201412_c220322.nc'

which doesn't include the time period we need.

We need a ubc_file_path starting in 1850. I see over files in the same directory but I have no clue there are ok to use.

Please advice. Thanks.

— Reply to this email directly, view it on GitHub https://github.com/NCAR/amwg_dev/issues/543#issuecomment-2143501315, or unsubscribe https://github.com/notifications/unsubscribe-auth/ADWZPO3AKEVR26TMPHZOAOTZFHXL7AVCNFSM6AAAAABITT6IU2VHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDCNBTGUYDCMZRGU . You are receiving this because you were mentioned.Message ID: @.***>

fvitt commented 5 months ago

@cecilehannay If all you need is H2O at the upper boundary then use:

 /glade/campaign/cesm/cesmdata/inputdata/atm/cam/chem/ubc/b.e21.BWHIST.f09_g17.CMIP6-historical-WACCM.ensAvg123.cam.h0zm.H2O.185001-201412_c230509cdf5.nc
cecilehannay commented 5 months ago

@fvitt and klindsay28 I will try to restart with these settings when the machine comes back up.

cecilehannay commented 5 months ago

@fvitt: I started the historical with teh file you recommended:

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

The model is still crashing, and I think the issue is that the file starts on 18500201. As a result, the model doesn't find the needed data to start on 18500101.

Could you add a month to your file so we can run the historical data? Thanks a bunch.

fvitt commented 5 months ago

@cecilehannay

Try this UBC file:

/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

Here I prepended 1850 as 1849 to the file.

wwieder commented 5 months ago

@olyson and @slevis-lmwg, just a heads up that there are historical cases from 092 and 098 that are running. When they are completed can you run diagnostics on these simulations? 092 is pretty far along (1960s), maybe we can compare it to a CESM2-hist case next week?

olyson commented 5 months ago

Yes, I have 092hist compared to LE2-001, through 1976, here:

https://webext.cgd.ucar.edu/BLTHIST/b.e23_alpha17f.BLTHIST.ne30_t232.092/lnd/b.e23_alpha17f.BLTHIST.ne30_t232.092_1957_1976-b.e21.BHISTcmip6.f09_g17.LE2-1001.001.1995-2014/setsIndex.html

The climatology is 1957-1976 for 092hist and 1995-2015 for LE2-001, so the trend plots are probably more useful at this point. I'll update this once 092hist is complete. 098hist is at 1871.

wwieder commented 5 months ago

You're amazing @olyson. thanks