[Thread Prev][Thread Next][Index]

SAMPLEXY funtion agian



 Dear Ansley,

 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=20:25,L=1],xlon,ylat)
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



Ansley Manke <ansley.b.manke@noaa.gov> wrote:
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


Do you Yahoo!?
vote.yahoo.com - Register online to vote today!
[Thread Prev][Thread Next][Index]

Dept of Commerce / NOAA / OAR / PMEL / TMAP

Contact Us | Privacy Policy | Disclaimer | Accessibility Statement