[2016年 9月 20日 火曜日 16:59:45 JST]
[~/2016.PolarLow/ERA-I/PV.gradT]
[am@aofd165]
$ srcdump.sh pv.gradt.up.gs plev.ctl sfc.ctl
aofd165.bio.mie-u.ac.jp
/work1/am/2016.PolarLow/ERA-I/PV.gradT
Tue Sep 20 16:59:56 JST 2016
======================
pv.gradt.up.gs
======================
'open sfc.ctl'
'open plev.ctl'
#'q ctlinfo 1'
#say result
header='ERA-I.Barents110116-110126'
plev1=925
plev2=400
varname='PV.' plev1 '.' plev2 '_GradT' plev1 '.UP'
outdir='Fig.'varname
'!mkdir -p 'outdir
#'set mproj latlon'
'set mproj nps'
'set lat 50 85'
'set lon -50 50'
i=1
while (i <= 34 )
'set t ' i
'cc'
'set mpdset mres'
'set frame off'
'q dims'
line=sublin(result,5)
datetime=subwrd(line,6)
'set strsiz 0.2 0.25'
'set string 1 c 5'
'draw string 4.5 8.5 'datetime
'mul 1 4 1 1 -yini 1'
plev=plev1
'set z 1'
var2='pv.2(lev=' plev ')*1.E6'
'set gxout shade2'
levs2='0.4 0.6 0.8 1.0 1.2 1.4'
kind2='white->lightskyblue->moccasin->salmon'
'color -levs ' levs2 ' -kind ' kind2 ' -xcbar 1 1.1 0.5 2.7 -fh 0.12 -fw 0.11 -edge triangle -line on'
'draw string 0.9 2.8 PVU'
'draw string 1.8 0.6 PV'plev
'set clab off'
'set grads off'
'd maskout(' var2 ',' sp.1 '-' plev '*100)'
'set gxout contour'
'set clab off'
'set cthick 1'
'set ccolor 0'
'set clevs 'levs2
'set grads off'
'd maskout(' var2 ',' sp.1 '-' plev '*100)'
'set gxout contour'
'set z 1'
levs1='960 965 970 975 980 985 990 995 1000 1005 1010 1015 1020'
'set clevs 'levs1
'set clab off'
'set cthick 5'
'set ccolor 0'
'set grads off'
'd ' msl.1 '/100'
'color -gxout contour -kind indigo->plum -levs ' levs1 ' -xcbar 7.5 7.6 2.3 6.5 -fh 0.12 -fw 0.11 -edge triangle -line on'
'draw string 7.5 6.6 hPa'
'set cthick 2'
'set grid on 5 15'
'set grads off'
'd ' msl.1 '/100'
'circlat 10 55 0.1'
'circlon 10 0.1'
'mul 1 4 1 2 -yini 0.1'
plev=plev1
km100=100000 ;* 100 km -> m
'set z 1'
pt'='t.2 '(lev=' plev1 ')*pow(1000/' plev1 ',0.2857)'
pi=3.14159265359
dtr=pi'/'180
a=6.371e6
dx '=' a '*cos(' dtr '*' lat ')*' dtr '*cdiff(' lon ',' x')'
dy '=' a '*' dtr '*cdiff(' lat ',' y ')'
dtdx '=cdiff(' pt ',' x ')/' dx
dtdy '=cdiff(' pt ',' y ')/' dy
grad '=mag(' dtdx ',' dtdy ')*'km100
'set gxout shade2'
kind1='white->skyblue->yellow->tomato'
'color 1 5 1 -kind ' kind1 ' -xcbar 1 1.1 3.3 5.5 -fh 0.12 -fw 0.11 -edge triangle -line on'
'draw string 0.9 5.6 K/100km'
# 'set strsiz 0.2 0.22'
'draw string 1.8 3.4 PT&`37`0PT'plev1
'set clab off'
'set grads off'
'd maskout(' grad ',' sp.1 '-' plev1 '*100)'
'set gxout contour'
'set clab off'
'set cthick 6'
'set ccolor 0'
'set cstyle 1'
'set cint 4'
'set grads off'
'd maskout(' pt '-273,' sp.1 '-' plev1 '*100)'
'set gxout contour'
'set clab on'
'set cthick 2'
'set ccolor 1'
'set cstyle 1'
'set cint 4'
'set grads off'
'd maskout(' pt '-273,' sp.1 '-' plev1 '*100)'
# 'set gxout contour'
# 'set z 1'
# levs1='960 965 970 975 980 985 990 995 1000 1005 1010 1015 1020'
# 'set clevs 'levs1
# 'set clab off'
# 'set cthick 5'
# 'set ccolor 0'
# 'set grads off'
# 'd ' msl.1 '/100'
# 'color -gxout contour -kind indigo->plum -levs ' levs1 ' -xcbar 7.5 7.6 2.3 6.5 -fh 0.12 -fw 0.11 -edge triangle -line on'
# 'draw string 7.5 6.6 hPa'
# 'set cthick 2'
# 'set grads off'
# 'd ' msl.1 '/100'
'circlat 20 55 0.1'
'circlon 10 0.1'
plev=plev2
'mul 1 4 1 3 -yini 0.1'
var2='pv.2(lev=' plev ')*1.E6'
'set gxout shade2'
'set grads off'
levs2='1 2 3 4 5 6'
kind2='white->lightsteelblue->greenyellow->gold->navajowhite->red'
'color -levs ' levs2 ' -kind ' kind2 ' -xcbar 1 1.1 6.0 8.2 -fh 0.12 -fw 0.11 -edge triangle -line on'
'draw string 0.9 8.3 PVU'
'draw string 1.8 6.1 PV'plev
'set clab off'
'd maskout(' var2 ',' sp.1 '-' plev '*100)'
'set gxout contour'
'set clab off'
'set cthick 1'
'set ccolor 0'
'set grads off'
'set cint 1'
'd maskout(' var2 ',' sp.1 '-' plev '*100)'
'set gxout contour'
'set z 1'
levs1='960 965 970 975 980 985 990 995 1000 1005 1010 1015 1020'
'set clevs 'levs1
'set clab off'
'set cthick 10'
'set ccolor 0'
'set grads off'
'd ' msl.1 '/100'
'color -gxout contour -kind indigo->plum -levs ' levs1 ' -xcbar 7.5 7.6 2.3 6.5 -fh 0.12 -fw 0.11 -edge triangle -line on'
'draw string 7.5 6.6 hPa'
'set cthick 2'
'set grads off'
'd ' msl.1 '/100'
'set gxout contour'
'set lev 'plev
levs1='0.1'
'set cint 'levs1
'set cthick 10'
'set ccolor 0'
'set grads off'
'd ' z.2 '/9.8/1000' ;*z.2' (lev=' plev ')/9.8/1000'
'set gxout contour'
'set clab on'
'set lev 'plev
'set cint 'levs1
'set cthick 4'
'set ccolor 4'
'set grid on 5 15'
'd ' z.2 '/9.8/1000' ;*z.2' (lev=' plev ')/9.8/1000'
'circlat 20 55 0.1'
'circlon 10 0.1'
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 'Output: ' figfile
'gxprint ' figfile
i=i+1
endwhile
'allclose'
exit
#'quit'
----------------------
End of pv.gradt.up.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
----------------------