シェルスクリプトはこちら

#/bin/csh!
gmtset HEADER_FONT 4
###########################################
set file=sample8
set fileps=sample8.ps
set resol=1
set range=138/145/35/45
set scale=m0.5
set fontsize=25
set fonttipe=7
set l=1/30
set g=2/2
set futosa=0.05
set tatenaga=0.2
set yokonaga=0.2
set type=-Sv$futosa/$tatenaga/$yokonaga
#############################################
pstext -JX5/5 -R-10/10/-10/10 -K -Y9 -X2 -P << END >! $fileps
-5 -7 40 0 $fonttipe 5 Sample8
END
####################
pscoast -J$scale -R$range -Ba1g1/a1g1:."":neWS -G155/220/0 -S220/220/255 -W1 -Dl -O -K -Y-7 -X0>> $fileps
####################
psxy rensyu1.dat -J$scale -R$range $type -W8 -O -K -G255/0/0 >> $fileps
pstext -JX -R-5/10/-5/10 -X-2 -O -K -P << END >> $fileps
-3 0 $fontsize 90 $fonttipe 5 Latitude(Deg)
END
pstext -JX -R-5/10/-10/10 -X2 -O -P -Y-2 << END >> $fileps
-4 -5 $fontsize 0 $fonttipe 5 Longitude(Deg)
END
gs $fileps
#


#! /bin/tcsh -f
gmtset HEADER_FONT 4
###########################################
set file=vector
set fileps=$file.ps
set resol=0.05
set range=135/145/42/50
set rangemask=135/139/42/45
set scale=M15c
set futosa=0.04
set tatenaga=0.15
set yokonaga=0.15
set type=-Sv$futosa/$tatenaga/$yokonaga
set gfile=test1.grd
#############################################
pscoast -J$scale -R$range -Ba1g1/a1g1:."Quickscat":neWS -G155/220/0 -S220/220/255 -W0.1 -Dh -K -P >! $fileps
psxy test.dat -J$scale -R$range $type -W0.1 -G255/100/100 -P -O -K >> $fileps
awk '{print $1" "$2" "$3}' wind.dat |\
blockmedian -I$resol -R$range |\
surface -G$gfile -I$resol -R$range -T1.0
awk '{print $1" "$2" "$4}' test.dat |\
blockmedian -I$resol -R$range |\
psmask -R$range -I0.9 -J$scale -B2 -P -O -K >> $fileps
grdcontour $gfile -C1 -L0/50 -J$scale -R$range -A1f15 -O -K -P >> $fileps
psmask -C -O >> $fileps
gs $fileps
もどります