$ srcdump.sh QV.V.975.gs
------------------------------
List of the following files:
------------------------------
QV.V.975.gs
------------------------------
Machine info
------------------------------
localhost
/work2/am/2017.KYUSHU-HOKUBU.HEAVY.RAIN/FNL.QV.V.975
Sun Feb 11 12:46:21 JST 2018
======================
QV.V.975.gs
======================
'reinit'
ctl='fnlgrib2.ctl'
'open 'ctl
say
say 'OPEN 'ctl
say
n=1
ne=1
plev=975
while (n<=ne)
'cc'
'set t 'n
'set lev 'plev
'set lat 15 50'
'set lon 90 170'
#'set lon 40 220'
# https://sites.google.com/site/afcanalysis/home/oanalysis/msm/msm_snapshot
'tc=TMPprs-273.15' ;*tc [degC]
'es=611.0/100.*pow(10, (7.5*'tc')/(237.3+'tc'))' ;*Teten's formula ;* [hPa]
'e=RHprs/100.*es' ;* [hPa]
'ws=621.97*(es/(' plev '-es))' ;*[g/kg]
'w=(RHprs/100*ws)' ;*[g/kg]
'd w'
'q dims'
line=sublin(result,5)
datetime=subwrd(line,6)
hh=substr(datetime,1,2)
dd=substr(datetime,4,2)
mmm=substr(datetime,6,3)
yyyy=substr(datetime,9,4)
timeutc=hh':00UTC' dd mmm yyyy
say timeutc
'draw string 6.7 7.65 'timeutc
# Header
'set strsiz 0.08 0.1'
'set string 1 l 2'
'draw string 0.1 8.35 /work2/am/2017.KYUSHU-HOKUBU.HEAVY.RAIN/FNL.QV.V.975'
'draw string 0.1 8.23 Qv.V.975.gs'
'set strsiz 0.10 0.12'
'draw string 0.1 7.90 NCEP FNL'
if(mmm='JAN');mm='01';endif
if(mmm='FEB');mm='02';endif
if(mmm='MAR');mm='03';endif
if(mmm='APR');mm='04';endif
if(mmm='MAY');mm='05';endif
if(mmm='JUN');mm='06';endif
if(mmm='JUL');mm='07';endif
if(mmm='AUG');mm='08';endif
if(mmm='SEP');mm='09';endif
if(mmm='OCT');mm='10';endif
if(mmm='NOV');mm='11';endif
if(mmm='DEC');mm='12';endif
fig='FNL_QV.V.'level'.'yyyy''mm''dd'_'hh'.eps'
'gxprint 'fig
say fig
say
n=n+1
endwhile ;*n
'allclose'
----------------------
End of QV.V.975.gs
----------------------