[gradsusr] About the calculation of standardized anomalies, Thanks.

[2008 duan]

Dear Everyone,



        Does anyone know how calculate the standard anomalies? [image:
Inline image 1]

        I downloaded the Reanalysis data ---- daily data of geopotential
data between 1981 to 2012 from
http://www.esrl.noaa.gov/psd/cgi-bin/db_search/DBListFiles.pl?did=32&tid=36240&vid=663.
It contains 31 files,

        First, I combined the data between 1981 to 2010 using  Climate Data
Operators (CDO) software. Then I calculated the mean  and standard
deviation using CDO.  Then I calculated sdσ using Grads software. sdσ = (X-
µ)/σ. And the results are as follows.
[image: Inline image 4]
[image: Inline image 3]
[image: Inline image 2][image: Inline image 1]
 *The above figures are standardized anomalies, and the order is 14 July,
13 July, 12 July and 11 July 2012. I am thinking there are some wrongs,
Because I didnot considering the seasonal about 21-day. But I donot know
how to center the 21-day. *
*
*
* So, I want to ask you how to center the 21-day. And my grads script is as
follows.*

'reinit'
'set grid off'
'set grads off'

*** for Windows
'sdfopen D:\Writing\disaster_management\hgt.20120711.nc'
'sdfopen D:\Writing\disaster_management\daily_mean.nc'
'sdfopen D:\Writing\disaster_management\hgt_std_1981_2010.nc'

*** for Linux
*'sdfopen /media/Lexar/hgt.20120711.nc'
*'sdfopen /media/Lexar/daily_mean.nc'

*** for output figure (type: eps)
*'enable print D:\Writing\disaster_management\hgtgap.gmf'
*'enable print dif_hgt.eps'
'set parea 1 9 1 8'
'set lat 20 50'
'set xlint 5'
'set lon 110 160'
'set ylint 5'



*'set map 1 1 1'
*'draw map'
*'set gxout shaded'
*'set ylopts 1 4 0.14'
*'set xlopts 1 4 0.14'

*'set cthick 5'



'set cint 0.5'
'set clskip 1'
'set rbcols auto'
'set strmden 1'
'set clevs  -4 -3 -2 -1 0 1 2 3 4 '
'set xlopts 1 4 0.15'
'set ylopts 1 4 0.15'

*'set define 1'
*'set t 180 660'
*'set define 2't
*'set t 180 660'
*'set define 3'
*'set t 180 660'
*'set define 4'
*'set t 1'

'define stdanomalies=(hgt.1(time=00Z11JUL2012)-hgt.2(time=00Z31DEC2010))/hgt.3(time=00Z31DEC2010)'


'set gxout shaded'
'd stdanomalies'
'cbarn 1 0.5 9.3 4.5'
*'run makebg.gs'

*** for output figure (type: eps)
*'print'
*'disable print'

*** for outut figure (type: png)
'printim D:\Writing\disaster_management\std_anomalies_20120711.png white'


 Thank you very very much!!!

 Best regards,

Duan




-- 
Weili Duan
Disaster Prevention Research Institute (DPRI), Kyoto University, Uji, Kyoto
611-0011, Japan
TEL: +81-774-38-4131 FAX: +81-774-38-4130
E-mail: duan.scut.cn at gmail.com
           duanwl1986 at flood.dpri.kyoto-u.ac.jp<dyuan.weiri.88r at st.kyoto-u.ac.jp><dyuan.weiri.88r at st.kyoto-u.ac.jp>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20130111/361d949b/attachment-0003.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: image/png
Size: 7564 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20130111/361d949b/attachment-0015.png 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: image/png
Size: 7590 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20130111/361d949b/attachment-0016.png 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: image/png
Size: 7650 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20130111/361d949b/attachment-0017.png 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: image/png
Size: 70160 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20130111/361d949b/attachment-0018.png 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: image/png
Size: 7563 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20130111/361d949b/attachment-0019.png