<div dir="ltr"><div dir="ltr"><div dir="ltr">Hi Lyndz,<div>1) You can usually safely ignore these error messages. I believe they usually occur due to timeseries (or individual gridpoints) having standard deviations of 0 or missing.</div><div>2) At line 50 you have this:</div><div> prob = ttest(xAve,xVar,xN, yAve,yVar,yN, iflag, False)</div><div>and then at line 58 you have this:</div><div> prob = prob(0,:)<br></div><div><br></div><div>The easiest solution is to rename prob to say prob2 at line 58, and then change prob to prob2 at lines 61 and 81. Alternatively, you can also overwrite the previous iteration of prob like this:</div><div>prob := prob(0,:)</div><div><br></div><div>Renaming the array is the better solution from a safe-coding perspective.</div><div><br></div><div>Hope that helps. If you have any further questions let ncl-talk know.</div><div>Adam</div><div><br></div><div><br></div></div></div></div><br><div class="gmail_quote"><div dir="ltr">On Fri, Dec 28, 2018 at 4:41 AM Lyndz <<a href="mailto:olagueralyndonmark429@gmail.com">olagueralyndonmark429@gmail.com</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"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr" class="gmail-m_-8301752397231960052gmail_signature"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><b>Dear NCL-experts,</b></div><div dir="ltr"><br></div><div dir="ltr">I want to perform a <i>t</i>-test similar to example 4 in the NCL webpage: ( <a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/ttest.shtml" target="_blank">https://www.ncl.ucar.edu/Document/Functions/Built-in/ttest.shtml</a>), but in this case, </div><div>I will be using a different equivalent sample size (degrees of freedom) following Davis (1976) (below image explains the details and also attached in this email):</div><div><br></div><div><div><img src="cid:ii_jq7ueetq0" alt="Screenshot from 2018-12-28 18-33-00.png" width="356" height="217" style="margin-right: 0px;"><br></div></div><div><br></div><div><div><i><b>Sources:</b></i></div><div><div><i>Davis, R. E., 1976: Predictability of sea surface temperature and sea level pressure anomalies over the North Pacific Ocean. J. Phys. Oceanogr., 6, 249–266.</i></div></div></div><div><br></div><div>Below are the links to my data sets. I also attached here the script. I want to save the output as a netcdf file with the test statistics and pvalues.</div><div><br></div><div><b>Link to File 1:</b> <a href="https://www.dropbox.com/s/cmjpmvcw9i4qxo0/rhum_P45P50_1979-1993_filt_850hPa.nc?dl=0" target="_blank">https://www.dropbox.com/s/cmjpmvcw9i4qxo0/rhum_P45P50_1979-1993_filt_850hPa.nc?dl=0</a></div><div><b>Link to File 2:</b> <a href="https://www.dropbox.com/s/o8dtkt2yznlgvql/rhum_P45P50_1994-2008_filt_850hPa.nc?dl=0" target="_blank">https://www.dropbox.com/s/o8dtkt2yznlgvql/rhum_P45P50_1994-2008_filt_850hPa.nc?dl=0</a></div><div><div><br></div></div><div><b><font color="#ff0000">PROBLEMS</font></b>:</div><div>1. I got the following warning message: </div><div><div> <i><font color="#ff0000"> warning:ttest: encountered 10512 cases where s1 and/or s2 were less than 2. Output set to missing values in these cases</font></i></div></div><div>I am not sure if this should be ignored or this indicates a mistake in my code. </div><div><br></div><div>2. I also got the following error message after running the script: </div><div><div><i> <font color="#ff0000">fatal:Number of dimensions on right hand side do not match number of dimension in left hand side</font></i></div><div><i><font color="#ff0000"> fatal:["Execute.c":8640]:Execute: Error occurred at or near line 58 in file ttest_filtered.ncl</font></i></div></div><div>Any idea how to solve this?</div><div><br></div><div>I would like to ask for your suggestions on how to do this correctly in NCL.</div><div>I'll appreciate any help.</div><div><br></div><div>Sincerely,</div><div><br></div><div><b>Lyndz</b></div><div><b><br></b></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></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><br>
</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>