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
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
|
/*
Copyright 2006, University Corporation for Atmospheric Research. See
the COPYRIGHT file for copying and redistribution conditions.
$Id: fort-lib.c,v 1.15 2009/02/13 15:58:00 ed Exp $
*/
/*
Extracted from of fort-lib.c. Used to supply four required netcdf4
functions used in Fortran2003 interfaces. These need to be external
and I don't want to mangle fort-lib.c just to define these four functions
Version 1.0, April, 2009 - based on netcdf-4.0.1 source
Version 2.0, April, 2010 - based on netcdf-4.1.1 source
Version 3.0, Jan. 2016 - added functions to return number of
groups, types, and compound field
dim_sizes. Needed to support move to
allocatable arrays for starts, strides, etc
modified by: Richard Weed Ph.D.
Center for Advanced Vehicular Systems
Mississippi State University
rweed@cavs.msstate.edu
*/
#include <config.h>
#include <stddef.h> /* for NULL */
#include <errno.h>
#include <string.h>
#include <stdlib.h>
#include <stdio.h>
#include "netcdf.h"
#ifdef USE_NETCDF4
/* These appear to only be defined in netcdf-4*/
/* Get the varids for a fortran function (i.e. add 1 to each
* varid.) */
extern int
nc_inq_varids_f(int ncid, int *nvars, int *fvarids)
{
int *varids, nvars1;
int i, ret = NC_NOERR;
/* Get the information from the C library. */
if ((ret = nc_inq_varids(ncid, &nvars1, NULL)))
return ret;
if (!(varids = malloc(nvars1 * sizeof(int))))
return NC_ENOMEM;
if ((ret = nc_inq_varids(ncid, NULL, varids)))
goto exit;
/* Add one to each, for fortran. */
for (i = 0; i < nvars1; i++)
fvarids[i] = varids[i] + 1;
/* Tell the user how many there are. */
if (nvars)
*nvars = nvars1;
exit:
free(varids);
return ret;
}
/* Get the dimids for a fortran function (i.e. add 1 to each
* dimid.) */
extern int
nc_inq_dimids_f(int ncid, int *ndims, int *fdimids, int parent)
{
int *dimids, ndims1;
int i, ret = NC_NOERR;
/* Get the information from the C library. */
if ((ret = nc_inq_dimids(ncid, &ndims1, NULL, parent)))
return ret;
if (!(dimids = malloc(ndims1 * sizeof(int))))
return NC_ENOMEM;
if ((ret = nc_inq_dimids(ncid, NULL, dimids, parent)))
goto exit;
/* Add one to each, for fortran. */
for (i = 0; i < ndims1; i++)
fdimids[i] = dimids[i] + 1;
/* Tell the user how many there are. */
if (ndims)
*ndims = ndims1;
exit:
free(dimids);
return ret;
}
/* Swap the dim sizes for fortran. */
extern int
nc_insert_array_compound_f(int ncid, int typeid, char *name,
size_t offset, nc_type field_typeid,
int ndims, int *dim_sizesp)
{
int *dim_sizes_f;
int i, ret;
if (ndims <= 0)
return NC_EINVAL;
/* Allocate some storage to hold ids. */
if (!(dim_sizes_f = malloc(ndims * sizeof(int))))
return NC_ENOMEM;
/* Create a backwards list of dimension sizes. */
for (i = 0; i < ndims; i++)
dim_sizes_f[i] = dim_sizesp[ndims - i - 1];
/* Call with backwards list. */
ret = nc_insert_array_compound(ncid, typeid, name, offset, field_typeid,
ndims, dim_sizes_f);
/* Clean up. */
free(dim_sizes_f);
return ret;
}
extern int
nc_inq_compound_field_ndims(int ncid, nc_type xtype, int fieldid, int *ndims)
{
int ret;
/* Find out how many dims. */
if ((ret = nc_inq_compound_field(ncid, xtype, fieldid, NULL, NULL,
NULL, ndims, NULL)))
return ret;
return NC_NOERR;
}
extern int
nc_inq_compound_field_f(int ncid, nc_type xtype, int fieldid, char *name,
size_t *offsetp, nc_type *field_typeidp, int *ndimsp,
int *dim_sizesp)
{
int ndims;
int ret;
/* Find out how many dims. */
if ((ret = nc_inq_compound_field(ncid, xtype, fieldid, NULL, NULL,
NULL, &ndims, NULL)))
return ret;
/* Call the function. */
if ((ret = nc_inq_compound_field(ncid, xtype, fieldid, name, offsetp,
field_typeidp, ndimsp, dim_sizesp)))
return ret;
/* Swap the order of the dimsizes. */
if (ndims)
{
int *f, *b, temp;
for (f = dim_sizesp, b = &dim_sizesp[ndims - 1]; f < b; f++, b--)
{
temp = *f;
*f = *b;
*b = temp;
}
}
return NC_NOERR;
}
extern int
nc_inq_numgrps(int ncid, int *numgrps)
{
int ret;
if ((ret = nc_inq_grps(ncid, numgrps, NULL)))
return ret;
}
extern int
nc_inq_numtypes(int ncid, int *numtypes)
{
int ret;
if ((ret = nc_inq_typeids(ncid, numtypes, NULL)))
return ret;
}
#endif /*USE_NETCDF4*/
/*
add a dummy nc_rename_grp function if it is not supported. This is include
here so we can build/test with netCDF < version 4.3.1 without
*/
#ifndef NC_HAVE_RENAME_GRP
extern int
nc_rename_grp(int ncid, const char *name)
{
printf("\n*** Warning - nc_rename_grp not supported in this netCDF version\n");
printf("*** Update your netCDF C libraries to version 4.3.1 or higher\n");
return NC_ENOGRP;
}
#endif
|
