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