psmask 地図上にマスクをかけてコンターを引く

データの無いところにはマスクをかけてコンターを引く。

シェルスクリプトはこちら
#! /bin/tcsh -f
set file=output
set file1=sd.dat
set psfile=$file.ps
set gfile=$file.grd
set range=125/160/30/55
set scale=X16/16
set resol=1
###
blockmedian $file -R$range -I$resol -V | surface -G$gfile -R$range -I$resol -T0.5 -V
pscoast -J$scale -R$range -Df -G0/255/0 -S200/200/255 -W0.01p/255/255/255 -Y9 -X4 -P -K -Bf5/f5:."1990-2000 nov":nWeS >! $psfile
psbasemap -R$range -J$scale -Lx5/1/0/200k -O -K >> $psfile
psmask $file -J$scale -R$range -I$resol -K -O >> $psfile
grdcontour $gfile -J$scale -R$range -C1 -A1f10 -W1 -V -P -O -K >> $psfile
psmask -C -O -K >> $psfile
blockmedian $file1 -R$range -I$resol -V | surface -G$gfile -R$range -I$resol -T0.5 -V
psmask $file1 -J$scale -R$range -I$resol -K -O >> $psfile
grdcontour $gfile -J$scale -R$range -C0.5 -A0.5f8 -W1/255/0/0 -V -P -O -K >> $psfile
psmask -C -O >> $psfile
gs $psfile

データはoutput ,sd.datともに経度、緯度、データの順。
出来上がりの図
モドル

Last modified: Sat May 11 18:49:55 JST 2002