The software uses innovation files to produce rank histograms and plots of skill (obs-fcst_mean) against spread (sd). One plot of each type is produced per cycle. Currently, the input data is obtained from innovation files.
The IDL source code can be found at
Note that a new program
The code uses library routines written by Andy Heaps (NCAS, University of
and these are accesssed by adding these two lines into your
export IDL_PATH=+/usr/local/itt . /data/nwp1/frbr/TIDLWorkspace/andy_lib/guide4.sh
The user must supply information as follows:
The program will then loop over the specified innovation files and generate (in the working directory) files with names such as N_groupname_variablename_hist.ps, where N is a simple counter, starting at 1 and incrementing by 1.
111 ; start cycle 112 ; end cycle /export/carrot/raid1/brugge/fcst_anal/ ; innovations file directory Innovations_qwq ; XXX - innovations filename XXXcccYYY .txt ; YYY - where ccc is the cycle number 3 ; selected group (eg 2 Surface, 3 Aircraft, 4 Satwind, etc) 2 ; selected type (eg 2 Temperature or theta, 4 u, 5 v, etc) -999.9 ; selected level, or -999.9 for all levels
where the records should been changed to reflect the input data and output requirements of the user. The text to the right of the ';' symbols can be kept as a reminder of the meaning of the data records.
The idl code can now been run in a batch job or, if run interactively, will loop through the specified cycles.
The figures show the output from two cycles-worth of aircraft temperature observations.
The ability for the software to use precipitation data from a non-innovation file.