<div dir="ltr"><div class="gmail_extra"><div><div dir="ltr">I have a few comments on your script.<br><div><div style="margin-left:80px"><span style="font-size:small"><font color="#38761d"><b></b></font></span></div><div style="margin-left:80px"><span style="font-size:small"><font color="#38761d"><b></b></font></span></div><font size="1"><br></font></div></div></div>
<br><div class="gmail_quote">On 30 October 2014 02:43, Ipshita Majhi <span dir="ltr"><<a href="mailto:ipmajhi@alaska.edu" target="_blank">ipmajhi@alaska.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div dir="ltr"><div>Dear all,<br><br></div>I want to calculate running correlation for All India Rainfall and I am not able to use run_cor because the time dimension does not make sense for me. I wrote a loop but I get the same value for all the period. There is no error in the code output. I am not sure where i am going wrong. It will be great if someone can help<br><br><div><div><br><br>;This is to get a 15 and 20 year running correlation for JJAS<br><br>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"<br>load "$NCARG_ROOT/lib/ncarg/nclscripts/contrib/run_cor.ncl"<br><br><br>ncol = 12<br>nrows=133<br><br>;Reading in the data from 1871 to 2012<br><br>air=readAsciiTable("~/Documents/PhD_June_2015/Data_AIR/air_1871_2012.txt", ncol, "integer",2)<br><br>;Extracting Indian summer monsoon months<br><br>JJAS=air(:,7:10)<br></div></div></div></blockquote><div>1. Why 7:10 is the data starting third column ? Any way not related to your question.</div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div dir="ltr"><div><div><br>JJAS_avg=dim_avg_n(JJAS,1) <br>;print(dimsizes(JJAS_avg))<br><br>years=air(:,1)<br></div></div></div></blockquote><div>2. Do you mean second column of data? </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div dir="ltr"><div><div><br><br>;This is running correlation for 15 years<br><br> yrStrt = 1871 <br> yrLast = 2012<br><br> nmos = 12<br> nyrs = yrLast-yrStrt+1<br> n_15 = nyrs/15 <br><br> running_correlation_15=new(n_15,"float")<br><br> do nyr=0,(15*n_15)-1,15<br> nyrStrt = nyr ; clarity<br> nyrLast = nyrStrt+14<br> running_correlation_15(nyr/15) = escorc(years,JJAS_avg )<br></div></div></div></blockquote><div>3. What is variable "years" ? is it some meaningful index ?</div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div dir="ltr"><div><div> <br> end do<br></div></div></div></blockquote><div>4. Any way this loop is not doing the 15 year window running correlation,</div><div><br></div><div>Suggestions</div><div>1. You need n_15=nyrs-15 to store 15 rear running correlation output.</div><div>2. You have to specify window while passing the data to the escorc function like escorc(index(y1:y15),JJAS_avg(y1:y15))</div><div><span class=""><br></span></div><div>3. Within the loop you can increment y1=y1+1 and y15=y15+1 for running correlation, or you can fix the time step increment as you wish. In the above loop you are using a 15 year time step by using do nyr=0,(15*n_15)-1,15. </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div dir="ltr"><div><div>print(running_correlation_15)<br><br><br></div></div></div>
<br>_______________________________________________<br>
ncl-talk mailing list<br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div></div>