Re: int2p_n_Wrap interpolation encountering missing value

From: Dennis Shea <shea_at_nyahnyahspammersnyahnyah>
Date: Sun Oct 27 2013 - 09:00:04 MDT

[1] Not enough information.
     Looking at the arrays, this looks like atmospheric data.
     If so, presumably the missing data are over (say) the Himalayasdata
would have to be extrapolated,
     Rocky Mountains, etc. In this case if your iso level is (say) 925,
     the

[2] Please do not send the same message to ncl-talk multiple times.

On 10/27/13 8:49 AM, Kai Zhang wrote:
> Dear All,
>
> I am using int2p_n_Wrap to interpolate pressure levels to iso levels
> (potential temperature). I have the following code:
>
> x = int2p_n_Wrap(theta,xin,iso_lev,0,1) --------------------
> ***************
> ; theta(time,plev,lat,lon) is potential temperature
> ; xin(time,plev,lat,lon) is the variable I want to interpolate to iso levels
> ;iso_lev is a one-dimensional levels containing iso levels I want to
> interpolate into.
>
> Here I have a problem, since both theta(time,plev,lat,lon) and xin
> (time,plev,lat,lon) have missing values, the output
> x(time,iso_lev,lat,lon) would be missing values.
>
> I know the example of one dimension is easy to interpolation by using
> ind(.not.ismissing) to miss out the missing values, but how to deal with
> four dimension dataset?? If I use three do loops, it will cost several days
> to runout the results. (ntime=612, nlat=96, nlon=144)
>
> So do you know if there is a way to miss out the missing values in the
> datasets??
>
> Thanks very much for your help! This is really important in my research!
> Without interpolation, I cannot see the results I need.
>
> Thanks again!!!
>
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Sun Oct 27 09:00:09 2013

This archive was generated by hypermail 2.1.8 : Fri Nov 01 2013 - 08:58:14 MDT