[netcdfgroup] bug mixing dap reads and netcdf4 writes

• To: support-netcdf@xxxxxxxxxxxxxxxx
• Subject: [netcdfgroup] bug mixing dap reads and netcdf4 writes
• From: Jeff Whitaker <jeffrey.s.whitaker@xxxxxxxx>
• Date: Wed, 13 Mar 2013 09:19:12 -0600

Here's a python program that triggers the error:

 import netCDF4
 
nc=netCDF4.Dataset('http://www.smast.umassd.edu:8080/thredds/dodsC/fvcom/hindcasts/wave_gom3'
  <http://www.smast.umassd.edu:8080/thredds/dodsC/fvcom/hindcasts/wave_gom3%27>)
 nco=netCDF4.Dataset('dummy.nc', 'w', format='NETCDF4')
 
nc2=netCDF4.Dataset('http://www.smast.umassd.edu:8080/thredds/dodsC/fvcom/hindcasts/30yr_gom3'
  <http://www.smast.umassd.edu:8080/thredds/dodsC/fvcom/hindcasts/30yr_gom3%27>)

Traceback (most recent call last):
  File "issue170.py", line 4, in <module>

    raise RuntimeError((<char *>nc_strerror(ierr)).decode('ascii'))
RuntimeError: NetCDF: Invalid argument

and here's the C counterpart

#include <stdlib.h>
#include <netcdf.h>
#define handle_error(e) {printf("Error: %s\n", nc_strerror(e)); exit(2);}

int main(int *argc, char **argv)
{

    char flnmo[] = "dummy.nc";
    int status;                        /* error status */
    int ncid, ncid2, ncido;            /* netCDF ID */

    status = nc_set_default_format(NC_NETCDF4, NULL);
    if (status != NC_NOERR) handle_error(status);

    status = nc_open(flnm, NC_NOWRITE, &ncid);
    if (status != NC_NOERR) handle_error(status);

    status = nc_open(flnmo, NC_WRITE | NC_CLOBBER, &ncid);
    if (status != NC_NOERR) handle_error(status);

    status = nc_open(flnm2, NC_NOWRITE, &ncid);
    if (status != NC_NOERR) handle_error(status);

    return ( 0 );
}

which produces

Error: NetCDF: Invalid argument

Regards, Jeff

• Follow-Ups:
  • Re: [netcdfgroup] bug mixing dap reads and netcdf4 writes
    • From: Russ Rew