USAGE:
          drawcoast @ ins-file [:target]


 PURPOSE:
          For producing and to visually inspect model boundaries.
          This program includes a GMT plot feature.
          You may consider program prpbound as a plot-less alternative.

 INS-FILE:

1. Open file block
   Open-string          I/O  File content                                  
   
   21 < flm-file       - I - M-flags
   22 < flz-file       - I - Z-flags and area/projection parameters
   31 B coast-xy       - O - cartesian coordinates, index-space
   32 B gmt-xy         - O - (optional) GMT coastline, index-space
   33 B coast-xy-prj   - O -
(optional) in cartesian-km or geographic space
      Q
(GMT in cartesian-km still missing as an option)
 2. Namelist &param

   Parameter      type       meaning                                      [default]

   include_flags  char*8     The flags to pass the array scan.             ['LaAS']
   qmarkland      logical    If .false., leave land areas blank.           [.true.]

 
  qkm            logical    Transform from index to kilometers, file 33   [.true.]
   qsphere        logical   
Transform from index to sphere, file 33      [.false.]

   qgmt           logical    Do the plot stage                            
[.true.]        

 
  If  qgmt:

 
 
fnps           char*128   postscript output file, default is derived
                             from the coast-xy file                          
                             Environment parameters will be resolved with
                             replace_envars_l                                 [' ']


   keep           char*3     '-K'   keeps the ps-file open.
                             File .gmtcommands4 is found in the directory
                             of script launch.      

   coastpath      char*128   Path to the coastline database. Defaults to
                             where OTEQ/PROG/s/glcorfgs.f defaults to         [' ']

   gmtopt         char*128   Options for global_coasts_gmt (in glcorfgs.f)
                                                     ['-A1000/1/4 -Df -I1 -I2 -I3']    
 
   plotsize       char*32    The `-JX´ -option in GMT psxy, eventually
                             with the portrait option, like '-JX6/8 -P'
                             If an aspect ratio of 1 is desired, the y-value
                             can be replaced with a `%´-sign; the program
                             will compute the value.
The unit (usually inch)  
                             depends on .gmtdefaults4
in the directory
                             of launch.                                  ['-JX9/6']

 EXAMPLES:

  (1) Only active and passive boundaries, only the input to a prospective psxy command

     setenv MDIR BALTIC
     drawcoast @ BALTIC/drawcoast.ins :XY

  
with instructions
XY>
21 < ${MDIR:[.]}/FLMH.DAT

22 < ${MDIR:[.]}/FLZ.DAT
31 B ${MDIR:[.]}/coast.dat
   Q
 &param
 qgmt=.false.
 include_flags='LAa', qmarkland=.false.
 &end

 (2) The same as (1), however creating the plot

    setenv MDIR BALTIC
    drawcoast @ BALTIC/drawcoast.ins :PLOT
    cd BALTIC/plot
    ps2png -m -d 144x144 -o ~/www/OTEQ/BALTIC/ apbounds.ps
    #
    # see http://holt.oso.chalmers.se/hgs/OTEQ/BALTIC/apbounds.png


   with instructions
PLOT>
21 < ${MDIR:[.]}/FLMH.DAT
22 < ${MDIR:[.]}/FLZ.DAT
31 B ${MDIR:[.]}/coast.dat
33 B ${MDIR:[.]}/coast-km.dat
   Q
 &param
 plotsize='-JX6.5/% -P'
 include_flags='LAa', qmarkland=.false.
 fnps='${MDIR:[.]}/plot/apbounds.ps'
 &end

 (3) Doing all

    setenv MDIR BALTIC
   
drawcoast @ BALTIC/drawcoast.ins
    cd BALTIC/plot
    ps2png -m -d 144x144 -o ~/www/OTEQ/BALTIC/ map.ps 

    #
    # see
http://holt.oso.chalmers.se/hgs/OTEQ/BALTIC/map.png

   with instructions
21 < ${MDIR:[.]}/FLMH.DAT
22 < ${MDIR:[.]}/FLZ.DAT
31 B ${MDIR:[.]}/coast.dat
32 B ${MDIR:[.]}/gmtcoast.dat
   Q
 &param
 plotsize='-JX6.5/% -P'
 fnps='${MDIR:[.]}/plot/map.ps'
 &end