00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038
00039
00040
00041
00042
00043
00044
00045
00046
00047
00048
00049
00050
00051
00052
00053
00054
00055
00056
00057
00058
00059
00060
00061
00062 #include <io.h>
00063
00065 int
00066 read_netcdf_xy(double **x, double **y, int *nx, int *ny, char *xname, char *yname, char *dimxname, char *dimyname, char *filename) {
00081 int istat;
00082
00083 size_t dimval;
00084
00085 int ncinid;
00086 int yinid;
00087 int xinid;
00088 nc_type vartype;
00089 int varndims;
00090 int vardimids[NC_MAX_VAR_DIMS];
00091 int xdiminid;
00092 int ydiminid;
00093
00094 size_t start[3];
00095 size_t count[3];
00096
00097 int ndims;
00098
00099
00100
00101
00102 printf("%s: Reading info from NetCDF input file %s\n", __FILE__, filename);
00103 istat = nc_open(filename, NC_NOWRITE, &ncinid);
00104 if (istat != NC_NOERR) handle_netcdf_error(istat, __FILE__, __LINE__);
00105
00107 ndims = 1;
00108
00109
00110 istat = nc_inq_dimid(ncinid, dimyname, &ydiminid);
00111 if (istat != NC_NOERR) handle_netcdf_error(istat, __FILE__, __LINE__);
00112 istat = nc_inq_dimlen(ncinid, ydiminid, &dimval);
00113 if (istat != NC_NOERR) handle_netcdf_error(istat, __FILE__, __LINE__);
00114 *ny = (int) dimval;
00115
00116 istat = nc_inq_dimid(ncinid, dimxname, &xdiminid);
00117 if (istat != NC_NOERR) handle_netcdf_error(istat, __FILE__, __LINE__);
00118 istat = nc_inq_dimlen(ncinid, xdiminid, &dimval);
00119 if (istat != NC_NOERR) handle_netcdf_error(istat, __FILE__, __LINE__);
00120 *nx = (int) dimval;
00121
00122
00123 istat = nc_inq_varid(ncinid, yname, &yinid);
00124 if (istat != NC_NOERR) handle_netcdf_error(istat, __FILE__, __LINE__);
00125
00126 istat = nc_inq_var(ncinid, yinid, (char *) NULL, &vartype, &varndims, vardimids, (int *) NULL);
00127 if (istat != NC_NOERR) handle_netcdf_error(istat, __FILE__, __LINE__);
00128 if (varndims != ndims) {
00129 (void) fprintf(stderr, "Error NetCDF type and/or dimensions %d != %d.\n", varndims, ndims);
00130 istat = ncclose(ncinid);
00131 if (istat != NC_NOERR) handle_netcdf_error(istat, __FILE__, __LINE__);
00132 return -1;
00133 }
00134
00135
00136 istat = nc_inq_varid(ncinid, xname, &xinid);
00137 if (istat != NC_NOERR) handle_netcdf_error(istat, __FILE__, __LINE__);
00138
00139 istat = nc_inq_var(ncinid, xinid, (char *) NULL, &vartype, &varndims, vardimids, (int *) NULL);
00140 if (istat != NC_NOERR) handle_netcdf_error(istat, __FILE__, __LINE__);
00141 if (varndims != ndims) {
00142 (void) fprintf(stderr, "Error NetCDF type and/or dimensions %d != %d.\n", varndims, ndims);
00143 istat = ncclose(ncinid);
00144 if (istat != NC_NOERR) handle_netcdf_error(istat, __FILE__, __LINE__);
00145 return -1;
00146 }
00147
00149
00150 start[0] = 0;
00151 start[1] = 0;
00152 start[2] = 0;
00153 count[0] = (size_t) (*ny);
00154 count[1] = 0;
00155 count[2] = 0;
00156 (*y) = (double *) malloc((*ny) * sizeof(double));
00157 if ((*y) == NULL) alloc_error(__FILE__, __LINE__);
00158
00159
00160 istat = nc_get_vara_double(ncinid, yinid, start, count, *y);
00161 if (istat != NC_NOERR) handle_netcdf_error(istat, __FILE__, __LINE__);
00162
00163
00164 start[0] = 0;
00165 start[1] = 0;
00166 start[2] = 0;
00167 count[0] = (size_t) (*nx);
00168 count[1] = 0;
00169 count[2] = 0;
00170 (*x) = (double *) malloc((*nx) * sizeof(double));
00171 if ((*x) == NULL) alloc_error(__FILE__, __LINE__);
00172
00173
00174 istat = nc_get_vara_double(ncinid, xinid, start, count, *x);
00175 if (istat != NC_NOERR) handle_netcdf_error(istat, __FILE__, __LINE__);
00176
00177 istat = ncclose(ncinid);
00178 if (istat != NC_NOERR) handle_netcdf_error(istat, __FILE__, __LINE__);
00179
00180
00181 return 0;
00182 }