Changeset 579

Timestamp:

05/21/09 15:01:36 (16 years ago)

Author:

Mathieu Morlighem

Message:

some fixing

File:

• issm/trunk/src/c/parallel/thermal.cpp (modified) (3 diffs)

issm/trunk/src/c/parallel/thermal.cpp

@@ -120 +120 @@
                 nsteps=(int)(ndt/dt);
                 time=(double*)xmalloc((nsteps+1)*sizeof(double));
+                time[0]=0;
 
                 /*allocate t_g and m_g arrays: */
@@ -136 +137 @@
 
                 for(i=0;i<nsteps;i++){
-                        if(debug)_printf_("time step: %i/%i\n",n,nsteps);
-                        time[i]=n*dt;
+                        if(debug)_printf_("time step: %i/%i\n",i+1,nsteps);
+                        time[i+1]=(i+1)*dt;
                         
                         if(debug)_printf_("computing temperatures:\n");
@@ -146 +147 @@
                         inputs->Add("temperature",t_g[i+1],1,numberofnodes);
                         inputs->Add("melting_offset",melting_offset);
-                        diagnostic_core_linear(&m_g[0],&femmodels[1],inputs,MeltingAnalysisEnum(),TransientAnalysisEnum());
+                        diagnostic_core_linear(&m_g[i+1],&femmodels[1],inputs,MeltingAnalysisEnum(),TransientAnalysisEnum());
                 }
         }