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 ¶m
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
¶m
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
¶m
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
¶m
plotsize='-JX6.5/% -P'
fnps='${MDIR:[.]}/plot/map.ps'
&end