[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