Skip to content

Commit

Permalink
Ignore FV3 restart checksum. Set stochastic parameters. Errors in goc…
Browse files Browse the repository at this point in the history 
…art input data
  • Loading branch information
@NickSzapiro-NOAA
NickSzapiro-NOAA committed Aug 12, 2024
1 parent 17f243b commit e6d4338
Show file tree
Hide file tree
Showing 4 changed files with 53 additions and 49 deletions.
7 changes: 4 additions & 3 deletions tests/fv3_conf/cpld_control_run.IN
View file
Original file line number Diff line number Diff line change
Expand Up @@ -169,13 +169,14 @@ else
#cp ${FV3_IC}/coupler.res ./INPUT
#coupler.res fv_core.res.nc fv_core.res.tile[1-6].nc fv_srf_wnd.res.tile[1-6].nc fv_tracer.res.tile[1-6].nc phy_data.tile[1-6].nc fv3_increment.nc
cp ${MOM_IC}/MOM*.nc ./INPUT
cp ${MOM_IC}/mom6_increment.nc ./INPUT
#cp ${MOM_IC}/ORAS5.mx025.ic.nc ./INPUT
cp ${MOM_IC}/interpolate_zgrid_26L.nc ./INPUT
cp ${ICE_IC}/cice*nc ./INPUT
ls -1 ./INPUT/cice*nc>ice.restart_file
cp ${FV3_IC}/ufs.cpld.cpl*nc .
ls -1 ufs.cpld.cpl*nc>rpointer.cpl
cp ${WAV_IC}/restart.ww3 .
cp ${WAV_IC}/restart.ww3 ufs.cpld.ww3.r.${RESTART_FILE_SUFFIX_SECS}
else
# C48mx500 warm start from ICs
cp ${FV3_IC}/*.nc ./INPUT
Expand Down Expand Up @@ -223,7 +224,7 @@ if [ $IAER = 2011 ]; then
cp @[INPUTDATA_ROOT]/FV3_input_data_INCCN_aeroclim/aer_data/LUTS/optics_DU.v15_3.dat optics_DU.dat
cp @[INPUTDATA_ROOT]/FV3_input_data_INCCN_aeroclim/aer_data/LUTS/optics_SS.v3_3.dat optics_SS.dat
cp @[INPUTDATA_ROOT]/FV3_input_data_INCCN_aeroclim/aer_data/LUTS/optics_SU.v1_3.dat optics_SU.dat
DIR_OPT=/scratch1/NCEPDEV/nems/Nick.Szapiro/tasks/input_data/gefs.v13/CHM/ExtData/ExtData/monochromatic/
DIR_OPT=/scratch1/NCEPDEV/nems/Nick.Szapiro/tasks/input_data/gefs.v13/CHM/ExtData/monochromatic/
cp ${DIR_OPT}/optics_BC.v1_3.nc optics_BC.nc
cp ${DIR_OPT}/optics_OC.v1_3.nc optics_OC.nc
cp ${DIR_OPT}/optics_DU.v15_3.nc optics_DU.nc
Expand All @@ -247,7 +248,7 @@ fi
#prognostic aerosols
if [ $CPLCHM = .true. ]; then
if [[ $BMIC == .true. ]]; then
ln -sf /scratch1/NCEPDEV/nems/Nick.Szapiro/tasks/input_data/gefs.v13/CHM/* .
ln -sf /scratch1/NCEPDEV/nems/Nick.Szapiro/tasks/input_data/gefs.v13/CHM/ExtData .
else
ln -sf @[INPUTDATA_ROOT]/GOCART/p8c_5d/ExtData .
fi
Expand Down
3 changes: 2 additions & 1 deletion tests/parm/global_control.nml.IN
View file
Original file line number Diff line number Diff line change
Expand Up @@ -94,6 +94,7 @@
agrid_vel_rst = .true.
read_increment = @[READ_INCREMENT]
res_latlon_dynamics = @[RES_LATLON_DYNAMICS]
ignore_rst_cksum = .true.
/

&external_ic_nml
Expand Down Expand Up @@ -398,7 +399,7 @@
SHUM_LSCALE=500000,
SHUMINT=3600,
ISEED_SHUM=1,
SPPT=0.56,0.28,0.14,0.056,0.028
SPPT=@[SPPT]
SPPT_TAU=2.16E4,2.592E5,2.592E6,7.776E6,3.1536E7
SPPT_LSCALE=500.E3,1000.E3,2000.E3,2000.E3,2000.E3
SPPT_LOGIT=.TRUE.,
Expand Down
86 changes: 43 additions & 43 deletions tests/parm/gocart/AERO_ExtData.rc
View file
Original file line number Diff line number Diff line change
Expand Up @@ -9,74 +9,74 @@ TROPP 'Pa' Y N - 0.0 1.0

#====== Dust Imports =================================================
# FENGSHA input files. Note: regridding should be N or E - Use files with _FillValue != NaN
DU_CLAY '1' Y E - none none clayfrac FENGSHA_2022_NESDIS_inputs_10km_v3.2.nc
DU_SAND '1' Y E - none none sandfrac FENGSHA_2022_NESDIS_inputs_10km_v3.2.nc
DU_CLAY '1' Y E - none none clayfrac ExtData/FENGSHA_2022_NESDIS_inputs_10km_v3.2.nc
DU_SAND '1' Y E - none none sandfrac ExtData/FENGSHA_2022_NESDIS_inputs_10km_v3.2.nc
DU_SILT '1' Y E - none none siltfrac /dev/null
DU_SSM '1' Y E - none none sep FENGSHA_2022_NESDIS_inputs_10km_v3.2.nc
DU_RDRAG '1' Y E %y4-%m2-%d2t12:00:00 none none albedo_drag FENGSHA_2022_NESDIS_inputs_10km_v3.2.nc
DU_UTHRES '1' Y E - none none uthres FENGSHA_2022_NESDIS_inputs_10km_v3.2.nc
DU_SSM '1' Y E - none none sep ExtData/FENGSHA_2022_NESDIS_inputs_10km_v3.2.nc
DU_RDRAG '1' Y E %y4-%m2-%d2t12:00:00 none none albedo_drag ExtData/FENGSHA_2022_NESDIS_inputs_10km_v3.2.nc
DU_UTHRES '1' Y E - none none uthres ExtData/FENGSHA_2022_NESDIS_inputs_10km_v3.2.nc

#====== Sulfate Sources =================================================
# Anthropogenic (BF & FF) emissions -- allowed to input as two layers
SU_ANTHROL1 NA N Y %y4-%m2-%d2t12:00:00 none none SO2 CEDS.2019.emis.20210325.nc
SU_ANTHROL2 NA N Y %y4-%m2-%d2t12:00:00 none none SO2_elev CEDS.2019.emis.20210325.nc
SU_ANTHROL1 NA N Y %y4-%m2-%d2t12:00:00 none none SO2 ExtData/CEDS.2019.emis.20210325.nc
SU_ANTHROL2 NA N Y %y4-%m2-%d2t12:00:00 none none SO2_elev ExtData/CEDS.2019.emis.20210325.nc

# Ship emissions
SU_SHIPSO2 NA N Y %y4-%m2-%d2t12:00:00 none none SO2_ship CEDS.2019.emis.20210325.nc
SU_SHIPSO4 NA N Y %y4-%m2-%d2t12:00:00 none none SO4_ship CEDS.2019.emis.20210325.nc
SU_SHIPSO2 NA N Y %y4-%m2-%d2t12:00:00 none none SO2_ship ExtData/CEDS.2019.emis.20210325.nc
SU_SHIPSO4 NA N Y %y4-%m2-%d2t12:00:00 none none SO4_ship ExtData/CEDS.2019.emis.20210325.nc

# Aircraft fuel consumption
SU_AIRCRAFT NA Y Y %y4-%m2-%d2t12:00:00 none none none /dev/null

# DMS concentration
SU_DMSO NA Y Y %y4-%m2-%d2t12:00:00 none none conc DMSclim_sfcconcentration.x360_y181_t12.Lana2011.nc4
SU_DMSO NA Y Y %y4-%m2-%d2t12:00:00 none none conc ExtData/DMSclim_sfcconcentration.x360_y181_t12.Lana2011.nc4

# Aviation emissions during the three phases of flight
SU_AVIATION_LTO NA Y Y %y4-%m2-%d2t12:00:00 none none so2_aviation htap-v2.2.emis_so2.aviation_lto.x3600_y1800_t12.2010.nc4
SU_AVIATION_CDS NA Y Y %y4-%m2-%d2t12:00:00 none none so2_aviation htap-v2.2.emis_so2.aviation_cds.x3600_y1800_t12.2010.nc4
SU_AVIATION_CRS NA Y Y %y4-%m2-%d2t12:00:00 none none so2_aviation htap-v2.2.emis_so2.aviation_crs.x3600_y1800_t12.2010.nc4
SU_AVIATION_LTO NA Y Y %y4-%m2-%d2t12:00:00 none none so2_aviation ExtData/htap-v2.2.emis_so2.aviation_lto.x3600_y1800_t12.2010.nc4
SU_AVIATION_CDS NA Y Y %y4-%m2-%d2t12:00:00 none none so2_aviation ExtData/htap-v2.2.emis_so2.aviation_cds.x3600_y1800_t12.2010.nc4
SU_AVIATION_CRS NA Y Y %y4-%m2-%d2t12:00:00 none none so2_aviation ExtData/htap-v2.2.emis_so2.aviation_crs.x3600_y1800_t12.2010.nc4

# H2O2, OH and NO3 mixing ratios
# --------------------------------------------------------------
# If using 64 levels please replace this section with the correct values (ie replace 127 with 64)

SU_H2O2 NA N Y %y4-%m2-%d2t12:00:00 none none h2o2 A2_ACCMIP_gmic_MERRA_oh_h2o2_no3.x144_y91_z127_t14.2021.nc
SU_OH NA N Y %y4-%m2-%d2t12:00:00 none none oh A2_ACCMIP_gmic_MERRA_oh_h2o2_no3.x144_y91_z127_t14.2021.nc
SU_NO3 NA N Y %y4-%m2-%d2t12:00:00 none none no3 A2_ACCMIP_gmic_MERRA_oh_h2o2_no3.x144_y91_z127_t14.2021.nc
SU_H2O2 NA N Y %y4-%m2-%d2t12:00:00 none none h2o2 ExtData/A2_ACCMIP_gmic_MERRA_oh_h2o2_no3.x144_y91_z127_t14.2021.nc
SU_OH NA N Y %y4-%m2-%d2t12:00:00 none none oh ExtData/A2_ACCMIP_gmic_MERRA_oh_h2o2_no3.x144_y91_z127_t14.2021.nc
SU_NO3 NA N Y %y4-%m2-%d2t12:00:00 none none no3 ExtData/A2_ACCMIP_gmic_MERRA_oh_h2o2_no3.x144_y91_z127_t14.2021.nc
#---------------------------------------------------------------

# Production of SO2 from OCS oxidation
pSO2_OCS NA Y Y %y4-%m2-%d2t12:00:00 none none biofuel /dev/null

#SU_regionMask NA N v - none none REGION_MASK ARCTAS.region_mask.x540_y361.2008.nc
#SU_regionMask NA N v - none none REGION_MASK ExtData/ARCTAS.region_mask.x540_y361.2008.nc

#=========== Carbonaceous aerosol sources ===========================================
# ORGANIC CARBON
# ---------------

# # VOCs - OFFLINE MEGAN BIOG
OC_ISOPRENE NA N Y %y4-%m2-%d2t12:00:00 none none isoprene MEGAN.OFFLINE.BIOVOC.2021.emis.20210325.nc
OC_LIMO NA N Y %y4-%m2-%d2t12:00:00 none none limo MEGAN.OFFLINE.BIOVOC.2021.emis.20210325.nc
OC_MTPA NA N Y %y4-%m2-%d2t12:00:00 none none mtpa MEGAN.OFFLINE.BIOVOC.2021.emis.20210325.nc
OC_MTPO NA N Y %y4-%m2-%d2t12:00:00 none none mtpo MEGAN.OFFLINE.BIOVOC.2021.emis.20210325.nc
OC_ISOPRENE NA N Y %y4-%m2-%d2t12:00:00 none none isoprene ExtData/MEGAN.OFFLINE.BIOVOC.2021.emis.20210325.nc
OC_LIMO NA N Y %y4-%m2-%d2t12:00:00 none none limo ExtData/MEGAN.OFFLINE.BIOVOC.2021.emis.20210325.nc
OC_MTPA NA N Y %y4-%m2-%d2t12:00:00 none none mtpa ExtData/MEGAN.OFFLINE.BIOVOC.2021.emis.20210325.nc
OC_MTPO NA N Y %y4-%m2-%d2t12:00:00 none none mtpo ExtData/MEGAN.OFFLINE.BIOVOC.2021.emis.20210325.nc

# Biofuel Source -- Included in AeroCom anthropogenic emissions
OC_BIOFUEL NA Y Y %y4-%m2-%d2t12:00:00 none none biofuel /dev/null

# Anthropogenic (BF & FF) emissions -- allowed to input as two layers
OC_ANTEOC1 NA N Y %y4-%m2-%d2t12:00:00 none none OC CEDS.2019.emis.20210325.nc
OC_ANTEOC2 NA N Y %y4-%m2-%d2t12:00:00 none none OC_elev CEDS.2019.emis.20210325.nc
OC_ANTEOC1 NA N Y %y4-%m2-%d2t12:00:00 none none OC ExtData/CEDS.2019.emis.20210325.nc
OC_ANTEOC2 NA N Y %y4-%m2-%d2t12:00:00 none none OC_elev ExtData/CEDS.2019.emis.20210325.nc

# EDGAR based ship emissions
OC_SHIP NA N Y %y4-%m2-%d2t12:00:00 none none OC_ship CEDS.2019.emis.20210325.nc
OC_SHIP NA N Y %y4-%m2-%d2t12:00:00 none none OC_ship ExtData/CEDS.2019.emis.20210325.nc

# Aircraft fuel consumption
OC_AIRCRAFT NA N Y %y4-%m2-%d2t12:00:00 none none oc_aviation /dev/null

# Aviation emissions during the three phases of flight
OC_AVIATION_LTO NA Y Y %y4-%m2-%d2t12:00:00 none none oc_aviation htap-v2.2.emis_oc.aviation_lto.x3600_y1800_t12.2010.nc4
OC_AVIATION_CDS NA Y Y %y4-%m2-%d2t12:00:00 none none oc_aviation htap-v2.2.emis_oc.aviation_cds.x3600_y1800_t12.2010.nc4
OC_AVIATION_CRS NA Y Y %y4-%m2-%d2t12:00:00 none none oc_aviation htap-v2.2.emis_oc.aviation_crs.x3600_y1800_t12.2010.nc4
OC_AVIATION_LTO NA Y Y %y4-%m2-%d2t12:00:00 none none oc_aviation ExtData/htap-v2.2.emis_oc.aviation_lto.x3600_y1800_t12.2010.nc4
OC_AVIATION_CDS NA Y Y %y4-%m2-%d2t12:00:00 none none oc_aviation ExtData/htap-v2.2.emis_oc.aviation_cds.x3600_y1800_t12.2010.nc4
OC_AVIATION_CRS NA Y Y %y4-%m2-%d2t12:00:00 none none oc_aviation ExtData/htap-v2.2.emis_oc.aviation_crs.x3600_y1800_t12.2010.nc4

# SOA production
pSOA_ANTHRO_VOC NA Y Y %y4-%m2-%d2t12:00:00 none none biofuel /dev/null
Expand All @@ -89,19 +89,19 @@ pSOA_ANTHRO_VOC NA Y Y %y4-%m2-%d2t12:00:00 none none biofuel /dev/null
BC_BIOFUEL NA Y Y %y4-%m2-%d2t12:00:00 none none biofuel /dev/null

# Anthropogenic (BF & FF) emissions -- allowed to input as two layers
BC_ANTEBC1 NA N Y %y4-%m2-%d2t12:00:00 none none BC CEDS.2019.emis.20210325.nc
BC_ANTEBC2 NA N Y %y4-%m2-%d2t12:00:00 none none BC_elev CEDS.2019.emis.20210325.nc
BC_ANTEBC1 NA N Y %y4-%m2-%d2t12:00:00 none none BC ExtData/CEDS.2019.emis.20210325.nc
BC_ANTEBC2 NA N Y %y4-%m2-%d2t12:00:00 none none BC_elev ExtData/CEDS.2019.emis.20210325.nc

# EDGAR based ship emissions
BC_SHIP NA N Y %y4-%m2-%d2t12:00:00 none none BC_ship CEDS.2019.emis.20210325.nc
BC_SHIP NA N Y %y4-%m2-%d2t12:00:00 none none BC_ship ExtData/CEDS.2019.emis.20210325.nc

# Aircraft fuel consumption
BC_AIRCRAFT NA N Y %y4-%m2-%d2t12:00:00 none none bc_aviation /dev/null

# Aviation emissions during the LTO, SDC and CRS phases of flight
BC_AVIATION_LTO NA Y Y %y4-%m2-%d2t12:00:00 none none bc_aviation htap-v2.2.emis_bc.aviation_lto.x3600_y1800_t12.2010.nc4
BC_AVIATION_CDS NA Y Y %y4-%m2-%d2t12:00:00 none none bc_aviation htap-v2.2.emis_bc.aviation_cds.x3600_y1800_t12.2010.nc4
BC_AVIATION_CRS NA Y Y %y4-%m2-%d2t12:00:00 none none bc_aviation htap-v2.2.emis_bc.aviation_crs.x3600_y1800_t12.2010.nc4
BC_AVIATION_LTO NA Y Y %y4-%m2-%d2t12:00:00 none none bc_aviation ExtData/htap-v2.2.emis_bc.aviation_lto.x3600_y1800_t12.2010.nc4
BC_AVIATION_CDS NA Y Y %y4-%m2-%d2t12:00:00 none none bc_aviation ExtData/htap-v2.2.emis_bc.aviation_cds.x3600_y1800_t12.2010.nc4
BC_AVIATION_CRS NA Y Y %y4-%m2-%d2t12:00:00 none none bc_aviation ExtData/htap-v2.2.emis_bc.aviation_crs.x3600_y1800_t12.2010.nc4

#============================================================================================================
# BROWN CARBON
Expand Down Expand Up @@ -134,19 +134,19 @@ BRC_AVIATION_CRS NA Y Y %y4-%m2-%d2t12:00:00 none none oc_aviation /dev/null
pSOA_BIOB_VOC NA Y Y %y4-%m2-%d2t12:00:00 none none biofuel /dev/null

# # ======= Nitrate Sources ========
# EMI_NH3_AG 'kg m-2 s-1' N Y %y4-%m2-%d2T12:00:00 none none NH3_ag CEDS.2019.emis.20210325.nc
# EMI_NH3_AG 'kg m-2 s-1' N Y %y4-%m2-%d2T12:00:00 none none NH3_ag ExtData/CEDS.2019.emis.20210325.nc
# EMI_NH3_EN 'kg m-2 s-1' N Y %y4-%m2-%d2T12:00:00 none none NH3_en /dev/null
# EMI_NH3_IN 'kg m-2 s-1' N Y %y4-%m2-%d2T12:00:00 none none NH3_in CEDS.2019.emis.20210325.nc
# EMI_NH3_RE 'kg m-2 s-1' N Y %y4-%m2-%d2T12:00:00 none none NH3_re CEDS.2019.emis.20210325.nc
# EMI_NH3_TR 'kg m-2 s-1' N Y %y4-%m2-%d2T12:00:00 none none NH3_tr CEDS.2019.emis.20210325.nc
# EMI_NH3_OC 'kg m-2 s-1' Y Y %y4-%m2-%d2T12:00:00 none none emiss_ocn GEIA.emis_NH3.ocean.x576_y361.t12.20080715_12z.nc4
# EMI_NH3_IN 'kg m-2 s-1' N Y %y4-%m2-%d2T12:00:00 none none NH3_in ExtData/CEDS.2019.emis.20210325.nc
# EMI_NH3_RE 'kg m-2 s-1' N Y %y4-%m2-%d2T12:00:00 none none NH3_re ExtData/CEDS.2019.emis.20210325.nc
# EMI_NH3_TR 'kg m-2 s-1' N Y %y4-%m2-%d2T12:00:00 none none NH3_tr ExtData/CEDS.2019.emis.20210325.nc
# EMI_NH3_OC 'kg m-2 s-1' Y Y %y4-%m2-%d2T12:00:00 none none emiss_ocn ExtData/GEIA.emis_NH3.ocean.x576_y361.t12.20080715_12z.nc4

# # --------------------------------------------------------------
# # If using 64 levels please replace this section with the correct values (ie replace 127 with 64)
# NITRATE_HNO3 'mol mol-1' Y N %y4-%m2-%d2T12:00:00 none 0.20 hno3 GMI.vmr_HNO3.x144_y91.t12.2006.nc4
# NITRATE_HNO3 'mol mol-1' Y N %y4-%m2-%d2T12:00:00 none 0.20 hno3 ExtData/GMI.vmr_HNO3.x144_y91.t12.2006.nc4
# # --------------------------------------------------------------
# NI_regionMask NA Y V - none none REGION_MASK ARCTAS.region_mask.x540_y361.2008.nc
SU_BIOMASS NA N Y %y4-%m2-%d2t12:00:00 none 0.7778 SO2 QFED_Blended_20210325_50days.nc
OC_BIOMASS NA N Y %y4-%m2-%d2t12:00:00 none 0.7778 OC QFED_Blended_20210325_50days.nc
BC_BIOMASS NA N Y %y4-%m2-%d2t12:00:00 none 0.7778 BC QFED_Blended_20210325_50days.nc
# NI_regionMask NA Y V - none none REGION_MASK ExtData/ARCTAS.region_mask.x540_y361.2008.nc
SU_BIOMASS NA N Y %y4-%m2-%d2t12:00:00 none 0.7778 SO2 ExtData/QFED_Blended_20210325_50days.nc
OC_BIOMASS NA N Y %y4-%m2-%d2t12:00:00 none 0.7778 OC ExtData/QFED_Blended_20210325_50days.nc
BC_BIOMASS NA N Y %y4-%m2-%d2t12:00:00 none 0.7778 BC ExtData/QFED_Blended_20210325_50days.nc
%%
6 changes: 4 additions & 2 deletions tests/tests/cpld_bmark_p8
View file
Original file line number Diff line number Diff line change
Expand Up @@ -147,6 +147,8 @@ export MOM6_DIAG_COORD_DEF_Z_FILE=interpolate_zgrid_26L.nc
export DO_SPPT=.true.
export DO_SHUM=.false.
export DO_SKEB=.true.
export SKEB="0.8,-999,-999,-999,-999"
export SPPT="0.56,0.28,0.14,0.056,0.028"
export READ_INCREMENT=.true.
export RES_LATLON_DYNAMICS="'fv3_increment.nc'"
export PERT_CLDS=.true.
Expand All @@ -162,8 +164,8 @@ export ODA_VINC_VAR="'v_pert'"
# MOM6 stochastics
export DO_OCN_SPPT=True
export PERT_EPBL=True
export OCN_SPPT=1.0
export EPBL=1.0
export OCN_SPPT="0.8,0.4,0.2,0.08,0.04"
export EPBL="0.8,0.4,0.2,0.08,0.04"

#Warm start
export RUNTYPE=continue
Expand Down

0 comments on commit e6d4338

Please sign in to comment.