<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&amp;tid=36240&amp;vid=663" target="_blank">http://www.esrl.noaa.gov/psd/cgi-bin/db_search/DBListFiles.pl?did=32&amp;tid=36240&amp;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:&quot;Times New Roman&quot;,&quot;serif&quot;">µ)/</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>&#39;</font><font>reinit</font><font>&#39;
&#39;</font><font>set</font><font> </font><font>grid</font><font> </font><font>off</font><font>&#39;
&#39;</font><font>set</font><font> </font><font>grads</font><font> </font><font>off</font><font>&#39;</font><font>

*** for Windows</font><font>
&#39;</font><font>sdfopen</font><font> </font><font>D</font><font>:\Writing\disaster_management\hgt.</font><font>20120711.</font><font>nc&#39; 
&#39;</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>&#39;
&#39;</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>&#39;</font><font>

*** for Linux
*&#39;</font><font>sdfopen</font><font> /media/Lexar/hgt.</font><font>20120711.</font><font>nc&#39;
*&#39;</font><font>sdfopen</font><font> /media/Lexar/<a href="http://daily_mean.nc" target="_blank">daily_mean.nc</a>&#39;

*** for output figure (type: eps)
*&#39;</font><font>enable</font><font> </font><font>print</font><font> </font><font>D</font><font>:\Writing\disaster_management\hgtgap.gmf&#39;
*&#39;</font><font>enable</font><font> </font><font>print</font><font> dif_hgt.eps&#39;</font><font>
&#39;</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>&#39;
&#39;</font><font>set</font><font> </font><font>lat</font><font> </font><font>20</font><font> </font><font>50</font><font>&#39;
&#39;</font><font>set</font><font> </font><font>xlint</font><font> </font><font>5</font><font>&#39;
&#39;</font><font>set</font><font> </font><font>lon</font><font> </font><font>110</font><font> </font><font>160</font><font>&#39;
&#39;</font><font>set</font><font> </font><font>ylint</font><font> </font><font>5</font><font>&#39;</font><font>



*&#39;</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>&#39;
*&#39;</font><font>draw</font><font> </font><font>map</font><font>&#39;
*&#39;</font><font>set</font><font> </font><font>gxout</font><font> </font><font>shaded</font><font>&#39;
*&#39;</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>&#39;
*&#39;</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>&#39;

*&#39;</font><font>set</font><font> </font><font>cthick</font><font> </font><font>5</font><font>&#39;</font><font>



&#39;</font><font>set</font><font> </font><font>cint</font><font> </font><font>0.5</font><font>&#39;
&#39;</font><font>set</font><font> </font><font>clskip</font><font> </font><font>1</font><font>&#39;
&#39;</font><font>set</font><font> </font><font>rbcols</font><font> </font><font>auto</font><font>&#39;
&#39;</font><font>set</font><font> </font><font>strmden</font><font> </font><font>1</font><font>&#39;
&#39;</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> &#39;
&#39;</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>&#39;
&#39;</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>&#39;</font><font>

*&#39;</font><font>set</font><font> </font><font>define</font><font> </font><font>1</font><font>&#39;
*&#39;</font><font>set</font><font> </font><font>t</font><font> </font><font>180</font><font> </font><font>660</font><font>&#39;
*&#39;</font><font>set</font><font> </font><font>define</font><font> </font><font>2</font><font>&#39;</font><font>t</font><font>
*&#39;</font><font>set</font><font> </font><font>t</font><font> </font><font>180</font><font> </font><font>660</font><font>&#39;
*&#39;</font><font>set</font><font> </font><font>define</font><font> </font><font>3</font><font>&#39;
*&#39;</font><font>set</font><font> </font><font>t</font><font> </font><font>180</font><font> </font><font>660</font><font>&#39;
*&#39;</font><font>set</font><font> </font><font>define</font><font> </font><font>4</font><font>&#39;
*&#39;</font><font>set</font><font> </font><font>t</font><font> </font><font>1</font><font>&#39;</font><font>

&#39;</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)&#39;


&#39;</font><font>set</font><font> </font><font>gxout</font><font> </font><font>shaded</font><font>&#39;
&#39;</font><font>d</font><font> stdanomalies&#39;
&#39;cbarn </font><font>1</font><font> </font><font>0.5</font><font> </font><font>9.3</font><font> </font><font>4.5</font><font>&#39;</font><font>
*&#39;</font><font>run</font><font> <a href="http://makebg.gs" target="_blank">makebg.gs</a>&#39;

*** for output figure (type: eps)
*&#39;</font><font>print</font><font>&#39;
*&#39;</font><font>disable</font><font> </font><font>print</font><font>&#39;

*** for outut figure (type: png)</font><font>
&#39;</font><font>printim</font><font> </font><font>D</font><font>:\Writing\disaster_management\std_anomalies_20120711.png white&#39;</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>