MSMデータの平面図を作図する

MSM.snapshot.run.sh

#!/bin/sh

exe=MSM.snapshot.gs

if [ ! -f $exe ]; then

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

exit 1

fi

# Parameters

input=0615.nc

prefix=MandA2012Leg3

lev=950 # Vertical 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

it=1

nt=8

while [ $it -le $nt ]; do

out=${prefix}.$(basename $input .nc).${var}.${lev}hPa.$(printf %02d $it).eps

grads -bcp "$exe -input $input $opt -t1 $it -lev $lev -var ${var} \

-ci $ci -out $out -q"

it=$(expr $it + 1 )

done

exit 0

MSM.snapshot.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'

'set rgb 100 153 76 0'

*

* 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 = '-t1' ); t1 = 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

'set t 't1

#'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=lev

say pha

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 'lonw' 'lone

'set lat 'lats' 'latn

'set mpdset hires'

'set map 1 1 7'

'set xlopts 1 5 0.15'

'set xlint 2'

'set ylopts 1 5 0.15'

'set ylint 2'

### ept

'set gxout shade2'

'color 326 350 2 -kind lightcyan->cyan->mediumspringgreen->lawngreen->moccasin->orange->orangered->red'

'd 'var

'cbarn'

'set gxout contour'

'set ccolor 0'

'set cthick 5'

'set clab off'

'set cint 2'

'd 'var

### vector

'set gxout vector'

'set ccolor 0'

'set cthick 20'

'd skip(u,8,8);skip(v,8,8)'

'set ccolor 1'

'set cthick 5'

'd skip(u,8,8);skip(v,8,8)'

'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.15'

'set clevs 340 348'

'd 'var

### Z

#'set gxout contour'

#'set cint '20

#'zsmth=smth9(z)'

#'set ccolor 0'

#'set cthick 20'

#'set clab off'

#'d 'zsmth

#'set cint '20

#'set clab on'

#'set cthick 5'

#'set ccolor 1'

#'set clopts 1 2 0.1'

#'d 'zsmth

'draw title 'datetime' 'pha' hPa'

'print 'out

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

return