Latitude-height cross sections using samplexy

["Gustafson, William I" <william.gustafson@xxxxxxx>]

Dear Ferreters,

I am trying to construct a vertical cross section with data plotted from one file on curvilinear coordinates with data from a rectilinear file overlaid on top. I was able to get this to work for longitude-height cross sections by subsampling a series of points along a line in the curvilinear file and then using samplexy to pull the corresponding points out of the rectilinear file as follows:

let jpt = 15
fill samplexy(data[d=rect],longitude[j=`jpt`,d=curv],latitude[j=`jpt`,d=curv]), longitude[j=`jpt`,d=curv], z[gz=p_rect]
contour/over data[j=`jpt`,d=curv], longitude[j=`jpt`,d=curv], z[gz=p_curv]

However, I cannot get the similar notation to work for a latitude-height cross section:

let ipt = 15
fill samplexy(data[d=rect],longitude[i=`ipt`,d=curv],latitude[i=`ipt`,d=curv]), latitude[i=`ipt`,d=curv], z[gz=p_rect]
contour/over data[i=`ipt`,d=curv], latitude[i=`ipt`,d=curv], z[gz=p_curv]

I get the error:

 **ERROR: inconsistent sizes of data regions: X coord field lies
          in different plane from data to be plotted
Presumably, this is because the "X" points given by latitude[i=`ipt`,d=curv] are not on the X-axis but on the Y-axis. 

So, my question is, how do I get this to work for the latitude-height cross sections? I have not been able to figure out how to transfer data between axes.

Thanks,
       Bill

_____________________________
Dr. William I. Gustafson Jr.
Pacific Northwest National Laboratory
P.O. Box 999, MSIN K9-30
Richland, WA 99352, USA
(509)372-6110
http://www.pnl.gov/atmos_sciences