Changeset 21548

Timestamp:

02/10/17 11:32:39 (8 years ago)

Author:

Mathieu Morlighem

Message:

NEW: enable dynamic spcs for enthalpy model (required for PDD scheme)

File:

• issm/trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp (modified) (1 diff)

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

@@ -40 +40 @@
                 }
         }
-        IoModelToConstraintsx(constraints,iomodel,spcvector,M,N,EnthalpyAnalysisEnum,finiteelement,0);
+
+        /*Specific case for PDD, we want the constaints to be updated by the PDD scheme itself*/
+        bool isdynamic = false;
+        if (iomodel->solution_enum==TransientSolutionEnum){
+                int smb_model;
+                iomodel->FindConstant(&smb_model,"md.smb.model");
+                if(smb_model==SMBpddEnum)     isdynamic=true;
+                if(smb_model==SMBd18opddEnum) isdynamic=true;
+        }
+
+        if(isdynamic){
+                IoModelToDynamicConstraintsx(constraints,iomodel,"md.thermal.spctemperature",ThermalAnalysisEnum,finiteelement);
+        }
+        else{
+                IoModelToConstraintsx(constraints,iomodel,spcvector,M,N,EnthalpyAnalysisEnum,finiteelement);
+        }
 
         /*Free ressources:*/