<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_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">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">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">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>