- Timestamp:
- 02/10/16 14:35:03 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
XIOS/trunk/src/transformation/domain_algorithm_interpolate.cpp
r809 r821 15 15 #include "netcdf.hpp" 16 16 #include "mapper.hpp" 17 #include "mpi_tag.hpp" 17 18 18 19 namespace xios { … … 95 96 double step=(domainSrc_->lon_end-domainSrc_->lon_start)/domainSrc_->ni_glo ; 96 97 for(int i=0; i<domainSrc_->ni_glo; ++i) lon_g(i)=domainSrc_->lon_start+i*step ; 97 step=(domainSrc_->lat_end-domainSrc_->lat_start)/domainSrc_->nj_glo ; 98 step=(domainSrc_->lat_end-domainSrc_->lat_start)/domainSrc_->nj_glo ; 98 99 for(int i=0; i<domainSrc_->ni_glo; ++i) lat_g(i)=domainSrc_->lat_start+i*step ; 99 100 } … … 160 161 double step=(domainDest_->lon_end-domainDest_->lon_start)/domainDest_->ni_glo ; 161 162 for(int i=0; i<domainDest_->ni_glo; ++i) lon_g(i)=domainDest_->lon_start+i*step ; 162 step=(domainDest_->lat_end-domainDest_->lat_start)/domainDest_->nj_glo ; 163 step=(domainDest_->lat_end-domainDest_->lat_start)/domainDest_->nj_glo ; 163 164 for(int i=0; i<domainDest_->ni_glo; ++i) lat_g(i)=domainDest_->lat_start+i*step ; 164 165 } … … 167 168 if (std::abs(poleValue - std::abs(lat_g(domainDest_->nj_glo-1))) < NumTraits<double>::epsilon()) isSouthPole = true; 168 169 169 170 171 170 171 172 172 173 if (isNorthPole && (0 == domainDest_->jbegin.getValue())) 173 174 { … … 469 470 MPI_INT, 470 471 itMap->first, 471 7,472 MPI_DOMAIN_INTERPOLATION_DEST_INDEX, 472 473 client->intraComm, 473 474 &sendRequest.back()); … … 477 478 MPI_INT, 478 479 itMap->first, 479 8,480 MPI_DOMAIN_INTERPOLATION_SRC_INDEX, 480 481 client->intraComm, 481 482 &sendRequest.back()); … … 485 486 MPI_DOUBLE, 486 487 itMap->first, 487 9,488 MPI_DOMAIN_INTERPOLATION_WEIGHT, 488 489 client->intraComm, 489 490 &sendRequest.back()); … … 504 505 MPI_INT, 505 506 MPI_ANY_SOURCE, 506 7,507 MPI_DOMAIN_INTERPOLATION_DEST_INDEX, 507 508 client->intraComm, 508 509 &recvStatus); … … 516 517 MPI_INT, 517 518 clientSrcRank, 518 8,519 MPI_DOMAIN_INTERPOLATION_SRC_INDEX, 519 520 client->intraComm, 520 521 &recvStatus); … … 524 525 MPI_DOUBLE, 525 526 clientSrcRank, 526 9,527 MPI_DOMAIN_INTERPOLATION_WEIGHT, 527 528 client->intraComm, 528 529 &recvStatus);
Note: See TracChangeset
for help on using the changeset viewer.