Changeset 2521

Timestamp:

10/23/09 12:10:51 (16 years ago)

Author:

seroussi

Message:

added extrusion of pressure if there is an existing pressure

File:

• issm/trunk/src/m/classes/public/extrude.m (modified) (1 diff)

issm/trunk/src/m/classes/public/extrude.m

@@ -211 +211 @@
 md.gridonwater=project3d(md,md.gridonwater,'node');
 
+%Put lithostatic pressure is there is an existing pressure
+if ~isnan(md.pressure),
+        md.pressure=md.g*md.rho_ice*(md.surface-md.z);
+end
+
 %special for thermal modeling:
 md.melting=project3d(md,md.melting,'node',1);