I have used the function of Samplexy to Subsample gridded fields onto profile coordinates as follow
---------------------------------------------
yes? let my_lon = lon1[d=my_grid.dat]
yes? let my_lat = lat1[d=my_grid.dat]
yes? LET x_var = samplexy(temp ,my_lon,my_lat)
---------------------------------------------------------------
and I can display the vertical section well along the axis's number.
When I try to use latitude as x-axis as follow
-------------------------------------------------------
let order = SORTI(my_lat)
let section = SAMPLEI(x_var, order)
DEFINE AXIS/Y/NAME=yax_sect/FROM_DATA/UNITS=degrees SAMPLEI(my_lat, order)
LET ysection = RESHAPE(section,Y[gy=yax_sect]+Z)
contour ysection
---------------------------------
**ERROR: unknown defining grid: GY=YAX_SECT --> g=user-or-pseudo-var not allowed ________________________________________
That is the problem.
Thanks
Xin
Xin Jin Inst.of Geophysics and Planetary Physics University of California Los Angeles 5845 Slichter Hall Los Angeles, CA 90095-4996, USA phone: (310) 825 4526 fax : (310) 206 3051 email: xjin@igpp.ucla.edu