Changeset 955

Timestamp:

06/12/09 14:42:39 (16 years ago)

Author:

Mathieu Morlighem

Message:

fixed load conservation

File:

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

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

@@ -72 +72 @@
                 gridset=m_t.gridset;
                 inputs=struct('pressure',pressure,'temperature',soln(n).t_g,'velocity',velocity,'dt',md.dt);
-                [soln(n+1).t_g m_t.loads_t  melting_offset]=thermal_core(m_t,analysis_type,inputs);
+                [soln(n+1).t_g m_t.loads  melting_offset]=thermal_core(m_t,analysis_type,inputs);
 
                 %Call core melting computation
@@ -79 +79 @@
                 inputs=struct('pressure',pressure,'temperature',soln(n+1).t_g,'melting_offset',melting_offset,'dt',md.dt);
                 soln(n+1).m_g=melting_core(m_m,analysis_type,inputs);
-                any(soln(n+1).m_g)
                 
         end