Convert Hybrid Sigma to Pressure Coords and Output new Vars to NetCDF

Hello NCL Talk,

I am attempting to convert a climate model's hybrid sigma coordinates to
pressure coordinates for several variables: zonal wind, meridional wind,
specific humidity, streamfunction, and velocity potential. I have done this
in the code below using the function, *vinth2p*. Now though, I want to
output to NetCDF these variables and their associated coordinate variables
on the new p levels. I am having trouble outputting the variables because
of the new pressure level array (which I call *lev_p*). I get a
"fatal:Execute: Error occurred at or near line 131 in file." The original
hybrid sigma levels were defined by *lev, *so I tried commenting this out
but the same error results. Does anyone know how to fix this? My code is
attached. Thanks.

-Zac

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