OISST日平均値の作図

実行例

データのダウンロード

$ sh get.oisst.sh

$ ls -lh sst.day.mean.2013.v2.nc

-rw-rw-r-- 1 am am 453M 11月 17 19:41 sst.day.mean.2013.v2.nc

平面図の作図

$ grads -bcp "oisst.gs -d 20130201 -q"

スクリプト

ダウンロード用シェルスクリプト: get.oisst.sh

#!/bin/bash

# Link to ftp site

# http://www.esrl.noaa.gov/psd/cgi-bin/db_search/DBListFiles.pl?did=132&tid=44137&vid=2423

#

wget -nc ftp://ftp.cdc.noaa.gov/Datasets/noaa.oisst.v2.highres/sst.day.mean.2013.v2.nc

作図用GrADSスクリプト: oisst.gs

function oisst( args )

# Default values

in='sst.day.mean.2013.v2.nc'

yyyymmdd=20130101

lonw=114

lone=144

lats=-14

latn=4

# Decode options

i = 1

while( 1 )

arg = subwrd( args, i )

i = i + 1;

if( arg = '' ); break; endif

while( 1 )

if( arg='-i'); in=subwrd(args,i);i=i+1;break;endif

if( arg='-d' ); yyyymmdd=subwrd(args,i) ;i=i+1; break; endif

if( arg='-q'); quitflag=true; break; endif

say 'Syntax error : arg= 'arg

return

endwhile

endwhile

'sdfopen 'in

yyyy=substr(yyyymmdd,1,4)

mm=substr(yyyymmdd,5,2)

dd=substr(yyyymmdd,7,2)

if(mm=01);mmm='JAN';endif

if(mm=02);mmm='FEB';endif

if(mm=03);mmm='MAR';endif

if(mm=04);mmm='APR';endif

if(mm=05);mmm='MAY';endif

if(mm=06);mmm='JUN';endif

if(mm=07);mmm='JUL';endif

if(mm=08);mmm='AUG';endif

if(mm=09);mmm='SEP';endif

if(mm=10);mmm='OCt';endif

if(mm=11);mmm='NOV';endif

if(mm=12);mmm='DEC';endif

datetime=dd''mmm''yyyy

'set time 'datetime

'set lon 'lonw' 'lone

'set lat 'lats' 'latn

'q dims'

say result

'set parea 1 7.5 1 10'

'set mpdset mres'

'set ylopts 1 5 0.2'

'set xlopts 1 5 0.2'

'set gxout shade2'

'color -kind mediumturquoise->greenyellow->khaki->orange->red 27 31 0.5'

#'rgbset2'

#'set clevs 27 27.5 28 28.5 29 29.5 30 30.5 31'

'd sst'

'draw title OISST 'datetime

'cbarn 0.9 1 8 5.5'

'set strsiz 0.15 0.2'

'draw string 8.1 7.8 `3.`0C'

'set gxout contour'

'set xlopts 1 0 0'

'set xlopts 1 0 0'

'set cthick 8'

'set ccolor 0'

'set clab off'

'd sst'

'set xlopts 1 0 0'

'set ylopts 1 0 0'

'set cthick 4'

'set ccolor 1'

'set clab on'

'set clopts 1 4 0.12'

'd sst'

eps='OISST.'yyyymmdd'.eps'

'print 'eps

png='OISST.'yyyymmdd'.png'

'!convert 'eps' 'png

if (quitflag = 'true')

quit

endif

return