------------------------------
List of the following files:
------------------------------
aprho.vs.oa3.run.sh
oaflux3_e.ctl
AphroJP_JST_DPREC.ctl
AphroJP_V1207_DPREC.ctl
------------------------------
Machine info
------------------------------
localhost
/work2/am/2017.Heavy.Rain/APHRO.vs.OAFLUX3
Mon Oct 16 12:01:07 JST 2017
======================
aprho.vs.oa3.run.sh
======================
#!/bin/bash
iys=1985
iye=2017 #2015
iy=$iys
export LANG=C
host=$(hostname)
cwd=$(pwd)
timestamp=$(date)
ctl=AphroJP_V1207_DPREC.ctl
ctl2=oaflux3_e.ctl
gs=$(basename $0 .run.sh).gs
cat<<EOF>$gs
*
*
*
function gs( args )
*
* Decode options
*
i = 1
while( 1 )
arg = subwrd( args, i )
i = i + 1;
if( arg = '' ); break; endif
while( 1 )
if( arg = '-ctl' ); ctl=subwrd(args,i);i=i+1; break;endif
if( arg = '-ctl2' ); ctl2=subwrd(args,i);i=i+1; break;endif
if( arg = '-date'); ddmmmyyyy=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
cmy2mmd=10.0/365.0 ;* cm/yr => mm/d
'cc'
'set mpdset hires'
xmax = 4
ymax = 2
xwid = 9.0/xmax
ywid = 5.5/ymax
nmap = 1
ymap = 1
while (ymap <= ymax)
xmap = 1
while (xmap <= xmax)
xs = 0.8 + (xwid+0.40)*(xmap-1)
xe = xs + xwid
ye = 7.3 - (ywid+0.70)*(ymap-1)
ys = ye - ywid
'set xlopts 1 2 0.0'
'set xlopts 1 2 0.0'
if (ymap = ymax)
'set xlopts 1 2 0.1'
else
'set xlopts 1 2 0.0'
'set xlopts 1 2 0.1'
endif
if (xmap = 1)
'set ylopts 1 2 0.1'
else
'set ylopts 1 2 0.0'
endif
if (nmap = 2)
'set ylopts 1 2 0.1'
else
'set ylopts 1 2 0.1'
endif
'set vpage 0.0 11.0 0.0 8.5'
'set parea 'xs ' 'xe' 'ys' 'ye
'set grads off'
'set font 4'
hh='00'
time0=hh'Z'ddmmmyyyy
#----- Plotting 1st panel -----
if (nmap = 1)
'open 'ctl
'set lon 128.5 133.5'
'set lat 30.0 35.0'
xlevs='129 131 133'
'set xlevs 'xlevs
'set ylint 1'
'set time 'time0
'q ctlinfo'
say result
'q dims'
say result
'set gxout shade2'
levs='10 20 50 100 200 300'
kind='white->blue->orange->red'
xcbar='1.5 9.5 0.58 0.7 -fw 0.1 -fh 0.12 -line on -edge circle'
if ( nmap = 1 )
'color -levs ' levs ' -kind ' kind ' -xcbar ' xcbar
'set string 1 l 2 0'
'set strsiz 0.10 0.12'
xx=9.6
yy=0.64
'draw string 'xx' 'yy' R [mm/d] '
endif
'set grads off'
'd precip'
'q gxinfo'
line=sublin(result,3)
xl=subwrd(line,4)
xr=subwrd(line,6)
line=sublin(result,4)
yb=subwrd(line,4)
yt=subwrd(line,6)
*
xx = xl+0.0
yy = yt+0.15
'set string 1 l 2 0'
'set strsiz 0.10 0.12'
'draw string 'xx' 'yy' R 'ddmmmyyyy
'allclose'
endif
#----- Plotting 1st panel -----
#----- Plotting other panels -----
if ( nmap >= 2 )
'open 'ctl2
'set time 'time0
'q dims'
tline=sublin(result,5)
t0=subwrd(tline,9)
#say t0
# Go backwards in time
t=t0-nmap+2
'set t 't
'q dims'
tline=sublin(result,5)
timenow=subwrd(tline,6)
daynow=substr(timenow,4,12)
'set lon 100. 140.0'
'set lat 5.0 40.0'
xlevs='100 110 120 130'
'set xlevs 'xlevs
'set ylint 10'
'set gxout shade2'
levs='5 10 15 20 25 30 40 60 80'
kind='white->blue->green->yellow->red' ;*'white->blue->red'
xcbar='1.5 9.5 0.18 0.3 -fw 0.1 -fh 0.12 -line on -edge circle'
if ( nmap = 2 )
'color -levs ' levs ' -kind ' kind ' -xcbar ' xcbar
'set string 1 l 2 0'
'set strsiz 0.10 0.12'
xx=9.6
yy=0.24
'draw string 'xx' 'yy' E [mm/d] '
else
'color -levs ' levs ' -kind ' kind
endif
'set grads off'
'd evapr*'cmy2mmd
'q gxinfo'
line=sublin(result,3)
xl=subwrd(line,4)
xr=subwrd(line,6)
line=sublin(result,4)
yb=subwrd(line,4)
yt=subwrd(line,6)
*
xx = xl+0.0
yy = yt+0.15
'set string 1 l 2 0'
'set strsiz 0.10 0.12'
'draw string 'xx' 'yy' E 'daynow
'allclose'
endif
'set parea off'
'set vpage off'
*
if (nmap = 22); break; endif
nmap = nmap + 1
xmap = xmap + 1
endwhile ;* xmap
ymap = ymap + 1
endwhile :* ymap
# Header
'set strsiz 0.08 0.1'
'set string 1 l 2'
'draw string 0.1 8.35 ${timestamp}'
'draw string 0.1 8.23 ${host}'
'draw string 0.1 8.11 ${cwd}'
'draw string 0.1 7.99 ${gs}'
'draw string 0.1 7.87 ${ctl}'
'draw string 0.1 7.75 ${ctl2}'
'gxprint 'out
say
say 'OUTPUT: 'out
say
if (quitopt="yes")then
quit
endif
exit
EOF
figdir="Fig"
if [ ! -d $figdir ]; then
mkdir -vp $figdir
fi
while [ $iy -le $iye ]; do
if [ $iy -le 2011 ]; then
ctl=AphroJP_V1207_DPREC.ctl
else
ctl=AphroJP_JST_DPREC.ctl
fi
#echo $iy $ctl
m=6
me=7
while [ $m -le $me ]; do
d=1
if [ $m -eq 7 ]; then
de=31
else
de=30
fi
while [ $d -le $de ]; do
yyyy=${iy}
mm=$(printf %02d $m)
dd=$(printf %02d $d)
if [ $m -eq 1 ]; then mmm="Jan"; fi
if [ $m -eq 2 ]; then mmm="Feb"; fi
if [ $m -eq 3 ]; then mmm="Mar"; fi
if [ $m -eq 4 ]; then mmm="Apr"; fi
if [ $m -eq 5 ]; then mmm="May"; fi
if [ $m -eq 6 ]; then mmm="Jun"; fi
if [ $m -eq 7 ]; then mmm="Jul"; fi
if [ $m -eq 8 ]; then mmm="Aug"; fi
if [ $m -eq 9 ]; then mmm="Sep"; fi
if [ $m -eq 10 ]; then mmm="Oct"; fi
if [ $m -eq 11 ]; then mmm="Nov"; fi
if [ $m -eq 12 ]; then mmm="Dec"; fi
ddmmmyyyy=${dd}${mmm}${yyyy}
ext="eps"
out="${figdir}/APHRO.vs.OA3_${yyyy}${mm}${dd}.${ext}"
opt="-ctl ${ctl} -ctl2 ${ctl2} -date ${ddmmmyyyy} -out ${out} -q"
grads -bcl "${gs} ${opt}"
d=$(expr $d + 1)
done
m=$(expr $m + 1)
done
iy=$(expr $iy + 1)
done
----------------------
End of aprho.vs.oa3.run.sh
----------------------
======================
oaflux3_e.ctl
======================
DSET /misc/raid112MANDA/KakenA.H28/OAFLUX3/evapr_oaflux_%y4.nc
OPTIONS TEMPLATE
DTYPE netcdf
TITLE WHOI Objectively Analyzed air-sea Fluxes Project
UNDEF 32766
XDEF 360 linear 0.5 1.0
YDEF 180 linear -89.5 1.0
ZDEF 1 levels 1
TDEF 42000 LINEAR 01JAN1985 1dy
VARS 2
evapr=>evapr 1 t, y, x daily mean evaporation rate cm yr-1
err=>err 1 t, y, x estimated error of analysed evaporation rate cm yr-1
ENDVARS
----------------------
End of oaflux3_e.ctl
----------------------
======================
AphroJP_JST_DPREC.ctl
======================
DSET /misc/raid112MANDA/KakenA.H28/APHRODITE/V1707/AphroJP_JST_DPREC.%y4
UNDEF -99.9
OPTIONS little_endian template
TITLE APHRO_JP V1207 daily precipitation
XDEF 460 LINEAR 123.025 0.05
YDEF 440 LINEAR 24.025 0.05
ZDEF 01 LEVELS 1
TDEF 42000 LINEAR 01JAN2009 1dy
VARS 3
precip 1 00 Daily precipitation (mm/day)
nobs 1 00 Number of stations in each grid [1 or 2]
eod 1 00 End of the day (UTC) = 15 UTC
ENDVARS
----------------------
End of AphroJP_JST_DPREC.ctl
----------------------
======================
AphroJP_V1207_DPREC.ctl
======================
DSET /misc/raid112MANDA/KakenA.H28/APHRODITE/V1207/AphroJP_V1207_DPREC.%y4
UNDEF -99.9
OPTIONS little_endian template
TITLE APHRO_JP V1207 daily precipitation
XDEF 460 LINEAR 123.025 0.05
YDEF 440 LINEAR 24.025 0.05
ZDEF 01 LEVELS 1
TDEF 42000 LINEAR 01JAN1900 1dy
VARS 1
precip 1 00 daily precipitation (mm/day)
ENDVARS
----------------------
End of AphroJP_V1207_DPREC.ctl
----------------------