a13x06 |
03-13-2010 09:27 AM |
Generic Mapping Tools (GMT) text help
Bellow is a script I'm running to produce several diagrams but I just can't get the pstext command to produce any text:(
Code:
#!/bin/csh
#Script "plotdata" for the assessment
#
#Setting GMT environmental variables
gmtset MEASURE_UNIT cm
gmtset PAPER_MEDIA a4+
gmtset PAGE_ORIENTATION portrait
# setting script environmental variables
set range = 341/345/27/31
set scale = 4i
set int = 1m
set search = 5k
#producing grid files
nearneighbor bathy_a_gmt.txt -R$range -I$int -S$search -Gbathy.grd -V
nearneighbor gravity_a_gmt.txt -R$range -I$int -S$search -Ggravity.grd -V
#making colour palettes
makecpt -Crainbow -T-80/316/20 > gravity.cpt
makecpt -Crainbow -T-2529/4745/500 > bathy.cpt
#2D model
grdimage gravity.grd -R$range -JM$scale -Cgravity.cpt -Y3 -X5 -K -V > plot.ps
pscoast -R$range -JM$scale -G200 -W -O -K -V >> plot.ps
grdcontour bathy.grd -R$range -JM$scale -L-2530/4746 -C500 -A1500 -O -K -V >> plot.ps
psbasemap -R$range -JM$scale -Ba2f2g2 -K -O -V >> plot.ps
#3D Model
grdview bathy.grd -Y13.5 -X-4 -JM5i -JZ2i -Qs -E135/10 -B2 -W -Cbathy.cpt -O -K >> plot.ps
#Adding scales
psscale -D3.5i/0i/7.2i/0.125ih -Cbathy.cpt -B/:"Bathymetry(m)": -O -K -V >> plot.ps
psscale -Y-15 -D3.5i/0i/7i/0.125ih -Cgravity.cpt -B/:"Gravity": -O -K -V >> plot.ps
#Adding Text
pstext -R$range -JM$scale -G1/0/0 -X5 -Y10 -V -O << END >>! plot.ps
1 1 20 1 4 BL GMT
END
#Producing Final Plot
gv plot.ps
Any help appreciated!!
Alex
|