Re: Puzzling **ERROR: inconsistent sizes of data regions

[hengkek@xxxxxxxxxx]

Thanks to Ansley and Vinu (who provided a similar solution).

	let var_diff=ev[d=1,i=1:133]-ev[d=2,i=1:133]

The above line did take care of it.
Based on the 'show grid' cmd as appended below, I conclude that the
problem was that the 1st dataset's axes were considered irregular, while
the 2nd, regular.

I.e. the 2nd dataset had the "point_spacing=even" netcdf attribute for
the axes while the 1st didn't.  After I inserted the said attribute in
the 2nd dataset, all was well again.  I presume the problem will also be
solved if I removed the said attribute from the 1st dataset.

Thanks for all the tips.

regards
-chk

-----------------------------------
yes? show grid ev[d=1]
    GRID GPR4
 name       axis              # pts   start                end
 LON       LONGITUDE          133mr   99E                  121E
 LAT       LATITUDE           199 r   9S                   24N
 normal    Z
 TIME      TIME                21 r   01-MAY-2001 00:00    06-MAY-2001
00:00
yes? show grid ev[d=2]
    GRID GKU4
 name       axis              # pts   start                end
 LON1      LONGITUDE          133mi   99E                  121E
 LAT1      LATITUDE           199 i   9S                   24N
 normal    Z
 TIME      TIME                21 r   01-MAY-2001 00:00    06-MAY-2001
00:00
-----------------------------------


-------- Original Message --------
Subject: Re: Puzzling **ERROR: inconsistent sizes of data regions
Date: Wed, 30 Jul 2003 08:16:57 -0700
From: Ansley Manke <Ansley.B.Manke@noaa.gov>
To: hengkek@nus.edu.sg
CC: ferret_users@ferret.wrc.noaa.gov

Hello,
This command might take care of it,

let var_diff=ev[d=1,i=1:133]-ev[d=2,i=1:133]

but you may still want to investigate why there seems to
be a discrepancy. Here are a few more commands that may
help you find out what is happening.  If you do determine the
cause, let us know.

  yes? show grid aaa
  yes? show grid bbb

The above commands will give you the names of the axes
in question.  Then use SHOW AXIS, naming the X axis of
each variable, or with /X or /I to list some of the coordinate
values.

  yes? show axis axisname

Are there any differences in the listings for your two
X axes?  The SHOW AXIS section in the Users Guide
lists the meaning of the output of this comand.  For instance,
is one of the axes a modulo axis, and the other not?  You
can make a more detailed comparison with the command

yes? show axis/x axisname   ! lists the X coordinates

You might also try this command at the start of your script

  yes? set mode diagnostic

There is an FAQ which discusses the output of this mode,
and it may help you see what is happening with your
variables.  The FAQ can be found at
http://www.ferret.noaa.gov/Ferret/FAQ/other/mode_diagnostic.html

Ansley Manke

hengkek@nus.edu.sg wrote:

> Greetings
>
> I'm getting a puzzling error:
>
>  *** NOTE: Ambiguous coordinates on X axis: EV[D=1]-EV[D=2]
>  **ERROR: inconsistent sizes of data regions: X axis
>           EV[D=2] has 134 points (I=000:133)
>           expression has 133 points (I=001:133)
>
> Both my datasets have 133 points and are identical in terms of grid and
> axes, but I'm puzzled as to how it becomes 134 in the 2nd data set as
> reported by the error above.
>
> Any clues how to correct this?
>
> My paraphrased session is appended.
> Thanks for any answers.
>
> -chk
> Physical Oceanography
> Singapore
>
> -------- Start my session ------------------
>  ! NOAA/PMEL TMAP
>  ! FERRET v5.51
>  ! Solaris 5.8 - 02/28/03
>  ! 30-Jul-03 16:25
> yes? use "aaa.nc"
> yes? use "bbb.nc"
> yes? let var_diff=ev[d=1]-ev[d=2]
>  !-> DEFINE VARIABLE var_diff=ev[d=1]-ev[d=2]
> yes? fill/l=2 var_diff
>  *** NOTE: Ambiguous coordinates on X axis: EV[D=1]-EV[D=2]
>  **ERROR: inconsistent sizes of data regions: X axis
>           EV[D=2] has 134 points (I=000:133)
>           expression has 133 points (I=001:133)
> yes? show dat
>      currently SET data sets:
>     1> aaa.nc
>  name     title                       I         J         K        L
> ...[deleted]...
>  EV      surface elevation          1:133     1:199     ...       1:21
>
>     2> bbb.nc  (default)
>  name     title                       I         J         K        L
> ...[deleted]...
>  EV      surface elevation          1:133     1:199     ...      1:21
> -------- end my session ------------------