Changeset 15556

Timestamp:

07/23/13 11:32:34 (12 years ago)

Author:

seroussi

Message:

FIX: finally found the problem with stokes front

File:

• issm/trunk-jpl/src/c/classes/Elements/Penta.cpp (modified) (3 diffs)

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

@@ -8216 +8216 @@
         IssmDouble zmax=xyz_list[0][2]; for(int i=1;i<6;i++) if(xyz_list[i][2]>zmax) zmax=xyz_list[i][2];
         IssmDouble zmin=xyz_list[0][2]; for(int i=1;i<6;i++) if(xyz_list[i][2]<zmin) zmin=xyz_list[i][2];
-        if(zmax>0 && zmin<0) gauss=new GaussPenta(area_coordinates,3,10); //refined in vertical because of the sea level discontinuity
+        if(zmax>0 && zmin<0) gauss=new GaussPenta(area_coordinates,3,30); //refined in vertical because of the sea level discontinuity
         else                 gauss=new GaussPenta(area_coordinates,3,3);
 
@@ -8234 +8234 @@
                         pe->values[4*i+1]+= pressure*Jdet*gauss->weight*normal[1]*basis[i];
                         pe->values[4*i+2]+= pressure*Jdet*gauss->weight*normal[2]*basis[i];
+                        pe->values[4*i+3]+= 0;
                 }
         }
@@ -8239 +8240 @@
         /*Transform coordinate system*/
         TransformLoadVectorCoord(pe,nodes,numnodes,XYZPEnum);
-
 
         /*Clean up and return*/