fixed dirac

NOAA-EMC · Jun 21, 2023 · d765863 · d765863
1 parent cfb5fb4
commit d765863
Show file tree

Hide file tree

Showing 2 changed files with 12 additions and 4 deletions.
diff --git a/test/soca/gw/CMakeLists.txt b/test/soca/gw/CMakeLists.txt
@@ -37,7 +37,7 @@ add_test(NAME test_gdasapp_soca_run_clean
 # Test JGDAS_GLOBAL_OCEAN_ANALYSIS_*
 set(jjob_list "JGDAS_GLOBAL_OCEAN_ANALYSIS_PREP"
               "JGDAS_GLOBAL_OCEAN_ANALYSIS_BMAT"
-#              "JGDAS_GLOBAL_OCEAN_ANALYSIS_BMAT_VRFY"
+              "JGDAS_GLOBAL_OCEAN_ANALYSIS_BMAT_VRFY"
               "JGDAS_GLOBAL_OCEAN_ANALYSIS_RUN"
               "JGDAS_GLOBAL_OCEAN_ANALYSIS_CHKPT"
               "JGDAS_GLOBAL_OCEAN_ANALYSIS_POST"

diff --git a/ush/ufsda/dirac_yaml.py b/ush/ufsda/dirac_yaml.py
@@ -34,9 +34,17 @@ def var2dirac(args):
     # Overwrite variables
     state_vars = ['cicen', 'hicen', 'hsnon', 'socn', 'tocn', 'ssh', 'hocn', 'mld', 'layer_depth']
     diracconfig['initial condition']['state variables'] = state_vars
-    for sub_bmat in range(len(diracconfig['background error']['components'])):
-        diracconfig['background error']['components'][sub_bmat]['covariance']['linear variable change']['input variables'] = state_vars
-        diracconfig['background error']['components'][sub_bmat]['covariance']['linear variable change']['output variables'] = state_vars
+
+    # str to integer
+    nmem=int(diracconfig['background error']['components'][1]['covariance']['members from template']['nmembers'] )
+    diracconfig['background error']['components'][1]['covariance']['members from template']['nmembers'] = nmem
+
+    # Set date to the middle of the window
+    window_middle = diracconfig['background error']['components'][1]['weight']['date']
+    diracconfig['initial condition']['date'] = window_middle
+
+    # Remove static B
+    diracconfig['background error'] = diracconfig['background error']['components'][1]['covariance']
 
     # Generate impulse indices
     ds = xr.open_dataset(grid)