Changeset 893

Timestamp:

06/11/09 08:47:26 (15 years ago)

Author:

seroussi

Message:

transient 3d serial fixed

File:

• issm/trunk/src/m/solutions/cielo/transient3d.m (modified) (3 diffs)

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

@@ -78 +78 @@
         inputs=add(inputs,'bed',solution(n).b_g,'doublevec',1,m_p.parameters.numberofnodes);
         inputs=add(inputs,'velocity',solution(n).u_g,'doublevec',3,m_p.parameters.numberofnodes);
-        %inputs=add(inputs,'velocity',solution(n).u_g,'doublevec',m_p.parameters.numberofdofspernode,m_p.parameters.numberofnodes);
         inputs=add(inputs,'pressure',solution(n).p_g,'doublevec',1,m_p.parameters.numberofnodes);
         inputs=add(inputs,'temperature',solution(n).t_g,'doublevec',1,m_t.parameters.numberofnodes);
@@ -85 +84 @@
         displaystring(md.debug,'\n%s',['    computing temperatures...']);
         [solution(n+1).t_g m_t.loads melting_offset]=thermal_core_nonlinear(m_t,inputs,'thermal','transient');
+        inputs=add(inputs,'temperature',solution(n+1).t_g,'doublevec',1,m_t.parameters.numberofnodes);
         
         displaystring(md.debug,'\n%s',['    computing melting...']);
@@ -142 +142 @@
         results(i).vz=solution(i).u_g(3:3:end)*yts;
         results(i).vel=sqrt(solution(i).u_g(1:3:end).^2+solution(i).u_g(2:3:end).^2+solution(i).u_g(3:3:end).^2)*yts;
+        results(i).pressure=solution(i).p_g;
         results(i).bed=solution(i).b_g;
         results(i).surface=solution(i).s_g;