Changeset 18696

Timestamp:

10/28/14 10:18:54 (10 years ago)

Author:

seroussi

Message:

CHG: no need to change enthalpy at base to compute basal melting

File:

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

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

@@ -1003 +1003 @@
                                 for(i=0;i<3;i++) vec_heatflux[i]=0.;
                                 break;
-                        case 1:
-                                // cold, wet base: keep at pressure melting point 
-                        case 2:
-                                // temperate, thin refreezing base: release spc
-
-                        case 3:
-                                // temperate, thin melting base: set spc
-                                enthalpies[vertexdown]=element->PureIceEnthalpy(pressures[vertexdown]); 
+                        case 1: case 2: case 3: 
+                                // case 1 : cold, wet base: keep at pressure melting point 
+                                // case 2: temperate, thin refreezing base: release spc
+                                // case 3: temperate, thin melting base: set spc
                                 enthalpy_input->GetInputDerivativeValue(&d1enthalpy[0],xyz_list,gauss);
                                 for(i=0;i<3;i++) vec_heatflux[i]=-kappa*d1enthalpy[i];