[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

961212: Trouble read a netcdf file on the cray

This archive contains answers to questions sent to Unidata support through mid-2025. Note that the archive is no longer being updated. We provide the archive for reference; many of the answers presented here remain technically correct, even if somewhat outdated. For the most up-to-date information on the use of NSF Unidata software and data services, please consult the Software Documentation first.

  • 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>