Add gsibec as a staticb option for fv3jedi_var.x (#166)

build.sh

@@ -60,7 +60,7 @@ case ${BUILD_TARGET} in
     module purge
     module use $dir_root/modulefiles
     module load GDAS/$BUILD_TARGET
-    CMAKE_OPTS+=" -DMPIEXEC_EXECUTABLE=$MPIEXEC_EXEC -DMPIEXEC_NUMPROC_FLAG=$MPIEXEC_NPROC"
+    CMAKE_OPTS+=" -DMPIEXEC_EXECUTABLE=$MPIEXEC_EXEC -DMPIEXEC_NUMPROC_FLAG=$MPIEXEC_NPROC -DBUILD_GSIBEC=ON"
     module list
     set -e
     ;;

parm/atm/berror/staticb_gsibec.yaml

@@ -0,0 +1,38 @@
+background error:
+  covariance model: hybrid
+  components:
+  - covariance:
+      covariance model: SABER
+      saber blocks:
+      - saber block name: gsi covariance
+        saber central block: true
+        input variables: &control_vars [eastward_wind,northward_wind,air_temperature,surface_pressure,
+                                 specific_humidity,mole_fraction_of_ozone_in_air]
+        output variables: *control_vars
+        gsi akbk: &akbk $(fv3jedi_fix_dir)/akbk.nc4
+        gsi error covariance file: &gsiberr $(fv3jedi_staticb_dir)/gsi-coeffs-gfs-global-l127x192y96.nc4
+        gsi berror namelist file: &gsibnml $(fv3jedi_staticb_dir)/gfs_gsi_global.nml
+        processor layout x direction: 3
+        processor layout y direction: 2
+        debugging mode: false
+      - saber block name: gsi interpolation to model grid
+        input variables: *control_vars
+        output variables: *control_vars
+        gsi akbk: *akbk
+        gsi error covariance file: *gsiberr
+        gsi berror namelist file: *gsibnml
+        processor layout x direction: 3
+        processor layout y direction: 2
+        debugging mode: false
+      linear variable change:
+        linear variable change name: Control2Analysis
+        input variables: *control_vars
+        output variables: &3dvars_anal [ua,va,t,ps,sphum,o3mr]
+    weight:
+      value: 1.00
+  - covariance:
+      covariance model: FV3JEDI-ID
+      date: $(window_begin)
+    weight:
+      value: 0.0
+

ush/examples/run_jedi_exe/3dvar_hera.yaml

@@ -16,6 +16,7 @@ executable options:
   dump: gdas
   case: C96
   case_anl: C48
+  staticb_type: bump
   dohybvar: false
   levs: 128
   interp_method: barycentric

ush/run_jedi_exe.py

@@ -73,6 +73,7 @@ def run_jedi_exe(yamlconfig):
         single_exec = True
         var_config = {
             'BERROR_YAML': executable_subconfig.get('berror_yaml', './'),
+            'STATICB_TYPE': executable_subconfig.get('staticb_type','gsibec'),
             'OBS_YAML_DIR': executable_subconfig['obs_yaml_dir'],
             'OBS_LIST': executable_subconfig['obs_list'],
             'atm': executable_subconfig.get('atm', False),

ush/ufsda/stage.py

@@ -32,9 +32,17 @@ def gdas_fix(input_fix_dir, working_dir, config):
     else:
         case_anl = config['CASE_ANL']
     layers = int(config['LEVS'])-1
-    # link static B files
-    ufsda.disk_utils.symlink(os.path.join(input_fix_dir, 'bump', case_anl),
+
+    # figure out staticb source
+    if config['STATICB_TYPE'] == 'bump':
+        staticb_source = 'bump'
+    else:
+        staticb_source = 'gsibec'
+
+    # link staticb
+    ufsda.disk_utils.symlink(os.path.join(input_fix_dir, staticb_source, case_anl),
                              config['fv3jedi_staticb_dir'])
+
     # link akbk file
     ufsda.disk_utils.symlink(os.path.join(input_fix_dir, 'fv3jedi',
                                           'fv3files', f"akbk{layers}.nc4"),