pv_isentropic

[2017年 1月 20日 金曜日 17:07:14 JST]

[~/2016.PolarLow/ERA-I/PV.Isentropic]

[am@aofd165]

$ ls

ERA-I.PV.Isentropic.110118-26.nc.txt colorcb.gs tlev.ctl

Fig/ isen.gs

PV.Isentropic.gs sfc.ctl

[2017年 1月 20日 金曜日 17:08:00 JST]

[~/2016.PolarLow/ERA-I/PV.Isentropic]

[am@aofd165]

$ srcdump.sh PV.Isentropic.gs sfc.ctl tlev.ctl

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

List of the following files:

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

PV.Isentropic.gs

sfc.ctl

tlev.ctl

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

Machine info

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

aofd165.bio.mie-u.ac.jp

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

Fri Jan 20 17:08:13 JST 2017

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

PV.Isentropic.gs

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

#

# PV distribution on an isentropic surface (constant potential temperature surface)

#

'open sfc.ctl'

'open tlev.ctl'

'q ctlinfo 1'

say result

'q ctlinfo 2'

say result

header='ERA-I.Barents110116-110126'

tlev2=300 ;# Potential temperature

varname='PV'tlev2'K'

outdir='Fig'

'!mkdir -p 'outdir

'set mproj nps'

'set mpvals -38 25 65 87'

if (tlev2 = 300)

cminv=1

cmaxv=6

cintv=1

endif

varlev=cminv' 'cmaxv' 'cintv

xmax = 2

ymax = 2

xwid = 6.0/xmax

ywid = 7.0/ymax

nmap = 1

;ymap = 1

sp = 0

ep = 8

np=sp ;* page (or day)

i=sp*xmax*ymax+1

while (np <= ep )

'cc'

'set mpdset mres'

ymap=1

while (ymap <= ymax)

xmap = 1

while (xmap <= xmax)

'set t ' i

xs = 0.8 + (xwid+0.2)*(xmap-1)

xe = xs + xwid

ye = 10.0 - (ywid+0.10)*(ymap-1)

ys = ye - ywid

#say xs' 'xe' 'ys' 'ye

'set vpage 0.0 8.5 0.0 11.0'

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

kind='white->blue->dodgerblue->skyblue->yellow->darkorange->red->darkred'

#----- Plotting each panel -----

'color ' varlev ' -kind ' kind ' -xcbar 4.0 6.5 2.6 2.7 -fh 0.12 -fw 0.11 -fstep 1 -edge circle -line on'

'set lev 'tlev2

pv2'=pv.2*1.E6'

'set z 1'

'set gxout shade2'

'set clab off'

'set frame off'

'set grads off'

'd 'pv2

'set gxout contour'

'set cthick 5'

'set ccolor 0'

'set clab off'

'set cint 1'

'set frame off'

'set grads off'

'set map 1 1 1'

'd 'pv2

'set gxout contour'

'set cthick 7'

'set ccolor 0'

'set clab off'

'set cint 4'

'set frame off'

'set grads off'

'set map 1 1 1'

'd 'msl.1'/'100

'set cthick 2'

'set ccolor 1'

'set cint 4'

'set clskip 2'

'set clopts 1 2 0.07'

'set clab on'

'set frame on'

'set grads off'

'set map 1 1 1'

'd 'msl.1'/'100

'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)

datetimeutc=hh'00UTC'dd''mmm''yyyy

say np' 'i' 'datetime

'q gxinfo'

line=sublin(result,3)

xl=subwrd(line,4)

xr=subwrd(line,6)

line=sublin(result,4)

yb=subwrd(line,4)

yt=subwrd(line,6)

*

xx = xl+xwid/2

yy = yt+0.2

'set string 1 c 4 0'

'set strsiz 0.13 0.15'

'draw string 'xx' 'yy' 'datetimeutc ;* For AMS format

#----- Plotting each panel -----

i=i+1

'set string 1 c 4 0'

'set strsiz 0.13 0.15'

'draw string 4.25 10.5 'header

'set string 1 c 4 0'

'set strsiz 0.10 0.12'

'draw string 6.8 2.65 [PVU]'

'draw string 5.25 2.85 'varname

'set parea off'

'set vpage off'

if (nmap = 22); break; endif

xmap = xmap + 1

endwhile ;* xmap

ymap = ymap + 1

endwhile ;* ymap

num=''

if (np<10)

num='00'np

endif

if (np>=10 & np<100)

num='0'np

endif

if (np>100)

num=np

endif

yyyymmmdd=yyyy''mmm''dd

figfile=outdir '/' header '_' varname '_' num '_'yyyymmmdd'.eps'

say 'Output: ' figfile

'gxprint ' figfile

np=np+1

endwhile ;* np

'allclose'

exit

#'quit'

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

End of PV.Isentropic.gs

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

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

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 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

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

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

tlev.ctl

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

dset /work4/data/ERA-Interim/ERA-I.PV.Isentropic.110116-26.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 3 levels 275 300 330

tdef 44 linear 00Z16JAN2011 360mn

vars 1

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

endvars

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

End of tlev.ctl

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

$ srcdump.sh /work4/data/ERA-Interim/get_ERA-I.Isentropic.sh

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

List of the following files:

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

/work4/data/ERA-Interim/get_ERA-I.Isentropic.sh

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

Machine info

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

aofd165.bio.mie-u.ac.jp

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

Fri Jan 20 17:06:30 JST 2017

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

/work4/data/ERA-Interim/get_ERA-I.Isentropic.sh

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

#!/bin/bash

url="http://stream.ecmwf.int/data/atls06/data/data01/scratch/"

nc="_grib2netcdf-atls06-95e2cf679cd58ee9b4db4dd119a05a8d-QEOeTg.nc"

wget $nc

if [ $? -ne 0 ]; then

Error in $0 : wget failed.

exit 1

fi

new="ERA-I.PV.Isentropic.110116-26.nc"

echo

mv -v $nc $new

echo

text=$(basename $new .nc).txt

cat <<END>$text

Final request

Stream:Atmospheric model

Area:90.0°N 180.0°W 30.0°N 180.0°E

Type:Analysis

Level:275, 300, 330

Dataset:interim_daily

Step:0

Version:1

Type of level:Potential temperature

Time:00:00:00, 06:00:00, 12:00:00, 18:00:00

Date:20110116 to 20110126

Grid:0.75° x 0.75°

Parameter:Potential vorticity

Class:ERA Interim

The status of the request is: complete

Download (9.8MB)

Request output:

END

ncdump -h $new >> $text

ncdump -c $new >> $text

echo

echo "netCDF file: $new"

echo "text file: $text"

echo

exit 0

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

End of /work4/data/ERA-Interim/get_ERA-I.Isentropic.sh

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

$ srcdump.sh /work4/data/ERA-Interim/ERA-I.PV.Isentropic.110116-26.txt

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

List of the following files:

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

/work4/data/ERA-Interim/ERA-I.PV.Isentropic.110116-26.txt

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

Machine info

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

aofd165.bio.mie-u.ac.jp

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

Fri Jan 20 17:07:14 JST 2017

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

/work4/data/ERA-Interim/ERA-I.PV.Isentropic.110116-26.txt

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

Final request

Stream:Atmospheric model

Area:90.0°N 180.0°W 30.0°N 180.0°E

Type:Analysis

Level:275, 300, 330

Dataset:interim_daily

Step:0

Version:1

Type of level:Potential temperature

Time:00:00:00, 06:00:00, 12:00:00, 18:00:00

Date:20110116 to 20110126

Grid:0.75° x 0.75°

Parameter:Potential vorticity

Class:ERA Interim

The status of the request is: complete

Download (9.8MB)

Request output:

netcdf ERA-I.PV.Isentropic.110116-26 {

dimensions:

longitude = 480 ;

latitude = 81 ;

level = 3 ;

time = UNLIMITED ; // (44 currently)

variables:

float longitude(longitude) ;

longitude:units = "degrees_east" ;

longitude:long_name = "longitude" ;

float latitude(latitude) ;

latitude:units = "degrees_north" ;

latitude:long_name = "latitude" ;

int level(level) ;

level:long_name = "levelist" ;

int time(time) ;

time:units = "hours since 1900-01-01 00:00:0.0" ;

time:long_name = "time" ;

time:calendar = "gregorian" ;

short pv(time, level, latitude, longitude) ;

pv:scale_factor = 2.71909973766549e-09 ;

pv:add_offset = 2.8075222107064e-05 ;

pv:_FillValue = -32767s ;

pv:missing_value = -32767s ;

pv:units = "K m**2 kg**-1 s**-1" ;

pv:long_name = "Potential vorticity" ;

// global attributes:

:Conventions = "CF-1.6" ;

:history = "2017-01-20 02:17:45 GMT by grib_to_netcdf-2.0.2: grib_to_netcdf /data/data01/scratch/_mars-atls03-95e2cf679cd58ee9b4db4dd119a05a8d-ldf6DE.grib -o /data/data01/scratch/_grib2netcdf-atls06-95e2cf679cd58ee9b4db4dd119a05a8d-QEOeTg.nc -utime" ;

}

netcdf ERA-I.PV.Isentropic.110116-26 {

dimensions:

longitude = 480 ;

latitude = 81 ;

level = 3 ;

time = UNLIMITED ; // (44 currently)

variables:

float longitude(longitude) ;

longitude:units = "degrees_east" ;

longitude:long_name = "longitude" ;

float latitude(latitude) ;

latitude:units = "degrees_north" ;

latitude:long_name = "latitude" ;

int level(level) ;

level:long_name = "levelist" ;

int time(time) ;

time:units = "hours since 1900-01-01 00:00:0.0" ;

time:long_name = "time" ;

time:calendar = "gregorian" ;

short pv(time, level, latitude, longitude) ;

pv:scale_factor = 2.71909973766549e-09 ;

pv:add_offset = 2.8075222107064e-05 ;

pv:_FillValue = -32767s ;

pv:missing_value = -32767s ;

pv:units = "K m**2 kg**-1 s**-1" ;

pv:long_name = "Potential vorticity" ;

// global attributes:

:Conventions = "CF-1.6" ;

:history = "2017-01-20 02:17:45 GMT by grib_to_netcdf-2.0.2: grib_to_netcdf /data/data01/scratch/_mars-atls03-95e2cf679cd58ee9b4db4dd119a05a8d-ldf6DE.grib -o /data/data01/scratch/_grib2netcdf-atls06-95e2cf679cd58ee9b4db4dd119a05a8d-QEOeTg.nc -utime" ;

data:

longitude = -180, -179.25, -178.5, -177.75, -177, -176.25, -175.5, -174.75,

-174, -173.25, -172.5, -171.75, -171, -170.25, -169.5, -168.75, -168,

-167.25, -166.5, -165.75, -165, -164.25, -163.5, -162.75, -162, -161.25,

-160.5, -159.75, -159, -158.25, -157.5, -156.75, -156, -155.25, -154.5,

-153.75, -153, -152.25, -151.5, -150.75, -150, -149.25, -148.5, -147.75,

-147, -146.25, -145.5, -144.75, -144, -143.25, -142.5, -141.75, -141,

-140.25, -139.5, -138.75, -138, -137.25, -136.5, -135.75, -135, -134.25,

-133.5, -132.75, -132, -131.25, -130.5, -129.75, -129, -128.25, -127.5,

-126.75, -126, -125.25, -124.5, -123.75, -123, -122.25, -121.5, -120.75,

-120, -119.25, -118.5, -117.75, -117, -116.25, -115.5, -114.75, -114,

-113.25, -112.5, -111.75, -111, -110.25, -109.5, -108.75, -108, -107.25,

-106.5, -105.75, -105, -104.25, -103.5, -102.75, -102, -101.25, -100.5,

-99.75, -99, -98.25, -97.5, -96.75, -96, -95.25, -94.5, -93.75, -93,

-92.25, -91.5, -90.75, -90, -89.25, -88.5, -87.75, -87, -86.25, -85.5,

-84.75, -84, -83.25, -82.5, -81.75, -81, -80.25, -79.5, -78.75, -78,

-77.25, -76.5, -75.75, -75, -74.25, -73.5, -72.75, -72, -71.25, -70.5,

-69.75, -69, -68.25, -67.5, -66.75, -66, -65.25, -64.5, -63.75, -63,

-62.25, -61.5, -60.75, -60, -59.25, -58.5, -57.75, -57, -56.25, -55.5,

-54.75, -54, -53.25, -52.5, -51.75, -51, -50.25, -49.5, -48.75, -48,

-47.25, -46.5, -45.75, -45, -44.25, -43.5, -42.75, -42, -41.25, -40.5,

-39.75, -39, -38.25, -37.5, -36.75, -36, -35.25, -34.5, -33.75, -33,

-32.25, -31.5, -30.75, -30, -29.25, -28.5, -27.75, -27, -26.25, -25.5,

-24.75, -24, -23.25, -22.5, -21.75, -21, -20.25, -19.5, -18.75, -18,

-17.25, -16.5, -15.75, -15, -14.25, -13.5, -12.75, -12, -11.25, -10.5,

-9.75, -9, -8.25, -7.5, -6.75, -6, -5.25, -4.5, -3.75, -3, -2.25, -1.5,

-0.75, 0, 0.75, 1.5, 2.25, 3, 3.75, 4.5, 5.25, 6, 6.75, 7.5, 8.25, 9,

9.75, 10.5, 11.25, 12, 12.75, 13.5, 14.25, 15, 15.75, 16.5, 17.25, 18,

18.75, 19.5, 20.25, 21, 21.75, 22.5, 23.25, 24, 24.75, 25.5, 26.25, 27,

27.75, 28.5, 29.25, 30, 30.75, 31.5, 32.25, 33, 33.75, 34.5, 35.25, 36,

36.75, 37.5, 38.25, 39, 39.75, 40.5, 41.25, 42, 42.75, 43.5, 44.25, 45,

45.75, 46.5, 47.25, 48, 48.75, 49.5, 50.25, 51, 51.75, 52.5, 53.25, 54,

54.75, 55.5, 56.25, 57, 57.75, 58.5, 59.25, 60, 60.75, 61.5, 62.25, 63,

63.75, 64.5, 65.25, 66, 66.75, 67.5, 68.25, 69, 69.75, 70.5, 71.25, 72,

72.75, 73.5, 74.25, 75, 75.75, 76.5, 77.25, 78, 78.75, 79.5, 80.25, 81,

81.75, 82.5, 83.25, 84, 84.75, 85.5, 86.25, 87, 87.75, 88.5, 89.25, 90,

90.75, 91.5, 92.25, 93, 93.75, 94.5, 95.25, 96, 96.75, 97.5, 98.25, 99,

99.75, 100.5, 101.25, 102, 102.75, 103.5, 104.25, 105, 105.75, 106.5,

107.25, 108, 108.75, 109.5, 110.25, 111, 111.75, 112.5, 113.25, 114,

114.75, 115.5, 116.25, 117, 117.75, 118.5, 119.25, 120, 120.75, 121.5,

122.25, 123, 123.75, 124.5, 125.25, 126, 126.75, 127.5, 128.25, 129,

129.75, 130.5, 131.25, 132, 132.75, 133.5, 134.25, 135, 135.75, 136.5,

137.25, 138, 138.75, 139.5, 140.25, 141, 141.75, 142.5, 143.25, 144,

144.75, 145.5, 146.25, 147, 147.75, 148.5, 149.25, 150, 150.75, 151.5,

152.25, 153, 153.75, 154.5, 155.25, 156, 156.75, 157.5, 158.25, 159,

159.75, 160.5, 161.25, 162, 162.75, 163.5, 164.25, 165, 165.75, 166.5,

167.25, 168, 168.75, 169.5, 170.25, 171, 171.75, 172.5, 173.25, 174,

174.75, 175.5, 176.25, 177, 177.75, 178.5, 179.25 ;

latitude = 90, 89.25, 88.5, 87.75, 87, 86.25, 85.5, 84.75, 84, 83.25, 82.5,

81.75, 81, 80.25, 79.5, 78.75, 78, 77.25, 76.5, 75.75, 75, 74.25, 73.5,

72.75, 72, 71.25, 70.5, 69.75, 69, 68.25, 67.5, 66.75, 66, 65.25, 64.5,

63.75, 63, 62.25, 61.5, 60.75, 60, 59.25, 58.5, 57.75, 57, 56.25, 55.5,

54.75, 54, 53.25, 52.5, 51.75, 51, 50.25, 49.5, 48.75, 48, 47.25, 46.5,

45.75, 45, 44.25, 43.5, 42.75, 42, 41.25, 40.5, 39.75, 39, 38.25, 37.5,

36.75, 36, 35.25, 34.5, 33.75, 33, 32.25, 31.5, 30.75, 30 ;

level = 275, 300, 330 ;

time = 973368, 973374, 973380, 973386, 973392, 973398, 973404, 973410,

973416, 973422, 973428, 973434, 973440, 973446, 973452, 973458, 973464,

973470, 973476, 973482, 973488, 973494, 973500, 973506, 973512, 973518,

973524, 973530, 973536, 973542, 973548, 973554, 973560, 973566, 973572,

973578, 973584, 973590, 973596, 973602, 973608, 973614, 973620, 973626 ;

}

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

End of /work4/data/ERA-Interim/ERA-I.PV.Isentropic.110116-26.txt

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