Interpolation of WRF output to isentropic surfaces / "cape_2d" calc crashing using wrf_user_getvar on standard WRF output

Hi all,

I have two very quick questions.

First, what is the most efficient way to vertically interpolate WRF
output data to isentropic surfaces? Is it to first interpolate the
wrfout data on model levels to pressure levels using wrf_interp_3d_z,
and then interp to isentropic surfaces using int2p_n?

Secondly, regarding CAPE. I have a standard wrfout file, and am trying
to extract CAPE using wrf_user_getvar(file,"cape_2d",-1). The routine
crashes, telling me:

capecalc3d: kmax got too big. kmax= 152

I know this comes from the rip_cape.f fortran routine, and refers to
the vertical levels of my output file somehow (I have output on 150
vertical levels). But I have no idea what to do with this error, and
more to the point, why the NCL function is crashing with a basic wrf
file?

Cheers,

Adam Morgan
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Nov 19 16:00:42 2012