[2016年 11月 11日 金曜日 18:09:12 JST]
[~/TEACHING/Tsu.Observatory.161115]
[am@aofd165]
$ cd /work4/data
$ mkdir CMAP
$ cd CMAP
$ wget ftp://ftp.cdc.noaa.gov/Datasets/cmap/std/precip.mon.ltm.nc
$ cd /work1/am/TEACHING/Tsu.Observatory.161115
$ grads -p
ga-> sdfopen /work4/data/CMAP/precip.mon.ltm.nc
ga-> q ctlinfo
dset /work4/data/CMAP/precip.mon.ltm.nc
title CPC Merged Analysis of Precipitation (excludes NCEP Reanalysis)
undef -9.96921e+36
dtype netcdf
xdef 144 linear 1.25 2.5
ydef 72 linear -88.75 2.5
zdef 1 linear 0 1
tdef 12 linear 00Z01JAN0001 1mo
vars 2
precip=>precip 0 t,y,x Long Term Mean Average Monthly Rate of Precipitation
valid_yr_count=>valid_yr_count 0 t,y,x count of non-missing values used in mean
endvars
ga-> cobesst.cmap.gs
Open cmapltm.ctl
white->blue
clevs= 5 10
ccols= 16 17 18
-> xcbar 0 7 5 5.2 -levcol 16 5 17 10 18
Open cobesst2ltm.ctl
Output: cobesst2-cmap_annual_mean.eps
closing 0 files
$ srcdump.sh cmapltm.ctl cobesst2ltm.ctl cobesst.cmap.gs
------------------------------
List of the following files:
------------------------------
cmapltm.ctl
cobesst2ltm.ctl
cobesst.cmap.gs
------------------------------
Machine info
------------------------------
aofd165.bio.mie-u.ac.jp
/work1/am/TEACHING/Tsu.Observatory.161115
Fri Nov 11 20:01:18 JST 2016
======================
cmapltm.ctl
======================
dset /work4/data/CMAP/precip.mon.ltm.nc
title CPC Merged Analysis of Precipitation (excludes NCEP Reanalysis)
undef -9.96921e+36
dtype netcdf
xdef 144 linear 1.25 2.5
ydef 72 linear -88.75 2.5
zdef 1 linear 0 1
tdef 12 linear 00Z01JAN0001 1mo
vars 2
precip=>precip 0 t,y,x Long Term Mean Average Monthly Rate of Precipitation
valid_yr_count=>valid_yr_count 0 t,y,x count of non-missing values used in mean
endvars
----------------------
End of cmapltm.ctl
----------------------
======================
cobesst2ltm.ctl
======================
dset /work4/data/COBE-SST/sst.mon.mean.nc
title created 12/2013 from data provided by JRA
undef 1e+20
dtype netcdf
xdef 360 linear 0.5 1
ydef 180 linear -89.5 1
zdef 1 linear 0 1
tdef 1992 linear 00Z01JAN1850 1mo
vars 1
sst=>sst 0 t,y,x Monthly Means of Global Sea Surface Temperature
endvars
----------------------
End of cobesst2ltm.ctl
----------------------
======================
cobesst.cmap.gs
======================
# 格子が一致する範囲がある場合、その範囲を指定すればよい。
#
# 解像度が異なる場合、演算はできない。
# 同時に開いていては重ね描きもできないが、
# 一つずつ描いてcloseすれば重ねられる。
ctl1='cmapltm.ctl'
ctl2='cobesst2ltm.ctl'
say
say 'Open 'ctl1
'open 'ctl1
say
'cc'
'set parea 1 8 6 9'
'set mproj scaled'
'set lat -25 25'
'set grid off'
'set xlab off'
'set ylab off'
'color -levs 5 10 -gxout shaded -kind white->blue -xcbar 0 7 5 5.2'
'set gxout shade2'
'd ave(precip,t=1,t=12)'
'close 1'
say
say 'Open 'ctl2
'open 'ctl2
say
'set parea 1 8 6 9'
'set mproj scaled'
'set lat -25 25'
'set gxout contour'
'set ccolor 0'
'set cthick 10'
'set clevs 27 28 29'
'd ave(sst,t=1,t=12)'
'set ccolor 2'
'set cthick 5'
'set clevs 27 28 29'
'set xlab on'
'set ylab on'
'd ave(sst,t=1,t=12)'
'basemap L'
'close 1'
fig='cobesst2-cmap_annual_mean.eps'
say
say 'Output: 'fig
say
'gxprint 'fig
'allclose'
exit
----------------------
End of cobesst.cmap.gs
----------------------