Changeset 335

Timestamp:

05/11/09 11:20:15 (16 years ago)

Author:

seroussi

Message:

problem of sign with melting

File:

• issm/trunk/src/c/objects/Tria.cpp (modified) (2 diffs)

issm/trunk/src/c/objects/Tria.cpp

@@ -2346 +2346 @@
                 /*Get bed slope: */
                 GetParameterDerivativeValue(&slope[0], &b[0],&xyz_list[0][0], gauss_l1l2l3);
-                dbdx=-slope[0];
-                dbdy=-slope[1];
+                dbdx=slope[0];
+                dbdy=slope[1];
 
                 /* Get Jacobian determinant: */
@@ -2358 +2358 @@
                 /*Build gaussian vector: */
                 for(i=0;i<numgrids;i++){
-                        pe_g_gaussian[i]=Jdet*gauss_weight*(vx*dbdx+vy*dbdy-meltingvalue)*L[i];
+                        pe_g_gaussian[i]=-Jdet*gauss_weight*(vx*dbdx+vy*dbdy-meltingvalue)*L[i];
                 }