Ignore:
Timestamp:
05/12/20 11:52:13 (2 years ago)
Author:
ymipsl
Message:

XIOS coupling branch
Some updates.

First coupling test is beginning to work...

YM

File:
1 edited

Legend:

Unmodified
Added
Removed
  • XIOS/dev/dev_ym/XIOS_COUPLING/src/transformation/domain_algorithm_interpolate.cpp

    r1787 r1875  
    171171      step=(domainSrc_->lat_end-domainSrc_->lat_start)/domainSrc_->nj_glo ; 
    172172      for (int i=0; i<domainSrc_->ni_glo; ++i) lat_g(i)=domainSrc_->lat_start+i*step ; 
     173    } 
     174    else if (!domainSrc_->lonvalue.isEmpty() && !domainSrc_->latvalue.isEmpty())  // data can come possibly from coupler so only lonvalue and latvalue are defined ? 
     175    { 
     176      CArray<double,1> lon(niSrc), lat(njSrc) ; 
     177      for(int i=0; i< niSrc;++i) lon(i) = domainSrc_->lonvalue(i) ; 
     178      for(int j=0, n=0; j<njSrc; ++j, n+=niSrc) lat(j) = domainSrc_->latvalue(n) ; 
     179      domainSrc_->AllgatherRectilinearLonLat(lon, lat, lon_g, lat_g) ; 
    173180    } 
    174181    else ERROR("void CDomainAlgorithmInterpolate::computeRemap()",<<"Cannot compute bounds for rectilinear domain") ; 
Note: See TracChangeset for help on using the changeset viewer.