<div dir="ltr"><div><div><div><div>Hi Dale,<br><br></div>I'm not really sure, and I'm not expert at this by any stretch, but it looks like the two variables, "P1" and "p", don't spatially overlap. For P1, the latitude range is:<br><br> lat_0: [18.42513..41.10219]<br><br></div>But for p, its<br><br> lat: [52.3799..73.3307]<br><br></div>Again, I really don't know for certain, just speculating...<br><br></div>Rick<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Nov 13, 2017 at 11:29 PM, dale zuri <span dir="ltr"><<a href="mailto:dalezuri@gmail.com" target="_blank">dalezuri@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>Hi,<br></div>Thanks.<br>Variable: P1<br>Type: float<br>Total Size: 8704800 bytes<br> 2176200 values<br>Number of Dimensions: 3<br>Dimensions and sizes: [forecast_time0 | 1209] x [lat_0 | 25] x [lon_0 | 72]<br>Coordinates: <br> forecast_time0: [6..7254]<br> lat_0: [18.42513..41.10219]<br> lon_0: [157.4998..224.0622]<br>Number Of Attributes: 11<br> initial_time : 05/02/2012 (18:00)<br> level : 2<br> level_type : Specified height level above ground (m)<br> parameter_template_discipline_<wbr>category_number : ( 0, 0, 0, 0 )<br> parameter_discipline_and_<wbr>category : Meteorological products, Temperature<br> grid_type : Gaussian latitude/longitude<br> _FillValue : 1e+20<br> units : K<br> long_name : Temperature<br> production_status : Operational products<br> center : US National Weather Service - NCEP (WMC)<br>(0) Temperature: min=273.771 max=302.267<br><br>Variable: p<br>Type: float<br>Total Size: 667584 bytes<br> 166896 values<br>Number of Dimensions: 3<br>Dimensions and sizes: [time | 366] x [lat | 12] x [lon | 38]<br>Coordinates: <br> time: [1297320..1306080]<br> lat: [52.3799..73.3307]<br> lon: [168.75..238.125]<br>Number Of Attributes: 15<br> _FillValue : -9.96921e+36<br> valid_range : ( -100, 400 )<br> actual_range : ( -69.19995, 237.2 )<br> missing_value : -9.96921e+36<br> parent_stat : Individual Obs<br> statistic : Mean<br> level_desc : Surface<br> dataset : NCEP Reanalysis Daily Averages<br> var_desc : Net Longwave Radiation Flux<br> GRIB_name : NLWRS<br> GRIB_id : 112<br> least_significant_digit : 0<br> precision : 1<br> units : W/m^2<br> long_name : mean Daily Net Longwave Radiation Flux at Surface<br><br>Variable: psa_gauss<br>Type: float<br>Total Size: 2205216 bytes<br> 551304 values<br>Number of Dimensions: 3<br>Dimensions and sizes: [forecast_time0 | 1209] x [lat | 12] x [lon | 38]<br>Coordinates: <br> forecast_time0: [6..7254]<br> lat: [52.3799..73.3307]<br> lon: [168.75..238.125]<br>Number Of Attributes: 11<br> _FillValue : 1e+20<br> initial_time : 05/02/2012 (18:00)<br> level : 2<br> level_type : Specified height level above ground (m)<br> parameter_template_discipline_<wbr>category_number : ( 0, 0, 0, 0 )<br> parameter_discipline_and_<wbr>category : Meteorological products, Temperature<br> grid_type : Gaussian latitude/longitude<br> units : K<br> long_name : Temperature<br> production_status : Operational products<br> center : US National Weather Service - NCEP (WMC)<br>(0) Temperature: min=1e+20 max=1e+20<br><br></div>Thanks<br></div><div class="HOEnZb"><div class="h5"><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><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="m_-2956002676927791929h5">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="m_-2956002676927791929h5"><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>
</div></div></blockquote></div><br></div>