[ncl-talk] Sub: linin2

[Rick Brownrigg]

Hi Dale,

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:

    lat_0: [18.42513..41.10219]

But for p, its

    lat: [52.3799..73.3307]

Again, I really don't know for certain, just speculating...

Rick

On Mon, Nov 13, 2017 at 11:29 PM, dale zuri <dalezuri at gmail.com> wrote:

> 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/0e7899e8/attachment.html>