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

#/bin/csh!
gmtset HEADER_FONT 4
###########################################0
set afile=Trend-start-Greenland
set efile=Trend-start-Labrador
set bfile=Trend-start-Okhotsk
set dfile=Trend-start-Bering
set cfile=Trend-start-Barents
set fileps=Trend-start.ps
set crange=0/65/65/80
set drange=160/210/50/70
set arange=315/360/55/85
set erange=260/315/40/80
set brange=135/165/40/65
set resol=1
set ascale=X2.7/2.7
set bscale=X2.7/2.7
set cscale=X2.7/2.7
set dscale=X2.7/2.7
set escale=X2.7/2.7
set fontsize=25
set fonttipe=7
set deg=90
#############################################1
pstext -JX5/5 -R-5/10/-5/10 -K -Y10.1 -X2.5 -P << END >! $fileps
-5 -3 40 0 $fonttipe 5 Trend-Start(1980-95)
END
blockmedian $afile -I$resol -R$arange > $afile.i
surface $afile.i -G$afile.grd -I$resol -R$arange -T0.5
grdimage $afile.grd -J$ascale -R$arange -Ctest.cpt -K -P -Y-2.6 -X-1.8 -O >> $fileps
#grdcontour -C1 -A1f9 -S3/3 $afile.grd -J$ascale -R$arange -K -P -O >> $fileps
pscoast -R$arange -J$ascale -Ba10g5/a10g5:."":neWS -G200 -W3 -Dl -A5000 -K -O >> $fileps
pstext -JX5/5 -R-5/10/-5/10 -X-2.3 -Y0 -O -K << END >> $fileps
0.2 -1.5 $fontsize $deg $fonttipe 5 Latitude
END
pstext -JX5/5 -R-5/10/-5/10 -X1.6 -Y-2.1 -O -K<< END >> $fileps -3 0 $fontsize 0 $fonttipe 5 Longitude(Greenland)
END
#############################################2
blockmedian $bfile -I$resol -R$brange > $bfile.i
surface $bfile.i -G$bfile.grd -I$resol -R$brange -T0.5
grdimage $bfile.grd -J$bscale -R$brange -Ctest.cpt -K -P -Y-1.2 -X0.7 -O >> $fileps
#grdcontour -C1 -A1f9 -S3/3 $bfile.grd -J$bscale -R$brange -K -P -O >> $fileps
pscoast -R$brange -J$bscale -Ba10g5/a10g5:."":neWS -G200 -W3 -Dl -A5000 -K -O>> $fileps
pstext -JX5/5 -R-5/10/-5/10 -X-2.3 -Y0 -O -K << END >> $fileps
0.5 -1 $fontsize $deg $fonttipe 5 Latitude
END
pstext -JX5/5 -R-5/10/-5/10 -X1.6 -Y-2.1 -O -K << END >> $fileps
-3 0 $fontsize 0 $fonttipe 5 Longitude(Okhotsk)
END
#############################################3
blockmedian $cfile -I$resol -R$crange > $cfile.i
surface $cfile.i -G$cfile.grd -I$resol -R$crange -T0.5
grdimage $cfile.grd -J$cscale -R$crange -Ctest.cpt -K -P -Y-1.3 -X0.7 -O >> $fileps
#grdcontour -C1 -A1f9 -S3/3 $cfile.grd -J$cscale -R$crange -K -P -O >> $fileps
pscoast -R$crange -J$cscale -Ba10g5/a10g5:."":neWS -G200 -W3 -Dl -A5000 -K -O >> $fileps
pstext -JX5/5 -R-5/10/-5/10 -X-2.3 -Y0 -O -K << END >> $fileps
0.7 -1 $fontsize $deg $fonttipe 5 Latitude
END
pstext -JX5/5 -R-5/10/-5/10 -X1.6 -Y-2.1 -O -K << END >> $fileps
-3 0 $fontsize 0 $fonttipe 5 Longitude(Barents)
END
#############################################4
blockmedian $dfile -I$resol -R$drange > $dfile.i
surface $dfile.i -G$dfile.grd -I$resol -R$drange -T0.5
grdimage $dfile.grd -J$dscale -R$drange -Ctest.cpt -K -P -Y7.2 -X4.4 -O >> $fileps
#grdcontour -C1 -A1f9 -S3/3 $dfile.grd -J$dscale -R$drange -K -P -O >> $fileps
pscoast -R$drange -J$dscale -Ba10g5/a10g5:."":neWS -G200 -W3 -Dl -A5000 -K -O>> $fileps
pstext -JX5/5 -R-5/10/-5/10 -X-2.3 -Y0 -O -K << END >> $fileps
0.5 -1 $fontsize $deg $fonttipe 5 Latitude
END
pstext -JX5/5 -R-5/10/-5/10 -X1.6 -Y-2.1 -O -K << END >> $fileps
-3 0 $fontsize 0 $fonttipe 5 Longitude(Bering)
END
#############################################5
blockmedian $efile -I$resol -R$erange > $efile.i
surface $efile.i -G$efile.grd -I$resol -R$erange -T0.5
grdimage $efile.grd -J$escale -R$erange -Ctest.cpt -K -P -Y-2 -X0.7 -O >> $fileps
#grdcontour -C1 -A1f9 -S3/3 $efile.grd -J$escale -R$erange -K -P -O >> $fileps
pscoast -R$erange -J$escale -Ba10g5/a10g5:."":neWS -G200 -W3 -Dl -A5000 -K -O >> $fileps
pstext -JX5/5 -R-5/10/-5/10 -X-2.3 -Y0 -O -K << END >> $fileps
0.5 -1.5 $fontsize $deg $fonttipe 5 Latitude
END
pstext -JX5/5 -R-5/10/-5/10 -X1.6 -Y-2.1 -O -K<< END >> $fileps
-3 0 $fontsize 0 $fonttipe 5 Longitude(Labrador)
END
psscale -Ctest.cpt -D5/-2/5/0.5p -O -Y6.5 -X-1.3 >> $fileps
gs $fileps
モドリマスカ？