[ncl-talk] Sub: linin2
Rick Brownrigg
brownrig at ucar.edu
Mon Nov 13 22:31:36 MST 2017
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/e9028c54/attachment.html>
More information about the ncl-talk
mailing list