aofd165.bio.mie-u.ac.jp
/work1/am/2016.Atmospheric.Data.Ornithology
Tue Jun 14 14:23:44 JST 2016
======================
JRA-55_T.snapshot.run.sh
======================
#!/bin/sh
exe=JRA-55_T.snapshot.gs
if [ ! -f $exe ]; then
echo Error in $0 : No such file, $exe
exit 1
fi
# Parameters
dir=/work4/data/JRA-55
infle="$dir/anl_surf.2015010100.ctl"
prefix=JRA-55_T
outdir="Fig"
#lev=1000 #950 # Vertical Level
lonw=120
lone=150
lats=22.4
latn=47.6
ci=5 #Contour Interval
opt="-lonw $lonw -lone $lone -lats $lats -latn $latn -ci $ci"
var=temp #ept #theta
time='00Z01JAN2015'
out=${prefix}_${time}.SFC.eps
grads -bcp "$exe -input $infle \
$opt -time $time -var ${var} \
-ci $ci -outdir $outdir -out $out -q"
echo
echo Output: $out
echo
exit 0
----------------------
End of JRA-55_T.snapshot.run.sh
----------------------
======================
JRA-55_T.snapshot.gs
======================
function jra55( args )
*
* Default values
*
infle='anl_surf.2015010100.ctl'
lonw=122.0 ; lone=134 ;*lone=150
lats=24.5 ; latn=35 ;* latn=47.5
out='ncep.snapshot.test.eps'
quitopt='no'
lev=1000
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 = '-input');infle = 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 = '-time' ); time = subwrd(args,i) ;i=i+1; break;endif
if( arg = '-ci' ); ci = subwrd(args,i) ;i=i+1; break;endif
if( arg = '-outdir' ); outdir = 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
say 'Open 'infle; 'open 'infle
say
'q ctlinfo'
say result
#'set lev 'lev ;#/100
'set time 'time
'q dims 3'
say result
line=sublin(result,5)
datetime=subwrd(line,6)
#line=sublin(result,4)
#pdb=subwrd(line,6)
pha=lev
say pha
say
say 'lon 'lonw' 'lone
say 'lat 'lats' 'latn
say
'set lon 'lonw' 'lone
'set lat 'lats' 'latn
'cc'
xmin=1; xmax=7
ymin=2; ymax=9
'set parea ' xmin ' ' xmax ' ' ymin ' ' ymax
'set lon 'lonw' 'lone
'set lat 'lats' 'latn
'set mpdset hires'
'set map 1 1 2'
'set xlopts 1 5 0.2'
'set ylopts 1 5 0.2'
'set grid off'
### Color shade
'set gxout shade2'
#'w=-w/11.769' ;# 11.769=rho*g
wmin=-20 ;#-0.05
wmax=20 ;#0.05
wint=5;#0.01
'color ' wmin ' ' wmax ' ' wint ' -kind blue->white->red'
'tc=' TMP2m '-273'
'set xlab off'
'set ylab off'
'd 'tc ;#w
'xcbar 7.6 7.7 3 8'
### contour
'set gxout contour'
'set ccolor 0'
'set cthick 2'
'set cint 'ci
'set clab off'
'set xlab off'
'set ylab off'
'd 'tc
### vector
'set gxout vector'
'set ccolor 0'
'set cthick 7'
'set xlab off'
'set ylab off'
'set arrscl .5 50'
'd skip(UGRD10m,1,1);skip(VGRD10m,1,1)'
'set ccolor 1'
'set cthick 1'
'set xlab on'
'set ylab on'
'set xlint 5'
'set ylint 5'
'set arrscl .5 50'
'd skip(UGRD10m,1,1);skip(VGRD10m,1,1)'
# Text
'set strsiz 0.10 0.12'
'set string 1 c 5 0'
'draw string 5.6 1.85 m/s'
'draw title 'datetime' SFC'
'gxprint 'outdir '/' out
if ( quitopt = "yes" ); quit; endif
return
----------------------
End of JRA-55_T.snapshot.gs
----------------------
======================
anl_p125_ugrd.ctl
======================
dset /work4/data/JRA-55/anl_p125_ugrd.2016010100
index /work4/data/JRA-55/anl_p125_ugrd.2016010100.idx
undef 9.999E+20
title anl_p125_ugrd.2016010100
* produced by grib2ctl v0.9.13
* command line options: -ncep_rean -no_prs anl_p125_ugrd.2016010100
dtype grib 255
options yrev
ydef 145 linear -90.000000 1.25
xdef 288 linear 0.000000 1.250000
tdef 1 linear 00Z01jan2016 1mo
* z has 37 levels, for prs
zdef 37 levels
1000 975 950 925 900 875 850 825 800 775 750 700 650 600 550 500 450 400 350 300 250 225 200 175 150 125 100 70 50 30 20 10 7 5 3 2 1
vars 1
UGRD 37 33,100,0 ** (profile) u wind [m/s]
ENDVARS
----------------------
End of anl_p125_ugrd.ctl
----------------------
======================
anl_p125_vgrd.ctl
======================
dset /work4/data/JRA-55/anl_p125_vgrd.2016010100
index /work4/data/JRA-55/anl_p125_vgrd.2016010100.idx
undef 9.999E+20
title anl_p125_vgrd.2016010100
* produced by grib2ctl v0.9.13
* command line options: -ncep_rean -no_prs anl_p125_vgrd.2016010100
dtype grib 255
options yrev
ydef 145 linear -90.000000 1.25
xdef 288 linear 0.000000 1.250000
tdef 1 linear 00Z01jan2016 1mo
* z has 37 levels, for prs
zdef 37 levels
1000 975 950 925 900 875 850 825 800 775 750 700 650 600 550 500 450 400 350 300 250 225 200 175 150 125 100 70 50 30 20 10 7 5 3 2 1
vars 1
VGRD 37 34,100,0 ** (profile) v wind [m/s]
ENDVARS
----------------------
End of anl_p125_vgrd.ctl
----------------------
======================
anl_p125_vvel.ctl
======================
dset /work4/data/JRA-55/anl_p125_vvel.2016010100
index /work4/data/JRA-55/anl_p125_vvel.2016010100.idx
undef 9.999E+20
title anl_p125_vvel.2016010100
* produced by grib2ctl v0.9.13
* command line options: -ncep_rean -no_prs anl_p125_vvel.2016010100
dtype grib 255
options yrev
ydef 145 linear -90.000000 1.25
xdef 288 linear 0.000000 1.250000
tdef 1 linear 00Z01jan2016 1mo
* z has 37 levels, for prs
zdef 37 levels
1000 975 950 925 900 875 850 825 800 775 750 700 650 600 550 500 450 400 350 300 250 225 200 175 150 125 100 70 50 30 20 10 7 5 3 2 1
vars 1
VVEL 37 39,100,0 ** (profile) Pressure vertical velocity [Pa/s]
ENDVARS
----------------------
End of anl_p125_vvel.ctl
----------------------