Changeset 22778 for issm/trunk-jpl/src/c/modules/ModelProcessorx/CreateElementsVerticesAndMaterials.cpp

Timestamp:

05/15/18 13:42:52 (7 years ago)

Author:

seroussi

Message:

NEW: starting to interpolate data between ice and ocean grids

File:

• issm/trunk-jpl/src/c/modules/ModelProcessorx/CreateElementsVerticesAndMaterials.cpp (modified) (4 diffs)

issm/trunk-jpl/src/c/modules/ModelProcessorx/CreateElementsVerticesAndMaterials.cpp

@@ -16 +16 @@
         bool dakota_analysis;
         bool adolc_analysis;
+        bool isoceancoupling;
         int nnat,dummy;
         int* nature=NULL;
@@ -24 +25 @@
         iomodel->FindConstant(&materials_type,"md.materials.type");
         iomodel->FindConstant(&adolc_analysis,"md.autodiff.isautodiff");
+        iomodel->FindConstant(&isoceancoupling,"md.transient.isoceancoupling");
 
         /*Did we already create the elements? : */
@@ -232 +234 @@
         if (iomodel->domaintype == Domain3DsurfaceEnum) iomodel->FetchData(3,"md.mesh.lat","md.mesh.long","md.mesh.r");
         else iomodel->FetchDataToInput(elements,"md.mesh.scale_factor",MeshScaleFactorEnum,1.);
+        if (isoceancoupling) iomodel->FetchData(2,"md.mesh.lat","md.mesh.long");
         
         CreateNumberNodeToElementConnectivity(iomodel,solution_type);
@@ -242 +245 @@
         iomodel->DeleteData(6,"md.mesh.x","md.mesh.y","md.mesh.z","md.geometry.base","md.geometry.thickness","md.mask.ice_levelset");
         if (iomodel->domaintype == Domain3DsurfaceEnum) iomodel->DeleteData(3,"md.mesh.lat","md.mesh.long","md.mesh.r");
+        if (isoceancoupling) iomodel->DeleteData(2,"md.mesh.lat","md.mesh.long");
 }