<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_category_number : ( 0, 0, 0, 0 )<br> parameter_discipline_and_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_category_number : ( 0, 0, 0, 0 )<br> parameter_discipline_and_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="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>