source: issm/oecreview/Archive/16133-16554/ISSM-16356-16357.diff@ 16556

../trunk-jpl/src/c/classes/Elements/Penta.cpp

@@ -4874 +4874 @@
                 /*Force temperature between [Tpmp-50 Tpmp] to disable penalties*/
                 IssmDouble pressure[numdof];
                 GetInputListOnVertices(&pressure[0],PressureEnum);
-                if(values[i]>matpar->TMeltingPoint(pressure[i]))     values[i]=matpar->TMeltingPoint(pressure[i]);
-                if(values[i]<matpar->TMeltingPoint(pressure[i])-50.) values[i]=matpar->TMeltingPoint(pressure[i])-50.;
+                for(i=0;i<numdof;i++){
+                        if(values[i]>matpar->TMeltingPoint(pressure[i]))     values[i]=matpar->TMeltingPoint(pressure[i]);
+                        if(values[i]<matpar->TMeltingPoint(pressure[i])-50.) values[i]=matpar->TMeltingPoint(pressure[i])-50.;
+                }
         }