<div dir="ltr"><font face="arial, helvetica, sans-serif">Hi Andreas,</font><div><font face="arial, helvetica, sans-serif">As there are singular values for reg_coef and trend_slope for each start and end time, I think it would be best to simply create a 3D array dimensioned (start_date, end_date,2), with the last dimension (0) = trend_slope and (1) representing reg_<wbr>coeff. If you <i>must</i> have a 4-dimensional array then you should create an array that is dimensioned 22 x 22 x 2 x 2, but (:,:,1,0) and (:,:,0,1) will be missing. In the example below I show how to create both types of arrays:</font></div><div><font face="arial, helvetica, sans-serif"><br></font></div><div><font face="arial, helvetica, sans-serif">finarr = new((/22,22,2/),typeof(value_test))</font></div><div><font face="arial, helvetica, sans-serif">finarr!0 = "startyear"</font></div><div><font face="arial, helvetica, sans-serif">finarr&startyear = ispan(1980,2001,1)</font></div><div><font face="arial, helvetica, sans-serif">finarr!1 = "endyear"</font></div><div><font face="arial, helvetica, sans-serif">finarr&endyear = ispan(1989,2010,1)</font></div><div><div><font face="arial, helvetica, sans-serif">finarrZ = new((/22,22,2,2/),</font><span style="font-family:arial,helvetica,sans-serif">typeof(value_test))</span></div><div><font face="arial, helvetica, sans-serif">finarrZ!0 = "startyear"</font></div><div><font face="arial, helvetica, sans-serif">finarrZ&startyear = ispan(1980,2001,1)</font></div><div><font face="arial, helvetica, sans-serif">finarrZ!1 = "endyear"</font></div><div><font face="arial, helvetica, sans-serif">finarrZ&endyear = ispan(1989,2010,1)</font></div><div><font face="arial, helvetica, sans-serif"><br></font></div><div><font face="arial, helvetica, sans-serif">do i=1980,2001,1</font></div><div><font face="arial, helvetica, sans-serif"> do j=1989,2010,1</font></div><div><font face="arial, helvetica, sans-serif"> ntimes=ispan(i,j,1)</font></div><div><font face="arial, helvetica, sans-serif"> if (j-i .ge. 9 ) then</font></div><div><font face="arial, helvetica, sans-serif"> test_rc=regline_stats(ntimes,<wbr>value_test({i:j}))</font></div><div><font face="arial, helvetica, sans-serif"> print(test_rc)</font></div><div><font face="arial, helvetica, sans-serif"> start_date=i</font></div><div><font face="arial, helvetica, sans-serif"> end_date=j</font></div><div><font face="arial, helvetica, sans-serif"> finarr({i},{j},0) = test_rc@b*10 ; prints slope of the linear trend for the selected range of years</font></div><div><font face="arial, helvetica, sans-serif"> finarr({i},{j},1) = test_rc@b95 ; prints 2.5% and 97.5% regression coefficient confidence intervals </font></div><div><font face="arial, helvetica, sans-serif"> finarrZ(</font><span style="font-family:arial,helvetica,sans-serif">{i},{j},0,0) = test_rc@b*10</span><span style="font-family:arial,helvetica,sans-serif"> </span></div><div><span style="font-family:arial,helvetica,sans-serif"> </span><font face="arial, helvetica, sans-serif"> finarrZ(</font><span style="font-family:arial,helvetica,sans-serif">{i},{j},1,1) = </span><span style="font-family:arial,helvetica,sans-serif">test_rc@b95 </span></div><div><font face="arial, helvetica, sans-serif"> else </font></div><div><font face="arial, helvetica, sans-serif"> print(“this age range has not been selected, since its less than 9 years")</font></div><div><font face="arial, helvetica, sans-serif"> end if </font></div><div><font face="arial, helvetica, sans-serif"> delete(ntimes)</font></div><div><font face="arial, helvetica, sans-serif"> end do</font></div><div><font face="arial, helvetica, sans-serif">end do </font></div></div><div><br></div><div>Hope that helps. If you have any further questions please respond to the ncl-talk email list.</div><div>Adam</div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Jun 5, 2017 at 9:48 AM, Andreas Chrysanthou <span dir="ltr"><<a href="mailto:eeac@leeds.ac.uk" target="_blank">eeac@leeds.ac.uk</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div style="word-wrap:break-word">
<div id="m_5767769911056133050divtagdefaultwrapper" dir="ltr">
<p style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
Hi NCL users,</p>
<p style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
<br>
</p>
<p style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
I created a loop for calculating a linear trend sensitivity for different start and end dates. </p>
<p style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
My timeseries span over 1980-2010 and I want to calculate the linear trend of the timeseries for when the start - end date are more than >= 9years.
<br>
</p>
<p style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
(1980-1989, 1980-1990, 1980-1991,..., 1980-2010</p>
<p style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
1981-1990, 1981-1991, ..., 1981-2010</p>
<p style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
...</p>
<p style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
...</p>
<p style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
...</p>
<p style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
2000-2009, 2001-2010<br>
</p>
<p style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
2001-2010)</p>
<div><br>
</div>
<div><br>
</div>
<div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
I’ve created the loops that calculate those values but I need to save them in an array as 4 different columns (start_date, end_date,trend_slope,reg_<wbr>coeff)</div>
<div style="margin-top:0px;margin-bottom:0px;font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
<span style="font-size:12pt">My aim is to plot those as a a contour box plot with x axis the start date and y axis the end date, and the calculated (i’m gonna filter out those values with some criteria) trend values based on a colormap.</span></div>
</div>
<div style="margin-top:0px;margin-bottom:0px;font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
<span style="font-size:12pt"><br>
</span></div>
<div style="margin-top:0px;margin-bottom:0px"><font face="Calibri, Arial, Helvetica, sans-serif" size="3">Can you help me in order to save the values to that new array so I can plot it?</font></div>
<div style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
<br>
</div>
<div style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
A snippet of the code follows:</div>
<div style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
<br>
</div>
<div>
<div><span style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:medium">do i=1980,2001,1</span></div>
<div><font face="Calibri, Arial, Helvetica, sans-serif" size="3"> </font><span style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:medium">do j=1989,2010,1</span></div>
<div><span style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:medium"> ntimes=ispan(i,j,1)</span></div>
<div><span style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:medium">if (j-i .ge. 9 ) then</span></div>
<div><font face="Calibri, Arial, Helvetica, sans-serif" size="3">
</font><span style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:medium"> test_rc=regline_stats(ntimes,<wbr>value_test({i:j}))</span></div>
<div><font face="Calibri, Arial, Helvetica, sans-serif" size="3"> print(test_rc)</font></div>
<div><font face="Calibri, Arial, Helvetica, sans-serif" size="3"> start_date=i</font></div>
<div><font face="Calibri, Arial, Helvetica, sans-serif" size="3"> end_date=j</font></div>
<div><font face="Calibri, Arial, Helvetica, sans-serif" size="3"> trend_slope=test_rc@b*10 ; prints slope of the linear trend for the selected range of years</font></div>
<div><font face="Calibri, Arial, Helvetica, sans-serif" size="3"> reg_coeff=test_rc@b95 ; prints 2.5% and 97.5% regression coefficient confidence intervals </font></div>
<div><font face="Calibri, Arial, Helvetica, sans-serif" size="3">else print(“this age range has not been selected, since its less than 9 years")</font></div>
<div><font face="Calibri, Arial, Helvetica, sans-serif" size="3"> end if </font></div>
<div><font face="Calibri, Arial, Helvetica, sans-serif" size="3"> delete(ntimes)</font></div>
<div><span style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:medium">end do</span></div>
<div><span style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:medium">end do </span></div>
</div>
<div style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
<br>
</div>
<div style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
Cheers,</div>
<div style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
Andreas</div><span class="HOEnZb"><font color="#888888">
<div style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
<br>
</div>
<p style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
<br>
</p>
<p style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
<br>
</p>
<p style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
<br>
</p>
<p style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
<br>
</p>
<p style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
<br>
</p>
<p style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
<br>
</p>
<p style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
<br>
</p>
<div id="m_5767769911056133050Signature" style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
<div id="m_5767769911056133050divtagdefaultwrapper" dir="ltr" style="font-size:12pt;color:#000000;font-family:Calibri,Arial,Helvetica,sans-serif">
<font face="HelveticaNeue-Light" size="1"><span style="font-size:13px;font-weight:normal;text-transform:none">
<div>
<div style="white-space:normal">
<div style="white-space:normal">
<div style="white-space:normal">
<div style="white-space:normal">
<div>Andreas</div>
<br>
</div>
</div>
</div>
</div>
</div>
</span></font>
<p><br>
</p>
</div>
</div>
</font></span></div>
</div>
<br>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div><div><span><font color="#888888">Adam Phillips <br></font></span></div><span><font color="#888888">Associate Scientist, </font></span><span><font color="#888888">Climate and Global Dynamics Laboratory, NCAR<br></font></span></div></div><div><span><font color="#888888"><a href="http://www.cgd.ucar.edu/staff/asphilli/" target="_blank">www.cgd.ucar.edu/staff/asphilli/</a> </font></span><span><font color="#888888">303-497-1726 </font></span></div><span><font color="#888888"></font></span><div><div><span><font color="#888888"><br></font></span><div><span><font color="#888888"><a href="http://www.cgd.ucar.edu/staff/asphilli" target="_blank"></a></font></span></div></div></div></div></div></div></div></div></div></div></div>
</div>