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_var_generic_val(double *buf, info_field_struct *info_field, char *filename, char *varname, int index) {
00077 int istat;
00078
00079 int ncinid;
00080 int varinid;
00081 nc_type vartype_main;
00082 int varndims;
00083 int vardimids[NC_MAX_VAR_DIMS];
00084
00085 float valf;
00086 char *tmpstr = NULL;
00087 size_t t_len;
00088 size_t *idx;
00089
00090
00091 idx = (size_t *) malloc(sizeof(size_t));
00092 if (idx == NULL) alloc_error(__FILE__, __LINE__);
00093 idx[0] = (size_t) index;
00094
00095
00096 tmpstr = (char *) malloc(MAXPATH * sizeof(char));
00097 if (tmpstr == NULL) alloc_error(__FILE__, __LINE__);
00098
00099
00100
00101
00102 printf("%s: Opening for reading 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
00106 printf("%s: READ %s %s\n", __FILE__, varname, filename);
00107
00108
00109 istat = nc_inq_varid(ncinid, varname, &varinid);
00110 if (istat != NC_NOERR) handle_netcdf_error(istat, __FILE__, __LINE__);
00111
00112
00113 istat = nc_inq_var(ncinid, varinid, (char *) NULL, &vartype_main, &varndims, vardimids, (int *) NULL);
00114 if (istat != NC_NOERR) handle_netcdf_error(istat, __FILE__, __LINE__);
00115
00116
00117 if (info_field != NULL) {
00118
00119 if (vartype_main == NC_FLOAT) {
00120 istat = nc_get_att_float(ncinid, varinid, "missing_value", &valf);
00121 if (istat != NC_NOERR)
00122 info_field->fillvalue = -9999.0;
00123 else
00124 info_field->fillvalue = (double) valf;
00125 }
00126 else if (vartype_main == NC_DOUBLE) {
00127 istat = nc_get_att_double(ncinid, varinid, "missing_value", &(info_field->fillvalue));
00128 if (istat != NC_NOERR)
00129 info_field->fillvalue = -9999.0;
00130 }
00131
00132
00133 istat = nc_inq_attlen(ncinid, varinid, "units", &t_len);
00134 if (istat == NC_NOERR) {
00135 handle_netcdf_error(istat, __FILE__, __LINE__);
00136 istat = nc_get_att_text(ncinid, varinid, "units", tmpstr);
00137 if (istat == NC_NOERR) {
00138 if (tmpstr[t_len-1] != '\0')
00139 tmpstr[t_len] = '\0';
00140 info_field->units = strdup(tmpstr);
00141 }
00142 else
00143 info_field->units = strdup("unknown");
00144 }
00145 else
00146 info_field->units = strdup("unknown");
00147
00148
00149 istat = nc_inq_attlen(ncinid, varinid, "height", &t_len);
00150 if (istat == NC_NOERR) {
00151 handle_netcdf_error(istat, __FILE__, __LINE__);
00152 istat = nc_get_att_text(ncinid, varinid, "height", tmpstr);
00153 if (istat == NC_NOERR) {
00154 if (tmpstr[t_len-1] != '\0')
00155 tmpstr[t_len] = '\0';
00156 info_field->height = strdup(tmpstr);
00157 }
00158 else
00159 info_field->height = strdup("unknown");
00160 }
00161 else
00162 info_field->height = strdup("unknown");
00163
00164
00165 istat = nc_inq_attlen(ncinid, varinid, "long_name", &t_len);
00166 if (istat == NC_NOERR) {
00167 handle_netcdf_error(istat, __FILE__, __LINE__);
00168 istat = nc_get_att_text(ncinid, varinid, "long_name", tmpstr);
00169 if (istat == NC_NOERR) {
00170 if (tmpstr[t_len-1] != '\0')
00171 tmpstr[t_len] = '\0';
00172 info_field->long_name = strdup(tmpstr);
00173 }
00174 else
00175 info_field->long_name = strdup(varname);
00176 }
00177 else
00178 info_field->long_name = strdup(varname);
00179 }
00180
00181
00182 istat = nc_get_var1_double(ncinid, varinid, idx, buf);
00183 if (istat != NC_NOERR) handle_netcdf_error(istat, __FILE__, __LINE__);
00184
00185
00186 istat = ncclose(ncinid);
00187 if (istat != NC_NOERR) handle_netcdf_error(istat, __FILE__, __LINE__);
00188
00189
00190 (void) free(tmpstr);
00191 (void) free(idx);
00192
00193
00194 return 0;
00195 }