segmentation fault: gasdf.c:3237, index out of bound

[Ryo Furue]

Hi Jennifer,

Thank you for your help.

| The seg fault happens when you try to display the 3-dimensional
| variable tr with the Z dimension set to something other than 1. That
| is a bug in the sdfopen code in GrADS 1.9. That code is no longer in
| version 2.0, and I am unlikely to spend any time debugging it since
| the problem is readily fixed if you use a 'dtype netcdf' descriptor
| file:

Um, I'm a bit confused.  I guess you say two things:

1. The bug is fixed in 2.0 and you don't want to fix it in the 1.9
  series.

2. The bug can be circumvented in GrADS 1.9 by using a 'dtype netcdf'
  descriptor file.

Correct?

I certainly understand that you don't want to fix a bug
in an old version that is already fixed in a newer one.
But, creating a descriptor file for a netCDF file is like
adding English subtitles to an English-speaking movie!
(They are useful for people with hearing disability,
 but that's not the point here . . .)  Anyway, I can live
with the workaround.  I don't want to push this case.

|     So, unless GrADS bothers to transpose the temperature array
|
| GrADS bothers to do this when you have 'options zrev' in your
| descriptor file or 'positive down' as a netcdf attribute.

Aha!  Thanks for the explanation.  That's clear.  Probably GrADS
doesn't actually transpose arrays on memory (which would be some
extra computation, which isn't really necessary), but it presents
such a "view" to the user.  My confusion was because I'm not
familiar with that view.  I usually use Ferret, where array indices
never change whether the axis is positive down or not.
The attribute "positive down" affects only how Ferret plots
the field.

Thank you again,
Ryo