wrf_cape_2d

Hi,
I am attempting to plot the four stability parameters output by the function wrf_cape_2d.
Looking at the output, the LCL and LFC parameters all look reasonable but the CAPE and CIN parameters are all constant fields of -0.1 so I am hoping to get some guidance from the experts. I have repeated this problem with both version 6.0 and 5.2.
I've attached a .txt file showing a single WRF-model grid point's incoming arguments and the 4 output parameters for that specific point to show how the data looks coming in and the problematic CAPE/CIN fields coming out.
It seems all the units and ordering are fine for each argument (Pressure: Pa, Temp: K, Q: kg/kg, Z: m, Zsfc: m, and Psfc: hPa).I also do not show any missing values that could influence the output for the single point that I am reporting (see attached ncl.txt file for the actual values). The point falls in central Arizona and is for January 1st.

The following is my code snippet that I've cleaned up a bit:
     NCEPf = addfile (NCEP+NCEPc(nf), "r") NCEPP = wrf_user_getvar(NCEPf, "pres", 0) NCEPT = wrf_user_getvar(NCEPf, "tk", 0) NCEPQ = NCEPf->QVAPOR(0,:,:,:) NCEPZ = wrf_user_getvar(NCEPf, "geopt", 0) ;comes in as m2/s2 NCEPZ = NCEPZ / 9.81 ;reduce to units of m NCEPsfc= NCEPf->PSFC(0,:,:) / 100. ; comes in as Pa, reduce to units of hPa NCEPCape = wrf_cape_2d( NCEPP, NCEPT, NCEPQ, NCEPZ, HGT, NCEPsfc, True)

Thanks for any advice.
Jon

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk