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']

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

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

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