Plotting WRF 500mb temperatures

Hello again. Now I am successfully plotting the 500mb temperature of
the FNL data thanks to your help. Next I would like to plot WRF
temperatures at the 500mb level. I know this requires a little
interpolation because the vertical levels in WRF are not on pressure
levels. I also found the function I think I want to use,
wrf_interp_3d_z. Even further, the example with the documentation
seems to be exactly what I'm trying to do.

However, it isn't working for me. I copied and pasted the following
example directly without any luck. When plotting var2d, the output is
very strange. I can best describe it as some sort of rainbow prism
looking effect. Is there something I am missing? I see the example
data is from 2000, have the fields in the WRF output changed since
then? Thanks!

James

  ; Temperature
    t = a->T(time,:,:,:) ; perturbation potential temperature
(theta-300)
    theta = t + 300. ; potential temperature
    p = a->P(time,:,:,:)
    pb = a->PB(time,:,:,:)
    pf = p + pb ; full pressure
    tk = wrf_tk (pf, theta) ; temperature in Kelvin

  ; Height in the model
    ph = a->PH(time,:,:,:) ; perturbation geopotential
    phb = a->PHB(time,:,:,:) ; base-state geopotential
    tmp = (ph+phb)/9.81 ; On full (w) levels
    dimv = dimsizes(tmp)
    z = 0.5*(tmp(0:dimv(0)-2,:,:)+tmp(1:dimv(0)-1,:,:)) ; On half
(mass) levels

  ; Interpolate tk to 500hPa
    var2d = wrf_interp_3d_z(tk,z,500.)
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Jul 25 2007 - 17:42:18 MDT