[219] | 1 | #include "nc4_data_input.hpp" |
---|
| 2 | |
---|
[599] | 3 | #include "context.hpp" |
---|
| 4 | #include "context_server.hpp" |
---|
[775] | 5 | #include "context_client.hpp" |
---|
| 6 | #include "domain.hpp" |
---|
| 7 | #include "axis.hpp" |
---|
[967] | 8 | #include "scalar.hpp" |
---|
[599] | 9 | |
---|
[335] | 10 | namespace xios |
---|
[219] | 11 | { |
---|
[802] | 12 | CNc4DataInput::CNc4DataInput(const StdString& filename, MPI_Comm comm_file, bool multifile, bool isCollective /*= true*/, const StdString& timeCounterName /*= "time_counter"*/) |
---|
[599] | 13 | : SuperClass() |
---|
[802] | 14 | , SuperClassWriter(filename, &comm_file, multifile, timeCounterName) |
---|
[599] | 15 | , comm_file(comm_file) |
---|
| 16 | , filename(filename) |
---|
| 17 | , isCollective(isCollective) |
---|
[775] | 18 | , readMetaDataDomains_(), readValueDomains_() |
---|
| 19 | , readMetaDataAxis_(), readValueAxis_() |
---|
[967] | 20 | , readMetaDataScalar_(), readValueScalar_() |
---|
[599] | 21 | { |
---|
| 22 | SuperClass::type = multifile ? MULTI_FILE : ONE_FILE; |
---|
| 23 | } |
---|
[219] | 24 | |
---|
[599] | 25 | CNc4DataInput::~CNc4DataInput(void) |
---|
| 26 | { /* Nothing more to do */ } |
---|
[219] | 27 | |
---|
[599] | 28 | StdSize CNc4DataInput::getFieldNbRecords_(CField* field) |
---|
| 29 | { |
---|
[770] | 30 | StdString fieldId = field->getFieldOutputName(); |
---|
[599] | 31 | |
---|
| 32 | if (SuperClassWriter::isTemporal(fieldId)) |
---|
| 33 | { |
---|
[811] | 34 | // return SuperClassWriter::getDimensions(&fieldId)[SuperClassWriter::getUnlimitedDimensionName()]; |
---|
| 35 | return SuperClassWriter::getDimensions(&fieldId)[SuperClassWriter::getTimeCounterName()]; |
---|
[599] | 36 | } |
---|
| 37 | |
---|
| 38 | return 1; |
---|
| 39 | } |
---|
| 40 | |
---|
| 41 | void CNc4DataInput::readFieldData_(CField* field) |
---|
| 42 | { |
---|
| 43 | CContext* context = CContext::getCurrent(); |
---|
| 44 | CContextServer* server = context->server; |
---|
| 45 | |
---|
| 46 | CGrid* grid = field->grid; |
---|
| 47 | |
---|
| 48 | if (!grid->doGridHaveDataToWrite()) |
---|
| 49 | if (SuperClass::type==MULTI_FILE || !isCollective) return; |
---|
| 50 | |
---|
[770] | 51 | StdString fieldId = field->getFieldOutputName(); |
---|
[599] | 52 | |
---|
| 53 | CArray<double,1> fieldData(grid->getWrittenDataSize()); |
---|
| 54 | if (!field->default_value.isEmpty()) fieldData = field->default_value; |
---|
| 55 | |
---|
| 56 | switch (SuperClass::type) |
---|
| 57 | { |
---|
| 58 | case MULTI_FILE: |
---|
[850] | 59 | SuperClassWriter::getData(fieldData, fieldId, isCollective, (field->getNStep() - 1)%field->nstepMax ); |
---|
[599] | 60 | break; |
---|
| 61 | case ONE_FILE: |
---|
| 62 | { |
---|
[765] | 63 | /* |
---|
[599] | 64 | std::vector<int> nZoomBeginGlobal = grid->getDistributionServer()->getZoomBeginGlobal(); |
---|
| 65 | std::vector<int> nZoomBeginServer = grid->getDistributionServer()->getZoomBeginServer(); |
---|
| 66 | std::vector<int> nZoomSizeServer = grid->getDistributionServer()->getZoomSizeServer(); |
---|
| 67 | |
---|
| 68 | int ssize = nZoomBeginGlobal.size(); |
---|
| 69 | |
---|
| 70 | std::vector<StdSize> start(ssize); |
---|
| 71 | std::vector<StdSize> count(ssize); |
---|
| 72 | |
---|
| 73 | for (int i = 0; i < ssize; ++i) |
---|
| 74 | { |
---|
| 75 | start[i] = nZoomBeginServer[ssize - i - 1] - nZoomBeginGlobal[ssize - i - 1]; |
---|
| 76 | count[i] = nZoomSizeServer[ssize - i - 1]; |
---|
| 77 | } |
---|
[765] | 78 | */ |
---|
[599] | 79 | |
---|
[765] | 80 | std::vector<int> nZoomBeginGlobal = grid->getDistributionServer()->getZoomBeginGlobal(); |
---|
| 81 | std::vector<int> nZoomBeginServer = grid->getDistributionServer()->getZoomBeginServer(); |
---|
| 82 | std::vector<int> nZoomSizeServer = grid->getDistributionServer()->getZoomSizeServer(); |
---|
| 83 | |
---|
| 84 | std::vector<StdSize> start, count; |
---|
| 85 | |
---|
[887] | 86 | CArray<int,1> axisDomainOrder = grid->axis_domain_order; |
---|
[765] | 87 | std::vector<StdString> domainList = grid->getDomainList(); |
---|
| 88 | std::vector<StdString> axisList = grid->getAxisList(); |
---|
| 89 | int numElement = axisDomainOrder.numElements(); |
---|
| 90 | int idxDomain = domainList.size() - 1, idxAxis = axisList.size() - 1; |
---|
| 91 | int idx = nZoomBeginGlobal.size() - 1; |
---|
| 92 | |
---|
| 93 | start.reserve(nZoomBeginGlobal.size()); |
---|
| 94 | count.reserve(nZoomBeginGlobal.size()); |
---|
| 95 | |
---|
| 96 | for (int i = numElement - 1; i >= 0; --i) |
---|
| 97 | { |
---|
[887] | 98 | if (2 == axisDomainOrder(i)) |
---|
[765] | 99 | { |
---|
| 100 | CDomain* domain = CDomain::get(domainList[idxDomain]); |
---|
| 101 | if ((domain->type) != CDomain::type_attr::unstructured) |
---|
| 102 | { |
---|
| 103 | start.push_back(nZoomBeginServer[idx] - nZoomBeginGlobal[idx]); |
---|
| 104 | count.push_back(nZoomSizeServer[idx]); |
---|
| 105 | } |
---|
| 106 | --idx ; |
---|
| 107 | start.push_back(nZoomBeginServer[idx] - nZoomBeginGlobal[idx]); |
---|
| 108 | count.push_back(nZoomSizeServer[idx]); |
---|
| 109 | --idx ; |
---|
| 110 | --idxDomain; |
---|
| 111 | } |
---|
[887] | 112 | else if (1 == axisDomainOrder(i)) |
---|
[765] | 113 | { |
---|
| 114 | start.push_back(nZoomBeginServer[idx] - nZoomBeginGlobal[idx]); |
---|
| 115 | count.push_back(nZoomSizeServer[idx]); |
---|
| 116 | --idx; |
---|
[967] | 117 | } |
---|
| 118 | else |
---|
| 119 | { |
---|
| 120 | if (1 == axisDomainOrder.numElements()) |
---|
| 121 | { |
---|
| 122 | start.push_back(0); |
---|
| 123 | count.push_back(1); |
---|
| 124 | } |
---|
| 125 | --idx; |
---|
| 126 | } |
---|
[765] | 127 | } |
---|
[777] | 128 | |
---|
[850] | 129 | SuperClassWriter::getData(fieldData, fieldId, isCollective, (field->getNStep() - 1)%field->nstepMax, &start, &count); |
---|
[599] | 130 | break; |
---|
| 131 | } |
---|
| 132 | } |
---|
| 133 | |
---|
| 134 | field->inputField(fieldData); |
---|
| 135 | |
---|
| 136 | if (!field->scale_factor.isEmpty() || !field->add_offset.isEmpty()) |
---|
| 137 | { |
---|
| 138 | double scaleFactor = 1.0, addOffset = 0.0; |
---|
| 139 | if (!field->scale_factor.isEmpty()) scaleFactor = field->scale_factor; |
---|
| 140 | if (!field->add_offset.isEmpty()) addOffset = field->add_offset; |
---|
| 141 | field->invertScaleFactorAddOffset(scaleFactor, addOffset); |
---|
| 142 | } |
---|
| 143 | } |
---|
| 144 | |
---|
[775] | 145 | void CNc4DataInput::readFieldAttributes_(CField* field, bool readAttributeValues) |
---|
| 146 | { |
---|
[777] | 147 | StdString fieldId = field->getFieldOutputName(); |
---|
[775] | 148 | |
---|
| 149 | CGrid* grid = field->grid; |
---|
| 150 | |
---|
| 151 | std::vector<CDomain*> domainP = grid->getDomains(); |
---|
| 152 | std::vector<CAxis*> axisP = grid->getAxis(); |
---|
[967] | 153 | std::vector<CScalar*> scalarP = grid->getScalars(); |
---|
[775] | 154 | int gridDim = domainP.size() * 2 + axisP.size(); |
---|
| 155 | |
---|
[967] | 156 | // Nothing to do with scalar without timestep |
---|
| 157 | if ((0 == gridDim) && (!SuperClassWriter::isTemporal(fieldId))) |
---|
| 158 | return; |
---|
| 159 | |
---|
[775] | 160 | // Verify the compatibility of dimension of declared grid and real grid in file |
---|
| 161 | int realGridDim = 1; |
---|
[783] | 162 | bool isUnstructuredGrid = SuperClassWriter::isUnstructured(fieldId); |
---|
[775] | 163 | std::map<StdString, StdSize> dimSizeMap = SuperClassWriter::getDimensions(&fieldId); |
---|
[807] | 164 | std::list<StdString> dimList = SuperClassWriter::getDimensionsList(&fieldId); |
---|
[825] | 165 | |
---|
[775] | 166 | realGridDim = SuperClassWriter::isTemporal(fieldId) ? dimSizeMap.size() - 1 : dimSizeMap.size(); |
---|
[783] | 167 | if (isUnstructuredGrid) ++realGridDim; |
---|
[775] | 168 | |
---|
| 169 | if (gridDim != realGridDim) |
---|
| 170 | ERROR("CNc4DataInput::readFieldAttributes_(CField* field, bool readAttributeValues)", |
---|
| 171 | << "Field '" << fieldId << "' has incorrect dimension " << std::endl |
---|
| 172 | << "Verify dimension of grid defined by 'grid_ref' or 'domain_ref'/'axis_ref' and dimension of grid in read file."); |
---|
| 173 | |
---|
| 174 | // Remove unlimited dimension from the map, we dont need it anymore |
---|
[825] | 175 | if (SuperClassWriter::isTemporal(fieldId)) |
---|
[807] | 176 | { |
---|
| 177 | dimSizeMap.erase(SuperClassWriter::getUnlimitedDimensionName()); |
---|
| 178 | dimList.pop_front() ; // assume time dimension is first |
---|
| 179 | } |
---|
[825] | 180 | |
---|
[783] | 181 | std::list<std::pair<StdString, StdSize> > listDimSize; |
---|
[807] | 182 | /* |
---|
[783] | 183 | for (std::map<StdString, StdSize>::const_iterator itMap = dimSizeMap.begin(); itMap != dimSizeMap.end(); ++itMap) |
---|
| 184 | listDimSize.push_front(*itMap); |
---|
[807] | 185 | */ |
---|
| 186 | for (std::list<StdString>::const_iterator it = dimList.begin(); it != dimList.end(); ++it) |
---|
| 187 | listDimSize.push_front(*dimSizeMap.find(*it)); |
---|
[775] | 188 | |
---|
| 189 | // Now process domain and axis |
---|
[887] | 190 | CArray<int,1> axisDomainOrder = grid->axis_domain_order; |
---|
[967] | 191 | int numElement = domainP.size() + axisP.size() + scalarP.size(); |
---|
[775] | 192 | int elementPosition = 0; |
---|
[967] | 193 | int idxDomain = 0, idxAxis = 0, idxScalar = 0; |
---|
[775] | 194 | |
---|
| 195 | std::pair<std::set<StdString>::iterator,bool> it; |
---|
| 196 | for (int i = 0; i < numElement; ++i) |
---|
| 197 | { |
---|
[887] | 198 | if(2 == axisDomainOrder(i)) |
---|
[775] | 199 | { |
---|
| 200 | if (readAttributeValues) |
---|
| 201 | { |
---|
| 202 | it = readValueDomains_.insert(domainP[idxDomain]->getId()); |
---|
[783] | 203 | if (it.second) readDomainAttributeValueFromFile(domainP[idxDomain], listDimSize, elementPosition, fieldId); |
---|
[775] | 204 | } |
---|
| 205 | else |
---|
| 206 | { |
---|
| 207 | it = readMetaDataDomains_.insert(domainP[idxDomain]->getId()); |
---|
[783] | 208 | if (it.second) readDomainAttributesFromFile(domainP[idxDomain], listDimSize, elementPosition, fieldId); |
---|
[775] | 209 | } |
---|
| 210 | ++idxDomain; |
---|
[783] | 211 | if (isUnstructuredGrid) ++elementPosition; |
---|
| 212 | else elementPosition += 2; |
---|
[775] | 213 | } |
---|
[887] | 214 | else if (1 == axisDomainOrder(i)) |
---|
[775] | 215 | { |
---|
| 216 | if (readAttributeValues) |
---|
| 217 | { |
---|
| 218 | it = readValueAxis_.insert(axisP[idxAxis]->getId()); |
---|
[783] | 219 | if (it.second) readAxisAttributeValueFromFile(axisP[idxAxis], listDimSize, elementPosition, fieldId); |
---|
[775] | 220 | } |
---|
| 221 | else |
---|
| 222 | { |
---|
| 223 | it = readMetaDataAxis_.insert(axisP[idxAxis]->getId()); |
---|
[783] | 224 | if (it.second) readAxisAttributesFromFile(axisP[idxAxis], listDimSize, elementPosition, fieldId); |
---|
[775] | 225 | } |
---|
| 226 | ++idxAxis; |
---|
| 227 | ++elementPosition; |
---|
| 228 | } |
---|
[967] | 229 | else |
---|
| 230 | { |
---|
| 231 | if (readAttributeValues) |
---|
| 232 | { |
---|
| 233 | it = readValueScalar_.insert(scalarP[idxScalar]->getId()); |
---|
| 234 | if (it.second) readScalarAttributeValueFromFile(scalarP[idxScalar], listDimSize, elementPosition, fieldId); |
---|
| 235 | } |
---|
| 236 | else |
---|
| 237 | { |
---|
| 238 | it = readMetaDataScalar_.insert(scalarP[idxScalar]->getId()); |
---|
| 239 | if (it.second) readScalarAttributesFromFile(scalarP[idxScalar], listDimSize, elementPosition, fieldId); |
---|
| 240 | } |
---|
| 241 | ++idxScalar; |
---|
| 242 | ++elementPosition; |
---|
| 243 | } |
---|
[775] | 244 | } |
---|
| 245 | } |
---|
| 246 | |
---|
| 247 | /*! |
---|
| 248 | Read attributes of a domain from a file |
---|
| 249 | \param [in] domain domain whose attributes are read from the file |
---|
| 250 | \param [in] dimSizeMap Dimensions and and their corresponding names and size read from file |
---|
| 251 | \param [in] emelentPosition position of domain in grid |
---|
[782] | 252 | \param [in] fieldId id (or name) associated with the grid |
---|
[775] | 253 | */ |
---|
[783] | 254 | void CNc4DataInput::readDomainAttributeValueFromFile(CDomain* domain, std::list<std::pair<StdString, StdSize> >& dimSizeMap, |
---|
[775] | 255 | int elementPosition, const StdString& fieldId) |
---|
| 256 | { |
---|
| 257 | // There are some optional attributes of a domain to retrieve from file // + lon lat? |
---|
[783] | 258 | std::list<std::pair<StdString, StdSize> >::const_iterator itMapNi = dimSizeMap.begin(), itMapNj, |
---|
| 259 | iteMap = dimSizeMap.end(); |
---|
[775] | 260 | |
---|
[783] | 261 | for (int i = 0; i < elementPosition; ++i, ++itMapNi) {} |
---|
| 262 | itMapNj = itMapNi; ++itMapNj; |
---|
[775] | 263 | |
---|
[1158] | 264 | if ((CDomain::type_attr::rectilinear == domain->type)) |
---|
[775] | 265 | { |
---|
| 266 | // Ok, try to read some f.. attributes such as longitude and latitude |
---|
[825] | 267 | bool hasLat = SuperClassWriter::hasVariable(itMapNj->first); |
---|
| 268 | if (hasLat) |
---|
| 269 | { |
---|
| 270 | domain->latvalue_rectilinear_read_from_file.resize(itMapNj->second); |
---|
| 271 | std::vector<StdSize> nBeginLat(1, 0), nSizeLat(1, itMapNj->second); |
---|
| 272 | readFieldVariableValue(domain->latvalue_rectilinear_read_from_file, itMapNj->first, nBeginLat, nSizeLat, true); |
---|
| 273 | } |
---|
[775] | 274 | |
---|
[825] | 275 | bool hasLon = SuperClassWriter::hasVariable(itMapNi->first); |
---|
| 276 | if (hasLon) |
---|
| 277 | { |
---|
| 278 | domain->lonvalue_rectilinear_read_from_file.resize(itMapNi->second); |
---|
| 279 | std::vector<StdSize> nBeginLon(1, 0), nSizeLon(1, itMapNi->second); |
---|
| 280 | readFieldVariableValue(domain->lonvalue_rectilinear_read_from_file, itMapNi->first, nBeginLon, nSizeLon, true); |
---|
[1158] | 281 | } |
---|
[775] | 282 | } |
---|
[1158] | 283 | else if ((CDomain::type_attr::curvilinear == domain->type)) |
---|
[783] | 284 | { |
---|
| 285 | int ni = domain->ni; |
---|
[785] | 286 | int nj = domain->nj; |
---|
| 287 | std::vector<StdSize> nBeginLatLon(2), nSizeLatLon(2); |
---|
[1158] | 288 | nBeginLatLon[0] = 0; nBeginLatLon[1] = 0; |
---|
| 289 | nSizeLatLon[0] = domain->nj_glo.getValue(); nSizeLatLon[1] = domain->ni_glo.getValue(); |
---|
[783] | 290 | |
---|
| 291 | StdString latName = this->getLatCoordName(fieldId); |
---|
[1158] | 292 | if (SuperClassWriter::hasVariable(latName)) |
---|
[825] | 293 | { |
---|
[1158] | 294 | domain->latvalue_curvilinear_read_from_file.resize(domain->ni_glo,domain->nj_glo); |
---|
| 295 | readFieldVariableValue(domain->latvalue_curvilinear_read_from_file, latName, nBeginLatLon, nSizeLatLon); |
---|
[825] | 296 | } |
---|
[783] | 297 | StdString lonName = this->getLonCoordName(fieldId); |
---|
[1158] | 298 | if (SuperClassWriter::hasVariable(lonName)) |
---|
[825] | 299 | { |
---|
[1158] | 300 | domain->lonvalue_curvilinear_read_from_file.resize(domain->ni_glo,domain->nj_glo); |
---|
| 301 | readFieldVariableValue(domain->lonvalue_curvilinear_read_from_file, lonName, nBeginLatLon, nSizeLatLon); |
---|
[825] | 302 | } |
---|
[783] | 303 | |
---|
| 304 | StdString boundsLatName = this->getBoundsId(latName); |
---|
| 305 | StdString boundsLonName = this->getBoundsId(lonName); |
---|
| 306 | |
---|
| 307 | int nbVertex = this->getNbVertex(fieldId); |
---|
[1158] | 308 | if (!domain->nvertex.isEmpty() && (domain->nvertex != nbVertex)) |
---|
| 309 | { |
---|
| 310 | ERROR("void CNc4DataInput::readDomainAttributeValueFromFile(...)", |
---|
| 311 | << "The domain " << domain->getDomainOutputName() |
---|
| 312 | << " has nvertex read from file " << this->filename << " and nvertex provided from model" |
---|
| 313 | << " are not coherent. They should be the same." << std::endl |
---|
| 314 | << " nvertex read from file: "<< nbVertex |
---|
| 315 | << " nvertex from model: "<< domain->nvertex << std::endl); |
---|
| 316 | } |
---|
| 317 | |
---|
| 318 | if (SuperClassWriter::hasVariable(boundsLatName) || SuperClassWriter::hasVariable(boundsLonName)) |
---|
[825] | 319 | domain->nvertex.setValue(nbVertex); |
---|
[785] | 320 | std::vector<StdSize> nBeginBndsLatLon(3), nSizeBndsLatLon(3); |
---|
[1158] | 321 | nBeginBndsLatLon[0] = 0; nSizeBndsLatLon[0] = domain->nj_glo.getValue(); |
---|
| 322 | nBeginBndsLatLon[1] = 0; nSizeBndsLatLon[1] = domain->nj_glo.getValue(); |
---|
[785] | 323 | nBeginBndsLatLon[2] = 0; nSizeBndsLatLon[2] = nbVertex; |
---|
[783] | 324 | |
---|
[1158] | 325 | if (SuperClassWriter::hasVariable(boundsLatName)) |
---|
[825] | 326 | { |
---|
[1158] | 327 | domain->bounds_latvalue_curvilinear_read_from_file.resize(nbVertex,domain->ni_glo,domain->nj_glo); |
---|
| 328 | readFieldVariableValue(domain->bounds_latvalue_curvilinear_read_from_file, boundsLatName, nBeginBndsLatLon, nSizeBndsLatLon); |
---|
[825] | 329 | |
---|
| 330 | } |
---|
[1158] | 331 | if (SuperClassWriter::hasVariable(boundsLonName)) |
---|
[825] | 332 | { |
---|
[1158] | 333 | domain->bounds_lonvalue_curvilinear_read_from_file.resize(nbVertex,domain->ni_glo,domain->nj_glo); |
---|
| 334 | readFieldVariableValue(domain->bounds_lonvalue_curvilinear_read_from_file, boundsLonName, nBeginBndsLatLon, nSizeBndsLatLon); |
---|
| 335 | } |
---|
[783] | 336 | } |
---|
[825] | 337 | else if ((CDomain::type_attr::unstructured == domain->type))// || (this->isUnstructured(fieldId))) |
---|
[775] | 338 | { |
---|
[785] | 339 | std::vector<StdSize> nBeginLatLon(1,0), nSizeLatLon(1,0); |
---|
| 340 | nSizeLatLon[0] = domain->ni_glo.getValue(); |
---|
| 341 | CArray<double,1> globalLonLat(domain->ni_glo.getValue()); |
---|
| 342 | |
---|
[782] | 343 | StdString latName = this->getLatCoordName(fieldId); |
---|
[1158] | 344 | if (SuperClassWriter::hasVariable(latName)) |
---|
[825] | 345 | { |
---|
[1158] | 346 | domain->latvalue_unstructured_read_from_file.resize(domain->ni_glo); |
---|
| 347 | readFieldVariableValue(domain->latvalue_unstructured_read_from_file, latName, nBeginLatLon, nSizeLatLon); |
---|
[825] | 348 | } |
---|
[785] | 349 | |
---|
[782] | 350 | StdString lonName = this->getLonCoordName(fieldId); |
---|
[825] | 351 | if (SuperClassWriter::hasVariable(lonName)) //(0 != lonName.compare("")) |
---|
| 352 | { |
---|
[1158] | 353 | // readFieldVariableValue(globalLonLat, lonName, nBeginLatLon, nSizeLatLon); |
---|
| 354 | domain->lonvalue_unstructured_read_from_file.resize(domain->ni_glo); |
---|
| 355 | readFieldVariableValue(domain->lonvalue_unstructured_read_from_file, lonName, nBeginLatLon, nSizeLatLon); |
---|
[825] | 356 | } |
---|
[782] | 357 | |
---|
| 358 | StdString boundsLatName = this->getBoundsId(latName); |
---|
| 359 | StdString boundsLonName = this->getBoundsId(lonName); |
---|
| 360 | |
---|
[783] | 361 | int nbVertex = this->getNbVertex(fieldId); |
---|
[1158] | 362 | if (!domain->nvertex.isEmpty() && (domain->nvertex != nbVertex)) |
---|
| 363 | { |
---|
| 364 | ERROR("void CNc4DataInput::readDomainAttributeValueFromFile(...)", |
---|
| 365 | << "The domain " << domain->getDomainOutputName() |
---|
| 366 | << " has nvertex read from file " << this->filename << " and nvertex provided from model" |
---|
| 367 | << " are not coherent. They should be the same." << std::endl |
---|
| 368 | << " nvertex read from file: "<< nbVertex |
---|
| 369 | << " nvertex from model: "<< domain->nvertex << std::endl); |
---|
| 370 | } |
---|
| 371 | |
---|
| 372 | if (SuperClassWriter::hasVariable(boundsLatName) || SuperClassWriter::hasVariable(boundsLonName)) |
---|
[825] | 373 | domain->nvertex.setValue(nbVertex); |
---|
| 374 | |
---|
[785] | 375 | std::vector<StdSize> nBeginBndsLatLon(2), nSizeBndsLatLon(2); |
---|
| 376 | nBeginBndsLatLon[0] = 0; nSizeBndsLatLon[0] = domain->ni_glo.getValue(); |
---|
| 377 | nBeginBndsLatLon[1] = 0; nSizeBndsLatLon[1] = nbVertex; |
---|
[782] | 378 | |
---|
[1158] | 379 | if (SuperClassWriter::hasVariable(boundsLatName)) |
---|
[825] | 380 | { |
---|
[1158] | 381 | domain->bounds_latvalue_unstructured_read_from_file.resize(nSizeBndsLatLon[1], nSizeBndsLatLon[0]); |
---|
| 382 | readFieldVariableValue(domain->bounds_latvalue_unstructured_read_from_file, boundsLatName, nBeginBndsLatLon, nSizeBndsLatLon); |
---|
[825] | 383 | } |
---|
[785] | 384 | |
---|
[1158] | 385 | if (SuperClassWriter::hasVariable(boundsLonName)) |
---|
[825] | 386 | { |
---|
[1158] | 387 | domain->bounds_lonvalue_unstructured_read_from_file.resize(nSizeBndsLatLon[1], nSizeBndsLatLon[0]); |
---|
| 388 | readFieldVariableValue(domain->bounds_lonvalue_unstructured_read_from_file, boundsLonName, nBeginBndsLatLon, nSizeBndsLatLon); |
---|
| 389 | } |
---|
[775] | 390 | } |
---|
[1158] | 391 | domain->fillInLonLat(); |
---|
[775] | 392 | } |
---|
| 393 | |
---|
| 394 | /*! |
---|
[782] | 395 | Read attribute value of a domain from a file |
---|
[775] | 396 | \param [in] domain domain whose attributes are read from the file |
---|
| 397 | \param [in] dimSizeMap Dimensions and and their corresponding names and size read from file |
---|
| 398 | \param [in] emelentPosition position of domain in grid |
---|
[782] | 399 | \param [in] fieldId id (or name) associated with the grid |
---|
[775] | 400 | */ |
---|
[783] | 401 | void CNc4DataInput::readDomainAttributesFromFile(CDomain* domain, std::list<std::pair<StdString, StdSize> >& dimSizeMap, |
---|
[775] | 402 | int elementPosition, const StdString& fieldId) |
---|
| 403 | { |
---|
| 404 | // There are some mandatory attributes of a domain to retrieve from file |
---|
| 405 | // + ni_glo, nj_glo |
---|
[783] | 406 | std::list<std::pair<StdString, StdSize> >::const_iterator itMapNi = dimSizeMap.begin(), itMapNj, |
---|
| 407 | iteMap = dimSizeMap.end(); |
---|
| 408 | for (int i = 0; i < elementPosition; ++i, ++itMapNi) {} |
---|
| 409 | itMapNj = itMapNi; ++itMapNj; |
---|
[775] | 410 | |
---|
[782] | 411 | if (this->isRectilinear(fieldId) || this->isCurvilinear(fieldId)) |
---|
[775] | 412 | { |
---|
[1158] | 413 | if (!domain->nj_glo.isEmpty() && (domain->nj_glo != itMapNj->second)) |
---|
| 414 | { |
---|
| 415 | ERROR("void CNc4DataInput::readDomainAttributesFromFile(...)", |
---|
| 416 | << "The domain " << domain->getDomainOutputName() |
---|
| 417 | << " has nj_glo read from file " << this->filename << " and nj_glo provided from model" |
---|
| 418 | << " are not coherent. They should be the same." << std::endl |
---|
| 419 | << " nj_glo read from file: "<< itMapNj->second |
---|
| 420 | << " nj_glo from model: "<< domain->nj_glo << std::endl); |
---|
| 421 | } |
---|
[775] | 422 | domain->nj_glo.setValue(itMapNj->second); |
---|
[1158] | 423 | |
---|
| 424 | if (!domain->ni_glo.isEmpty() && (domain->ni_glo != itMapNi->second)) |
---|
| 425 | { |
---|
| 426 | ERROR("void CNc4DataInput::readDomainAttributesFromFile(...)", |
---|
| 427 | << "The domain " << domain->getDomainOutputName() |
---|
| 428 | << " has ni_glo read from file " << this->filename << " and ni_glo provided from model" |
---|
| 429 | << " are not coherent. They should be the same." << std::endl |
---|
| 430 | << " ni_glo read from file: "<< itMapNi->second |
---|
| 431 | << " ni_glo from model: "<< domain->ni_glo << std::endl); |
---|
| 432 | } |
---|
[783] | 433 | domain->ni_glo.setValue(itMapNi->second); |
---|
[775] | 434 | } |
---|
| 435 | else if (this->isUnstructured(fieldId)) |
---|
| 436 | { |
---|
[783] | 437 | domain->nj_glo.setValue(1); |
---|
[1158] | 438 | |
---|
| 439 | if (!domain->ni_glo.isEmpty() && (domain->ni_glo != itMapNi->second)) |
---|
| 440 | { |
---|
| 441 | ERROR("void CNc4DataInput::readDomainAttributesFromFile(...)", |
---|
| 442 | << "The domain " << domain->getDomainOutputName() |
---|
| 443 | << " has ni_glo read from file " << this->filename << " and ni_glo provided from model" |
---|
| 444 | << " are not coherent. They should be the same." << std::endl |
---|
| 445 | << " ni_glo read from file: "<< itMapNi->second |
---|
| 446 | << " ni_glo from model: "<< domain->ni_glo << std::endl); |
---|
| 447 | } |
---|
[783] | 448 | domain->ni_glo.setValue(itMapNi->second); |
---|
[775] | 449 | } |
---|
| 450 | } |
---|
| 451 | |
---|
| 452 | /*! |
---|
| 453 | Read attributes of an axis from a file |
---|
| 454 | \param [in] axis axis whose attributes are read from the file |
---|
| 455 | \param [in] dimSizeMap Dimensions and and their corresponding names and size read from file |
---|
| 456 | \param [in] emelentPosition position of axis in grid |
---|
[782] | 457 | \param [in] fieldId id (or name) associated with the grid |
---|
[775] | 458 | */ |
---|
[783] | 459 | void CNc4DataInput::readAxisAttributesFromFile(CAxis* axis, std::list<std::pair<StdString, StdSize> >& dimSizeMap, |
---|
[775] | 460 | int elementPosition, const StdString& fieldId) |
---|
| 461 | { |
---|
[783] | 462 | std::list<std::pair<StdString, StdSize> >::const_iterator itMapN = dimSizeMap.begin(), |
---|
| 463 | iteMap = dimSizeMap.end(); |
---|
[775] | 464 | for (int i = 0; i < elementPosition; ++i, ++itMapN) {} |
---|
[1158] | 465 | |
---|
| 466 | if (!axis->n_glo.isEmpty() && (axis->n_glo != itMapN->second)) |
---|
| 467 | { |
---|
| 468 | ERROR("void CNc4DataInput::readAxisAttributesFromFile(...)", |
---|
| 469 | << "The axis " << axis->getAxisOutputName() |
---|
| 470 | << " has n_glo read from file " << this->filename << " and n_glo provided from model" |
---|
| 471 | << " are not coherent. They should be the same." << std::endl |
---|
| 472 | << " n_glo read from file: "<< itMapN->second |
---|
| 473 | << " n_glo from model: "<< axis->n_glo << std::endl); |
---|
| 474 | } |
---|
[775] | 475 | axis->n_glo.setValue(itMapN->second); |
---|
| 476 | } |
---|
| 477 | |
---|
| 478 | /*! |
---|
[782] | 479 | Read attribute value of an axis from a file |
---|
[775] | 480 | \param [in] axis axis whose attributes are read from the file |
---|
| 481 | \param [in] dimSizeMap Dimensions and and their corresponding names and size read from file |
---|
| 482 | \param [in] emelentPosition position of axis in grid |
---|
[782] | 483 | \param [in] fieldId id (or name) associated with the grid |
---|
[775] | 484 | */ |
---|
[783] | 485 | void CNc4DataInput::readAxisAttributeValueFromFile(CAxis* axis, std::list<std::pair<StdString, StdSize> >& dimSizeMap, |
---|
[775] | 486 | int elementPosition, const StdString& fieldId) |
---|
| 487 | { |
---|
[783] | 488 | std::list<std::pair<StdString, StdSize> >::const_iterator itMapN = dimSizeMap.begin(), |
---|
| 489 | iteMap = dimSizeMap.end(); |
---|
[775] | 490 | for (int i = 0; i < elementPosition; ++i, ++itMapN) {} |
---|
| 491 | |
---|
| 492 | { // Read axis value |
---|
| 493 | std::vector<StdSize> nBegin(1, 0), nSize(1, itMapN->second); |
---|
| 494 | CArray<double,1> readAxisValue(itMapN->second); |
---|
| 495 | readFieldVariableValue(readAxisValue, itMapN->first, nBegin, nSize, true); |
---|
| 496 | int begin = 0, n = itMapN->second; |
---|
| 497 | if (!axis->begin.isEmpty()) begin = axis->begin.getValue(); |
---|
| 498 | if (!axis->n.isEmpty()) n = axis->n.getValue(); |
---|
| 499 | axis->value.resize(n); |
---|
| 500 | for (int i = 0; i < n; ++i) axis->value(i) = readAxisValue(begin + i); |
---|
| 501 | } |
---|
| 502 | } |
---|
| 503 | |
---|
[967] | 504 | /*! |
---|
| 505 | Read attributes of a scalar from a file |
---|
| 506 | \param [in] scalar scalar whose attributes are read from the file |
---|
| 507 | \param [in] dimSizeMap Dimensions and and their corresponding names and size read from file |
---|
| 508 | \param [in] emelentPosition position of scalar in grid |
---|
| 509 | \param [in] fieldId id (or name) associated with the grid |
---|
| 510 | */ |
---|
| 511 | void CNc4DataInput::readScalarAttributesFromFile(CScalar* scalar, std::list<std::pair<StdString, StdSize> >& dimSizeMap, |
---|
| 512 | int elementPosition, const StdString& fieldId) |
---|
| 513 | { |
---|
[1158] | 514 | /*Nothing to do */ |
---|
[967] | 515 | } |
---|
| 516 | |
---|
| 517 | /*! |
---|
| 518 | Read attribute value of an axis from a file |
---|
| 519 | \param [in] axis axis whose attributes are read from the file |
---|
| 520 | \param [in] dimSizeMap Dimensions and and their corresponding names and size read from file |
---|
| 521 | \param [in] emelentPosition position of axis in grid |
---|
| 522 | \param [in] fieldId id (or name) associated with the grid |
---|
| 523 | */ |
---|
| 524 | void CNc4DataInput::readScalarAttributeValueFromFile(CScalar* scalar, std::list<std::pair<StdString, StdSize> >& dimSizeMap, |
---|
| 525 | int elementPosition, const StdString& fieldId) |
---|
| 526 | { |
---|
[1158] | 527 | /*Nothing to do */ |
---|
[967] | 528 | } |
---|
| 529 | |
---|
[599] | 530 | void CNc4DataInput::closeFile_(void) |
---|
| 531 | { |
---|
| 532 | SuperClassWriter::close(); |
---|
| 533 | } |
---|
[335] | 534 | } // namespace xios |
---|