[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
961212: Trouble read a netcdf file on the cray
- Subject: 961212: Trouble read a netcdf file on the cray
- Date: Thu, 12 Dec 96 14:31:31 -0700
Tom,
>Date: Thu, 12 Dec 1996 11:03:24 -0700 (MST)
>From: address@hidden (Thomas Mayer)
>Organization: National Center for Atmospheric Research NCAR
>To: address@hidden (Steve Emmerson)
>Subject: Re: 961212: Trouble read a netcdf file on the cray
>Keywords: 199612102207.AA21583
In the above message you wrote:
> I understand. Note the code does work on our SUN. I bought it first
> to Jeff Kuehn in the consulting office. When he couldn't find any
> problem in the code, we both stated to wonder if there wasn't a
> problem with the netcdf FORTRAN interface
I don't think the problem lies with the FORTRAN interface. I just built
a modified version of your `netcdf.f' program (I added a printing of
variable `ihgt500' after reading it in and then an immediate call to
exit(). Here's the first page of output:
shavano$ cf77 -c netcdf.f
shavano$ cf77 -o netcdf netcdf.o ../libsrc/libnetcdf.a
shavano$ ./netcdf
netCDF file = z500
NETID,VARID 0 1
MISSING,FILL_VALUE,ADD_OFFSET,SCALE_FACTOR 32766 -32767 32066.00000
1.00000
DIMS(I), DIMNAME(I), DIMSZ(I) 1 lon 144
DIMS(I), DIMNAME(I), DIMSZ(I) 2 lat 73
DIMS(I), DIMNAME(I), DIMSZ(I) 3 level 12
DIMS(I), DIMNAME(I), DIMSZ(I) 4 time 33
I,NMDIMS,START(I),COUNT(I) 1 4 1 144
I,NMDIMS,START(I),COUNT(I) 2 4 1 73
I,NMDIMS,START(I),COUNT(I) 3 4 1 12
I,NMDIMS,START(I),COUNT(I) 4 4 1 33
ihgt500 = 164*-26811, 5*-26810, 3*-26809, 2*-26808, 3*-26807, 2*-26806, 2
*-26805, 2*-26804, 2*-26803, 2*-26802, 2*-26801,
2*-26800, 2*-26799, 2*-26798, 3*-26797, 4*-26796, 8*-26795, 2*-26794, 15
*-26795, 3*-26796, 16*-26795, 5*-26796, 3*-26797,
3*-26798, 2*-26799, 2*-26800, 2*-26801, -26802, 2*-26803, 2*-26804, 2*-2
6805, 2*-26806, 2*-26807, 2*-26808, 3*-26809,
4*-26810, 5*-26811, -26787, -26786, 2*-26785, -26784, 2*-26783, 2*-26782
, 9*-26781, 3*-26782, 10*-26783, 2*-26782,
2*-26781, 3*-26780, 2*-26779, 13*-26778, 3*-26779, 2*-26780, 2*-26781, -
26782, 2*-26783, -26784, -26785, -26787, -26788,
-26790, -26791, -26793, -26794, -26796, -26798, -26800, -26801, -26803,
-26805, -26806, -26807, 2*-26808, 3*-26809,
2*-26808, -26807, -26805, -26804, -26802, -26800, -26799, -26797, -2679
5, -26793, -26790, -26788, -26787, -26785,
-26783, -26781, -26780, -26778, -26777, -26776, -26775, -26774, -26773,
2*-26772, 5*-26771, 2*-26772, 2*-26773, -26774,
-26775, -26776, -26777, -26779, -26780, -26781, -26782, -26783, -26784,
-26785, -26786, -26787, 2*-26788, 7*-26789,
2*-26788, -26787, -26775, -26773, -26770, -26768, -26765, -26762, -2675
9, -26756, -26754, -26751, -26749, -26747,
-26746, -26745, 3*-26744, 2*-26745, -26746, -26747, -26749, -26750, -26
751, 2*-26752, 2*-26753, 4*-26754, 4*-26753,
2*-26754, 2*-26755, -26757, -26758, -26759, -26761, -26762, -26764, -26
766, -26768, -26770, -26771, -26773, -26774,
-26775, -26777, -26778, 2*-26779, -26780, -26781, -26782, -26784, -2678
5, -26787, -26790, -26792, -26796, -26799,
-26804, -26809, -26814, -26820, -26826, -26833, -26839, -26846, -26852,
-26858, -26864, -26868, -26872, -26875,
2*-26877, -26876, -26875, -26872, -26868, -26863, -26857, -26851, -2684
4, -26837, -26830, -26823, -26816, -26809,
-26803, -26797, -26791, -26786, -26781, -26777, -26772, -26768, -26765,
-26761, -26757, -26754, -26751, -26748, -26745
-26743, -26741, -26739, 2*-26737, 2*-26736, -26737, -26738, -26740, -26
742, -26744, -26747, -26750, -26753, -26756,
-26759, -26762, -26765, -26767, -26770, -26772, -26774, -26776, -26777,
As you can see, the values look Ok. It therefore appears that the
problem lies elsewhere in the program.
Please let me know if this helps.
--------
Steve Emmerson <address@hidden>