<div dir="ltr">Hi Lyndz,<div>latGauWgt will only work with gaussian latitudes (and is therefore returning 72 weights instead of 73), and you have a regular 2.5 degree spaced latitude grid. There are other ways to account for the shrinking size of grid boxes as you approach the poles. You could weight by the cosine of the latitude:</div><div><div>pi = 4.*atan(1.0)</div><div>rad = (pi/180.)</div><div>coslat = cos(xtmp&lat*rad) ; input array = x</div></div><div><br></div><div><div>xN = wgt_areaave (xEqv, coslat, 1., 0) </div><div>yN = wgt_areaave (yEqv, coslat, 1., 0)</div></div><div><br></div><div>See example one on the wgt_areaave page here for a good summary of weighting options:</div><div><a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/wgt_areaave.shtml">https://www.ncl.ucar.edu/Document/Functions/Built-in/wgt_areaave.shtml</a><br></div><div><br></div><div>Hope that helps!</div><div>Adam</div><div> </div></div><div class="gmail_extra"><br><div class="gmail_quote">On Sun, Apr 29, 2018 at 2:23 AM, Lyndz <span dir="ltr"><<a href="mailto:olagueralyndonmark429@gmail.com" target="_blank">olagueralyndonmark429@gmail.com</a>></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><div class="m_3131801359802042935gmail_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"><div>Dear NCL-experts,</div><div><br></div><div>I am trying to perform an F-test. Here are the files:</div><div><br></div><div><b>Link to File 1: </b></div><div><a href="https://drive.google.com/file/d/1BCKb7UdM50ckVqIvhnSbO9tkyMK0DmiF/view?usp=sharing" target="_blank">https://drive.google.com/file/<wbr>d/<wbr>1BCKb7UdM50ckVqIvhnSbO9tkyMK0D<wbr>miF/view?usp=sharing</a><br></div><div><br></div><div><b>Link to File 2:</b></div><div><a href="https://drive.google.com/file/d/1YtYd4hTIUC3tupaAm4oc69xhilkP9qcj/view?usp=sharing" target="_blank">https://drive.google.com/file/<wbr>d/<wbr>1YtYd4hTIUC3tupaAm4oc69xhilkP9<wbr>qcj/view?usp=sharing</a><br></div><div><br></div><div>Here's my script:</div><div><br></div><div><div>a = addfile("P45P50_bpass_mean_<wbr><a href="http://anom_1979-1993.nc">anom_1979-1993.nc</a>","r")</div><div>fa = a->xBPF</div><div><br></div><div>b = addfile("P45P50_bpass_mean_<wbr><a href="http://anom_1994-2008.nc">anom_1994-2008.nc</a>","r")</div><div>fb = b->xBPF</div><div><br></div><div>xtmp = fa(lat|:,lon|:,time|:)</div><div>ttmp = fb(lat|:,lon|:,time|:)</div></div><div><br></div><div><div>xVar = dim_variance (xtmp)</div><div>yVar = dim_variance (ytmp)</div><div><br></div><div>sigr = 0.05</div><div>xEqv = equiv_sample_size (xtmp, sigr,0)</div><div>yEqv = equiv_sample_size (ytmp, sigr,0)</div><div><br></div><div>nlat = 73</div><div>gwt = latGauWgt(nlat, "lat", "gaussian weights", "") ; 72 points</div></div><div><br></div><div><div>xN = wgt_areaave (xEqv, gwt, 1., 0) </div><div>yN = wgt_areaave (yEqv, gwt, 1., 0)</div><div><br></div><div>sigf = 0.10</div><div>prob = ftest(xVar,xN, yVar,yN, False)</div></div><div><br></div><div><br></div><div>Here's the <b>printVarSummary</b> of fa:</div><div><br></div><div><div>Variable: fa</div><div>Type: float</div><div>Total Size: 630720 bytes</div><div> 157680 values</div><div>Number of Dimensions: 3</div><div>Dimensions and sizes:<span style="white-space:pre-wrap"> </span>[time | 15] x [lat | 73] x [lon | 144]</div><div>Coordinates: </div><div> time: [1574700..1697436]</div><div> lat: [90..-90]</div><div> lon: [ 0..357.5]</div><div>Number Of Attributes: 5</div><div> long_name :<span style="white-space:pre-wrap"> </span>Lanczos Bpassfilt</div><div> units :<span style="white-space:pre-wrap"> </span>W m-2</div><div> _FillValue :<span style="white-space:pre-wrap"> </span>32766</div><div> missing_value :<span style="white-space:pre-wrap"> </span>32766</div><div> wgt_runave_op_ncl :<span style="white-space:pre-wrap"> </span>wgt_runave_n</div></div><div><br></div><div><br></div><div>After the <b>gwt</b> I encounter the following errors:</div><div><br></div><div><div>fatal:wgt_areaave: wgty must be a scalar or a 1-dimensional vector the same size as the second-to-the-last dimension of x</div><div>fatal:["Execute.c":8638]:<wbr>Execute: Error occurred at or near line 19</div></div><div><br></div><div><br></div><div>Any suggestion on how to do this correctly in NCL?</div><div><br></div><div>I'll appreciate any help.</div><div><br></div><div>Sincerely,</div><div><br></div><div>Lyndz</div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div></div></div></div></div></div></div></div></div></div></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>