[ncl-talk] Sub: linin2

[dale zuri]

Hi,
Thanks.
Variable: P1
Type: float
Total Size: 8704800 bytes
            2176200 values
Number of Dimensions: 3
Dimensions and sizes:   [forecast_time0 | 1209] x [lat_0 | 25] x [lon_0 |
72]
Coordinates:
            forecast_time0: [6..7254]
            lat_0: [18.42513..41.10219]
            lon_0: [157.4998..224.0622]
Number Of Attributes: 11
  initial_time :        05/02/2012 (18:00)
  level :        2
  level_type :  Specified height level above ground (m)
  parameter_template_discipline_category_number :       ( 0, 0, 0, 0 )
  parameter_discipline_and_category :   Meteorological products, Temperature
  grid_type :   Gaussian latitude/longitude
  _FillValue :  1e+20
  units :       K
  long_name :   Temperature
  production_status :   Operational products
  center :      US National Weather Service - NCEP (WMC)
(0)     Temperature: min=273.771   max=302.267

Variable: p
Type: float
Total Size: 667584 bytes
            166896 values
Number of Dimensions: 3
Dimensions and sizes:   [time | 366] x [lat | 12] x [lon | 38]
Coordinates:
            time: [1297320..1306080]
            lat: [52.3799..73.3307]
            lon: [168.75..238.125]
Number Of Attributes: 15
  _FillValue :  -9.96921e+36
  valid_range : ( -100, 400 )
  actual_range :        ( -69.19995, 237.2 )
  missing_value :       -9.96921e+36
  parent_stat : Individual Obs
  statistic :   Mean
  level_desc :  Surface
  dataset :     NCEP Reanalysis Daily Averages
  var_desc :    Net Longwave Radiation Flux
  GRIB_name :   NLWRS
  GRIB_id :     112
  least_significant_digit :     0
  precision :   1
  units :       W/m^2
  long_name :   mean Daily Net Longwave Radiation Flux at Surface

Variable: psa_gauss
Type: float
Total Size: 2205216 bytes
            551304 values
Number of Dimensions: 3
Dimensions and sizes:   [forecast_time0 | 1209] x [lat | 12] x [lon | 38]
Coordinates:
            forecast_time0: [6..7254]
            lat: [52.3799..73.3307]
            lon: [168.75..238.125]
Number Of Attributes: 11
  _FillValue :  1e+20
  initial_time :        05/02/2012 (18:00)
  level :        2
  level_type :  Specified height level above ground (m)
  parameter_template_discipline_category_number :       ( 0, 0, 0, 0 )
  parameter_discipline_and_category :   Meteorological products, Temperature
  grid_type :   Gaussian latitude/longitude
  units :       K
  long_name :   Temperature
  production_status :   Operational products
  center :      US National Weather Service - NCEP (WMC)
(0)     Temperature: min=1e+20   max=1e+20

Thanks

On Mon, Nov 13, 2017 at 8:31 PM, Rick Brownrigg <brownrig at ucar.edu> wrote:

> Hi Dale,
>
> 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:
>
>  p = pp(:,{51:75},{168:239})
>
> Rick
>
>
> On Mon, Nov 13, 2017 at 9:33 PM, dale zuri <dalezuri at gmail.com> wrote:
>
>> Hi all,
>> Please help me understand this issue.
>> I didn't get any error. But the interpolated output are missing numbers.
>>  iz = addfile("tmp2m.2012050218.time.nc","r")
>>   print(iz)
>>   printVarSummary(iz)
>>
>> P1=iz->TMP_P0_L103_GGA0(:,75:51,168:239)
>>   printVarSummary(P1)
>>   printMinMax(P1,0)
>>   ds = addfile("nlwrs.sfc.gauss.1948.nc","r")
>>   pp=ds->nlwrs(:,::-1,:)
>>   p = pp(:,{51:75},{168:239})
>>
>>   printVarSummary(p)
>> ;****************************************************
>>  psa_gauss = linint2_Wrap(P1&lon_0,P1&lat_0,P1, False, p&lon,p&lat, 0)
>>   printVarSummary(psa_gauss)
>>   printMinMax(psa_gauss,0)
>>
>>
>> ;;;;
>> (0)    Temperature: min=1e+20   max=1e+20
>>
>> Thanks
>>
>>
>>
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20171113/88ed6c5e/attachment.html>