<br><div class="gmail_quote">Dear Everyone,<div><br></div><div> </div>
<div><br></div><div> Does anyone know how calculate the standard anomalies? <img src="cid:ii_13c2780599239c2f" alt="Inline image 1"></div><div><br></div><div>
I downloaded the Reanalysis data ---- daily data of geopotential data between 1981 to 2012 from <a href="http://www.esrl.noaa.gov/psd/cgi-bin/db_search/DBListFiles.pl?did=32&tid=36240&vid=663" target="_blank">http://www.esrl.noaa.gov/psd/cgi-bin/db_search/DBListFiles.pl?did=32&tid=36240&vid=663</a>. It contains 31 files, </div>
<div><br></div><div> 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-<span style="font-size:11.0pt;line-height:150%;font-family:"Times New Roman","serif"">µ)/</span>σ. And the results are as follows.</div>
<div><img src="cid:ii_13c27749c2804a84" alt="Inline image 4"></div><div><img src="cid:ii_13c27746e5805665" alt="Inline image 3"></div><div><img src="cid:ii_13c27742fb03eb04" alt="Inline image 2"><img src="cid:ii_13c2773d67be983e" alt="Inline image 1"><br>
</div><div> <b><font size="4">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. </font></b></div>
<div><b><font size="4"><br></font></b></div><div><b><font size="4"> So, I want to ask you how to center the 21-day. And my grads script is as follows.</font></b></div><div><pre><font style="font-family:Consolas,monospace;font-size:17px;line-height:17px"><font>'</font><font>reinit</font><font>'
'</font><font>set</font><font> </font><font>grid</font><font> </font><font>off</font><font>'
'</font><font>set</font><font> </font><font>grads</font><font> </font><font>off</font><font>'</font><font>
*** for Windows</font><font>
'</font><font>sdfopen</font><font> </font><font>D</font><font>:\Writing\disaster_management\hgt.</font><font>20120711.</font><font>nc'
'</font><font>sdfopen</font><font> </font><font>D</font><font>:\Writing\disaster_management\<a href="http://daily_mean.nc" target="_blank">daily_mean.nc</a>'
'</font><font>sdfopen</font><font> </font><font>D</font><font>:\Writing\disaster_management\<a href="http://hgt_std_1981_2010.nc" target="_blank">hgt_std_1981_2010.nc</a>'</font><font>
*** for Linux
*'</font><font>sdfopen</font><font> /media/Lexar/hgt.</font><font>20120711.</font><font>nc'
*'</font><font>sdfopen</font><font> /media/Lexar/<a href="http://daily_mean.nc" target="_blank">daily_mean.nc</a>'
*** for output figure (type: eps)
*'</font><font>enable</font><font> </font><font>print</font><font> </font><font>D</font><font>:\Writing\disaster_management\hgtgap.gmf'
*'</font><font>enable</font><font> </font><font>print</font><font> dif_hgt.eps'</font><font>
'</font><font>set</font><font> </font><font>parea</font><font> </font><font>1</font><font> </font><font>9</font><font> </font><font>1</font><font> </font><font>8</font><font>'
'</font><font>set</font><font> </font><font>lat</font><font> </font><font>20</font><font> </font><font>50</font><font>'
'</font><font>set</font><font> </font><font>xlint</font><font> </font><font>5</font><font>'
'</font><font>set</font><font> </font><font>lon</font><font> </font><font>110</font><font> </font><font>160</font><font>'
'</font><font>set</font><font> </font><font>ylint</font><font> </font><font>5</font><font>'</font><font>
*'</font><font>set</font><font> </font><font>map</font><font> </font><font>1</font><font> </font><font>1</font><font> </font><font>1</font><font>'
*'</font><font>draw</font><font> </font><font>map</font><font>'
*'</font><font>set</font><font> </font><font>gxout</font><font> </font><font>shaded</font><font>'
*'</font><font>set</font><font> </font><font>ylopts</font><font> </font><font>1</font><font> </font><font>4</font><font> </font><font>0.14</font><font>'
*'</font><font>set</font><font> </font><font>xlopts</font><font> </font><font>1</font><font> </font><font>4</font><font> </font><font>0.14</font><font>'
*'</font><font>set</font><font> </font><font>cthick</font><font> </font><font>5</font><font>'</font><font>
'</font><font>set</font><font> </font><font>cint</font><font> </font><font>0.5</font><font>'
'</font><font>set</font><font> </font><font>clskip</font><font> </font><font>1</font><font>'
'</font><font>set</font><font> </font><font>rbcols</font><font> </font><font>auto</font><font>'
'</font><font>set</font><font> </font><font>strmden</font><font> </font><font>1</font><font>'
'</font><font>set</font><font> </font><font>clevs</font><font> -</font><font>4</font><font> -</font><font>3</font><font> -</font><font>2</font><font> -</font><font>1</font><font> </font><font>0</font><font> </font><font>1</font><font> </font><font>2</font><font> </font><font>3</font><font> </font><font>4</font><font> '
'</font><font>set</font><font> </font><font>xlopts</font><font> </font><font>1</font><font> </font><font>4</font><font> </font><font>0.15</font><font>'
'</font><font>set</font><font> </font><font>ylopts</font><font> </font><font>1</font><font> </font><font>4</font><font> </font><font>0.15</font><font>'</font><font>
*'</font><font>set</font><font> </font><font>define</font><font> </font><font>1</font><font>'
*'</font><font>set</font><font> </font><font>t</font><font> </font><font>180</font><font> </font><font>660</font><font>'
*'</font><font>set</font><font> </font><font>define</font><font> </font><font>2</font><font>'</font><font>t</font><font>
*'</font><font>set</font><font> </font><font>t</font><font> </font><font>180</font><font> </font><font>660</font><font>'
*'</font><font>set</font><font> </font><font>define</font><font> </font><font>3</font><font>'
*'</font><font>set</font><font> </font><font>t</font><font> </font><font>180</font><font> </font><font>660</font><font>'
*'</font><font>set</font><font> </font><font>define</font><font> </font><font>4</font><font>'
*'</font><font>set</font><font> </font><font>t</font><font> </font><font>1</font><font>'</font><font>
'</font><font>define</font><font> stdanomalies=(hgt.</font><font>1</font><font>(</font><font>time</font><font>=00Z11JUL2012)-hgt.</font><font>2</font><font>(</font><font>time</font><font>=00Z31DEC2010))/hgt.</font><font>3</font><font>(</font><font>time</font><font>=00Z31DEC2010)'
'</font><font>set</font><font> </font><font>gxout</font><font> </font><font>shaded</font><font>'
'</font><font>d</font><font> stdanomalies'
'cbarn </font><font>1</font><font> </font><font>0.5</font><font> </font><font>9.3</font><font> </font><font>4.5</font><font>'</font><font>
*'</font><font>run</font><font> <a href="http://makebg.gs" target="_blank">makebg.gs</a>'
*** for output figure (type: eps)
*'</font><font>print</font><font>'
*'</font><font>disable</font><font> </font><font>print</font><font>'
*** for outut figure (type: png)</font><font>
'</font><font>printim</font><font> </font><font>D</font><font>:\Writing\disaster_management\std_anomalies_20120711.png white'</font></font>
</pre></div><div><br></div><div> Thank you very very much!!!</div><div><br></div><div> Best regards,</div><div><br></div><div>Duan</div><span class="HOEnZb"><font color="#888888"><div><div><br></div></div></font></span></div>
<br><br clear="all"><div><br></div>-- <br>Weili Duan<br>Disaster Prevention Research Institute (DPRI), Kyoto University, Uji, Kyoto 611-0011, Japan<br>TEL: +81-774-38-4131 FAX: +81-774-38-4130<br>E-mail: <a href="mailto:duan.scut.cn@gmail.com" target="_blank">duan.scut.cn@gmail.com</a><br>
<a href="mailto:dyuan.weiri.88r@st.kyoto-u.ac.jp" target="_blank">duanwl1986@flood.dpri.kyoto-u.ac.jp</a><a href="mailto:dyuan.weiri.88r@st.kyoto-u.ac.jp" target="_blank"></a>