USAGE:
          sol2xyz [options] [-Oaddopts -Oaddopts ...] file

 PURPOSE:
          In the context of OTEQ solutions or Global Tide grids,
          prepares ascii-data for GMT (OTEQ/bin/solplot)
          Prints Z-array, amplitude or phase in x-y-z format.
          Files are expected in $dir, see -d option.

 Notation: #number $string compulsory with emphasis
 All file names can be written after a blankspace following the option

 OPTIONS:
          RETRIEVAL AND TIDE VALUES:

 -T$tide      - two-character symbol                                          [M2]
 -V#version   - OTEQ version number for retrieval.
                -1 for the first array in the file.                            [1]
 -A#region    - #W/#E/#S/#N                                                  [all]
 -t$target    - include node types                                   [all:  AaOoS]
 -s#scale     - amplitude scaling factor                                     [1.0]

 -p[#cut]     - write phase component,                                      [amps]
                      cut phase circle at#cut degrees                        [0.0]
 -Z[J]        - write complex values, J for conjugate                       [amps]
 -r#min,#max  - output values only within specified range                    [all]
 -q           - query array size                               [print grid values]

          PROCESSING:

 -agrad       - sum of east and north successive differences
 -egrad -ngrad
              - Not implemented yet
 -P#s,#n      - Virtually extend the grid scope to south and north.
                In global grids: complement array rows to cover
                the sphere completely.                                       [0,0]
 -X#w,#e      - Virtually extend the grid scope to west and east
                In global grids: complement array columns to cover
                the sphere completely; shift the column where array
                scan begins (modulus m and modulus 360. is always applied)
 -R[F]        - Reorient tide grid to north-south (F also the flag grid)
                Calls ZREVO (IREVO)                                         [S->N]
                For a rotation of array longitudes 

          COORDINATES IN OUTPUT:

 -C           - cartesian, units km, regional model,
 -CG[R][:#c]  - write cartesian global grid (degrees)
                R for reverse N-S orientation                               [S->N]
                Cut the longitudes at #c                                     [0.0]
             
 -S           - stereographic (deg), regional model,                  [grid units]
 -SG[R][:#c]  - Global models: Dummy projection on the sphere
                R for reverse N-S orientation                               [S->N]
                Cut the longitudes at #c                                   [180.0]
 -c#lon,#lat  - origin for dummy projection                              [0.0,0.0]
        
          FILES:

 -E           - use the Endian of the old HP-UX machine when
                reading the complex array

 -d$dir       - input file directory (add "/")                                [./]

 -fNONE       - don't include node flags in processing = the default for -f
 -fMAKE[:#mrs[,#acc]]
              - create flags from tide grid. Dry = (||z|-#mrs) < acc
                #mrs : critical value or symbolic value for dry nodes,       [9.0]
                #acc : detection accuracy                                    [0.1]
             
 -fSAME[,-A]  - read flags from tide file, ,-A to skip the area parameters
 -f[[$t]:#m,]$file     
              - use flag file                                           [FLZU.DAT]
                grid type $t (M or Z)                           [tide's grid type]
                rotate around y-axis by #m grid units                          [0]

 -a{+|-}      - get (+) or don't get (-) area from flag file   [context dependent]

 -F$file      - output flag file

 -z[$file]    - include o-o-area tide from file                      [OOTIDES.DAT]                                       
 -o$file      - output to named file                                      [STDOUT]
                (You'll have to pipe STDOUT through  | fgrep -v '>' )
 -oNONE       - no output; for testing input files

          ADDITIONAL OPTIONS: (must be coded last, but before the file name)

 -O$opts      - additional options into OPT parameter
                of print_zm subroutine, cumulative. See next paragraph.
        L:#dl - shift longitudes (really!) by dl                             [0.0]
           +P - print a protocol message '<PrinZm>>>...'
            G - global grid (circular in x-direction,
                area centre (doubtful defaults!)
                resolution derived from X-dimension)
    O:#x0,#y0 - distance of coordinate origin from grid origin

The options that the PRINT_ZM subroutine understands are:
 C G          - set by -C
 S[:#cc,]     - set by -S
 R            - set by -R
 Z[J]         - set by -Z

The following must be set explicitly with sol2xyz's option -O $code, :
 +P           - print a protocol message '<PrinZm>>>...'
 O:#x0,#y0    - centre point in global grid, longitude and latitude      [0.0,0.0]
  L:#cl,       - Longitude coordinate offset from grid index X=1              [0.0]
                (some of our global otep-grids are North-to-south oriented
                and start at longitude -180.)
 

 EXAMPLES     - see plot/solplot . Try:
      cd ~/Otide
      sol2xyz -fNONE -TM2 -V-1 -CG -A0/5/0/5 -E lepr-otepM2.dat

 for a simple check of data retrieval, try
     
      sol2xyz -dLEPR/ -fNONE -TM2 -V-1 -CG -R -P1,1 -A0/360/-90/90 \
              -o LEPR/lepr-otepM2.xyz lepr-otepM2.dat
      makecpt ... >! amp.cpt
      xyz2grd -Gtmp.grd LEPR/lepr-otepM2.xyz -R0/360/0/90 -I0.5/0.5
      grdview tmp.grd -Ba30f10/a30f10 -Camp.cpt -S3 -Ts -JX10/5 >! plot/lepr-otepM2.ps

.bye