Changeset 16748 for issm/trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp

Timestamp:

11/13/13 17:12:08 (11 years ago)

Author:

seroussi

Message:

BUG: fixed enthalpy analysis solution

File:

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

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

@@ -237 +237 @@
                                 break;
                         case PatersonEnum:
-                                for(i=0;i<numnodes;i++) B[i]=Paterson(values[i]);
+                                for(i=0;i<numnodes;i++) B[i]=Paterson(temperature[i]);
                                 element->AddMaterialInput(MaterialsRheologyBEnum,&B[0],P1Enum);
                                 break;
                         case ArrheniusEnum:
                                 element->GetVerticesCoordinates(&xyz_list);
-                                for(i=0;i<numnodes;i++) B[i]=Arrhenius(values[i],surface[i]-xyz_list[i*3+2],element->GetMaterialParameter(MaterialsRheologyNEnum));
+                                for(i=0;i<numnodes;i++) B[i]=Arrhenius(temperature[i],surface[i]-xyz_list[i*3+2],element->GetMaterialParameter(MaterialsRheologyNEnum));
                                 element->AddMaterialInput(MaterialsRheologyBEnum,&B[0],P1Enum);
                                 break;
                         case LliboutryDuvalEnum:
                                 //for(i=0;i<numnodes;i++) B[i]=LliboutryDuval(values[i],pressure[i],material->GetN(),matpar->GetBeta(),matpar->GetReferenceTemperature(),matpar->GetHeatCapacity(),matpar->GetLatentHeat()); for(i=0;i<numnodes;i++) B[i]=Paterson(values[i]); 
-                                for(i=0;i<numnodes;i++) B[i]=Paterson(values[i]);
+                                for(i=0;i<numnodes;i++) B[i]=Paterson(temperature[i]);
                                 element->AddMaterialInput(MaterialsRheologyBEnum,&B[0],P1Enum); 
                                 break;