<div dir="ltr"><div><div 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"><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">https://drive.google.com/file/d/1BCKb7UdM50ckVqIvhnSbO9tkyMK0DmiF/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">https://drive.google.com/file/d/1YtYd4hTIUC3tupaAm4oc69xhilkP9qcj/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_anom_1979-1993.nc","r")</div><div>fa = a->xBPF</div><div><br></div><div>b  = addfile("P45P50_bpass_mean_anom_1994-2008.nc","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">   </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">      </span>Lanczos Bpassfilt</div><div>  units :<span style="white-space:pre">   </span>W m-2</div><div>  _FillValue :<span style="white-space:pre">  </span>32766</div><div>  missing_value :<span style="white-space:pre">       </span>32766</div><div>  wgt_runave_op_ncl :<span style="white-space:pre">   </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]: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>