1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
|
#include "config.h"
#include <stdlib.h>
#include <stdio.h>
#include <string.h>
#include <netcdf.h>
#define FILE_NAME "http://nomads.ncdc.noaa.gov:80/dods/SEAWINDS/clm/uvclm95to05"
#define ERRCODE 2
#define ERR(e) {printf("Error: %s\n", nc_strerror(e)); exit(ERRCODE);}
#undef DEBUG
int
main()
{
int ncid, varid;
int retval,i;
size_t indx;
double timevals[2][12];
double expected[2][2] = {{14.0, 14.0},{45.0, 45.0}};
memset((void*)timevals,0,sizeof(timevals));
if ((retval = nc_open(FILE_NAME, NC_NOWRITE, &ncid)))
ERR(retval);
if ((retval = nc_inq_varid(ncid, "time", &varid)))
ERR(retval);
for(indx=0;indx<2;indx++) {
for(i=0;i<2;i++) {
if ((retval = nc_get_var1_double(ncid, varid, &indx, &timevals[indx][i])))
ERR(retval);
#ifdef DEBUG
printf("expected[%d][%d]=%g timevals[%d][%d]=%g\n",
indx,i,expected[indx][i],indx,i,timevals[indx][i]);
#endif
if(expected[indx][i] != timevals[indx][i]) {
printf("*** FAIL: expected[%d][%d]=%g timevals[%d][%d]=%g\n",
indx,i,expected[indx][i],indx,i,timevals[indx][i]);
exit(1);
}
}
}
if ((retval = nc_close(ncid)))
ERR(retval);
printf("*** PASS\n");
return 0;
}
|
