Trying to use EP5 restart files. Errors in gocart

NickSzapiro-NOAA · Jun 13, 2024 · d6c416e · d6c416e
1 parent 949aa06
commit d6c416e
Show file tree

Hide file tree

Showing 8 changed files with 691 additions and 63 deletions.
diff --git a/tests/fv3_conf/cpld_control_run.IN b/tests/fv3_conf/cpld_control_run.IN
@@ -18,9 +18,13 @@ ICERES=${OCNRES:0:1}.${OCNRES:1}
 V2_SFC_FILE=${V2_SFC_FILE:-false}
 
 if [[ $BMIC == .true. ]]; then
-  FV3_IC=@[INPUTDATA_ROOT_BMIC]/${SYEAR}${SMONTH}${SDAY}${SHOUR}/@[FV3BMIC]/@[ATMRES]_L@[NPZ]/INPUT
-  MOM_IC=@[INPUTDATA_ROOT_BMIC]/${SYEAR}${SMONTH}${SDAY}${SHOUR}/mom6_da
-  ICE_IC=@[INPUTDATA_ROOT_BMIC]/${SYEAR}${SMONTH}${SDAY}${SHOUR}/cpc
+  #FV3_IC=@[INPUTDATA_ROOT_BMIC]/${SYEAR}${SMONTH}${SDAY}${SHOUR}/@[FV3BMIC]/@[ATMRES]_L@[NPZ]/INPUT
+  #MOM_IC=@[INPUTDATA_ROOT_BMIC]/${SYEAR}${SMONTH}${SDAY}${SHOUR}/mom6_da
+  #ICE_IC=@[INPUTDATA_ROOT_BMIC]/${SYEAR}${SMONTH}${SDAY}${SHOUR}/cpc
+  FV3_IC=@[INPUTDATA_ROOT_BMIC]
+  MOM_IC=@[INPUTDATA_ROOT_BMIC]
+  ICE_IC=@[INPUTDATA_ROOT_BMIC]
+  WAV_IC=@[INPUTDATA_ROOT_BMIC]
 else
   if [[ $ATMRES == C48 ]]; then
      FV3_IC=@[INPUTDATA_ROOT]/${FV3_DIR}/INPUT_L127_mx${OCNRES}/${SYEAR}${SMONTH}${SDAY}${SHOUR}
@@ -154,11 +158,21 @@ else
      fi
 
   else
-     # C48mx500 warm start from ICs
-     cp    ${FV3_IC}/*.nc ./INPUT
-     cp    ${MOM_IC}/MOM.res.nc ./INPUT
-     cp    ${ICE_IC}/* .
-     cp    @[INPUTDATA_ROOT]/CMEPS_IC/${ATMRES}mx${OCNRES}/${SYEAR}${SMONTH}${SDAY}${SHOUR}/* .
+     if [[ $BMIC == .true. ]]; then
+       cp ${FV3_IC}/gfs_ctrl.nc ./INPUT
+       cp ${FV3_IC}/gfs_data*.nc ./INPUT
+       cp ${FV3_IC}/sfc_data*.nc ./INPUT
+       cp ${MOM_IC}/MOM*.nc ./INPUT
+       cp ${ICE_IC}/cice*nc .
+       cp ${FV3_IC}/ufs.cpld.cpl*nc .
+       cp ${WAV_IC}/restart.ww3 .
+     else
+       # C48mx500 warm start from ICs
+       cp    ${FV3_IC}/*.nc ./INPUT
+       cp    ${MOM_IC}/MOM.res.nc ./INPUT
+       cp    ${ICE_IC}/* .
+       cp    @[INPUTDATA_ROOT]/CMEPS_IC/${ATMRES}mx${OCNRES}/${SYEAR}${SMONTH}${SDAY}${SHOUR}/* .
+     fi
   fi
 fi
 
@@ -180,6 +194,21 @@ if [ $IAER = 1011 ]; 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
+fi  
+
+#rad_c_gocart
+if [ $IAER = 2011 ]; then
+  cp    @[INPUTDATA_ROOT]/FV3_input_data_INCCN_aeroclim/aer_data/LUTS/optics_BC.v1_3.dat  optics_BC.dat
+  cp    @[INPUTDATA_ROOT]/FV3_input_data_INCCN_aeroclim/aer_data/LUTS/optics_OC.v1_3.dat  optics_OC.dat
+  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/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
+  cp ${DIR_OPT}/optics_SS.v3_3.nc  optics_SS.nc
+  cp ${DIR_OPT}/optics_SU.v1_3.nc  optics_SU.nc
 fi
 
 cp    @[INPUTDATA_ROOT]/FV3_input_data/ugwp_c384_tau.nc ./ugwp_limb_tau.nc
@@ -195,5 +224,9 @@ fi
 
 #prognostic aerosols
 if [ $CPLCHM = .true. ]; then
-  ln -sf  @[INPUTDATA_ROOT]/GOCART/p8c_5d/ExtData .
+  if [[ $BMIC == .true. ]]; then
+    ln -sf /scratch1/NCEPDEV/nems/Nick.Szapiro/tasks/input_data/gefs.v13/CHM/ExtData .
+  else
+    ln -sf  @[INPUTDATA_ROOT]/GOCART/p8c_5d/ExtData .
+  fi
 fi
diff --git a/tests/parm/diag_table/diag_table_gefsv13 b/tests/parm/diag_table/diag_table_gefsv13
diff --git a/tests/parm/field_table/field_table_gefsv13 b/tests/parm/field_table/field_table_gefsv13
@@ -0,0 +1,197 @@
+# added by FRE: sphum must be present in atmos
+# specific humidity for moist runs
+ "TRACER", "atmos_mod", "sphum"
+           "longname",     "specific humidity"
+           "units",        "kg/kg"
+       "profile_type", "fixed", "surface_value=3.e-6" /
+# prognostic cloud water mixing ratio
+ "TRACER", "atmos_mod", "liq_wat"
+           "longname",     "cloud water mixing ratio"
+           "units",        "kg/kg"
+       "profile_type", "fixed", "surface_value=1.e30" /
+# prognostic ice water mixing ratio
+ "TRACER", "atmos_mod", "ice_wat"
+           "longname",     "cloud ice mixing ratio"
+           "units",        "kg/kg"
+       "profile_type", "fixed", "surface_value=1.e30" /
+# prognostic rain water mixing ratio
+ "TRACER", "atmos_mod", "rainwat"
+           "longname",     "rain water mixing ratio"
+           "units",        "kg/kg"
+       "profile_type", "fixed", "surface_value=1.e30" /
+# prognostic snow water mixing ratio
+ "TRACER", "atmos_mod", "snowwat"
+           "longname",     "snow water mixing ratio"
+           "units",        "kg/kg"
+       "profile_type", "fixed", "surface_value=1.e30" /
+# prognostic Grau water mixing ratio
+ "TRACER", "atmos_mod", "graupel"
+           "longname",     "graupel mixing ratio"
+           "units",        "kg/kg"
+       "profile_type", "fixed", "surface_value=1.e30" /
+# prognostic cloud water number concentration - not for non-aerosol runs
+# "TRACER", "atmos_mod", "water_nc"
+#           "longname",     "cloud  liquid water  number concentration"
+#           "units",        "/kg"
+#       "profile_type", "fixed", "surface_value=0.0" /
+# prognostic cloud ice number concentration
+ "TRACER", "atmos_mod", "ice_nc"
+           "longname",     "cloud ice water number concentration"
+           "units",        "/kg"
+       "profile_type", "fixed", "surface_value=0.0" /
+# prognostic rain number concentration
+ "TRACER", "atmos_mod", "rain_nc"
+           "longname",     "rain number concentration"
+           "units",        "/kg"
+       "profile_type", "fixed", "surface_value=0.0" /
+# prognostic ozone mixing ratio tracer
+ "TRACER", "atmos_mod", "o3mr"
+           "longname",     "ozone mixing ratio"
+           "units",        "kg/kg"
+       "profile_type", "fixed", "surface_value=1.e30" /
+# water- and ice-friendly aerosols (Thompson) - not for non-aerosol runs
+# "TRACER", "atmos_mod", "liq_aero"
+#           "longname",     "water-friendly aerosol number concentration"
+#           "units",        "/kg"
+#       "profile_type", "fixed", "surface_value=0.0" /
+# "TRACER", "atmos_mod", "ice_aero"
+#           "longname",     "ice-friendly aerosol number concentration"
+#           "units",        "/kg"
+#       "profile_type", "fixed", "surface_value=0.0" /
+# prognostic subgrid scale turbulent kinetic energy
+ "TRACER", "atmos_mod", "sgs_tke"
+           "longname",     "subgrid scale turbulent kinetic energy"
+           "units",        "m2/s2"
+       "profile_type", "fixed", "surface_value=0.0" /
+# prognotsitc sigmab tracer
+ "TRACER", "atmos_mod", "sigmab"
+           "longname",     "sigma fraction"
+           "units",        "fraction"
+       "profile_type", "fixed", "surface_value=0.0" /
+# prognostic aerosols
+ "TRACER", "atmos_mod", "so2"
+           "longname",     "so2 mixing ratio"
+           "units",        "ppm"
+       "tracer_usage", "chemistry"
+       "profile_type", "fixed", "surface_value=0.0" /
+ "TRACER", "atmos_mod", "so4"
+           "longname",     "sulfate mixing ratio"
+           "units",        "ug/kg"
+       "tracer_usage", "chemistry"
+       "profile_type", "fixed", "surface_value=0.0" /
+ "TRACER", "atmos_mod", "dms"
+           "longname",     "DMS mixing ratio"
+           "units",        "ppm"
+       "tracer_usage", "chemistry"
+       "profile_type", "fixed", "surface_value=0.0" /
+ "TRACER", "atmos_mod", "msa"
+           "longname",     "msa mixing ratio"
+           "units",        "ppm"
+       "tracer_usage", "chemistry"
+       "profile_type", "fixed", "surface_value=0.0" /
+ "TRACER", "atmos_mod", "bc1"
+           "longname",     "hydrophobic black carbon mixing ratio"
+           "units",        "ug/kg"
+       "tracer_usage", "chemistry"
+       "profile_type", "fixed", "surface_value=0.0" /
+ "TRACER", "atmos_mod", "bc2"
+           "longname",     "hydrophillic black carbon mixing ratio"
+           "units",        "ug/kg"
+       "tracer_usage", "chemistry"
+       "profile_type", "fixed", "surface_value=0.0" /
+ "TRACER", "atmos_mod", "oc1"
+           "longname",     "hydrophobic organic carbon mixing ratio"
+           "units",        "ug/kg"
+       "tracer_usage", "chemistry"
+       "profile_type", "fixed", "surface_value=0.0" /
+ "TRACER", "atmos_mod", "oc2"
+           "longname",     "hydrophillic organic carbon mixing ratio"
+           "units",        "ug/kg"
+       "tracer_usage", "chemistry"
+       "profile_type", "fixed", "surface_value=0.0" /
+ "TRACER", "atmos_mod", "dust1"
+           "longname",     "fine dust1 mixing ratio"
+           "units",        "ug/kg"
+       "tracer_usage", "chemistry"
+       "profile_type", "fixed", "surface_value=0.0" /
+ "TRACER", "atmos_mod", "dust2"
+           "longname",     "fine dust2 mixing ratio"
+           "units",        "ug/kg"
+       "tracer_usage", "chemistry"
+       "profile_type", "fixed", "surface_value=0.0" /
+ "TRACER", "atmos_mod", "dust3"
+           "longname",     "coarse dust3 mixing ratio"
+           "units",        "ug/kg"
+       "tracer_usage", "chemistry"
+       "profile_type", "fixed", "surface_value=0.0" /
+ "TRACER", "atmos_mod", "dust4"
+           "longname",     "coarse dust4 mixing ratio"
+           "units",        "ug/kg"
+       "tracer_usage", "chemistry"
+       "profile_type", "fixed", "surface_value=0.0" /
+ "TRACER", "atmos_mod", "dust5"
+           "longname",     "coarse dust5 mixing ratio"
+           "units",        "ug/kg"
+       "tracer_usage", "chemistry"
+       "profile_type", "fixed", "surface_value=0.0" /
+ "TRACER", "atmos_mod", "seas1"
+           "longname",     "seasalt1 mixing ratio"
+           "units",        "ug/kg"
+       "tracer_usage", "chemistry"
+       "profile_type", "fixed", "surface_value=0.0" /
+   "TRACER", "atmos_mod", "seas2"
+           "longname",     "seasalt2 mixing ratio"
+           "units",        "ug/kg"
+       "tracer_usage", "chemistry"
+       "profile_type", "fixed", "surface_value=0.0" /
+ "TRACER", "atmos_mod", "seas3"
+           "longname",     "seasalt3 mixing ratio"
+           "units",        "ug/kg"
+       "tracer_usage", "chemistry"
+       "profile_type", "fixed", "surface_value=0.0" /
+ "TRACER", "atmos_mod", "seas4"
+           "longname",     "seasalt4 mixing ratio"
+           "units",        "ug/kg"
+       "tracer_usage", "chemistry"
+       "profile_type", "fixed", "surface_value=0.0" /
+ "TRACER", "atmos_mod", "seas5"
+           "longname",     "seasalt5 mixing ratio"
+           "units",        "ug/kg"
+       "tracer_usage", "chemistry"
+       "profile_type", "fixed", "surface_value=0.0" /
+# "TRACER", "atmos_mod", "nh3"
+#           "longname",     "primary NH3 mixing ratio"
+#           "units",        "ug/kg"
+#       "tracer_usage", "chemistry"
+#       "profile_type", "fixed", "surface_value=0.0" /
+# "TRACER", "atmos_mod", "nh4a"
+#           "longname",     "primary NH4a mixing ratio"
+#           "units",        "ug/kg"
+#       "tracer_usage", "chemistry"
+#       "profile_type", "fixed", "surface_value=0.0" /
+# "TRACER", "atmos_mod", "no3an1"
+#           "longname",     "primary NO3an1 mixing ratio"
+#           "units",        "ug/kg"
+#       "tracer_usage", "chemistry"
+#       "profile_type", "fixed", "surface_value=0.0" /
+# "TRACER", "atmos_mod", "no3an2"
+#           "longname",     "primary NO3an2 mixing ratio"
+#           "units",        "ug/kg"
+#       "tracer_usage", "chemistry"
+#       "profile_type", "fixed", "surface_value=0.0" /
+# "TRACER", "atmos_mod", "no3an3"
+#           "longname",     "primary NO3an3 mixing ratio"
+#           "units",        "ug/kg"
+#       "tracer_usage", "chemistry"
+#       "profile_type", "fixed", "surface_value=0.0" /
+# diagnostic PM tracers
+ "TRACER", "atmos_mod", "pm25"
+           "longname",     "primary PM25 mixing ratio"
+           "units",        "ug/m3"
+       "tracer_usage", "chemistry", "type=diagnostic"
+       "profile_type", "fixed", "surface_value=0.0" /
+ "TRACER", "atmos_mod", "pm10"
+           "longname",     "primary PM10 mixing ratio"
+           "units",        "ug/m3"
+       "tracer_usage", "chemistry", "type=diagnostic"
+       "profile_type", "fixed", "surface_value=0.0" /
diff --git a/tests/parm/global_control.nml.IN b/tests/parm/global_control.nml.IN
@@ -272,7 +272,7 @@
 @[HIDE_UGWPV1] knob_ugwp_lzmin    = 0.75e3
 @[HIDE_UGWPV1] knob_ugwp_lzstar   = 2.0e3
 @[HIDE_UGWPV1] knob_ugwp_taumin   = 0.25e-3
-@[HIDE_UGWPV1] knob_ugwp_tauamp   = 3.0e-3
+@[HIDE_UGWPV1] knob_ugwp_tauamp   = 0.8e-3
 @[HIDE_UGWPV1] knob_ugwp_lhmet    = 200.0e3
 @[HIDE_UGWPV1] knob_ugwp_orosolv  = 'pss-1986'
 /

diff --git a/tests/parm/model_configure.IN b/tests/parm/model_configure.IN
@@ -8,7 +8,9 @@ nhours_fcst:             @[FHMAX]
 fhrot:                   @[FHROT]
 
 dt_atmos:                @[DT_ATMOS]
+calendar:                'julian'
 restart_interval:        @[RESTART_INTERVAL]
+output_1st_tstep_rst:    .false.
 
 quilting:                @[QUILTING]
 quilting_restart:        @[QUILTING_RESTART]
@@ -25,6 +27,7 @@ output_grid:             @[OUTPUT_GRID]
 output_file:             @[OUTPUT_FILE]
 zstandard_level:         @[ZSTANDARD_LEVEL]
 ideflate:                @[IDEFLATE]
+nbits:                   14
 quantize_mode:           'quantize_bitround'
 quantize_nsd:            @[QUANTIZE_NSD]
 ichunk2d:                @[ICHUNK2D]

diff --git a/tests/rt.conf b/tests/rt.conf
@@ -38,7 +38,8 @@ RUN | cpld_control_sfs                                  | - noaacloud
 COMPILE | s2swa_32bit_pdlib_debug | intel | -DAPP=S2SWA -D32BIT=ON -DCCPP_SUITES=FV3_GFS_v17_coupled_p8_ugwpv1 -DPDLIB=ON -DDEBUG=ON | - noaacloud jet | fv3 |
 RUN | cpld_debug_gfsv17                                 | - noaacloud jet                      | baseline |
 
-COMPILE | s2swa | intel | -DAPP=S2SWA -DCCPP_SUITES=FV3_GFS_v17_coupled_p8 | | fv3 |
+#COMPILE | s2swa | intel | -DAPP=S2SWA -DCCPP_SUITES=FV3_GFS_v17_coupled_p8 | | fv3 |
+COMPILE | s2swa | intel | -DAPP=S2SWA -DCCPP_SUITES=FV3_GFS_v17_coupled_p8_ugwpv1 | | fv3 |
 RUN | cpld_control_p8                                   | - noaacloud                          | baseline |
 RUN | cpld_control_p8.v2.sfc                            | - noaacloud                          | baseline |
 RUN | cpld_restart_p8                                   | - noaacloud                          |          | cpld_control_p8

diff --git a/tests/rt.sh b/tests/rt.sh
@@ -1020,8 +1020,10 @@ if [[ "${CREATE_BASELINE}" == false ]] ; then
 fi
 
 INPUTDATA_ROOT=${INPUTDATA_ROOT:-${DISKNM}/NEMSfv3gfs/input-data-20240501}
-INPUTDATA_ROOT_WW3=${INPUTDATA_ROOT}/WW3_input_data_20240214
-INPUTDATA_ROOT_BMIC=${INPUTDATA_ROOT_BMIC:-${DISKNM}/NEMSfv3gfs/BM_IC-20220207}
+#INPUTDATA_ROOT_WW3=${INPUTDATA_ROOT}/WW3_input_data_20240214
+INPUTDATA_ROOT_WW3=/scratch1/NCEPDEV/nems/Nick.Szapiro/tasks/updateToEP5/WW3_inputdata
+#INPUTDATA_ROOT_BMIC=${INPUTDATA_ROOT_BMIC:-${DISKNM}/NEMSfv3gfs/BM_IC-20220207}
+INPUTDATA_ROOT_BMIC=/scratch1/NCEPDEV/nems/Nick.Szapiro/tasks/input_data/gefs.v13/IC
 
 shift $((OPTIND-1))
 if [[ $# -gt 1 ]]; then