渦位の水平分布

aofd165.bio.mie-u.ac.jp

/work1/am/2016.PolarLow/ERA-I/PV/Horizontal

Wed Jun 15 14:40:51 JST 2016

======================

pv.h.ctl

======================

dset /work4/data/ERA-Interim/ERA-I_110118-110125_Barents_Plev_90W-90E.nc

title

undef -9.99e+33

dtype netcdf

UNPACK scale_factor add_offset

xdef 241 linear -90 0.75

ydef 81 linear 30 0.75

zdef 29 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 70

50

tdef 24 linear 00Z18JAN2011 360mn

vars 6

pv=>pv 29 t,z,y,x Potential vorticity

z=>z 29 t,z,y,x Geopotential

t=>t 29 t,z,y,x Temperature

w=>w 29 t,z,y,x Vertical velocity

u=>u 29 t,z,y,x U component of wind

v=>v 29 t,z,y,x V component of wind

endvars

----------------------

End of pv.h.ctl

----------------------

======================

pv.h.gs

======================

#

'open pv.h.ctl'

header='ERA-I_110118-110125_Barents'

varname='PV'

outdir='Fig'

'mkdir -p 'outdir

'set mproj nps'

'set lat 30 90'

'set lon -90 90'

i=1

while (i <= 24)

'set t ' i

var=pv

'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=950

'set lev 'plev

'sclvar=' var '* 1E6'

'set gxout shade2'

'set grads off'

'color -levs 10 11 12 13 14 15 16 17 -kind white->lightsteelblue->greenyellow->gold->navajowhite->red -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 2 0.7 'plev '[hPa]'

'set clab off'

'd 'sclvar

'circlat 20 90 0.1'

'circlon 10 0.1'

'mul 1 4 1 2 -yini 0.1'

plev=850

'set lev 'plev

'sclvar=' var '* 1E6'

'set gxout shade2'

'set grads off'

'color -levs 3 4 5 6 7 8 9 10 -kind white->lightsteelblue->greenyellow->gold->navajowhite->red -xcbar 1 1.1 3.3 5.5 -fh 0.12 -fw 0.11 -edge triangle -line on'

'draw string 0.9 5.6 PVU'

'draw string 2 3.5 'plev '[hPa]'

'set clab off'

'd 'sclvar

'circlat 20 90 0.1'

'circlon 10 0.1'

'mul 1 4 1 3 -yini 0.1'

plev=300

'set lev 'plev

'sclvar=' var '* 1E6'

'set gxout shade2'

'set grads off'

'color -levs 2 4 6 8 10 12 14 16 -kind white->lightsteelblue->greenyellow->gold->navajowhite->red -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 2 6.2 'plev '[hPa]'

'set clab off'

'd 'sclvar

'circlat 20 90 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 '.png'

say 'Output: ' figfile

'gxprint ' figfile

i=i+1

endwhile

'allclose'

#'quit'

----------------------

End of pv.h.gs

----------------------