Save the correct field for ssh_u (#622)

* save the correct field for ssh_u ... and stop 0ing it out ... * fixed typo
NOAA-EMC · Sep 13, 2023 · 516a5bd · 516a5bd
1 parent 2678ea1
commit 516a5bd
Show file tree

Hide file tree

Showing 4 changed files with 9 additions and 7 deletions.
diff --git a/parm/soca/berror/saber_blocks.yaml b/parm/soca/berror/saber_blocks.yaml
@@ -32,21 +32,23 @@ components:
         ice_filename: 'ice.bkgerr_stddev.incr.{{ATM_WINDOW_BEGIN}}.nc'
         date: '{{ATM_WINDOW_MIDDLE}}'
         t_min: 0.1
-        t_max: 10.0
+        t_max: 5.0
         s_min: 0.1
-        s_max: 10.0
+        s_max: 1.0
         ssh_min: 0.0 # std ssh=0 => ssh balance applied as
-        ssh_max: 0.0 #              strong constraint
+        ssh_max: 1.0 #              strong constraint
         cicen_min: 0.1
         cicen_max: 0.5
         hicen_min: 0.0
         hicen_max: 0.0
         standard deviation: true
 
       - linear variable change name: BalanceSOCA
+        ksshts:
+          nlayers: 10
 
   weight:
-    value: 1.0
+    value: 0.25
 - covariance:
     covariance model: ensemble
     members from template:

diff --git a/parm/soca/berror/soca_ensweights.yaml b/parm/soca/berror/soca_ensweights.yaml
@@ -18,7 +18,7 @@ background:
 weights:
   # Need to provide weights^2 when reading from file
   ice: 0.0025   #  5% of original variance
-  ocean: 0.01   # 10%      "       "
+  ocean: 0.0625   # 25%      "       "
 
 output:
   datadir: ./

diff --git a/scripts/exgdas_global_marine_analysis_vrfy.py b/scripts/exgdas_global_marine_analysis_vrfy.py
@@ -48,7 +48,7 @@
 config = plotConfig(grid_file=grid_file,
                     data_file=data_file,
                     variables_horiz={'ave_ssh': [-1, 1]},
-                    colormap='gist_ncar',
+                    colormap='seismic',
                     comout=os.path.join(comout, 'vrfy', 'recentering_error'))
 recErrPlotter = statePlotter(config)
 recErrPlotter.plot()

diff --git a/utils/soca/gdas_ens_handler.h b/utils/soca/gdas_ens_handler.h
@@ -231,7 +231,7 @@ namespace gdasapp {
       oops::Log::info() << "mean non-steric ssh: " << sshMean << std::endl;
       oops::Log::info() << "std non-steric ssh: " << sshNonStericStd << std::endl;
       eckit::LocalConfiguration nonStericSshOutputConfig(fullConfig, "ssh output.unbalanced");
-      sshStd.write(nonStericSshOutputConfig);
+      sshNonStericStd.write(nonStericSshOutputConfig);
 
       // Compute filtered ensemble moments
       ensMean.zero();