SLP, 層厚, Z500, T500, Z300, PV300

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

List of the following files:

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

map.panel.2.gs

map.panel.2.run.sh

plev.ctl

sfc.ctl

color.map.panel.2.gs

color.map.panel.2.run.sh

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

Machine info

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

localhost

/work2/am/2016.PolarLow/ERA-I/map.panel

Thu Nov 30 11:30:03 JST 2017

スクリプト（白黒）

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

map.panel.2.gs

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

#

# SFC 500 300

# SLP Z500-Z1000 Z T Z PV

function profile( args )

'reinit'

# red grn blu

'set rgb 16 50 50 50'

*

* Decode options

*

i = 1

while( 1 )

arg = subwrd( args, i )

i = i + 1;

if( arg = '' ); break; endif

while( 1 )

if( arg = '-cwd' ); cwd= subwrd(args,i);i=i+1; break;endif

if( arg = '-host' ); host= subwrd(args,i);i=i+1; break;endif

if( arg = '-tstamp' ); tstamp= subwrd(args,i);i=i+1; break;endif

if( arg = '-gs' ); gs= subwrd(args,i);i=i+1; break;endif

if( arg = '-q' ); quitopt=yes ;i=i+1; break; endif

say 'Syntax error : arg= 'arg

return

endwhile

endwhile

say

say 'datetime = ' datetime

say

'open sfc.ctl'

'open plev.ctl'

'q ctlinfo 1'

say result

'q ctlinfo 2'

say result

header='ERA-I.Barents110116-110126'

plev1=500

plev2=300

varname='SFC.P+Thickness_'plev1'.Z+T_'plev2'.Z+PV'

outdir='Fig'

'!mkdir -p 'outdir

'set mproj nps'

'set mpvals -50 45 60 85'

#'set mpvals -38 12 62 85'

i=13

nt=38 ;*40

while (i <= nt )

'cc'

'set vpage 0.0 8.5 0.0 11.0'

'set grads off'

xmax = 3

ymax = 4

xwid = 7.5/xmax

ywid = 6.0/ymax ;*9.0/ymax

nmap = 1

ymap = 1

while (ymap <= ymax)

xmap = 1

while (xmap <= xmax)

'set t ' i

'set vpage 0.0 8.5 0.0 11.0'

xs = 0.5 + (xwid+0.10)*(xmap-1)

xe = xs + xwid

ye = 10 - (ywid+0.30)*(ymap-1)

ys = ye - ywid

'set parea 'xs ' 'xe' 'ys' 'ye

'set mpdset mres'

'set poli off'

if (xmap = 1)

'set gxout contour'

'set dfile 1'

'set z 1'

var=msl.1 '/100'

'set cmin 965'

'set cmax 1030'

'set cint 5'

'set grads off'

'set frame off'

'set clab off'

'set ccolor 0'

'set cthick 3'

'set mpdraw on'

'set frame on'

'set map 16 1 1'

'd 'var

'set clab forced'

'set cint 5'

'set clskip 2'

'set clopts 1 2 0.05'

'set ccolor 1'

'set cthick 1'

'set grid off'

'set mpdraw off'

'd 'var

'set dfile 2'

'set lev 500'

'var1=z/98.0' ;*GPH in deca meters

'set lev 1000'

'var2=z/98.0' ;*GPH in deca meters

'set gxout contour'

'set cmin 484'

'set cmax 548'

'set cint 4'

'set cstyle 3'

'set cthick 2'

'set ccolor 1'

'set clab forced'

'set clskip 2'

'set grads off'

'set frame off'

'set grid on'

'd 'var1'-'var2

endif

if (xmap = 2)

'set lev 500'

var=t.2 '-273'

'set gxout shade2'

'set clevs -80 -40 0'

'set ccols 15 15 0'

'set mpdraw on'

'set map 16 1 1'

'd 'var

'set gxout contour'

'set ccolor 1'

'set cthick 1'

'set cmin -40'

'set cmax 10'

'set cint 2'

'set clab forced'

'set clskip 4'

'set clopts 1 1 0.05'

'set grads off'

'set frame off'

'd 'var

var=z.2 '/98.0' ;*GPH in deca meters

'set gxout contour'

'set grads off'

'set frame off'

'set map 16 1 1'

'set cmin 450'

'set cmax 600'

'set cint 10'

'set clab off'

'set ccolor 0'

'set cthick 3'

'set grads off'

'set frame off'

'd 'var

'set cint 10'

'set ccolor 1'

'set cthick 1'

'set clab forced'

'set clskip 2'

'set clopts 1 2 0.05'

'set grads off'

'set frame on'

'set mpdraw off'

'd 'var

endif

if (xmap = 3)

'set lev 300'

var=pv.2 '*1.E6'

'set gxout shade2'

'set clevs 0 6 15'

'set ccols 0 0 15'

'set mpdraw on'

'set map 16 1 1'

'd 'var

'set gxout contour'

'set cmin 1'

'set cmax 15'

'set cint 2'

'set ccolor 1'

'set cthick 1'

'set clab forced'

'set clskip 1'

'set cstyle 3'

'set clopts 1 2 0.05'

'set grads off'

'set mpdraw on'

'set map 16 1 1'

'set frame off'

'd 'var

#if (xmap = 3 & ymap = 4)

#xcopt='-xcbar 6.2 8.2 2 2.1 -edge circle -line on'

#'color 'copt' 'xcopt

#endif

var=z.2 '/98.0' ;*GPH in deca meters

'set gxout contour'

'set grads off'

'set frame off'

'set cmin 700'

'set cmax 900'

'set cint 10'

'set clab off'

'set ccolor 0'

'set cthick 5'

'set grads off'

'set frame off'

'set frame on'

'set mpdraw off'

'd 'var

'set ccolor 1'

'set cthick 1'

'set clab forced'

'set cint 10'

'set clskip 2'

'set clopts 1 2 0.05'

'set grads off'

'set frame off'

'set mpdraw off'

'd 'var

endif

if (nmap = 22); break; endif

nmap = nmap + 1

xmap = xmap + 1

endwhile ;* xmap

'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'UTC' dd mmm yyyy

'set strsiz 0.11 0.13'

'set string 1 c 2 90'

say timeutc' Row='ymap

ye = 10 - (ywid+0.30)*(ymap-1)

yc = ys + ywid/2

'draw string 0.2 ' yc ' 'timeutc

'set string 1 c 2 0'

i=i+1

ymap = ymap + 1

endwhile :* ymap

'set vpage 0.0 8.5 0.0 11.0'

'set strsiz 0.11 0.13'

yc=10.2

xmap=1

xs = 0.5 + (xwid+0.10)*(xmap-1)

xc = xs + xwid/2

'draw string 'xc' ' yc' SLP Z1000-Z500'

xmap=2

xs = 0.5 + (xwid+0.10)*(xmap-1)

xc = xs + xwid/2

'draw string 'xc' ' yc' Z&T ' plev1'hPa'

xmap=3

xs = 0.5 + (xwid+0.10)*(xmap-1)

xc = xs + xwid/2

'draw string 'xc' ' yc' Z&PV ' plev2'hPa'

* Footer

'set strsiz 0.08 0.1'

'set string 1 l 2'

'draw string 0.1 1.35 'tstamp' 'host

'draw string 0.1 1.23 'cwd' 'gs

'set strsiz 0.10 0.12'

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 '/' gs '_' varname '_'dd''mmm''yyyy'.eps'

say 'Output: ' figfile

say

'gxprint ' figfile

endwhile ;*i

'allclose'

if (quitopt="yes")then

quit

endif

exit

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

End of map.panel.2.gs

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

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

map.panel.2.run.sh

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

#!/bin/bash

gs=$(basename $0 .run.sh).gs

export LANG=C

host=$(hostname)

cwd=$(pwd)

timestamp=$(date)

host=$(echo $host |tr ' ' '_')

cwd=$(echo $cwd |tr ' ' '_')

timestamp=$(echo $timestamp |tr ' ' '_')

script=$(basename $gs .gs)

grads -bcp "${gs} -tstamp $timestamp -host $host -cwd $cwd -gs $script"

exit 0

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

End of map.panel.2.run.sh

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

コントロールファイル

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

plev.ctl

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

dset /work4/data/ERA-Interim/PL1101/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/PL1101/ERA-I.Barents_110116-26_sfc.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 1 linear 0 1

tdef 44 linear 00Z16JAN2011 360mn

vars 13

ci=>ci 0 t,y,x Sea-ice cover

sst=>sst 0 t,y,x Sea surface temperature

p55.162=>p55.162 0 t,y,x Vertical integral of water vapour

p56.162=>p56.162 0 t,y,x Vertical integral of cloud liquid water

sp=>sp 0 t,y,x Surface pressure

tcwv=>tcwv 0 t,y,x Total column water vapour

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

lcc=>lcc 0 t,y,x Low cloud cover

hcc=>hcc 0 t,y,x High cloud cover

endvars

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

End of sfc.ctl

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

カラー

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

color.map.panel.2.gs

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

#

# SFC 500 300

# SLP Z500-Z1000 Z T Z PV

function profile( args )

'reinit'

# red grn blu

'set rgb 16 50 50 50'

*

* Decode options

*

i = 1

while( 1 )

arg = subwrd( args, i )

i = i + 1;

if( arg = '' ); break; endif

while( 1 )

if( arg = '-cwd' ); cwd= subwrd(args,i);i=i+1; break;endif

if( arg = '-host' ); host= subwrd(args,i);i=i+1; break;endif

if( arg = '-tstamp' ); tstamp= subwrd(args,i);i=i+1; break;endif

if( arg = '-gs' ); gs= subwrd(args,i);i=i+1; break;endif

if( arg = '-q' ); quitopt=yes ;i=i+1; break; endif

say 'Syntax error : arg= 'arg

return

endwhile

endwhile

say

say 'datetime = ' datetime

say

'open sfc.ctl'

'open plev.ctl'

'q ctlinfo 1'

say result

'q ctlinfo 2'

say result

header='ERA-I.Barents110116-110126'

plev1=500

plev2=300

varname='SFC.P+Thickness_'plev1'.Z+T_'plev2'.Z+PV'

outdir='Fig'

'!mkdir -p 'outdir

'set mproj nps'

'set mpvals -55 45 55 85'

#'set mpvals -38 12 62 85'

i=9

nt=34 ;*40

while (i <= nt )

'cc'

'set vpage 0.0 8.5 0.0 11.0'

'set grads off'

xmax = 3

ymax = 4

xwid = 7.5/xmax

ywid = 6.0/ymax ;*9.0/ymax

nmap = 1

ymap = 1

while (ymap <= ymax)

xmap = 1

while (xmap <= xmax)

'set t ' i

'set vpage 0.0 8.5 0.0 11.0'

xs = 0.5 + (xwid+0.10)*(xmap-1)

xe = xs + xwid

ye = 10 - (ywid+0.30)*(ymap-1)

ys = ye - ywid

'set parea 'xs ' 'xe' 'ys' 'ye

'set mpdset mres'

'set poli off'

if (xmap = 1)

'set dfile 2'

'set lev 500'

'var1=z/98.0' ;*GPH in deca meters

'set lev 1000'

'var2=z/98.0' ;*GPH in deca meters

'set gxout shade2'

'set cmin 484'

'set cmax 548'

'set cint 4'

'set grads off'

'set frame on'

'set mpdraw on'

'set map 16 1 1'

'd 'var1'-'var2

'xcbar 0.5 2.9 2.9 3 -line on -edge circle -fw 0.05 -fh 0.07 -ft 1'

'set gxout contour'

'set cmin 484'

'set cmax 548'

'set cint 4'

'set grads off'

'set frame off'

'set clab off'

'set ccolor 0'

'set mpdraw off'

'set cthick 1'

'd 'var1'-'var2

'set gxout contour'

'set dfile 1'

'set z 1'

var=msl.1 '/100'

'set cmin 965'

'set cmax 1030'

'set cint 5'

'set grads off'

'set frame off'

'set clab off'

'set ccolor 0'

'set mpdraw off'

'set cthick 3'

'd 'var

'set clab forced'

'set cint 5'

'set clskip 2'

'set clopts 1 2 0.05'

'set ccolor 1'

'set cthick 1'

'set mpdraw off'

'd 'var

endif

if (xmap = 2)

'set lev 500'

var=t.2 '-273'

'set gxout shade2'

'set clevs -80 -40 0'

'set ccols 15 15 0'

'set mpdraw on'

'set map 16 1 1'

'd 'var

'set gxout contour'

'set ccolor 1'

'set cthick 1'

'set cmin -60'

'set cmax 10'

'set cint 2'

'set clab forced'

'set clskip 4'

'set clopts 1 1 0.05'

'set grads off'

'set frame off'

'd 'var

var=z.2 '/98.0' ;*GPH in deca meters

'set gxout contour'

'set grads off'

'set frame off'

'set map 16 1 1'

'set cmin 450'

'set cmax 600'

'set cint 10'

'set clab off'

'set ccolor 0'

'set cthick 3'

'set grads off'

'set frame off'

'd 'var

'set cint 10'

'set ccolor 1'

'set cthick 1'

'set clab forced'

'set clskip 2'

'set clopts 1 2 0.05'

'set grads off'

'set frame on'

'set mpdraw off'

'd 'var

endif

if (xmap = 3)

'set lev 300'

var=pv.2 '*1.E6'

'set gxout shade2'

'set clevs 0 6 15'

'set ccols 0 0 15'

'set mpdraw on'

'set map 16 1 1'

'd 'var

'set gxout contour'

'set cmin 1'

'set cmax 15'

'set cint 2'

'set ccolor 1'

'set cthick 1'

'set clab forced'

'set clskip 1'

'set cstyle 3'

'set clopts 1 2 0.05'

'set grads off'

'set mpdraw on'

'set map 16 1 1'

'set frame off'

'd 'var

#if (xmap = 3 & ymap = 4)

#xcopt='-xcbar 6.2 8.2 2 2.1 -edge circle -line on'

#'color 'copt' 'xcopt

#endif

var=z.2 '/98.0' ;*GPH in deca meters

'set gxout contour'

'set grads off'

'set frame off'

'set cmin 700'

'set cmax 900'

'set cint 10'

'set clab off'

'set ccolor 0'

'set cthick 5'

'set grads off'

'set frame off'

'set frame on'

'set mpdraw off'

'd 'var

'set ccolor 1'

'set cthick 1'

'set clab forced'

'set cint 10'

'set clskip 2'

'set clopts 1 2 0.05'

'set grads off'

'set frame off'

'set mpdraw off'

'd 'var

endif

if (nmap = 22); break; endif

nmap = nmap + 1

xmap = xmap + 1

endwhile ;* xmap

'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'UTC' dd mmm yyyy

'set strsiz 0.11 0.13'

'set string 1 c 2 90'

say timeutc' Row='ymap

ye = 10 - (ywid+0.30)*(ymap-1)

yc = ys + ywid/2

'draw string 0.2 ' yc ' 'timeutc

'set string 1 c 2 0'

i=i+1

ymap = ymap + 1

endwhile :* ymap

'set vpage 0.0 8.5 0.0 11.0'

'set strsiz 0.11 0.13'

yc=10.2

xmap=1

xs = 0.5 + (xwid+0.10)*(xmap-1)

xc = xs + xwid/2

'draw string 'xc' ' yc' SLP Z1000-Z500'

xmap=2

xs = 0.5 + (xwid+0.10)*(xmap-1)

xc = xs + xwid/2

'draw string 'xc' ' yc' Z&T ' plev1'hPa'

xmap=3

xs = 0.5 + (xwid+0.10)*(xmap-1)

xc = xs + xwid/2

'draw string 'xc' ' yc' Z&PV ' plev2'hPa'

* Footer

'set strsiz 0.08 0.1'

'set string 1 l 2'

'draw string 0.1 1.35 'tstamp' 'host

'draw string 0.1 1.23 'cwd' 'gs

'set strsiz 0.10 0.12'

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 '/' gs '_' varname '_'dd''mmm''yyyy'.eps'

say 'Output: ' figfile

say

'gxprint ' figfile

endwhile ;*i

'allclose'

if (quitopt="yes")then

quit

endif

exit

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

End of color.map.panel.2.gs

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

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

color.map.panel.2.run.sh

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

#!/bin/bash

gs=$(basename $0 .run.sh).gs

export LANG=C

host=$(hostname)

cwd=$(pwd)

timestamp=$(date)

host=$(echo $host |tr ' ' '_')

cwd=$(echo $cwd |tr ' ' '_')

timestamp=$(echo $timestamp |tr ' ' '_')

script=$(basename $gs .gs)

echo ${gs} -tstamp $timestamp -host $host -cwd $cwd -gs $script

grads -bcp "${gs} -tstamp $timestamp -host $host -cwd $cwd -gs $script"

exit 0

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

End of color.map.panel.2.run.sh

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