Changeset 1614 for XIOS/trunk/extern/remap/src/libmapper.cpp

Timestamp:

11/26/18 10:16:06 (5 years ago)

Author:

ymipsl

Message:

Interpolation : enhancement : you have now the possibility to give computed area to the remaper for better global conservation.

YM

File:

• XIOS/trunk/extern/remap/src/libmapper.cpp (modified) (1 diff)

XIOS/trunk/extern/remap/src/libmapper.cpp

@@ -41 +41 @@
         assert(n_cell_dst >= 4);
         assert(1 <= order && order <= 2);
-
+  double* src_area=NULL ;
+  double* dst_area=NULL ;
   mapper = new Mapper(MPI_COMM_WORLD);
   mapper->setVerbosity(PROGRESS) ;
-  mapper->setSourceMesh(src_bounds_lon, src_bounds_lat, n_vert_per_cell_src, n_cell_src, src_pole ) ;
-  mapper->setTargetMesh(dst_bounds_lon, dst_bounds_lat, n_vert_per_cell_dst, n_cell_dst, dst_pole ) ;
+  mapper->setSourceMesh(src_bounds_lon, src_bounds_lat, src_area, n_vert_per_cell_src, n_cell_src, src_pole ) ;
+  mapper->setTargetMesh(dst_bounds_lon, dst_bounds_lat, dst_area, n_vert_per_cell_dst, n_cell_dst, dst_pole ) ;
 
 /*