[Opengrads-devel] errors on sdfopen with hdf

[Jennifer A. <jm...@co...>]

On Feb 2, 2008, at 3:21 PM, Arlindo da Silva wrote:
>
> 4) Of the 70 tests, gradsv2 is failing 16. They are related to 2  
> issues:
>
> a) The "sdfopen model.hdf" is not quite working: it reports the  
> file as having 1 timestep instead of 5. I verified that the RH  
> binary Jennifer posted has exactly the same problem, so it does not  
> appear to be my build.

Arlindo,
I have put some print statments in gasdf.c to see why this is happening.
The HDF routine SDdiminfo returns the following information on  
model.hdf:

dim_id,name,dimsize,dtype,ndatts = 327680 time 0 6 1
dim_id,name,dimsize,dtype,ndatts = 327681 longitude 72 6 2
dim_id,name,dimsize,dtype,ndatts = 327682 latitude 46 6 2
dim_id,name,dimsize,dtype,ndatts = 327683 levels 7 6 5

Looks like the file has a time dimension of size 0. I remember now  
that when I discovered this during development, I put in the  
following line:
	if (dimsize==0) dimsize=1; /* would you have a coordinate axis with  
size==0 */

The interesting thing is that when I run ncdump (from the HDF4.2r2  
library build) it gives me this:

dimensions:
         time = UNLIMITED ; // (5 currently)
         longitude = 72 ;
         latitude = 46 ;
         levels = 7 ;

--Jennifer