[v2.1.3-8.4] Changes in preparation for SFE

Author: chunhuazhou

fix/gsi_bec/gsiparm_regional.anl

@@ -4,6 +4,7 @@
  /
  &GRIDOPTS
    regional=.true.,fv3_regional=.false.,mpas_regional=.true.,
+   use_fv3_grid_spec=.false.,
    nlon=@GSIBEC_NLON@,nlat=@GSIBEC_NLAT@,nsig=@GSIBEC_NSIG@,
    rlat_start=@GSIBEC_LAT_START@,rlat_end=@GSIBEC_LAT_END@,
    rlon_start=@GSIBEC_LON_START@,rlon_end=@GSIBEC_LON_END@,
@@ -31,7 +32,7 @@ met_guess::
   tv       @GSIBEC_NSIG@       2         virtual_temperature  tv
   tsen     @GSIBEC_NSIG@       2         sensible_temperature tsen 
   q        @GSIBEC_NSIG@       2         specific_humidity    sphum
-  prse     @GSIBEC_NSIG@      -1         air_pressure         prse    
+  prsl     @GSIBEC_NSIG@      -1         air_pressure         prsl    
 #  oz       @GSIBEC_NSIG@       2         ozone                o3mr 
 #  delp     @GSIBEC_NSIG@      -1         fv3_del_pressure     delp    
 ::

parm/hrrrv5/namelist.init_atmosphere

@@ -54,6 +54,10 @@
 &physics
     config_tsk_seaice_threshold = 271.4
 /
+&lbcs
+    config_lbc_hydrometeors_rrfs = ${lbc_hydrometeors_rrfs}
+    config_lbc_hydrometeors_gfs = ${lbc_hydrometeors_gfs}
+/
 &io
     config_pio_num_iotasks = 0
     config_pio_stride = 1

parm/jedivar.yaml

@@ -158,6 +158,7 @@ cost function:
            engine:
              type: H5File
              obsfile: data/obs/ioda_mrms_refl.nc
+             missing file action: "warn"
          obsdataout:
            engine:
              type: H5File
@@ -3035,14 +3036,20 @@ cost function:
          simulated variables: [windEastward, windNorthward]
 
        obs operator:
-         name: SfcCorrected
-         correction scheme to use: GSL
-         variables:
-         - name: windEastward
-         - name: windNorthward
+           name: VertInterp
+           vertical coordinate: air_pressure
+           observation vertical coordinate: pressure
+           observation vertical coordinate group: MetaData
+           interpolation method: log-linear
+           #gsi geovals:
+           #  filename: "obsout/sfc_uv_geoval_2022052619.nc"
+           #  levels_are_top_down: False
+           variables:
+           - name: windEastward
+           - name: windNorthward
 
        linear obs operator:
-         name: Identity
+         name: VertInterp
 
        obs error:
          covariance model: diagonal

scripts/exrrfs_ic.sh

@@ -14,6 +14,8 @@ init_case=7
 start_time=$(date -d "${CDATE:0:8} ${CDATE:8:2}" +%Y-%m-%d_%H:%M:%S)
 end_time=${start_time}
 
+lbc_hydrometeors_rrfs=true
+lbc_hydrometeors_gfs=false
 if [[ "${prefix}" == "RAP" || "${prefix}" == "HRRR" ]]; then
   nfglevels=51
   nfgsoillevels=9

scripts/exrrfs_lbc.sh

@@ -29,6 +29,9 @@ start_time=$(date -d "${EDATE:0:8} ${EDATE:8:2}" +%Y-%m-%d_%H:%M:%S)
 EDATE=$(${NDATE} "${fhr_end}" "${CDATEin}")
 end_time=$(date -d "${EDATE:0:8} ${EDATE:8:2}" +%Y-%m-%d_%H:%M:%S)
 
+lbc_hydrometeors_rrfs=true
+lbc_hydrometeors_gfs=false
+
 if [[ "${prefix}" == "RAP" || "${prefix}" == "HRRR" ]]; then
   nfglevels=51
   nfgsoillevels=9
@@ -38,9 +41,13 @@ elif  [[ "${prefix}" == "RRFS" ]]; then
 elif  [[ "${prefix}" == "GFS" ]]; then
   nfglevels=58
   nfgsoillevels=4
+  lbc_hydrometeors_rrfs=false
+  lbc_hydrometeors_gfs=true
 elif  [[ "${prefix}" == "GEFS" ]]; then
   nfglevels=32
   nfgsoillevels=4
+  lbc_hydrometeors_rrfs=false
+  lbc_hydrometeors_gfs=false
 fi
 nsoillevels=${NSOIL_LEVELS}