$ srcdump.sh pv.theta.cross.section.zonal.gs plev.ctl sfc.ctl
aofd165.bio.mie-u.ac.jp
/work1/am/2016.PolarLow/ERA-I/PV.Theta.CrossSection.Zonal
Tue Sep 20 17:01:13 JST 2016
======================
pv.theta.cross.section.zonal.gs
======================
'reinit'
#'open plev.ctl'
'sdfopen /work4/data/ERA-Interim/ERA-I.Barents_110116-26_plev.nc'
'open sfc.ctl'
#'q ctlinfo 1'
#say 'ctlinfo 1'
#say result
#say
#'q ctlinfo 2'
#say 'ctlinfo 2'
#say result
header='ERA-I.Barents110116-110126'
time1='00Z21JAN2011'
lat1=70
lon1=-20
lon2=15
lev1=1000
lev2=300
if(lat1 < 0)
lat1o=-lat1'S'
else
lat1o=lat1'N'
endif
if(lon1 < 0)
lon1o=-lon1'W'
else
lon1o=lon1'E'
endif
if(lon2 < 0)
lon2o=-lon2'W'
else
lon2o=lon2'E'
endif
varname='PV.Theta.Cross.Section.Zonal'
outdir='Fig.'varname
'!mkdir -p 'outdir
'set parea 1 8 1 9'
'set grid off'
'set gxout contour'
'set time ' time1
'cc'
'set mpdset mres'
'set frame off'
'q dims'
line=sublin(result,5)
datetime=subwrd(line,6)
say datetime
line=sublin(result,5)
i=subwrd(line,9)
'set strsiz 0.2 0.25'
'set string 1 c 5'
'draw string 4.5 9.7 'datetime
'set strsiz 0.15 0.2'
'set string 1 c 5'
'draw string 4.5 9.3 'lat1o' 'lon1o'-'lon2o
'set lat 'lat1
'set lon 'lon1' 'lon2
'set lev 'lev1' 'lev2
'define spr=' sp.2 '(z=1)/100'
# 'set gxout shade2'
# 'set clab forced'
# 'd tmp'
# 'cbarn'
# exit
'set zlog on'
levs1='260 265 270 275 280 285 290 295 300 305 310 315 320 325 330'
'set clevs 'levs1
'set clab forced'
'set ccolor 1'
'set cthick 2'
'set cstyle 2'
'set xlab off'
'set ylab off'
'set grads off'
'd ' t.1 '*pow(1000/' lev ',0.2857)'
levs1='0.5 1 1.5 3 4 6 8'
levs2=' 2 '
'set clevs 'levs1
'set clopts 1 4 0.12'
'set clab forced'
'set xlab on'
'set ylab on'
'set xlopts 1 4 0.15'
'set ylopts 1 4 0.15'
'set frame on'
'set cthick 5'
'set cstyle 1'
'set grads off'
# 'd 'pv.1'*1.E6'
'd maskout(' pv.1 '*1.E6,' spr '-' lev ')'
'set clevs 'levs2
'set ccolor 1'
'set cthick 8'
'set cstyle 1'
'set grads off'
# 'd 'pv.1'*1.E6'
'd maskout(' pv.1 '*1.E6,' spr '-' lev ')'
'set gxout contour'
'set clevs 0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30'
'set ccolor 15'
'set cthick 10'
'set clab off'
'd -(spr -' lev ')'
'draw ylab P [hPa]'
num=''
if (i<10)
num='00'i
endif
if (i>=10 & i<100)
num='0'i
endif
if (i>100)
num=i
endif
figfile=outdir '/' header '_' varname '_' num '_' datetime '.png'
say
say 'Output: ' figfile
say
'gxprint ' figfile
'allclose'
exit
#'quit'
----------------------
End of pv.theta.cross.section.zonal.gs
----------------------
======================
plev.ctl
======================
dset /work4/data/ERA-Interim/ERA-I.Barents_110116-26_plev.nc
title
undef -9.99e+33
dtype netcdf
UNPACK scale_factor add_offset
OPTIONS yrev zrev
xdef 480 linear -180 0.75
ydef 81 linear 30 0.75
zdef 27 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
tdef 44 linear 00Z16JAN2011 360mn
vars 10
pv=>pv 27 t,z,y,x Potential vorticity
z=>z 27 t,z,y,x Geopotential
t=>t 27 t,z,y,x Temperature
q=>q 27 t,z,y,x Specific humidity
w=>w 27 t,z,y,x Vertical velocity
r=>r 27 t,z,y,x Relative humidity
clwc=>clwc 27 t,z,y,x Specific cloud liquid water content
ciwc=>ciwc 27 t,z,y,x Specific cloud ice water content
u=>u 27 t,z,y,x U component of wind
v=>v 27 t,z,y,x V component of wind
endvars
----------------------
End of plev.ctl
----------------------
======================
sfc.ctl
======================
dset /work4/data/ERA-Interim/ERA-I.Barents_110116-26_sfc.nc
title
undef -9.99e+33
dtype netcdf
UNPACK scale_factor add_offset
OPTIONS yrev
xdef 480 linear -180 0.75
ydef 81 linear 30 0.75
zdef 1 linear 0 1
tdef 44 linear 00Z16JAN2011 360mn
vars 8
ci=>ci 0 t,y,x Sea-ice cover
sst=>sst 0 t,y,x Sea surface temperature
sp=>sp 0 t,y,x Surface pressure
msl=>msl 0 t,y,x Mean sea level pressure
u10=>u10 0 t,y,x 10 metre U wind component
v10=>v10 0 t,y,x 10 metre V wind component
t2m=>t2m 0 t,y,x 2 metre temperature
d2m=>d2m 0 t,y,x 2 metre dewpoint temperature
endvars
----------------------
End of sfc.ctl
----------------------