$ srcdump.sh qa.sh qa.gs qa.ctl
aofd165.bio.mie-u.ac.jp
/work1/am/Specific.Humidity_NearSurface/J-OFURO3
Fri Jul 15 09:52:21 JST 2016
======================
qa.sh
======================
#!/bin/bash
yyyy=2008
mm=08
mmm=AUG
dd=25
date=${dd}${mmm}${yyyy}
dataset=J-OFURO3
datadir=/work4/data/${dataset}/Specific.Humidity
ctl=$(basename $0 .sh).ctl
gs=$(basename $0 .sh).gs
cat <<EOF>$ctl
dset ${datadir}/J-OFURO3_QA_MULTI_V0.8_DAILY_${yyyy}.nc
title
undef -9999
dtype netcdf
xdef 1440 linear 0.125 0.25
ydef 720 linear -89.875 0.25
zdef 1 linear 0 1
tdef 366 linear 00Z01JAN${yyyy} 1440mn
vars 1
QA=>qa 0 t,y,x J-OFURO3 Multi QA V0.8
endvars
EOF
outdir='Fig'
mkdir -p $outdir
cat <<EOF>$gs
'open '$ctl
'cc'
'set mpdset hires'
'set lon 115 150'
'set lat 20 45'
date=${date}
'set time ' date
'set gxout shade2'
'color 15 22 1 -kind white->grainbow -xcbar 2 7 1.2 1.4 -fh 0.12 -fw 0.11 -edge triangle -line on'
'd 'qa
'set gxout contour'
'set clevs 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29'
'set clab off'
'set cthick 10'
'set ccolor 0'
'd 'qa
'set clevs 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29'
'set clab on'
'set cthick 1'
'set ccolor 1'
'd 'qa
'q dims'
line=sublin(result,5)
datetime=subwrd(line,6)
'draw title ${dataset} 'date
figfile=${outdir} '/qa2m.' ${yyyy}''${mm}''${dd} '_${dataset}.eps'
say
say figfile
say
'gxprint 'figfile
'allclose'
'quit'
EOF
grads -bcp "$gs"
exit 0
----------------------
End of qa.sh
----------------------
======================
qa.gs
======================
'open 'qa.ctl
'cc'
'set mpdset hires'
'set lon 115 150'
'set lat 20 45'
date=25AUG2008
'set time ' date
'set gxout shade2'
'color 15 22 1 -kind white->grainbow -xcbar 2 7 1.2 1.4 -fh 0.12 -fw 0.11 -edge triangle -line on'
'd 'qa
'set gxout contour'
'set clevs 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29'
'set clab off'
'set cthick 10'
'set ccolor 0'
'd 'qa
'set clevs 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29'
'set clab on'
'set cthick 1'
'set ccolor 1'
'd 'qa
'q dims'
line=sublin(result,5)
datetime=subwrd(line,6)
'draw title J-OFURO3 'date
figfile=Fig '/qa2m.' 2008''08''25 '_J-OFURO3.eps'
say
say figfile
say
'gxprint 'figfile
'allclose'
'quit'
----------------------
End of qa.gs
----------------------
======================
qa.ctl
======================
dset /work4/data/J-OFURO3/Specific.Humidity/J-OFURO3_QA_MULTI_V0.8_DAILY_2008.nc
title
undef -9999
dtype netcdf
xdef 1440 linear 0.125 0.25
ydef 720 linear -89.875 0.25
zdef 1 linear 0 1
tdef 366 linear 00Z01JAN2008 1440mn
vars 1
QA=>qa 0 t,y,x J-OFURO3 Multi QA V0.8
endvars
----------------------
End of qa.ctl
----------------------