[netcdfgroup] 1995-era code

• To: "netcdfgroup@xxxxxxxxxxxxxxxx" <netcdfgroup@xxxxxxxxxxxxxxxx>
• Subject: [netcdfgroup] 1995-era code
• From: "Cronin, Catherine Kemble (LARC-B702)" <catherine.k.cronin@xxxxxxxx>
• Date: Fri, 10 Jul 2009 15:12:12 -0500

I have a C program that was written in the 1995 time frame that reads a netCDF 
file.  It uses "ncopen" and "ncinquire" calls, which I assume means it was from 
version 2.  I did not write this program and do not have any experience with 
netCDF files directly.  My SA installed the netCDF software (I believe it's the 
latest version), and the HDF4 software and I have been able to successfully 
compile the C program.  However, the ncopen call fails upon execution.

After a brief search, I found that I could check the libnetcdf.a file to see if 
the installation was backward compatible with the old version.  Which I did, 
and here are the results from that:

nm -A libnetcdf.a | grep ncopen

libnetcdf.a:v2i.o:00000052 R $XBK2SpSrz1VKG1G.ncopen.__func__
libnetcdf.a:v2i.o:00000f40 T ncopen
libnetcdf.a:fort-v2compat.o:         U ncopen

So, it appears that the library does contain the old functions.  Also, I can 
successfully run an ncdump on the file, so I don't believe the netCDF file is 
corrupted.

Does anyone have any other ideas on what I could check or why I may not be able 
to open the file using ncopen?  Any help would be greatly appreciated.

Cathy

/*----------------------------------------------------------------------------*/
#include <stdio.h>
#include <string.h>
#include <malloc.h>
#include <time.h>
/*----------------------------------------------------------------------------*/

#include "/usr/local/src/HDF4.2r3/hdf4/include/netcdf.h"

#define MAX_CHAR_STR_LEN 80

static char attname[MAX_CHAR_STR_LEN];
static char format_value[MAX_CHAR_STR_LEN];
static char output_format[MAX_CHAR_STR_LEN];

/*
   Find PrISM variables in HLFC netCDF file
*/
int netcdf2PrISM(netCDFfile, fp_out)
char *netCDFfile;
FILE *fp_out;
{

   int  ncid;                   /* netCDF id */
   int ndims;                  /* number of dimensions */
   int nvars;                  /* number of variables */
   int ngatts;    /* number of global attributes */
   int xdimid;                 /* id of unlimited dimension */

   char varname[MAX_NC_NAME];
   nc_type vartype;
   int varndims;
   int vardims[MAX_VAR_DIMS];
   int varnatts;
   long int start[1], count[1];
   void *value;

   int varid;
   int i;
   nc_type atttype;
   int attlen;
   void *attvalue;
   int default_format;


   char dimname[MAX_NC_NAME];
   long dimsize;

   strcpy (attname, "prism");
   /* customize netCDF error processing */
   /*dpr ncopts = NC_VERBOSE; dpr*/
   ncopts = 0;

   /* open netCDF file for reading */
   if ((ncid = ncopen (netCDFfile, NC_NOWRITE)) == -1)
    {
      fprintf (stderr,"HLFC netCDF file could not be opened\n");
      return(-1);
    }

   /* get information about netCDF file */
   if (ncinquire(ncid, &ndims, &nvars, &ngatts, &xdimid) == -1)
    {
      fprintf (stderr,"HLFC netCDF file ncinquire read error\n");
      return(-1);
    }

   /* loop thru the number of variables in netCDF file */
   for (varid=0; varid<nvars; varid++)
   {
      if (ncvarinq(ncid,varid,varname,&vartype,&varndims,vardims,&varnatts)
          == -1)
      {
         fprintf (stderr,"HLFC netCDF file ncvarinq read error\n");
         return(-1);
      }

.......................................................... It goes on further, 
but hopefully this is enough to get the idea....................

• Follow-Ups:
  • Re: [netcdfgroup] 1995-era code
    • From: Roy Mendelssohn