MSMデータの平面図を作図する（1日分をパネル表示）

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