Several bug fixes: (#558)

1) fix streak issue by change namelist to fv_sg_adj = 300 n_sponge = -1 sg_cutoff = 40000.0 2) use mesonet uselist in retros 3) change modules in run_enkfupdt.local.lua for running enkf after wcoss2 updates. Co-authored-by: ming hu <[email protected]>
NOAA-EMC · Nov 27, 2024 · dd62411 · dd62411
1 parent 4ad5e4f
commit dd62411
Show file tree

Hide file tree

Showing 4 changed files with 22 additions and 17 deletions.
diff --git a/modulefiles/tasks/wcoss2/run_enkfupdt.local.lua b/modulefiles/tasks/wcoss2/run_enkfupdt.local.lua
@@ -1,4 +1,4 @@
-load("python_srw")
+load(pathJoin("envvar", os.getenv("envvar_ver")))
 
 load(pathJoin("PrgEnv-intel", os.getenv("PrgEnv_intel_ver")))
 load(pathJoin("intel", os.getenv("intel_ver")))
@@ -9,14 +9,14 @@ load(pathJoin("cray-pals", os.getenv("cray_pals_ver")))
 prepend_path("MODULEPATH", os.getenv("modulepath_compiler"))
 prepend_path("MODULEPATH", os.getenv("modulepath_mpi"))
 
+load(pathJoin("jasper", os.getenv("jasper_ver")))
+load(pathJoin("zlib", os.getenv("zlib_ver")))
+load(pathJoin("libpng", os.getenv("libpng_ver")))
+load(pathJoin("hdf5", os.getenv("hdf5_ver")))
 load(pathJoin("netcdf", os.getenv("netcdf_ver")))
+
 load(pathJoin("bacio", os.getenv("bacio_ver")))
-load(pathJoin("w3emc", os.getenv("w3emc_ver")))
-load(pathJoin("sp", os.getenv("sp_ver")))
 load(pathJoin("ip", os.getenv("ip_ver")))
-load(pathJoin("nemsio", os.getenv("nemsio_ver")))
-
-load(pathJoin("udunits", os.getenv("udunits_ver")))
-load(pathJoin("gsl", os.getenv("gsl_ver")))
-load(pathJoin("nco", os.getenv("nco_ver")))
-
+load(pathJoin("sp", os.getenv("sp_ver")))
+load(pathJoin("w3emc", os.getenv("w3emc_ver")))
+load(pathJoin("w3nco", os.getenv("w3nco_ver")))
diff --git a/parm/input.nml.FV3 b/parm/input.nml.FV3
@@ -80,7 +80,7 @@
     fill = .true.
     full_zs_filter = .false.
     fv_debug = .false.
-    fv_sg_adj = 7200
+    fv_sg_adj = 300
     gfs_phil = .false.
     hord_dp = 6
     hord_mt = 6
@@ -98,7 +98,7 @@
     make_nh = .true.
     mountain = .false.
     n_split = 5
-    n_sponge = 65 
+    n_sponge = -1
     n_zs_filter = 0
     na_init = 1
     ncep_ic = .false.
@@ -119,7 +119,7 @@
     res_latlon_dynamics = 'fv3_increment.nc'
     reset_eta = .false.
     rf_cutoff = 20.e2
-    sg_cutoff = -1
+    sg_cutoff = 40000.0
     tau = 5.0
     use_hydro_pressure = .false.
     vtdm4 = 0.075

diff --git a/parm/input.nml.RRFSFW b/parm/input.nml.RRFSFW
@@ -82,7 +82,7 @@
     fill = .true.
     full_zs_filter = .false.
     fv_debug = .false.
-    fv_sg_adj = 7200
+    fv_sg_adj = 300
     gfs_phil = .false.
     hord_dp = 6
     hord_mt = 6
@@ -100,7 +100,7 @@
     make_nh = .true.
     mountain = .false.
     n_split = 5
-    n_sponge = 65
+    n_sponge = -1
     n_zs_filter = 0
     na_init = 1
     ncep_ic = .false.
@@ -123,7 +123,7 @@
     res_latlon_dynamics = 'fv3_increment.nc'
     reset_eta = .false.
     rf_cutoff = 2000.0
-    sg_cutoff = -1
+    sg_cutoff = 40000.0
     tau = 5.0
     use_hydro_pressure = .false.
     vtdm4 = 0.02

diff --git a/scripts/exrrfs_run_analysis.sh b/scripts/exrrfs_run_analysis.sh
@@ -775,7 +775,8 @@ fi
 
 # Get aircraft reject list
 for reject_list in "${AIRCRAFT_REJECT}/current_bad_aircraft.txt" \
-                   "${AIRCRAFT_REJECT}/${AIR_REJECT_FN}"
+                   "${AIRCRAFT_REJECT}/${AIR_REJECT_FN}" \
+                   "${FIX_GSI}/current_bad_aircraft.txt"
 do
   if [ -r $reject_list ]; then
     cp $reject_list current_bad_aircraft
@@ -791,7 +792,8 @@ fi
 gsd_sfcobs_uselist="gsd_sfcobs_uselist.txt"
 for use_list in "${SFCOBS_USELIST}/current_mesonet_uselist.txt" \
                 "${SFCOBS_USELIST}/${MESO_USELIST_FN}"      \
-                "${SFCOBS_USELIST}/gsd_sfcobs_uselist.txt"
+                "${SFCOBS_USELIST}/gsd_sfcobs_uselist.txt"  \
+                "${FIX_GSI}/gsd_sfcobs_uselist.txt"
 do 
   if [ -r $use_list ] ; then
     cp $use_list  $gsd_sfcobs_uselist
@@ -1041,6 +1043,7 @@ else
   cat fit_p1 fit_w1 fit_t1 fit_q1 fit_pw1 fit_rad1 fit_rw1 > $COMOUT/rrfs.t${HH}z.fits.tm00
   cat fort.208 fort.210 fort.211 fort.212 fort.213 fort.220 > $COMOUT/rrfs.t${HH}z.fits2.tm00
   cat fort.238 > $COMOUT/rrfs.t${HH}z.fits3.tm00
+  cp -L dbzobs.nc  $COMOUT/rrfs.mrms.${YYYYMMDDHH}.nc
 fi
 #
 #-----------------------------------------------------------------------
@@ -1185,6 +1188,8 @@ if [ "${DO_GSIDIAG_OFFLINE}" = "FALSE" ]; then
       cp ./satbias_pc.out ${satbias_dir}/rrfs.${spinup_or_prod_rrfs}.${YYYYMMDDHH}_satbias_pc
       cp ./satbias_out ${COMOUT}/rrfs.${spinup_or_prod_rrfs}.${YYYYMMDDHH}_satbias
       cp ./satbias_pc.out ${COMOUT}/rrfs.${spinup_or_prod_rrfs}.${YYYYMMDDHH}_satbias_pc
+      cp -L dbzobs.nc  $COMOUT/rrfs.mrms.${YYYYMMDDHH}.nc
+
     fi
   fi
 fi # run diag inline (with GSI)