<div dir="ltr">Hi Edwin,<div>I put a few print/printVarSummary statements into your script:</div><div> ind_store_end:=ind_store_start+92*tres-1<br> print(ind_extract_start+" "+ind_extract_end)<br> printVarSummary(var)</div><div> var_lag(ind_store_start:ind_store_end)=var(y,ind_extract_start:ind_extract_end)<br></div><div><br></div><div>Output:</div><div>(0) -20 347 <br><br>Variable: var<br>Type: float<br>Total Size: 182400 bytes<br> 45600 values<br>Number of Dimensions: 2<br>Dimensions and sizes: [evn | 8] x [time | 5700]<br>Coordinates: <br> evn: [1..8]<br> time: [396..14050]<br>Number Of Attributes: 4<br> long_name : EOF: Amplitude: Climate Hazards group InfraRed Precipitation with Stations<br> _FillValue : -9999<br> matrix : covariance<br> ts_mean : ( 1.154463e-10, 5.018545e-10, -3.159135e-10, 5.93204e-10, -1.381243e-09, -3.113667e-10, -1.764017e-10, -6.856012e-10 )<br></div><div><br></div><div>So your ind_extract_start value is -20, which is outside the valid index range of (0-5699) as is shown in the output of printVarSummary. My guess as to the problem is that you are asking for data early on in your record and you do not have 48 timesteps before your chosen day to calculate your lead/lag correlation. Outside of that recommendation I am afraid I cannot offer any further advice.</div><div>Good luck,<br>Adam</div><div><br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Mar 3, 2020 at 7:45 AM Edwin Marco via ncl-talk <<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><p class="MsoNormal" style="margin:0cm 0cm 8pt;line-height:107%;font-size:11pt;font-family:Calibri,sans-serif">Dear ncl users,</p><p class="MsoNormal" style="margin:0cm 0cm 8pt;line-height:107%;font-size:11pt;font-family:Calibri,sans-serif"><span style="font-size:11pt">I was wondering if could help me to find out the source of error I got in </span><span style="font-size:11pt">plotting correlation of principal components (PCs)</span><span style="font-size:11pt">. I have used filtered precipitation datasets to calculate
and plot eofs and generate associate PCs and write them into netcdf file (EOF_my.nc).</span></p>
<p class="MsoNormal" style="margin:0cm 0cm 8pt;line-height:107%;font-size:11pt;font-family:Calibri,sans-serif">Then using the written netcdf file (EOF_my.nc), I am getting
error in plotting correlation of the PCs.</p>
<p class="MsoNormal" style="margin:0cm 0cm 8pt;line-height:107%;font-size:11pt;font-family:Calibri,sans-serif"><b><span style="font-family:"Adobe Devanagari",serif;background:yellow">fatal: Subscript out of range, error in subscript #1</span></b></p>
<p class="MsoNormal" style="margin:0cm 0cm 8pt;line-height:107%;font-size:11pt;font-family:Calibri,sans-serif"><b><span style="font-family:"Adobe Devanagari",serif;background:yellow">fatal:["Execute.c":8637]:Execute: Error occurred at or near
line 98 in file plot_correlation.ncl</span></b><b><span style="font-family:"Adobe Devanagari",serif"></span></b></p>
<p class="MsoNormal" style="margin:0cm 0cm 8pt;line-height:107%;font-size:11pt;font-family:Calibri,sans-serif">I have done a thorough investigation of the script and suspect
line 75 to be a source of error of line 98, however, I have not been able to solve
the problem. </p>
<p class="MsoNormal" style="margin:0cm 0cm 8pt;line-height:107%;font-size:11pt;font-family:Calibri,sans-serif"><b>Line 75: </b></p>
<p class="MsoNormal" style="margin:0cm 0cm 8pt;line-height:107%;font-size:11pt;font-family:Calibri,sans-serif"><b><span style="font-family:"Adobe Devanagari",serif;background:silver">ind_mar_1_1981=ind(cd_cal(:,0) .eq. yrs .and. cd_cal(:,1) .eq. mns
.and. cd_cal(:,2) .eq. 1 .and. cd_cal(:,3) .eq. 0)</span></b><b><span style="font-family:"Adobe Devanagari",serif"></span></b></p>
<p class="MsoNormal" style="margin:0cm 0cm 8pt;line-height:107%;font-size:11pt;font-family:Calibri,sans-serif"> <span style="font-size:11pt">Please I need your help!!</span></p>
<p class="MsoNormal" style="margin:0cm 0cm 8pt;line-height:107%;font-size:11pt;font-family:Calibri,sans-serif">I attach the script and the <a href="http://eof_my.nc" target="_blank">eof_my.nc</a> datasets.</p><p class="MsoNormal" style="margin:0cm 0cm 8pt;line-height:107%;font-size:11pt;font-family:Calibri,sans-serif">Thanks</p><p class="MsoNormal" style="margin:0cm 0cm 8pt;line-height:107%;font-size:11pt;font-family:Calibri,sans-serif">Edwin</p></div>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">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/mailman/listinfo/ncl-talk</a></blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="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>