MSMデータから水平温度移流のホフメラー図を作成する

/work2/kunoki/to_manda_sensei/Tools/MSM.Temp.Adv.Hovmoller

実行用シェルスクリプト: MSM.hovmoller.run.sh

#!/bin/sh

exe=MSM.hovmoller.gs

if [ ! -f $exe ]; then

echo Error in $0 : No such file, $exe

exit 1

fi

# Parameters

input=0615.nc

prefix=MandA2012Leg3

z1=3 # Vertical Level

lonw=129.5

lone=130.0

opt="-input $input -lonw $lonw -lone $lone"

var=theta #ept

ci=1 #3 #Contour Interval

out=${prefix}.$(basename $input .nc).T.ADV.${var}.hovmoller.eps

grads -bcp "$exe $opt -z1 $z1 -var ${var} -ci $ci -out $out -q"

exit 0

GrADSスクリプト: MSM.hovmoller.gs

function msm( args )

*

* Default values

*

indir='/work2/kunoki/to_manda_sensei/Data/MSM/Original_data/Pressure'

input='0515.nc'

lonw=129.0; lone=129.5

lats=28 ; latn=35

out='msm.hovmoller.test.eps'

quitopt='no'

z1=1

t1=1; t2=8

ci=5

var='temp'

*

* Decode options

*

i = 1

while( 1 )

arg = subwrd( args, i )

i = i + 1;

if( arg = '' ); break; endif

while( 1 )

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

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

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

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

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

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

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

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

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

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

if( arg = '-q' ); quitopt=yes ;i=i+1; break; endif

say 'Syntax error : arg= 'arg

return

endwhile

endwhile

#'sdfopen 'indir'/'input

'open MSM.ctl'

#'q ctlinfo'

#say result

'set z 'z1

'set time 00Z15JUN2012 00Z16JUN2012'

#'q file'

#say result

'q dims'

say result

line=sublin(result,5)

datetime=subwrd(line,6)

line=sublin(result,4)

pdb=subwrd(line,6)

pha=pdb*100

say pha

# 'set lon 'lonw' 'lone

'define dtx = cdiff(temp,x)'

'define dty = cdiff(temp,y)'

'define dx = cdiff(lon,x)*3.1416/180'

'define dy = cdiff(lat,y)*3.1416/180'

hadv'=-1*( ('u'*'dtx')/(cos('lat'*3.1416/180)*'dx') + 'v'*'dty'/'dy' )/6.37e6*84600'

if (var = "ept" )

'es=6.1173*exp(((2.501*pow(10,6))/461.50)*(1/273.16 - 1/temp))'

'ws=621.97*(es/(' pha '-es))'

'w=(rh*ws)/(100*1000)'

'ept=(temp+((2.501*pow(10,6))/1004)*w)*pow((1000/' pha '),287/1004)'

endif

if (var = "theta" )

'theta=temp*pow((1000/' pha '),287/1004)'

endif

'set lon 129'

varave='ave(' var ',lon=' lonw ',lon=' lone ')'

#varave='ave(' var ',lon=129,lon=129.5)'

'set lat 'lats' 'latn

'set cint 'ci

'set parea 2 8 2 10.5'

#'set clopts 1 3 0.1'

'set cthick 5'

'set tlsupp year'

'set xlint 1'

'set xlopts 1 5 0.18'

'set ylopts 1 5 0.18'

'set gxout shaded'

'rgbset2'

'set clevs -100 -50 -20 -10 10 20 50 100'

'set clopts 1 4 0.2'

'd hadv'

'cbarn 1 0 4.5 1'

'set strsiz 0.1 0.15'

'draw string 8 1 [K/d]'

'set gxout contour'

'set ccolor 3'

'set clopts 3 3 0.15'

'd 'varave

'q w2xy 32.5414 02z15jun2012'

say result

x1= subwrd(result,3)

y1= subwrd(result,6)

'set line 1'

'draw mark 9 ' x1 ' 'y1' '0.25

x2=x1-0.25

'set string 1 c 5 0'

'set strsiz 0.2'

'draw string ' x2 ' ' y1 ' 1'

'q w2xy 32.3677 03z15jun2012'

say result

x1= subwrd(result,3)

y1= subwrd(result,6)

'set line 1'

'draw mark 9 ' x1 ' 'y1' '0.25

x2=x1-0.25

'set strsiz 0.2'

'draw string ' x2 ' ' y1 ' 2'

'q w2xy 32.1521 04z15jun2012'

say result

x1= subwrd(result,3)

y1= subwrd(result,6)

'set line 1'

'draw mark 9 ' x1 ' 'y1' '0.25

x2=x1-0.25

'set strsiz 0.2'

'draw string ' x2 ' ' y1 ' 3'

'q w2xy 32.0480 05z15jun2012'

say result

x1= subwrd(result,3)

y1= subwrd(result,6)

'set line 1'

'draw mark 9 ' x1 ' 'y1' '0.25

x2=x1-0.25

'set strsiz 0.2'

'draw string ' x2 ' ' y1 ' 4'

'q w2xy 31.9070 06z15jun2012'

say result

x1= subwrd(result,3)

y1= subwrd(result,6)

'set line 1'

'draw mark 9 ' x1 ' 'y1' '0.25

x2=x1-0.25

'set strsiz 0.2'

'draw string ' x2 ' ' y1 ' 5'

'q w2xy 31.7205 07z15jun2012'

say result

x1= subwrd(result,3)

y1= subwrd(result,6)

'set line 1'

'draw mark 9 ' x1 ' 'y1' '0.25

x2=x1-0.25

'set strsiz 0.2'

'draw string ' x2 ' ' y1 ' 6'

'q w2xy 31.4864 08z15jun2012'

say result

x1= subwrd(result,3)

y1= subwrd(result,6)

'set line 1'

'draw mark 9 ' x1 ' 'y1' '0.25

x2=x1-0.25

'set strsiz 0.2'

'draw string ' x2 ' ' y1 ' 7'

'q w2xy 31.3434 09z15jun2012'

say result

x1= subwrd(result,3)

y1= subwrd(result,6)

'set line 1'

'draw mark 9 ' x1 ' 'y1' '0.25

x2=x1-0.25

'set strsiz 0.2'

'draw string ' x2 ' ' y1 ' 8'

'q w2xy 31.1561 10z15jun2012'

say result

x1= subwrd(result,3)

y1= subwrd(result,6)

'set line 1'

'draw mark 9 ' x1 ' 'y1' '0.25

x2=x1-0.25

'set strsiz 0.2'

'draw string ' x2 ' ' y1 ' 9'

'q w2xy 30.9609 11z15jun2012'

say result

x1= subwrd(result,3)

y1= subwrd(result,6)

'set line 1'

'draw mark 9 ' x1 ' 'y1' '0.25

x2=x1-0.3

'set strsiz 0.2'

'draw string ' x2 ' ' y1 ' 10'

'q w2xy 30.7565 12z15jun2012'

say result

x1= subwrd(result,3)

y1= subwrd(result,6)

'set line 1'

'draw mark 9 ' x1 ' 'y1' '0.25

x2=x1-0.3

'set strsiz 0.2'

'draw string ' x2 ' ' y1 ' 11'

'q w2xy 30.5166 13z15jun2012'

say result

x1= subwrd(result,3)

y1= subwrd(result,6)

'set line 1'

'draw mark 9 ' x1 ' 'y1' '0.25

x2=x1-0.3

'set strsiz 0.2'

'draw string ' x2 ' ' y1 ' 12'

'q w2xy 30.3081 14z15jun2012'

say result

x1= subwrd(result,3)

y1= subwrd(result,6)

'set line 1'

'draw mark 9 ' x1 ' 'y1' '0.25

x2=x1-0.3

'set strsiz 0.2'

'draw string ' x2 ' ' y1 ' 13'

'q w2xy 30.2938 15z15jun2012'

say result

x1= subwrd(result,3)

y1= subwrd(result,6)

'set line 1'

'draw mark 9 ' x1 ' 'y1' '0.25

x2=x1-0.3

'set strsiz 0.2'

'draw string ' x2 ' ' y1 ' 14'

'q w2xy 30.1333 16z15jun2012'

say result

x1= subwrd(result,3)

y1= subwrd(result,6)

'set line 1'

'draw mark 9 ' x1 ' 'y1' '0.25

x2=x1-0.3

'set strsiz 0.2'

'draw string ' x2 ' ' y1 ' 15'

'q w2xy 29.9550 17z15jun2012'

say result

x1= subwrd(result,3)

y1= subwrd(result,6)

'set line 1'

'draw mark 9 ' x1 ' 'y1' '0.25

x2=x1-0.3

'set strsiz 0.2'

'draw string ' x2 ' ' y1 ' 16'

'q w2xy 29.7799 18z15jun2012'

say result

x1= subwrd(result,3)

y1= subwrd(result,6)

'set line 1'

'draw mark 9 ' x1 ' 'y1' '0.25

x2=x1-0.3

'set strsiz 0.2'

'draw string ' x2 ' ' y1 ' 17'

'q w2xy 29.3908 20z15jun2012'

say result

x1= subwrd(result,3)

y1= subwrd(result,6)

'set line 1'

'draw mark 9 ' x1 ' 'y1' '0.25

x2=x1-0.3

'set strsiz 0.2'

'draw string ' x2 ' ' y1 ' 18'

'q w2xy 28.9772 22z15jun2012'

say result

x1= subwrd(result,3)

y1= subwrd(result,6)

'set line 1'

'draw mark 9 ' x1 ' 'y1' '0.25

x2=x1-0.3

'set strsiz 0.2'

'draw string ' x2 ' ' y1 ' 19'

'draw title 'datetime' 'pha' hPa'

'print 'out

if ( quitopt = "yes" ); quit; endif

return