 
      CALL SEQUENCE...
       hdplot  in input output \
                   [ -t1 VAL1  -t2 VAL2  -t3    VAL3  -name VAL4 -lv1 VAL5 \ 
                     -lx VAL6  -ly VAL7  -scal  VAL8  -a    VAL9 -b  VAL10 \ 
                     -c  VAL11 -d  VAL12 -label VAL13 -kind VAL14 ]
 
      DESCRIPTION...
       HARMONIC DIAL PLOT
 
 
      AUTHOR  - R.Laprise, MAY 12/83.
      LAST REVISION: $Header: /data/armnraid1/cvs/RMNLIB/utils/rdiag/lspgm/hdplot.ptn,v 4.7 2010/05/06 16:24:37 dugas Exp $
 
 
      PURPOSE - PLOTS A POLAR PLOT (HARMONIC DIAL) OF ONE SELECTED WAVENUMBER
                OR COMPONENT FROM A TIME SERIES OF COMPLEX FIELDS USING THE 
                NCAR PLOT PACKAGE.
 
      NOTE    - PROGRAM READS DATA LABELLED WITH STANDARD 8 WORD LABEL.
                THE MAXIMUM NUMBER OF POINTS IS 1000 .
 
 
      INPUT FILE...
       IN = FILE CONTAINING TIME SERIES OF COMPLEX DATA.
 
 
      INPUT PARAMETERS...
       NT1,NT2       =  FIRST AND LAST STEP NUMBERS IN THE TIME SERIES
       (t1,t2)         (-1 OR BLANK INPUT DEFAULTS TO ALL)
       INT  (t3)     =  GRAPHING INTERVAL (IN TIME) (DEFAULTS TO 1)
       NAME (name)   =  NAME  OF FIELD TO BE PLOTTED (BLANK DEFAULTS TO NEXT)
       LVL  (lv1)    =  LEVEL OF FIELD TO BE PLOTTED (BLANK IMPLIES DON'T CARE)
       MW   (lx)     =  THE FOURIER WAVE NUMBER (DEFAULTS TO 1)
       NW   (ly)     = (IF KIND= 'SPEC'), THE ORDER OF THE SPHERICAL HARMONIC
                       (IF KIND= 'FOUR'), THE LATITUDE NUMBER (FROM S.P. TO N.P)
                       (IF KIND= 'CMPL'), SAME AS 'FOUR'
       SCAL (scal)   =  SCALE FACTOR FOR THE AMPLITUDES
                        SCALING IS AUTOMATIC IF XLO=0.
       XLO,XHI (a,b) =  X PLOT LIMITS
       YLO,YHI (c,d) =  Y PLOT LIMITS
       KIND          =  VERTICAL LEVEL TYPE. IT CAN ONLY BE USED WITH CMC/RPN STD
                        FILES AND THEN LV1 AND LV2 ARE ASSUMED TO CONTAIN DECODED
                        VALUES. FURTHERMORE, THIS PARAMETRE CAN ONLY BE SPECIFIED
                        IN COMMAND-LINE MODE. THE CURRENTLY SUPPORTED VALUES ARE:
                        'HAUTEUR'   , heights (m) with respect to sea level
                        'SIGMA'     , classical terrain following (0.0 -> 1.0)
                        'ETA'       , first version of the GEM ETA coordinate
                        'PRESSION'  , pressure (mb)
                        'ARBITRAIRE', arbitrary level (no units apply)
                        'HAUTEURS'  , heights (M) with respect to the surface
                        'HYBRIDE'   , hybrid GEM coordinate
                        'THETA'     , potential temperature coordinate
                        'GALCHEN'   , GalChen generalized coordinate
 
       LABEL   = 80 CHARACTER LABEL PRINTED ON THE PLOT
 
      READ(5,5010) NT1,NT2,INT,NAME,LVL,MW,NW,SCAL,XLO,XHI,YLO,YHI
      READ(5,5020) LABEL
 5010 FORMAT(10X,2I10,I5,1X,A4,I5,1X,2I2,5F6.0)
 5020 FORMAT(10A8)
 
 
     EXAMPLE OF INPUT CARDS...
 
    HDPLOT         0        96    1 TEMP  500  6 6    1.  -1.0  +1.0  -1
   HARMONIC DIAL OF TEMPERATURE WAVE
 
         0123456789012345678901234567890123456789012345678901234567890
         1         2         3         4         5         6         7
 
 
      EXIT CONDITIONS...
        0 NO PROBLEMS
        1 UNABLE TO READ FILE IN
        2 TOO MANY SAMPLE RECORDS
        3 DATA IN FILE IN IS NOT COMPLEX
        4 SPEC TYPE: ILLEGAL NW,MW VALUES
        5 SPEC TYPE: FIELD TO SMALL
        6 FOUR TYPE: FIELD TO SMALL
        7 CMPL TYPE: FIELD TO SMALL
