Changeset 549

Timestamp:

05/21/09 08:42:46 (16 years ago)

Author:

Mathieu Morlighem

Message:

save transient results in thermaltransient_results field

File:

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

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

@@ -44 +44 @@
                 t_g=m_t.parameters.t_g;
                 melting_g=m_t.parameters.melting_g;
-
                 nsteps=m_t.parameters.ndt/m_t.parameters.dt;
 
@@ -50 +49 @@
                 soln.t_g=t_g;
                 soln.melting_g=melting_g;
+                soln.time=0;
 
                 for n=1:nsteps, 
 
                         displaystring(md.debug,'\n%s%i/%i\n','time step: ',n,nsteps);
+                        soln(n+1).time=n*m_t.parameters.dt;
                         
                         displaystring(md.debug,'\n%s',['    computing temperatures...']);
@@ -67 +68 @@
                 
                 %Wrap up
-                solution_temperature=struct('temperature',{});
+                solution=struct('time',{},'temperature',{},'melting',{});
                 for n=1:nsteps+1,
-                        solution_temperature(n).temperature=soln(n).t_g;
+                        solution(n).temperature=soln(n).t_g;
+                        solution(n).melting=soln(n).melting_g*md.yts; %in m/year
+                        solution(n).time=soln(n).time/md.yts;         %in years
                 end
-                md.temperature=solution_temperature;
-                
-                solution_melting=struct('melting',{});
-                for n=1:nsteps+1,
-                        solution_melting(n).melting=soln(n).melting_g*md.yts; %in m/a
-                end
-                md.melting=solution_melting;
+                md.thermaltransient_results=solution;
         end