[Dart-dev] [4168] DART/trunk/diagnostics: The obs_diag program now has the following ' copies':

[nancy at ucar.edu]

Revision: 4168
Author:   thoar
Date:     2009-11-25 15:21:03 -0700 (Wed, 25 Nov 2009)
Log Message:
-----------
The obs_diag program now has the following 'copies':
Nposs
Nused
NbigQC
NbadIZ
NbadUV
NbadLV
rmse
bias
spread
totalspread
NbadDARTQC
observation
ens_mean
N_DARTqc_0
N_DARTqc_1
N_DARTqc_2
N_DARTqc_3
N_DARTqc_4
N_DARTqc_5
N_DARTqc_6
N_DARTqc_7

which can be plotted with the usual diagnostics/matlab routines.
So - now it is possible to not only plot the bias as a function of height,
you can now plot the (mean) observation value as a function of height.
Additionally, obs_diag and the associated plotting functions
also support observations on model levels (curious obs).

Modified Paths:
--------------
    DART/trunk/diagnostics/matlab/plot_bias_xxx_profile.m
    DART/trunk/diagnostics/matlab/plot_profile.m
    DART/trunk/diagnostics/matlab/plot_rmse_xxx_profile.m
    DART/trunk/diagnostics/threed_sphere/obs_diag.f90

-------------- next part --------------
Modified: DART/trunk/diagnostics/matlab/plot_bias_xxx_profile.m
===================================================================
--- DART/trunk/diagnostics/matlab/plot_bias_xxx_profile.m	2009-11-25 21:36:59 UTC (rev 4167)
+++ DART/trunk/diagnostics/matlab/plot_bias_xxx_profile.m	2009-11-25 22:21:03 UTC (rev 4168)
@@ -284,6 +284,7 @@
    h = legend(h1, str_bias_pr, str_bias_po, str_other_pr, str_other_po, ...
           'Location','East');
    legend(h,'boxoff')
+   set(h,'Interpreter','none')
 
    axlims = [plotdat.Xrange plotdat.Yrange];
    axis(axlims)
@@ -462,6 +463,10 @@
       ymax =  ceil(max(glommed));
    end
 
+   if (ymin == 0 && ymax == 0)
+       ymax = 1;
+   end
+   
    if (ymin == ymax)
      ymin = ymin - 0.1*ymin;
      ymax = ymax + 0.1*ymax;

Modified: DART/trunk/diagnostics/matlab/plot_profile.m
===================================================================
--- DART/trunk/diagnostics/matlab/plot_profile.m	2009-11-25 21:36:59 UTC (rev 4167)
+++ DART/trunk/diagnostics/matlab/plot_profile.m	2009-11-25 22:21:03 UTC (rev 4168)
@@ -269,6 +269,7 @@
    set(h1,'LineWidth',plotdat.linewidth);
    h = legend(h1,str_other_pr, str_other_po, 'Location', 'East');
    legend(h,'boxoff')
+   set(h,'Interpreter','none')
 
    axlims = [plotdat.Xrange plotdat.Yrange];
    axis(axlims)
@@ -452,6 +453,10 @@
       ymax =  ceil(max(glommed));
    end
 
+   if (ymin == 0 && ymax == 0)
+       ymax = 1;
+   end
+   
    if (ymin == ymax)
      ymin = ymin - 0.1*ymin;
      ymax = ymax + 0.1*ymax;

Modified: DART/trunk/diagnostics/matlab/plot_rmse_xxx_profile.m
===================================================================
--- DART/trunk/diagnostics/matlab/plot_rmse_xxx_profile.m	2009-11-25 21:36:59 UTC (rev 4167)
+++ DART/trunk/diagnostics/matlab/plot_rmse_xxx_profile.m	2009-11-25 22:21:03 UTC (rev 4168)
@@ -284,6 +284,7 @@
    h = legend(h1, str_rmse_pr, str_rmse_po, str_other_pr, str_other_po, ...
           'Location','East');
    legend(h,'boxoff')
+   set(h,'Interpreter','none')
 
    axlims = [plotdat.Xrange plotdat.Yrange];
    axis(axlims)
@@ -462,6 +463,10 @@
       ymax =  ceil(max(glommed));
    end
 
+   if (ymin == 0 && ymax == 0)
+       ymax = 1;
+   end
+   
    if (ymin == ymax)
      ymin = ymin - 0.1*ymin;
      ymax = ymax + 0.1*ymax;

Modified: DART/trunk/diagnostics/threed_sphere/obs_diag.f90
===================================================================
--- DART/trunk/diagnostics/threed_sphere/obs_diag.f90	2009-11-25 21:36:59 UTC (rev 4167)
+++ DART/trunk/diagnostics/threed_sphere/obs_diag.f90	2009-11-25 22:21:03 UTC (rev 4168)
@@ -24,7 +24,6 @@
 ! 'priorspred' should really be 'priorvar' since you have to accumulate variances 
 ! the math is correct as it is, but the variable names don't make it easy ...
 
-
 use        types_mod, only : r4, r8, digits12, MISSING_R8, MISSING_R4
 use obs_sequence_mod, only : read_obs_seq, obs_type, obs_sequence_type, get_first_obs, &
                              get_obs_from_key, get_obs_def, get_copy_meta_data, &