Changeset 894

Timestamp:

06/11/09 08:50:11 (16 years ago)

Author:

seroussi

Message:

problem with pressure in ice transient

File:

• issm/trunk/src/m/solutions/ice/transient3d.m (modified) (2 diffs)

issm/trunk/src/m/solutions/ice/transient3d.m

@@ -39 +39 @@
 if ~isnan(md.vy), solution.u_g(2:6:fem.m_dh.gridset.gsize)=md.vy/md.yts; end
 if ~isnan(md.vz), solution.u_g(3:6:fem.m_dh.gridset.gsize)=md.vz/md.yts; end
-solution.pressure=zeros(fem.m_dh.gridset.gsize,1); solution.pressure(1:6:fem.m_dh.gridset.gsize)=md.rho_ice*md.g*(md.surface-md.z); %lithostatic pressure to spin up
+solution.pressure=zeros(fem.m_dh.gridset.gsize,1); 
+if ~isnan(md.pressure), solution.pressure(1:6:end)=md.pressure; end
 solution.t_g=zeros(fem.m_t.gridset.gsize,1);
 if ~isempty(md.temperature), 
@@ -153 +154 @@
         solution2(i).vz=solution(i).u_g(indz)*md.yts;
         solution2(i).vel=sqrt(solution2(i).vx.^2+solution2(i).vy.^2+solution2(i).vz.^2);
-        solution2(i).pressure=solution(i).pressure(indx)/10^5;
+        solution2(i).pressure=solution(i).pressure(indx);
         solution2(i).temperature=solution(i).t_g(indx);
         solution2(i).melting=solution(i).melting_g(indx);