[Dart-users] DART observation-space rank histogram ready for use.

[Tim Hoar]

I extended obs_diag.f90 to also output rank histogram information if
there are any 'prior ensemble member xxx' copies in the obs_seq.final  
files.
(They exist if you specify input.nml:filter_nml:num_output_obs_members  
== <ens_size>.)

After you create the obs_diag_output.nc file, you can run

DART/diagnostics/matlab/plot_rank_histogram.m

to view the results. It has the 'expected' help file - available from  
the matlab
command line.

At present, ALL observations with DART QC 0,1,2,3, and 7 are used to  
create the
rank histogram. 7 is the outlier observation, so even these are being  
'counted'.
The observation error -as specified in each observation- is used to  
generate
the sampling error signal for the ensemble members' expected  
observations that
define the bin edges for the rank histogram.

obs_diag.f90 has also been updated to provide better input guidance  
for the
instance when the user has specified a time window that does not  
intersect with
the observations in the files specified.

Please

svn update DART/diagnostics

rebuild your obs_diag, and let me know how I can improve it. I have  
one idea already,
so stay tuned for a minor tweak to the matlab script on Monday.

Tim Hoar, Associate Scientist
National Center for Atmospheric Research
thoar at ucar.edu
303 497 1708

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/dart-users/attachments/20101008/26e86777/attachment.html