Changeset 6244

Timestamp:

10/12/10 09:01:33 (15 years ago)

Author:

seroussi

Message:

fixed basalfriction problem in thermal ice sheets

File:

• issm/trunk/src/c/objects/Elements/Tria.cpp (modified) (1 diff)

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

@@ -4511 +4511 @@
                 geothermalflux_input->GetParameterValue(&geothermalflux_value,gauss);
                 friction->GetAlpha2(&alpha2,gauss,VxEnum,VyEnum,VzEnum);
+                vx_input->GetParameterValue(&vx,gauss);
+                vy_input->GetParameterValue(&vy,gauss);
+                basalfriction=alpha2*(pow(vx,2.0)+pow(vy,2.0));
                 
                 scalar=gauss->weight*Jdet*(basalfriction+geothermalflux_value)/(heatcapacity*rho_ice);