Add offline tool to compute delp increments (NOAA-EMC#520)

This PR adds an offline python tool to produce delp increments from
increments of ps. It also adds the infrastructure to add delp as an
analysis variable to the fv3-jedi ctests, but that is currently turned
off since we want to use an offline tool to do this. There are also
various other consistency and formatting changes. It was previously
found that using ps as control variable and analysis variable did not
diagnostically update delp. Delp increments are needed for model
restart.

Author: delippi

parm/jcb-rdas/model/fv3/atmosphere/atmosphere_background_error_3denvar_mgbf.yaml.j2

@@ -19,13 +19,12 @@ members from template:
       soilm: {{soilm}}
       eastward_wind_at_surface: {{eastward_wind_at_surface}}
       northward_wind_at_surface: {{northward_wind_at_surface}}
-      equivalent_reflectivity_factor: {{equivalent_reflectivity_factor}}
       rain_water: {{rain_water}}
       snow_water: {{snow_water}}
       graupel: {{graupel}}
       upward_air_velocity: {{upward_air_velocity}}
+      equivalent_reflectivity_factor: {{equivalent_reflectivity_factor}}
     datapath: {{ atmosphere_background_ensemble_path }}
-    #filename is datetime templated: true
     filename_core: {{atmosphere_background_time_prefix}}{{filename_core}}
     filename_trcr: {{atmosphere_background_time_prefix}}{{filename_trcr}}
     filename_phys: {{atmosphere_background_time_prefix}}{{filename_phys}}

parm/jcb-rdas/model/fv3/atmosphere/atmosphere_background_error_3dvar_gsibec.yaml.j2

@@ -1,6 +1,7 @@
 covariance model: SABER
 saber central block:
   saber block name: gsi static covariance
+  active variables: {{control_variables}}
   read:
     gsi akbk: fv3_akbk
     gsi error covariance file: berror_stats

parm/jcb-rdas/model/fv3/atmosphere/atmosphere_background_error_hybrid3denvar_bump.yaml.j2

@@ -4,6 +4,7 @@ components:
     covariance model: SABER
     saber central block:
       saber block name: gsi static covariance
+      active variables: {{control_variables}}
       read:
         gsi akbk: fv3_akbk
         gsi error covariance file: berror_stats
@@ -54,6 +55,7 @@ components:
           northward_wind: {{northward_wind}}
           air_temperature: {{air_temperature}}
           air_pressure_thickness: {{air_pressure_thickness}}
+          air_pressure_at_surface: {{air_pressure_at_surface}}
           water_vapor_mixing_ratio_wrt_moist_air: {{water_vapor_mixing_ratio_wrt_moist_air}}
           cloud_liquid_ice: {{cloud_liquid_ice}}
           cloud_liquid_water: {{cloud_liquid_water}}
@@ -77,7 +79,7 @@ components:
         filename_sfcw: {{atmosphere_background_time_prefix}}{{filename_sfcw}}
         filename_cplr: {{atmosphere_background_time_prefix}}{{filename_cplr}}
       pattern: "%mem%"
-      nmembers: 30
+      nmembers: {{atmosphere_number_ensemble_members}}
       zero padding: 3
     localization:
       localization method: SABER

parm/jcb-rdas/model/fv3/atmosphere/atmosphere_background_error_hybrid3denvar_mgbf.yaml.j2

@@ -4,6 +4,7 @@ components:
     covariance model: SABER
     saber central block:
       saber block name: gsi static covariance
+      active variables: {{control_variables}}
       read:
         gsi akbk: fv3_akbk
         gsi error covariance file: berror_stats
@@ -54,6 +55,7 @@ components:
           northward_wind: {{northward_wind}}
           air_temperature: {{air_temperature}}
           air_pressure_thickness: {{air_pressure_thickness}}
+          air_pressure_at_surface: {{air_pressure_at_surface}}
           water_vapor_mixing_ratio_wrt_moist_air: {{water_vapor_mixing_ratio_wrt_moist_air}}
           cloud_liquid_ice: {{cloud_liquid_ice}}
           cloud_liquid_water: {{cloud_liquid_water}}
@@ -77,7 +79,7 @@ components:
         filename_sfcw: {{atmosphere_background_time_prefix}}{{filename_sfcw}}
         filename_cplr: {{atmosphere_background_time_prefix}}{{filename_cplr}}
       pattern: "%mem%"
-      nmembers: 30
+      nmembers: {{atmosphere_number_ensemble_members}}
       zero padding: 3
     localization:
       localization method: SABER
@@ -106,6 +108,6 @@ components:
          forward interpolator:
            local interpolator type: oops unstructured grid interpolator
          inverse interpolator:
-           local nterpolator type: oops unstructured grid interpolator
+           local interpolator type: oops unstructured grid interpolator
   weight:
     value: 0.85

parm/jcb-rdas/test/ci/jcb-rrfs_fv3jedi_2024052700_3denvar.yaml

@@ -35,15 +35,16 @@ analysis_variables:
 - northward_wind
 - air_temperature
 - water_vapor_mixing_ratio_wrt_moist_air
-- air_pressure_at_surface
 - ozone_mass_mixing_ratio
+#- air_pressure_thickness
 #- cloud_liquid_water
 #- cloud_liquid_ice
 #- rain_water
 #- snow_water
 #- graupel
 #- upward_air_velocity
 #- equivalent_reflectivity_factor
+- air_pressure_at_surface
 
 # Testing
 # -------

parm/jcb-rdas/test/ci/jcb-rrfs_fv3jedi_2024052700_3denvar_mgbf.yaml

@@ -35,15 +35,16 @@ analysis_variables:
 - northward_wind
 - air_temperature
 - water_vapor_mixing_ratio_wrt_moist_air
-- air_pressure_at_surface
 - ozone_mass_mixing_ratio
+#- air_pressure_thickness
 #- cloud_liquid_water
 #- cloud_liquid_ice
 #- rain_water
 #- snow_water
 #- graupel
 #- upward_air_velocity
 #- equivalent_reflectivity_factor
+- air_pressure_at_surface
 
 # Testing
 # -------
@@ -191,7 +192,7 @@ observations:
 - aircar_airTemperature_133
 - aircar_specificHumidity_133
 - aircar_winds_233
-#
+
 #- aircft_airTemperature_130
 #- aircft_airTemperature_131
 #- aircft_airTemperature_134

parm/jcb-rdas/test/ci/jcb-rrfs_fv3jedi_2024052700_3denvar_refl.yaml

@@ -35,15 +35,16 @@ analysis_variables:
 - northward_wind
 - air_temperature
 - water_vapor_mixing_ratio_wrt_moist_air
-- air_pressure_at_surface
 - ozone_mass_mixing_ratio
+#- air_pressure_thickness
 - cloud_liquid_water
 - cloud_liquid_ice
 - rain_water
 - snow_water
 - graupel
 - upward_air_velocity
 - equivalent_reflectivity_factor
+- air_pressure_at_surface
 
 # Testing
 # -------
@@ -190,7 +191,7 @@ observations:
 #- aircar_airTemperature_133
 #- aircar_specificHumidity_133
 #- aircar_winds_233
-#
+
 #- aircft_airTemperature_130
 #- aircft_airTemperature_131
 #- aircft_airTemperature_134

parm/jcb-rdas/test/ci/jcb-rrfs_fv3jedi_2024052700_3dvar.yaml

@@ -35,22 +35,22 @@ analysis_variables:
 - northward_wind
 - air_temperature
 - water_vapor_mixing_ratio_wrt_moist_air
-- air_pressure_at_surface
 - ozone_mass_mixing_ratio
+#- air_pressure_thickness
 #- cloud_liquid_water
 #- cloud_liquid_ice
 #- rain_water
 #- snow_water
 #- graupel
 #- upward_air_velocity
 #- equivalent_reflectivity_factor
+- air_pressure_at_surface
 
 control_variables:
 - eastward_wind
 - northward_wind
 - virtual_temperature
 - water_vapor_mixing_ratio_wrt_moist_air
-- air_pressure_at_surface
 - ozone_mass_mixing_ratio
 #- cloud_liquid_water
 #- cloud_liquid_ice
@@ -59,16 +59,17 @@ control_variables:
 #- graupel
 #- upward_air_velocity
 #- equivalent_reflectivity_factor
+- air_pressure_at_surface
 
 state_variables_to_inverse:
 - eastward_wind
 - northward_wind
 - air_temperature
 - water_vapor_mixing_ratio_wrt_moist_air
-- air_pressure_at_surface
-- air_pressure_thickness
-- geopotential_height_times_gravity_at_surface
 - ozone_mass_mixing_ratio
+- geopotential_height_times_gravity_at_surface
+- air_pressure_thickness
+- air_pressure_at_surface
 
 # Testing
 # -------

parm/jcb-rdas/test/ci/jcb-rrfs_fv3jedi_2024052700_3dvar_conv_surface.yaml

@@ -35,22 +35,22 @@ analysis_variables:
 - northward_wind
 - air_temperature
 - water_vapor_mixing_ratio_wrt_moist_air
-- air_pressure_at_surface
 - ozone_mass_mixing_ratio
+#- air_pressure_thickness
 #- cloud_liquid_water
 #- cloud_liquid_ice
 #- rain_water
 #- snow_water
 #- graupel
 #- upward_air_velocity
 #- equivalent_reflectivity_factor
+- air_pressure_at_surface
 
 control_variables:
 - eastward_wind
 - northward_wind
 - virtual_temperature
 - water_vapor_mixing_ratio_wrt_moist_air
-- air_pressure_at_surface
 - ozone_mass_mixing_ratio
 #- cloud_liquid_water
 #- cloud_liquid_ice
@@ -59,16 +59,17 @@ control_variables:
 #- graupel
 #- upward_air_velocity
 #- equivalent_reflectivity_factor
+- air_pressure_at_surface
 
 state_variables_to_inverse:
 - eastward_wind
 - northward_wind
 - air_temperature
 - water_vapor_mixing_ratio_wrt_moist_air
-- air_pressure_at_surface
-- air_pressure_thickness
-- geopotential_height_times_gravity_at_surface
 - ozone_mass_mixing_ratio
+- geopotential_height_times_gravity_at_surface
+- air_pressure_thickness
+- air_pressure_at_surface
 
 # Testing
 # -------

parm/jcb-rdas/test/ci/jcb-rrfs_fv3jedi_2024052700_3dvar_conv_upperair.yaml

@@ -35,22 +35,22 @@ analysis_variables:
 - northward_wind
 - air_temperature
 - water_vapor_mixing_ratio_wrt_moist_air
-- air_pressure_at_surface
 - ozone_mass_mixing_ratio
+#- air_pressure_thickness
 #- cloud_liquid_water
 #- cloud_liquid_ice
 #- rain_water
 #- snow_water
 #- graupel
 #- upward_air_velocity
 #- equivalent_reflectivity_factor
+- air_pressure_at_surface
 
 control_variables:
 - eastward_wind
 - northward_wind
 - virtual_temperature
 - water_vapor_mixing_ratio_wrt_moist_air
-- air_pressure_at_surface
 - ozone_mass_mixing_ratio
 #- cloud_liquid_water
 #- cloud_liquid_ice
@@ -59,16 +59,17 @@ control_variables:
 #- graupel
 #- upward_air_velocity
 #- equivalent_reflectivity_factor
+- air_pressure_at_surface
 
 state_variables_to_inverse:
 - eastward_wind
 - northward_wind
 - air_temperature
 - water_vapor_mixing_ratio_wrt_moist_air
-- air_pressure_at_surface
-- air_pressure_thickness
-- geopotential_height_times_gravity_at_surface
 - ozone_mass_mixing_ratio
+- geopotential_height_times_gravity_at_surface
+- air_pressure_thickness
+- air_pressure_at_surface
 
 # Testing
 # -------