Changeset 22506 for issm/trunk-jpl/src/c/analyses/MasstransportAnalysis.cpp

Timestamp:

03/07/18 10:19:07 (7 years ago)

Author:

Mathieu Morlighem

Message:

CHG: do not compute melt before marshalling

File:

• issm/trunk-jpl/src/c/analyses/MasstransportAnalysis.cpp (modified) (2 diffs)

issm/trunk-jpl/src/c/analyses/MasstransportAnalysis.cpp

@@ -156 +156 @@
         iomodel->FetchDataToInput(elements,"md.mask.groundedice_levelset",MaskGroundediceLevelsetEnum);
         iomodel->FetchDataToInput(elements,"md.basalforcings.groundedice_melting_rate",BasalforcingsGroundediceMeltingRateEnum);
-        iomodel->FetchDataToInput(elements,"md.basalforcings.floatingice_melting_rate",BasalforcingsFloatingiceMeltingRateEnum);
         iomodel->FetchDataToInput(elements,"md.initialization.vx",VxEnum);
         iomodel->FetchDataToInput(elements,"md.initialization.vy",VyEnum);
+
+        /*Get what we need for ocean-induced basal melting*/
+        int basalforcing_model;
+        iomodel->FindConstant(&basalforcing_model,"md.basalforcings.model");
+        switch(basalforcing_model){
+                case FloatingMeltRateEnum:
+                        iomodel->FetchDataToInput(elements,"md.basalforcings.floatingice_melting_rate",BasalforcingsFloatingiceMeltingRateEnum);
+                        break;
+                case LinearFloatingMeltRateEnum:
+                        break;
+                case MismipFloatingMeltRateEnum:
+                        break;
+                case MantlePlumeGeothermalFluxEnum:
+                        break;
+                case BasalforcingsPicoEnum:
+                        iomodel->FetchDataToInput(elements,"md.basalforcings.basin_id",BasalforcingsPicoBasinIdEnum);
+                        break;
+                default:
+                        _error_("Basal forcing model "<<EnumToStringx(basalforcing_model)<<" not supported yet");
+        }
 
         if(isoceancoupling){
@@ -177 +196 @@
                 iomodel->FetchDataToInput(elements,"md.mesh.vertexonbase",MeshVertexonbaseEnum);
                 iomodel->FetchDataToInput(elements,"md.mesh.vertexonsurface",MeshVertexonsurfaceEnum);
-        }
-
-        int basalforcing_model;
-        iomodel->FindConstant(&basalforcing_model,"md.basalforcings.model");
-        switch(basalforcing_model){
-                case BasalforcingsPicoEnum:
-                        iomodel->FetchDataToInput(elements,"md.basalforcings.basin_id",BasalforcingsPicoBasinIdEnum);
-                        break;
-                default: // Do nothing
-                        break;
         }