Re: ncl 5.1 vs hdf

All,

We are investigating this offline and will post back when we
have more information.

--Mary

On Mon, 9 Mar 2009, Sijikumar S wrote:

> Hello,
>
> After upgrading to version 5.1 my script which reads HDF files is
> crashing with following error.
>
>
> warning:Invalid stride: stride must be positive non-zero integer
> fatal:HDF: An error occurred while attempting to read variable (Longitude)
> from file
> (/home/dbs/AIRS/2003/AIRS.2003.01.01.L3.RetStd031.v5.0.14.0.G07294211513.hdf)
> fatal:Coordinate variables must have one dimension only
> fatal:No coordinate variable exists for dimension (lon) in variable (lon)
> Segmentation fault
>
>
> It was working fine with NCL 5.0
>
>
> I think the problem is with 2 dimensional Lat-Lon array reading, when
> Longitude is reading as it is, it works fine
>
>
> lon1=f[0]->Longitude
> print(lon1)
>
> Variable: lon1
> Type: float
> Total Size: 259200 bytes
> 64800 values
> Number of Dimensions: 2
> Dimensions and sizes: [YDim_location | 180] x [XDim_location | 360]
> Coordinates:
> Number Of Attributes: 2
> _FillValue : -9999
> hdf_name : Longitude
> (0,0) -179.5
> (0,1) -178.5
> (0,2) -177.5
> (0,3) -176.5
> (0,4) -175.5
> (0,5) -174.5
> (0,6) -173.5
>
>
> But when it reads as follows
> ncl 38> lon=f[0]->Longitude(0,:) gives
>
> warning:Invalid stride: stride must be positive non-zero integer
> fatal:Subscript out of range, error in subscript #1
> fatal:Memory allocation failure:[errno=12]
> fatal:Execute: Error occurred at or near line 37
>
> ncl 39> print(lon)
>
> Variable: lon
> Type: float
> Total Size: 720 bytes
> 180 values
> Number of Dimensions: 2
> Dimensions and sizes: [YDim_location | 180] x [YDim_location | 1]
> Coordinates:
> Number Of Attributes: 2
> _FillValue : -9999
> hdf_name : Longitude
> (0,0) -2.660818e-05
> (1,0) -2.660818e-05
> (2,0) 4.362793e-32
> (3,0) 4.362793e-32
> (4,0) 4.879321e-42
> (5,0) 4.362316e-32
> (6,0) 2.574143e-39
>
>
> NCL is from ncl_ncarg-5.1.0.Linux_i686_gcc412.tar.gz and uses
> ncl_static as ncl.
>
>
> Regards
> Siji
>
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Mar 09 2009 - 21:43:27 MDT