Re: Latitude-height cross sections using samplexy

[Ansley Manke <Ansley.B.Manke@xxxxxxxx>]

Hi William,
I have not worked through this, but it looks to me as if you can simply
apply the XSEQUENCE function to put your latitude points onto an
abstract X axis.

Ansley Manke

"Gustafson, William I" wrote:

> 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

--
Ansley Manke  Pacific Marine Environmental Laboratory  Seattle WA