<div dir="ltr"><div>Hi Rick,</div><div><br></div><div>begin<br>;;;;;;;;;;;;;;;;;ERA;;;;;;;;;;;;;;;;;;;;;;;;;;;<br>; iz = addfile("<a href="http://psa_mask_cfsv2.nc">psa_mask_cfsv2.nc</a>","r")<br> iz = addfile("<a href="http://tmp2m.2012050218.time.nc">tmp2m.2012050218.time.nc</a>","r")<br> print(iz)<br> printVarSummary(iz)<br><br>P1=iz->TMP_P0_L103_GGA0(:,15:40,180:255)<br> printVarSummary(P1)<br> printMinMax(P1,0)<br> ds = addfile("<a href="http://nlwrs.sfc.gauss.1948.nc">nlwrs.sfc.gauss.1948.nc</a>","r")<br> pp=ds->nlwrs(:,:,:)<br> p = pp(:,{76:51},{168:239})<br><br> printVarSummary(p)<br>;****************************************************<br> psa_gauss = linint2_Wrap(P1&lon_0,P1&lat_0,P1, False, p&lon,p&lat, 0)<br><br> printVarSummary(psa_gauss)<br> printMinMax(psa_gauss,0)<br>return<br><br></div><div><br></div><div><br></div><div>I changed the line. Then also the same issue of getting a missing number as output. <br></div><div>P1=iz->TMP_P0_L103_GGA0(:,15:40,180:255)</div><div> p = pp(:,{76:51},{168:239})</div><div><br></div><div>Variable: P1<br>Type: float<br>Total Size: 9555936 bytes<br> 2388984 values<br>Number of Dimensions: 3<br>Dimensions and sizes: [forecast_time0 | 1209] x [lat_0 | 26] x [lon_0 | 76]<br>Coordinates: <br> forecast_time0: [6..7254]<br> lat_0: [75.11737..51.49582]<br> lon_0: [168.7498..239.0622]</div><div><br></div><div>Variable: p<br>Type: float<br>Total Size: 723216 bytes<br> 180804 values<br>Number of Dimensions: 3<br>Dimensions and sizes: [time | 366] x [lat | 13] x [lon | 38]<br>Coordinates: <br> time: [1297320..1306080]<br> lat: [75.2351..52.3799]<br> lon: [168.75..238.125]<br><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Nov 13, 2017 at 8:31 PM, Rick Brownrigg <span dir="ltr"><<a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@ucar.edu</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">Hi Dale,<div><br></div><div>It would really help if you could show the results of those various printVarSummary()'s And include one for pp/p -- are you certain your coordinate subscripting is correct in this line:</div><span class=""><div><br></div><div><span style="font-size:12.8px"> p = pp(:,{51:75},{168:239})</span><br></div><div><span style="font-size:12.8px"><br></span></div></span><div><span style="font-size:12.8px">Rick</span></div><div><span style="font-size:12.8px"><br></span></div></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div class="h5">On Mon, Nov 13, 2017 at 9:33 PM, dale zuri <span dir="ltr"><<a href="mailto:dalezuri@gmail.com" target="_blank">dalezuri@gmail.com</a>></span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5"><div dir="ltr"><div><div><div>Hi all,<br></div>Please help me understand this issue.<br></div>I didn't get any error. But the interpolated output are missing numbers.<br> iz = addfile("<a href="http://tmp2m.2012050218.time.nc" target="_blank">tmp2m.2012050218.time<wbr>.nc</a>","r")<br> print(iz)<br> printVarSummary(iz)<br><br>P1=iz->TMP_P0_L103_GGA0(:,75:5<wbr>1,168:239)<br> printVarSummary(P1)<br> printMinMax(P1,0)<br> ds = addfile("<a href="http://nlwrs.sfc.gauss.1948.nc" target="_blank">nlwrs.sfc.gauss.1948.<wbr>nc</a>","r")<br> pp=ds->nlwrs(:,::-1,:)<br> p = pp(:,{51:75},{168:239})<br><br> printVarSummary(p)<br>;*****************************<wbr>***********************<br> psa_gauss = linint2_Wrap(P1&lon_0,P1&lat_0<wbr>,P1, False, p&lon,p&lat, 0)<br> printVarSummary(psa_gauss)<br> printMinMax(psa_gauss,0)<br><br><br>;;;;<br>(0) Temperature: min=1e+20 max=1e+20<br><br></div>Thanks <br><div><div><br><br></div></div></div>
<br></div></div>______________________________<wbr>_________________<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/mailma<wbr>n/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>
</blockquote></div><br></div>