Index: /issm/trunk/src/m/classes/public/extrude.m
===================================================================
--- /issm/trunk/src/m/classes/public/extrude.m	(revision 2520)
+++ /issm/trunk/src/m/classes/public/extrude.m	(revision 2521)
@@ -211,4 +211,9 @@
 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);