Add errtable comparison tool and align selected wind obs errors with GSI (NOAA-EMC#571)

* Adds a new `check_j2_errtable.py` utility for comparing JEDI YAML
`ObsErrorModelStepwiseLinear` pressure-dependent error tables against
`errtable.rrfs`.
* Updates selected SATWND YAMLs for types 245, 246, and 247 and
satellite identifiers 270/272 to align error handling and QC behavior
with GSI.
* Adds `use_error_inflation` Jinja guards around the SATWND error
inflation blocks so error inflation can be toggled consistently for
diagnostic experiments.
* Moves the SATWND `0.5` observation-error scaling earlier in the filter
sequence, alongside the other error inflation steps, rather than
applying it later after additional QC filters.
* Relaxes SATWND wind ObsValue bounds from `[-130, 130]` to `[-200,
200]`.
* Updates SATWND type 246 pressure bounds from `150-400 mb` to `150-300
mb`.
* Updates SATWND type 247 initial errors to use the large GSI-style
placeholder error values from `errtable.rrfs`.
* Updates `sfcshp_winds_284.yaml.j2` initial wind errors from `1.079` to
`1.5874`.
* Updates affected FV3-JEDI 3DVar reference test outputs.

Author: delippi

parm/jcb-rdas/observations/atmosphere/satwnd_winds_245_270.yaml.j2

@@ -140,6 +140,7 @@
                  xvals: [110000, 105000, 100000, 95000, 90000, 85000, 80000, 75000, 70000, 65000, 60000, 55000, 50000, 45000, 40000, 35000, 30000, 25000, 20000, 15000, 10000, 7500, 5000, 4000, 3000, 2000, 1000, 500, 400, 300, 200, 100, 0]
                  errors: [3.8, 3.8, 3.8, 3.8, 3.8, 3.8, 3.8, 3.8, 3.9, 3.9, 4.0, 4.0, 4.0, 4.0, 4.0, 4.0, 4.0, 4.0, 4.0, 4.0, 4.0, 4.0, 4.0, 4.0, 4.0, 4.0, 4.0, 4.0, 4.0, 4.0, 4.0, 4.0, 4.0]
 
+{% if use_error_inflation | default(true) %}  # BEGIN use_error_inflation
          # Error inflation (windEastward) based on pressure check (setupw.f90)
          - filter: Perform Action
            udescriptor: "error_inflation_pressure_check"
@@ -176,14 +177,30 @@
                  variable: windNorthward
                  inflation factor: 4.0
 
+         # Halve the ob error (based on read_satwnd.f90)
+         - filter: Perform Action
+           udescriptor: "ob_error_inflation"
+           filter variables:
+           - name: windEastward
+           - name: windNorthward
+           where:
+           - variable: ObsType/windEastward
+             is_in: 245
+           - variable: MetaData/satelliteIdentifier
+             is_in: 270
+           action:
+             name: inflate error
+             inflation factor: 0.5
+{% endif %}  # END use_error_inflation
+
          # Bounds Check (ObsValue)
          - filter: Bounds Check
            udescriptor: "bounds_check_ObsValue"
            filter variables:
            - name: windEastward
            - name: windNorthward
-           minvalue: -130
-           maxvalue: 130
+           minvalue: -200
+           maxvalue: 200
            where:
            - variable: ObsType/windEastward
              is_in: 245
@@ -332,21 +349,6 @@
            action:
              name: reject
 
-         # Halve the ob error (based on read_satwnd.f90)
-         - filter: Perform Action
-           udescriptor: "ob_error_inflation"
-           filter variables:
-           - name: windEastward
-           - name: windNorthward
-           where:
-           - variable: ObsType/windEastward
-             is_in: 245
-           - variable: MetaData/satelliteIdentifier
-             is_in: 270
-           action:
-             name: inflate error
-             inflation factor: 0.5
-
          # Reject when pressure is more than 50 mb above tropopause
          - filter: Difference Check
            udescriptor: "difference_check_tropopause_pressure"

parm/jcb-rdas/observations/atmosphere/satwnd_winds_245_272.yaml.j2

@@ -140,6 +140,7 @@
                  xvals: [110000, 105000, 100000, 95000, 90000, 85000, 80000, 75000, 70000, 65000, 60000, 55000, 50000, 45000, 40000, 35000, 30000, 25000, 20000, 15000, 10000, 7500, 5000, 4000, 3000, 2000, 1000, 500, 400, 300, 200, 100, 0]
                  errors: [3.8, 3.8, 3.8, 3.8, 3.8, 3.8, 3.8, 3.8, 3.9, 3.9, 4.0, 4.0, 4.0, 4.0, 4.0, 4.0, 4.0, 4.0, 4.0, 4.0, 4.0, 4.0, 4.0, 4.0, 4.0, 4.0, 4.0, 4.0, 4.0, 4.0, 4.0, 4.0, 4.0]
 
+{% if use_error_inflation | default(true) %}  # BEGIN use_error_inflation
          # Error inflation (windEastward) based on pressure check (setupw.f90)
          - filter: Perform Action
            udescriptor: "error_inflation_pressure_check"
@@ -176,14 +177,30 @@
                  variable: windNorthward
                  inflation factor: 4.0
 
+         # Halve the ob error (based on read_satwnd.f90)
+         - filter: Perform Action
+           udescriptor: "ob_error_inflation"
+           filter variables:
+           - name: windEastward
+           - name: windNorthward
+           where:
+           - variable: ObsType/windEastward
+             is_in: 245
+           - variable: MetaData/satelliteIdentifier
+             is_in: 272
+           action:
+             name: inflate error
+             inflation factor: 0.5
+{% endif %}  # END use_error_inflation
+
          # Bounds Check (ObsValue)
          - filter: Bounds Check
            udescriptor: "bounds_check_ObsValue"
            filter variables:
            - name: windEastward
            - name: windNorthward
-           minvalue: -130
-           maxvalue: 130
+           minvalue: -200
+           maxvalue: 200
            where:
            - variable: ObsType/windEastward
              is_in: 245
@@ -332,21 +349,6 @@
            action:
              name: reject
 
-         # Halve the ob error (based on read_satwnd.f90)
-         - filter: Perform Action
-           udescriptor: "ob_error_inflation"
-           filter variables:
-           - name: windEastward
-           - name: windNorthward
-           where:
-           - variable: ObsType/windEastward
-             is_in: 245
-           - variable: MetaData/satelliteIdentifier
-             is_in: 272
-           action:
-             name: inflate error
-             inflation factor: 0.5
-
          # Reject when pressure is more than 50 mb above tropopause
          - filter: Difference Check
            udescriptor: "difference_check_tropopause_pressure"

parm/jcb-rdas/observations/atmosphere/satwnd_winds_246_270.yaml.j2

@@ -140,6 +140,7 @@
                  xvals: [110000, 105000, 100000, 95000, 90000, 85000, 80000, 75000, 70000, 65000, 60000, 55000, 50000, 45000, 40000, 35000, 30000, 25000, 20000, 15000, 10000, 7500, 5000, 4000, 3000, 2000, 1000, 500, 400, 300, 200, 100, 0]
                  errors: [3.8, 3.8, 3.8, 3.8, 3.8, 3.8, 3.8, 3.8, 3.9, 3.9, 4.0, 4.0, 4.1, 5.0, 6.0, 6.3, 6.6, 7.0, 7.0, 7.0, 7.0, 7.0, 7.0, 7.0, 7.0, 7.0, 7.0, 7.0, 7.0, 7.0, 7.0, 7.0, 7.0]
 
+{% if use_error_inflation | default(true) %}  # BEGIN use_error_inflation
          # Error inflation (windEastward) based on pressure check (setupw.f90)
          - filter: Perform Action
            udescriptor: "error_inflation_pressure_check"
@@ -176,14 +177,30 @@
                  variable: windNorthward
                  inflation factor: 4.0
 
+         # Halve the ob error (based on read_satwnd.f90)
+         - filter: Perform Action
+           udescriptor: "ob_error_inflation"
+           filter variables:
+           - name: windEastward
+           - name: windNorthward
+           where:
+           - variable: ObsType/windEastward
+             is_in: 246
+           - variable: MetaData/satelliteIdentifier
+             is_in: 270
+           action:
+             name: inflate error
+             inflation factor: 0.5
+{% endif %}  # END use_error_inflation
+
          # Bounds Check (ObsValue)
          - filter: Bounds Check
            udescriptor: "bounds_check_ObsValue"
            filter variables:
            - name: windEastward
            - name: windNorthward
-           minvalue: -130
-           maxvalue: 130
+           minvalue: -200
+           maxvalue: 200
            where:
            - variable: ObsType/windEastward
              is_in: 246
@@ -227,7 +244,7 @@
            action:
              name: reject
 
-         # Reject obs with pressure < 150 mb and  > 400 mb
+         # Reject obs with pressure < 150 mb and  > 300 mb
          - filter: Bounds Check
            udescriptor: "bounds_check_pressure"
            filter variables:
@@ -236,7 +253,7 @@
            test variables:
            - name: MetaData/pressure
            minvalue: 15000.
-           maxvalue: 40000.
+           maxvalue: 30000.
            where:
            - variable: ObsType/windEastward
              is_in: 246
@@ -315,21 +332,6 @@
            action:
              name: reject
 
-         # Halve the ob error (based on read_satwnd.f90)
-         - filter: Perform Action
-           udescriptor: "ob_error_inflation"
-           filter variables:
-           - name: windEastward
-           - name: windNorthward
-           where:
-           - variable: ObsType/windEastward
-             is_in: 246
-           - variable: MetaData/satelliteIdentifier
-             is_in: 270
-           action:
-             name: inflate error
-             inflation factor: 0.5
-
          # Reject when pressure is more than 50 mb above tropopause
          - filter: Difference Check
            udescriptor: "difference_check_tropopause_pressure"

parm/jcb-rdas/observations/atmosphere/satwnd_winds_246_272.yaml.j2

@@ -140,6 +140,7 @@
                  xvals: [110000, 105000, 100000, 95000, 90000, 85000, 80000, 75000, 70000, 65000, 60000, 55000, 50000, 45000, 40000, 35000, 30000, 25000, 20000, 15000, 10000, 7500, 5000, 4000, 3000, 2000, 1000, 500, 400, 300, 200, 100, 0]
                  errors: [3.8, 3.8, 3.8, 3.8, 3.8, 3.8, 3.8, 3.8, 3.9, 3.9, 4.0, 4.0, 4.1, 5.0, 6.0, 6.3, 6.6, 7.0, 7.0, 7.0, 7.0, 7.0, 7.0, 7.0, 7.0, 7.0, 7.0, 7.0, 7.0, 7.0, 7.0, 7.0, 7.0]
 
+{% if use_error_inflation | default(true) %}  # BEGIN use_error_inflation
          # Error inflation (windEastward) based on pressure check (setupw.f90)
          - filter: Perform Action
            udescriptor: "error_inflation_pressure_check"
@@ -176,14 +177,30 @@
                  variable: windNorthward
                  inflation factor: 4.0
 
+         # Halve the ob error (based on read_satwnd.f90)
+         - filter: Perform Action
+           udescriptor: "ob_error_inflation"
+           filter variables:
+           - name: windEastward
+           - name: windNorthward
+           where:
+           - variable: ObsType/windEastward
+             is_in: 246
+           - variable: MetaData/satelliteIdentifier
+             is_in: 272
+           action:
+             name: inflate error
+             inflation factor: 0.5
+{% endif %}  # END use_error_inflation
+
          # Bounds Check (ObsValue)
          - filter: Bounds Check
            udescriptor: "bounds_check_ObsValue"
            filter variables:
            - name: windEastward
            - name: windNorthward
-           minvalue: -130
-           maxvalue: 130
+           minvalue: -200
+           maxvalue: 200
            where:
            - variable: ObsType/windEastward
              is_in: 246
@@ -227,7 +244,7 @@
            action:
              name: reject
 
-         # Reject obs with pressure < 150 mb and  > 400 mb
+         # Reject obs with pressure < 150 mb and  > 300 mb
          - filter: Bounds Check
            udescriptor: "bounds_check_pressure"
            filter variables:
@@ -236,7 +253,7 @@
            test variables:
            - name: MetaData/pressure
            minvalue: 15000.
-           maxvalue: 40000.
+           maxvalue: 30000.
            where:
            - variable: ObsType/windEastward
              is_in: 246
@@ -315,21 +332,6 @@
            action:
              name: reject
 
-         # Halve the ob error (based on read_satwnd.f90)
-         - filter: Perform Action
-           udescriptor: "ob_error_inflation"
-           filter variables:
-           - name: windEastward
-           - name: windNorthward
-           where:
-           - variable: ObsType/windEastward
-             is_in: 246
-           - variable: MetaData/satelliteIdentifier
-             is_in: 272
-           action:
-             name: inflate error
-             inflation factor: 0.5
-
          # Reject when pressure is more than 50 mb above tropopause
          - filter: Difference Check
            udescriptor: "difference_check_tropopause_pressure"

parm/jcb-rdas/observations/atmosphere/satwnd_winds_247_270.yaml.j2

@@ -138,8 +138,10 @@
                  xvar:
                    name: MetaData/pressure
                  xvals: [110000, 105000, 100000, 95000, 90000, 85000, 80000, 75000, 70000, 65000, 60000, 55000, 50000, 45000, 40000, 35000, 30000, 25000, 20000, 15000, 10000, 7500, 5000, 4000, 3000, 2000, 1000, 500, 400, 300, 200, 100, 0]
-                 errors: [3.8, 3.8, 3.8, 3.8, 3.8, 3.8, 3.8, 3.8, 3.9, 3.9, 4.0, 4.0, 4.1, 5.0, 6.0, 6.3, 6.6, 7.0, 7.0, 7.0, 7.0, 7.0, 7.0, 7.0, 7.0, 7.0, 7.0, 7.0, 7.0, 7.0, 7.0, 7.0, 7.0]
+                 errors: [1000000000.0, 1000000000.0, 1000000000.0, 1000000000.0, 1000000000.0, 1000000000.0, 1000000000.0, 1000000000.0, 1000000000.0, 1000000000.0, 1000000000.0, 1000000000.0, 1000000000.0, 1000000000.0, 1000000000.0, 1000000000.0, 1000000000.0, 1000000000.0, 1000000000.0, 1000000000.0, 1000000000.0, 1000000000.0, 1000000000.0, 1000000000.0, 1000000000.0, 1000000000.0, 1000000000.0, 1000000000.0, 1000000000.0, 1000000000.0, 1000000000.0, 1000000000.0, 1000000000.0]
+                 #errors: [3.8, 3.8, 3.8, 3.8, 3.8, 3.8, 3.8, 3.8, 3.9, 3.9, 4.0, 4.0, 4.1, 5.0, 6.0, 6.3, 6.6, 7.0, 7.0, 7.0, 7.0, 7.0, 7.0, 7.0, 7.0, 7.0, 7.0, 7.0, 7.0, 7.0, 7.0, 7.0, 7.0]
 
+{% if use_error_inflation | default(true) %}  # BEGIN use_error_inflation
          # Error inflation (windEastward) based on pressure check (setupw.f90)
          - filter: Perform Action
            udescriptor: "error_inflation_pressure_check"
@@ -176,14 +178,30 @@
                  variable: windNorthward
                  inflation factor: 4.0
 
+         # Halve the ob error (based on read_satwnd.f90)
+         - filter: Perform Action
+           udescriptor: "ob_error_inflation"
+           filter variables:
+           - name: windEastward
+           - name: windNorthward
+           where:
+           - variable: ObsType/windEastward
+             is_in: 247
+           - variable: MetaData/satelliteIdentifier
+             is_in: 270
+           action:
+             name: inflate error
+             inflation factor: 0.5
+{% endif %}  # END use_error_inflation
+
          # Bounds Check (ObsValue)
          - filter: Bounds Check
            udescriptor: "bounds_check_ObsValue"
            filter variables:
            - name: windEastward
            - name: windNorthward
-           minvalue: -130
-           maxvalue: 130
+           minvalue: -200
+           maxvalue: 200
            where:
            - variable: ObsType/windEastward
              is_in: 247
@@ -296,21 +314,6 @@
            action:
              name: reject
 
-         # Halve the ob error (based on read_satwnd.f90)
-         - filter: Perform Action
-           udescriptor: "ob_error_inflation"
-           filter variables:
-           - name: windEastward
-           - name: windNorthward
-           where:
-           - variable: ObsType/windEastward
-             is_in: 247
-           - variable: MetaData/satelliteIdentifier
-             is_in: 270
-           action:
-             name: inflate error
-             inflation factor: 0.5
-
          # Reject when pressure is more than 50 mb above tropopause
          - filter: Difference Check
            udescriptor: "difference_check_tropopause_pressure"