Changeset 14679

Timestamp:

04/19/13 17:38:03 (12 years ago)

Author:

bdef

Message:

CHG: putting real parameters in hydro matrix building

File:

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

issm/trunk-jpl/src/c/classes/objects/Elements/Tria.cpp

@@ -5894 +5894 @@
         this->parameters->FindParam(&dt,TimesteppingTimeStepEnum);
 
-        /*Compute Storing coefficient fro the above */
+        /*Compute Storing coefficient from the above */
 
         sediment_storing = rho_freshwater*gravity*sediment_porosity*sediment_thickness*
@@ -5908 +5908 @@
 
                 /*Diffusivity*/
-                D_scalar=sediment_porosity*gauss->weight*Jdet;
+                D_scalar=sediment_transmitivity*gauss->weight*Jdet;
                 if(reCast<bool,IssmDouble>(dt)) D_scalar=D_scalar*dt;
                 D[0][0]=D_scalar; D[0][1]=0.;
@@ -5921 +5921 @@
                 if(reCast<bool,IssmDouble>(dt)){
                         GetL(&L[0], &xyz_list[0][0], gauss,NDOF1);
-                        D_scalar=gauss->weight*Jdet;
+                        D_scalar=sediment_storing*gauss->weight*Jdet;
                         
                         TripleMultiply(&L[0],numdof,1,0,