[Dart-dev] DART two_experiments comparison

[Tim Hoar]

Just to put an example behind the email announcement of the Matlab  
script that
compares 'two' experments - I'll include the calling sequence and an  
example plot
that shows the results of two of Hui's ernesto experiments:

orient tall,
wysiwyg;
files = {'/fs/image/home/hliu/DART/models/wrf/work/ernesto/ctl/ 
obs_diag_output.nc',
          '/fs/image/home/hliu/DART/models/wrf/work/ernesto/airsQ/ 
obs_diag_output.nc'};
titles   = {'Control', 'airsQ'};
varnames = {'RADIOSONDE_U_WIND_COMPONENT', 'RADIOSONDE_TEMPERATURE'};
qtty     = 'rmse';     % rmse, spread, totalspread, bias, etc.
prpo     = 'both';     % prior & posterior

two_experiments_profile(files, titles, varnames, qtty, prpo)

Please make note that the 'region' string (in this case 'Northern  
Hemisphere')
comes from whatever was used when 'obs_diag' was run. It this case, I  
would
have declared something like:

&obs_diag_nml
    reg_names  = 'Ernesto',




Here's the enticement to keep reading. I called it 'two_experiments',  
but in reality,
it is N_experiments. As things get more cluttered, you can elect to  
plot only the
priors (or posteriors, if you want to incur the Wrath of Jeff). I put  
a lot of effort
into improving the reliability and labeling of the two axes (nobs, and  
'rmse').
Each experiment has its own nobs 'possible' and 'used' that are the  
same color as
the corresponding line.

I am in the process of putting a little watermark on the bottom with  
the filenames
that are being plotted. At present, I'm not automatically printing  
these things.
I might add an option ...

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-dev/attachments/20100405/7fa118ce/attachment-0001.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: WRF_example.png
Type: image/png
Size: 20148 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/dart-dev/attachments/20100405/7fa118ce/attachment-0001.png