[cig-commits] r8351 - mc/3D/CitcomCU/trunk/src

[tan2 at geodynamics.org]

Author: tan2
Date: 2007-11-28 17:40:04 -0800 (Wed, 28 Nov 2007)
New Revision: 8351

Modified:
   mc/3D/CitcomCU/trunk/src/Output.c
Log:
Write 3 components of velocity in velo files

Modified: mc/3D/CitcomCU/trunk/src/Output.c
===================================================================
--- mc/3D/CitcomCU/trunk/src/Output.c	2007-11-29 01:27:38 UTC (rev 8350)
+++ mc/3D/CitcomCU/trunk/src/Output.c	2007-11-29 01:40:04 UTC (rev 8351)
@@ -120,6 +120,13 @@
 		}
 		fclose(E->filed[10]);
 
+		sprintf(output_file, "%s.velo.%d.%d", E->control.data_file2, E->parallel.me, file_number);
+		E->filed[10] = fopen(output_file, "w");
+		fprintf(E->filed[10], "%6d %6d %.5e\n", E->lmesh.nno, E->advection.timesteps, E->monitor.elapsed_time);                                         
+		for(i = 1; i <= E->lmesh.nno; i++)
+			fprintf(E->filed[10], "%.6e %.6e %.6e\n", E->V[1][i], E->V[2][i], E->V[3][i]);
+		fclose(E->filed[10]);
+
 		if(E->parallel.me_loc[3] == E->parallel.nprocz - 1)
 		{
 			sprintf(output_file, "%s.th_t.%d.%d", E->control.data_file2, E->parallel.me, file_number);