Re: interpolate WRF output to isentropic levels

It is the exact same approach as outlined at
      http://www.ncl.ucar.edu/Applications/isent.shtml,%e2%80%8d

Even easier because wrf_user_getvar returns theta.

load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"

  a = addfile("...","r")
  time = -1
  theta= wrf_user_getvar(a,"theta",time)
  printVarSummary(theta)

  x = wrf_user_getvar(a,"tc" ,time)
  printVarSummary(x)
  printMinMax(x, 0)
                                                         ; isentropic
level(s)
  lvl = 320. ; (/300., 320., 340.=
/)

  lvl@description = "isentropic level"
  lvl@units = theta@units

  xlvl = int2p_n_Wrap (theta, x, lvl, 0, 1)
  printVarSummary(xlvl)
  print("xlvl: min="+min(xlvl)+" max="+max(xlvl))
-----
Please send WRF related questions to wrfhelp@ucar.edu
You can cc ncl-talk@ucar.edu




On Sun, Jul 6, 2014 at 8:35 PM, Li Qi <liqi123sh@qq.com> wrote:

> Dear all,
>
> I'd like to interpolate a variable in WRF output to a specified isentropi=
c
> level.
>
> GrADS has a function (i.e., isen(field,tgrid,pgrid,tlev)‍) to com=
plete
> the task, like
> * 1) Display vertical velocity field on 320K surface:
> *
> * "d "isen(w,t,PP,320)‍
>
> Is there any corresponding function in NCL?
>
> I've checked the examples and got this one:
> http://www.ncl.ucar.edu/Applications/isent.shtml,%e2%80%8d
>
> but I have no idea about the hybrid coef in WRF.
>
> Best,
>
> Li Qi
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>

Received on Mon Jul 07 02:00:42 2014