SAMPLEXY funtion works fine for the examples you mentioned but it did not work with my data. I used following commands which generate the following error:
yes? LET xlon=82+I[I=1:6]
yes? LET dely = 24./49
yes? LET ylat = 24.5 - dely*i[i=1:6] + dely
yes? PLOT/VS/LINE/SYM=27 xlon,ylat ! gives slanted line
yes? LET slanttk = samplexy(tk[x=80:90,y=
yes? FILL slanttk
Bailing out of external function "samplexy":
cannot handle MODULO axis that is not regular
**ERROR: : error in external function
I wonder what caused the error.
thanking you in advance
jagadish
ps : my dataset
yes? show data
currently SET data sets:
1> ./ATMJJAS-DAILYAVG.CDF (default)
name title I J K L
U Zonal Wind 1:63 1:66 1:18 1:122
V Meridional
Wind 1:63 1:66 1:18 1:122
TK Temperature 1:63 1:66 1:18 1:122
QD Mixing Ratio 1:63 1:66 1:18 1:122
QC Cloud Mixing Ratio 1:63 &n!
bsp;
1:66 1:18 1:122
MSE Moist Static Energy 1:63 1:66 1:18 1:122
RH Relative Humidity 1:63 1:66 1:18 1:122
HGT Geopotential Height 1:63 1:66 1:18 1:122
PS Surface
Pressure 1:63 1:66 ... 1:122
RT Total Precip 1:63 1:66 ... 1:122
TGRND Ground Temperature 1:63 1:66 ... 1:122
SMT Total Soil Water 1:63 1:66  !
;
... 1:122
RB Base Flow 1:63 1:66 ... 1:122
Hello Jagadish,
There is an example which does this kind of thing in the
Users Guide. Look up "SAMPLEXY, function definition".
The X axis of the plot that is made has just a count of points,
1 for the first point, 2 for the second point and so on. Here
is the example from the Users Guide, which shows how to define
a set of points along a line. You could also define them as lists as I
did in the previous email. This example has depth vs location; time
vs location would be analogous.
USE ocean_atlas_annual
LET xlon = 234.5 + I[I=1:50] ! define the slant line
LET dely = 24./49
LET ylat = 24.5 - dely*i[i=1:50] + dely
LET slantsalt = samplexy(salt[x=200:300,y=0:30],xlon,ylat)
FILL/LEVELS=(33.2,35.2,0.1)/VLIMITS=0:4000 slantsalt
Ansley