<div dir="ltr"><div><div>Thank you will keep that in mind<br></div>Best Regards<br></div>Ipshita<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Oct 30, 2014 at 11:45 AM, gibies george <span dir="ltr">&lt;<a href="mailto:gibiesgeorge@gmail.com" target="_blank">gibiesgeorge@gmail.com</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><p dir="ltr">Question 1 and 2 meant to remind you that NCL use FORTRAN indexing. 0, 1, 2, 3, ......  </p><div class="HOEnZb"><div class="h5">
<div class="gmail_quote">On Oct 30, 2014 8:38 PM, &quot;Ipshita Majhi&quot; &lt;<a href="mailto:ipmajhi@alaska.edu" target="_blank">ipmajhi@alaska.edu</a>&gt; wrote:<br type="attribution"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div>Thank you . Yes! I realised that my script was not doing the running correlation instead it was doing correlation for every 15 year time step.Answer to the questions are as follows:-<br><br><div>1. Why <span><span>7:10</span></span> is the data starting third column ? Any way not related to your question.<br></div><div><br></div><div>I am calculating average for summer months of monsoon June, July, August and September and that is the column 7:10<br></div><br>JJAS_avg=dim_avg_n(JJAS,1) <br>;print(dimsizes(JJAS_avg))<br><br>years=air(:,1)<br><br></div>2. Years is the first column from air which has information for year going from 1871 to 2012<br><br></div>I am incorporating the suggestion recommended.<br><br></div>Thank You<br>Best Regards<br>Ipshita<br><div><div><div><br><br></div></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Oct 30, 2014 at 2:48 AM, gibies george <span dir="ltr">&lt;<a href="mailto:gibies@tropmet.res.in" target="_blank">gibies@tropmet.res.in</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><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"><span>On 30 October 2014 02:43, Ipshita Majhi <span dir="ltr">&lt;<a href="mailto:ipmajhi@alaska.edu" target="_blank">ipmajhi@alaska.edu</a>&gt;</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 &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl&quot;<br>load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/contrib/run_cor.ncl&quot;<br><br><br>ncol = 12<br>nrows=133<br><br>;Reading in the data from 1871 to 2012<br><br>air=readAsciiTable(&quot;~/Documents/PhD_June_2015/Data_AIR/air_1871_2012.txt&quot;, ncol, &quot;integer&quot;,2)<br><br>;Extracting Indian summer monsoon months<br><br>JJAS=air(:,7:10)<br></div></div></div></blockquote></span><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><span><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,&quot;float&quot;)<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></span><div>3. What is variable &quot;years&quot; ? 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><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>
</blockquote></div><br></div>
</blockquote></div>
</div></div></blockquote></div><br></div>