[ncl-talk] Sub: linin2

Rick Brownrigg brownrig at ucar.edu
Tue Nov 14 07:36:06 MST 2017


Hi Dale,

Hard to say from a distance. I take you've made some basic plots of your
extracted variables so that you know what your data looks like and should
be?

Rick

On Tue, Nov 14, 2017 at 12:33 AM, dale zuri <dalezuri at gmail.com> wrote:

> Hi Rick,
>
> begin
> ;;;;;;;;;;;;;;;;;ERA;;;;;;;;;;;;;;;;;;;;;;;;;;;
> ;  iz = addfile("psa_mask_cfsv2.nc","r")
>   iz = addfile("tmp2m.2012050218.time.nc","r")
>   print(iz)
>   printVarSummary(iz)
>
> P1=iz->TMP_P0_L103_GGA0(:,15:40,180:255)
>   printVarSummary(P1)
>   printMinMax(P1,0)
>   ds = addfile("nlwrs.sfc.gauss.1948.nc","r")
>   pp=ds->nlwrs(:,:,:)
>   p = pp(:,{76:51},{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)
> return
>
>
>
> I changed the line. Then also the same issue of getting a missing number
> as output.
> P1=iz->TMP_P0_L103_GGA0(:,15:40,180:255)
>  p = pp(:,{76:51},{168:239})
>
> Variable: P1
> Type: float
> Total Size: 9555936 bytes
>             2388984 values
> Number of Dimensions: 3
> Dimensions and sizes:    [forecast_time0 | 1209] x [lat_0 | 26] x [lon_0 |
> 76]
> Coordinates:
>             forecast_time0: [6..7254]
>             lat_0: [75.11737..51.49582]
>             lon_0: [168.7498..239.0622]
>
> Variable: p
> Type: float
> Total Size: 723216 bytes
>             180804 values
> Number of Dimensions: 3
> Dimensions and sizes:    [time | 366] x [lat | 13] x [lon | 38]
> Coordinates:
>             time: [1297320..1306080]
>             lat: [75.2351..52.3799]
>             lon: [168.75..238.125]
>
>
> 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/20171114/ea27dc65/attachment.html>


More information about the ncl-talk mailing list