[ncl-talk] wrf_user_intrp3d problem

[Bill Ladwig]

The wrf_user_intrp3d does not do extrapolation below ground, so you're
seeing the fill values there. You need to use the wrf_user_vert_interp
function for this (
https://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_user_vert_interp.shtml).
Note however that if you use this function with pressure as your field
variable and height as your vertical coordinate, the values in the
extrapolated region will be in the wrong units of hPa (currently working on
a fix for this). In your case, where pressure is your vertical coordinate,
you should be fine.

Hope this helps,

Bill

On Tue, Oct 9, 2018 at 3:35 PM Daniel Vecellio <djvecellio at tamu.edu> wrote:

> All,
>
> I'm attempting to extrapolate heights down to 1000mb over my domain in
> order to calculate 1000-500mb thickness. However, at some of my higher
> elevations where pressure only gets to 800mb or so at the surface, the
> extrapolation in the wrf_user_intrp3d function becomes wonky, and I get
> values that are 9.9e+36 for 1000mb height.
>
> The line of code I use is: z_1000 =
> wrf_user_intrp3d(z,p,"h",1000.,0.,False) where z and p are height and
> pressure obtain from wrf_user_getvar
>
> Is there a way to provide reasonable numbers here? Given that we have
> perfectly fine operational thickness maps over the Rockies, it has to be
> possible, right?
>
> Thanks in advance,
> Dan
> _______________________________________________
> 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/20181009/1edbc6a2/attachment.html>