Changeset 21549

Timestamp:

02/10/17 11:49:37 (8 years ago)

Author:

Mathieu Morlighem

Message:

BUG: fixing enthalpy spc bug

Location:

issm/trunk-jpl/src/c

Files:

• analyses/EnthalpyAnalysis.cpp (modified) (1 diff)
• cores/transient_core.cpp (modified) (1 diff)

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

@@ -36 +36 @@
         _assert_(N>0); _assert_(M>=iomodel->numberofvertices);
         for(int i=0;i<iomodel->numberofvertices;i++){
-                for(int j=0;i<N;j++){
-                        spcvector[i*N+j] = heatcapacity*(spcvector[i]-referencetemperature);
+                for(int j=0;j<N;j++){
+                        spcvector[i*N+j] = heatcapacity*(spcvector[i*N+j]-referencetemperature);
                 }
         }

issm/trunk-jpl/src/c/cores/transient_core.cpp

@@ -86 +86 @@
                 if(isthermal && domaintype==Domain3DEnum){ 
                         if(issmb){
+                                bool isenthalpy;
+                                femmodel->parameters->FindParam(&isenthalpy,ThermalIsenthalpyEnum);
                                 femmodel->parameters->FindParam(&smb_model,SmbEnum);
-                                if(smb_model==SMBpddEnum) ResetBoundaryConditions(femmodel,ThermalAnalysisEnum);
-                                if(smb_model==SMBd18opddEnum) ResetBoundaryConditions(femmodel,ThermalAnalysisEnum);
+                                if(isenthalpy){
+                                        if(smb_model==SMBpddEnum)     ResetBoundaryConditions(femmodel,EnthalpyAnalysisEnum);
+                                        if(smb_model==SMBd18opddEnum) ResetBoundaryConditions(femmodel,EnthalpyAnalysisEnum);
+                                }
+                                else{
+                                        if(smb_model==SMBpddEnum)     ResetBoundaryConditions(femmodel,ThermalAnalysisEnum);
+                                        if(smb_model==SMBd18opddEnum) ResetBoundaryConditions(femmodel,ThermalAnalysisEnum);
+                                }
                         }
                         if(VerboseSolution()) _printf0_("   computing thermal regime\n");