[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