/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