MSM.snapshot.panel.run.sh
#!/bin/sh
exe=MSM.snapshot.panel.gs
if [ ! -f $exe ]; then
echo Error in $0 : No such file, $exe
exit 1
fi
# Parameters
input=0615.nc
prefix=MandA2012Leg3
lev=950 # Pressure Level
lonw=120
lone=132
lats=22.5
latn=35.0
opt=# "-lonw $lonw -lone $lonw -lats $lats -latn $latn"
var=ept #theta
ci=3 #Contour Interval
out=${prefix}.$(basename $input .nc).${var}.${lev}hPa.Panel.eps
echo "$exe -input $input $opt -lev $lev -var ${var} \
-ci $ci -out $out -q"
grads -bcp "$exe -input $input $opt -lev $lev -var ${var} \
-ci $ci -out $out -q"
exit 0
MSM.snapshot.panel.gs
function msm( args )
*
* Default values
*
indir='/work2/kunoki/to_manda_sensei/Data/MSM/Original_data/Pressure'
input='0522.nc'
lonw=122.0 ; lone=134 ;*lone=150
lats=24.5 ; latn=35 ;* latn=47.5
out='msm.snapshot.test.eps'
quitopt='no'
lev=950
ci=20
var='temp'
sclv=40
*
* 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 = '-lev' ); lev = 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
#'q ctlinfo'
#say result
'set lev 'lev/100
pha=lev
say pha
'set lon 'lonw' 'lone
'set lat 'lats' 'latn
######
t1=1
'mul 3 3 1 3 -yint 1'
'set t 't1
'q dims'
say result
line=sublin(result,5)
datetime=subwrd(line,6)
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 mpdset mres'
'set map 1 1 5'
'set xlopts 1 3 0.12'
'set xlint 4'
'set ylopts 1 3 0.12'
'set ylint 2'
### ept
'set gxout shade2'
'color 328 348 4 -kind lightcyan->cyan->mediumspringgreen->lawngreen->moccasin->orange->orangered->red'
'd 'var
#'cbarn'
'set gxout contour'
'set ccolor 0'
'set cthick 4'
'set clab off'
'set cint 2'
'd 'var
### vector
'set gxout vector'
'set arrscl 0.5 'sclv
'set arrlab off'
'set ccolor 0'
'set cthick 10'
'd skip(u,10,10);skip(v,10,10)'
'set ccolor 1'
'set cthick 3'
'set arrscl 0.5 'sclv
'd skip(u,10,10);skip(v,10,10)'
'set strsiz 0.10 0.15'
'set string 1 c 5 0'
# 'draw string 6.6 0.95 m/s'
### ept (coutour label)
'set gxout contour'
'set ccolor 0'
'set cthick 0'
'set clab on'
'set clopts 1 4 0.1'
'set clevs 340'
'd 'var
'set strsiz 0.15 0.15'
'set string 1 c 6'
'draw string 1.8 10.8 'datetime
######
t1=2
'mul 3 3 2 3'
'set t 't1
'q dims'
say result
line=sublin(result,5)
datetime=subwrd(line,6)
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 mpdset mres'
'set map 1 1 5'
'set xlopts 1 3 0.12'
'set xlint 4'
'set ylopts 1 3 0.12'
'set ylint 2'
### ept
'set gxout shade2'
'color 328 348 4 -kind lightcyan->cyan->mediumspringgreen->lawngreen->moccasin->orange->orangered->red'
'd 'var
#'cbarn'
'set gxout contour'
'set ccolor 0'
'set cthick 4'
'set clab off'
'set cint 2'
'd 'var
### vector
'set gxout vector'
'set arrscl 0.5 'sclv
'set arrlab off'
'set ccolor 0'
'set cthick 10'
'd skip(u,10,10);skip(v,10,10)'
'set ccolor 1'
'set cthick 3'
'set arrscl 0.5 'sclv
'd skip(u,10,10);skip(v,10,10)'
'set strsiz 0.10 0.15'
'set string 1 c 5 0'
# 'draw string 6.6 0.95 m/s'
### ept (coutour label)
'set gxout contour'
'set ccolor 0'
'set cthick 0'
'set clab on'
'set clopts 1 4 0.1'
'set clevs 340'
'd 'var
'set strsiz 0.15 0.15'
'set string 1 c 6'
'draw string 4.3 10.8 'datetime
######
t1=3
'mul 3 3 3 3'
'set t 't1
'q dims'
say result
line=sublin(result,5)
datetime=subwrd(line,6)
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 mpdset mres'
'set map 1 1 5'
'set xlopts 1 3 0.12'
'set xlint 4'
'set ylopts 1 3 0.12'
'set ylint 2'
### ept
'set gxout shade2'
'color 328 348 4 -kind lightcyan->cyan->mediumspringgreen->lawngreen->moccasin->orange->orangered->red'
'd 'var
#'cbarn'
'set gxout contour'
'set ccolor 0'
'set cthick 4'
'set clab off'
'set cint 2'
'd 'var
### vector
'set gxout vector'
'set arrscl 0.5 'sclv
'set arrlab off'
'set ccolor 0'
'set cthick 10'
'd skip(u,10,10);skip(v,10,10)'
'set ccolor 1'
'set cthick 3'
'set arrscl 0.5 'sclv
'd skip(u,10,10);skip(v,10,10)'
'set strsiz 0.10 0.15'
'set string 1 c 5 0'
# 'draw string 6.6 0.95 m/s'
### ept (coutour label)
'set gxout contour'
'set ccolor 0'
'set cthick 0'
'set clab on'
'set clopts 1 4 0.1'
'set clevs 340'
'd 'var
'set strsiz 0.15 0.15'
'set string 1 c 6'
'draw string 6.8 10.8 'datetime
######
t1=4
'mul 3 3 1 2 -yint 1'
'set t 't1
'q dims'
say result
line=sublin(result,5)
datetime=subwrd(line,6)
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 mpdset mres'
'set map 1 1 5'
'set xlopts 1 3 0.12'
'set xlint 4'
'set ylopts 1 3 0.12'
'set ylint 2'
### ept
'set gxout shade2'
'color 328 348 4 -kind lightcyan->cyan->mediumspringgreen->lawngreen->moccasin->orange->orangered->red'
'd 'var
#'cbarn'
'set gxout contour'
'set ccolor 0'
'set cthick 4'
'set clab off'
'set cint 2'
'd 'var
### vector
'set gxout vector'
'set arrscl 0.5 'sclv
'set arrlab off'
'set ccolor 0'
'set cthick 10'
'd skip(u,10,10);skip(v,10,10)'
'set ccolor 1'
'set cthick 3'
'set arrscl 0.5 'sclv
'd skip(u,10,10);skip(v,10,10)'
'set strsiz 0.10 0.15'
'set string 1 c 5 0'
# 'draw string 6.6 0.95 m/s'
### ept (coutour label)
'set gxout contour'
'set ccolor 0'
'set cthick 0'
'set clab on'
'set clopts 1 4 0.1'
'set clevs 340'
'd 'var
'set strsiz 0.15 0.15'
'set string 1 c 6'
'draw string 1.8 7.3 'datetime
######
t1=5
'mul 3 3 2 2'
'set t 't1
'q dims'
say result
line=sublin(result,5)
datetime=subwrd(line,6)
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 mpdset mres'
'set map 1 1 5'
'set xlopts 1 3 0.12'
'set xlint 4'
'set ylopts 1 3 0.12'
'set ylint 2'
### ept
'set gxout shade2'
'color 328 348 4 -kind lightcyan->cyan->mediumspringgreen->lawngreen->moccasin->orange->orangered->red'
'd 'var
#'cbarn'
'set gxout contour'
'set ccolor 0'
'set cthick 4'
'set clab off'
'set cint 2'
'd 'var
### vector
'set gxout vector'
'set arrscl 0.5 'sclv
'set arrlab off'
'set ccolor 0'
'set cthick 10'
'd skip(u,10,10);skip(v,10,10)'
'set ccolor 1'
'set cthick 3'
'set arrscl 0.5 'sclv
'd skip(u,10,10);skip(v,10,10)'
'set strsiz 0.10 0.15'
'set string 1 c 5 0'
# 'draw string 6.6 0.95 m/s'
### ept (coutour label)
'set gxout contour'
'set ccolor 0'
'set cthick 0'
'set clab on'
'set clopts 1 4 0.1'
'set clevs 340'
'd 'var
'set strsiz 0.15 0.15'
'set string 1 c 6'
'draw string 4.3 7.3 'datetime
######
t1=6
'mul 3 3 3 2'
'set t 't1
'q dims'
say result
line=sublin(result,5)
datetime=subwrd(line,6)
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 mpdset mres'
'set map 1 1 5'
'set xlopts 1 3 0.12'
'set xlint 4'
'set ylopts 1 3 0.12'
'set ylint 2'
### ept
'set gxout shade2'
'color 328 348 4 -kind lightcyan->cyan->mediumspringgreen->lawngreen->moccasin->orange->orangered->red'
'd 'var
#'cbarn'
'set gxout contour'
'set ccolor 0'
'set cthick 4'
'set clab off'
'set cint 2'
'd 'var
### vector
'set gxout vector'
'set arrscl 0.5 'sclv
'set arrlab off'
'set ccolor 0'
'set cthick 10'
'd skip(u,10,10);skip(v,10,10)'
'set ccolor 1'
'set cthick 3'
'set arrscl 0.5 'sclv
'd skip(u,10,10);skip(v,10,10)'
'set strsiz 0.10 0.15'
'set string 1 c 5 0'
# 'draw string 6.6 0.95 m/s'
### ept (coutour label)
'set gxout contour'
'set ccolor 0'
'set cthick 0'
'set clab on'
'set clopts 1 4 0.1'
'set clevs 340'
'd 'var
'set strsiz 0.15 0.15'
'set string 1 c 6'
'draw string 6.8 7.3 'datetime
######
t1=7
'mul 3 3 1 1'
'set t 't1
'q dims'
say result
line=sublin(result,5)
datetime=subwrd(line,6)
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 mpdset mres'
'set map 1 1 5'
'set xlopts 1 3 0.12'
'set xlint 4'
'set ylopts 1 3 0.12'
'set ylint 2'
### ept
'set gxout shade2'
'color 328 348 4 -kind lightcyan->cyan->mediumspringgreen->lawngreen->moccasin->orange->orangered->red'
'd 'var
#'cbarn'
'set gxout contour'
'set ccolor 0'
'set cthick 4'
'set clab off'
'set cint 2'
'd 'var
### vector
'set gxout vector'
'set arrscl 0.5 'sclv
'set arrlab off'
'set ccolor 0'
'set cthick 10'
'd skip(u,10,10);skip(v,10,10)'
'set ccolor 1'
'set cthick 3'
'set arrscl 0.5 'sclv
'd skip(u,10,10);skip(v,10,10)'
'set strsiz 0.10 0.15'
'set string 1 c 5 0'
# 'draw string 6.6 0.95 m/s'
### ept (coutour label)
'set gxout contour'
'set ccolor 0'
'set cthick 0'
'set clab on'
'set clopts 1 4 0.1'
'set clevs 340'
'd 'var
'set strsiz 0.15 0.15'
'set string 1 c 6'
'draw string 1.8 3.8 'datetime
######
t1=8
'mul 3 3 2 1'
'set t 't1
'q dims'
say result
line=sublin(result,5)
datetime=subwrd(line,6)
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 mpdset mres'
'set map 1 1 5'
'set xlopts 1 3 0.12'
'set xlint 4'
'set ylopts 1 3 0.12'
'set ylint 2'
### ept
'set gxout shade2'
'color 328 348 4 -kind lightcyan->cyan->mediumspringgreen->lawngreen->moccasin->orange->orangered->red'
'd 'var
'cbarn 0.75 1 6.5 2.4'
'set strsiz 0.15 0.15'
'set string 1 c 6'
'draw string 6.78 3.8 [K]'
'set gxout contour'
'set ccolor 0'
'set cthick 4'
'set clab off'
'set cint 2'
'd 'var
### vector
'set gxout vector'
'set arrscl 0.5 'sclv
'set arrlab on'
'set ccolor 0'
'set cthick 10'
'd skip(u,10,10);skip(v,10,10)'
'set ccolor 1'
'set cthick 3'
'set arrscl 0.5 'sclv
'd skip(u,10,10);skip(v,10,10)'
'set strsiz 0.10 0.15'
'set string 1 c 5 0'
'draw string 3.85 0.85 m/s'
### ept (coutour label)
'set gxout contour'
'set ccolor 0'
'set cthick 0'
'set clab on'
'set clopts 1 4 0.1'
'set clevs 340'
'd 'var
'set strsiz 0.15 0.15'
'set string 1 c 6'
'draw string 4.3 3.8 'datetime
'print 'out
if ( quitopt = "yes" ); quit; endif
return