These are the full instructions for converting the NHC wind probability to RPN format
and merging with the CMC 21-member dynamic ensemble. The procedure begins on whxlab with
the grib2 converter (requires linux) then on pollux which extracts the whxlab file and
does the merge using armnlib and r.diag functions. This script was set up in Nov 2009.

1. Copy the storm grib2 file of interest from the NHC ftp site at 
ftp://ftp.tpc.ncep.noaa.gov/wsp/ to the storm directory you've set up on the Z-drive
under "Knabb_Wind_Project"

2. Transfer that file to the /net/whxlab/home/fogartyc/Wind_Prob_Project directory on 
whxlab.

3. Run the converter script "NHC_to_RPN.sc" or "NHC_to_RPN_newdir.sc" if you want to recompile
the grib2 converter. The script creates "NHC_probs.fst".

4. Now log onto pollux system, go to directory WORK/GEM_extraction/GEM_ensemble_wind_project/past_cases
and run xrarc to extract the historical case of interest using "Ensemble_past_cases" directive file
within xrarc itself (also see "directives_for_xrarc.txt" in the "past_cases" directory)

5. Change the "date" and "heure" in the xrarc directives then execute. There should be 126 files 
that are being accessed.

6. After the download from xrarc, run the "get_ensemble_winds.sc" script.  Make sure to modify the
RUNSTAMP and HR variables. If this is a real-time or recent event, run "get_ensemble_winds.sc"
from the regular "/users/dor/armn/cfo/WORK/GEM_extraction/GEM_ensemble_wind_project" directory.

7. Copy the 120-hr, 34-kt wind probability file to the "past_cases" directory with name format
"${SNAME}_cmltv_prob_34kt_120" where $SNAME is set in "Ensemble_merger.sc".

8. Run "Ensemble_merger.sc".  This will generate the final output file "${SNAME}_merged_120HR_34kt.fst"