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

[tan2 at geodynamics.org]

Author: tan2
Date: 2007-11-08 15:33:04 -0800 (Thu, 08 Nov 2007)
New Revision: 8252

Modified:
   mc/3D/CitcomCU/trunk/src/Pan_problem_misc_functions.c
   mc/3D/CitcomCU/trunk/src/Viscosity_structures.c
   mc/3D/CitcomCU/trunk/src/input1
Log:
Disable system call since some mpi implementation doesn't support it.

Modified: mc/3D/CitcomCU/trunk/src/Pan_problem_misc_functions.c
===================================================================
--- mc/3D/CitcomCU/trunk/src/Pan_problem_misc_functions.c	2007-11-08 23:32:23 UTC (rev 8251)
+++ mc/3D/CitcomCU/trunk/src/Pan_problem_misc_functions.c	2007-11-08 23:33:04 UTC (rev 8252)
@@ -66,7 +66,10 @@
 	{
 		sprintf(unique_name, "%06d.%s-%s", E->control.PID, comment, name);
 		sprintf(command, "cp -f %s %s\n", name, unique_name);
+#if 0
+        /* disable copying file, since some MPI implementation doesn't support system call. */
 		system(command);
+#endif
 	}
 
 	return;

Modified: mc/3D/CitcomCU/trunk/src/Viscosity_structures.c
===================================================================
--- mc/3D/CitcomCU/trunk/src/Viscosity_structures.c	2007-11-08 23:32:23 UTC (rev 8251)
+++ mc/3D/CitcomCU/trunk/src/Viscosity_structures.c	2007-11-08 23:33:04 UTC (rev 8252)
@@ -345,7 +345,28 @@
 
 	}
 
+for(i = 1; i <= E->lmesh.nel; i++) {
 
+    /* get temperature field on vpts */
+        for(jj = 1; jj <= vpts; jj++) {
+                TT[jj] = 0.0;
+                        for(kk = 1; kk <= ends; kk++)
+                                     TT[jj] += E->T[E->ien[i].node[kk]] * E->N.vpt[GNVINDEX(kk,jj)];
+                                         }
+
+                                             /* write to the log files */
+                                                 for(jj = 1; jj <= vpts; jj++) {
+                                                         fprintf(E->fp, "%d %d %e %e\n", i, jj,
+                                                                         E->EVI[E->mesh.levmax][(i - 1) * vpts + jj],
+                                                                                         TT[jj]);
+                                                                                             }
+                                                                                              
+                                                                                                
+                                                                                                }
+
+
+
+
 /*
 	fprintf(E->fp,"aaa\n");
 	for(i=1;i<=nel;i++)

Modified: mc/3D/CitcomCU/trunk/src/input1
===================================================================
--- mc/3D/CitcomCU/trunk/src/input1	2007-11-08 23:32:23 UTC (rev 8251)
+++ mc/3D/CitcomCU/trunk/src/input1	2007-11-08 23:33:04 UTC (rev 8252)
@@ -13,7 +13,7 @@
 
 	stokes_flow_only=0	# 1: only solve velocity once; 0: time-dependent problem.
 	maxstep=20000 		# max time steps
-	storage_spacing=50	# write data every ...
+	storage_spacing=10	# write data every ...
 
 # 2. Geometry, Ra numbers, Internal heating, Thermochemical/Purely thermal convection
 
@@ -34,16 +34,16 @@
 	
 # 3. Grid And Multiprocessor Information
 
-	nprocx=2		# number of processors in x or theta dir	
+	nprocx=1		# number of processors in x or theta dir	
 	nprocz=1		# number of processors in z or r dir
-	nprocy=2		# number of processors in y or fi dir
+	nprocy=1		# number of processors in y or fi dir
 
-	nodex=33 nodez=33 nodey=33	# only relevant with conj-grad
+	nodex=5  nodez=5  nodey=5 	# only relevant with conj-grad
 
-	mgunitx=6		# multigrid base level in x or t 
-	mgunitz=6		# multigrid base level in z or r
-	mgunity=6		# multigrid base level in y or f
-	levels=4		# and how many times it gets doubled
+	mgunitx=2		# multigrid base level in x or t 
+	mgunitz=2		# multigrid base level in z or r
+	mgunity=2		# multigrid base level in y or f
+	levels=2		# and how many times it gets doubled
 
 # 4. Coordinate Information
 
@@ -208,7 +208,7 @@
 	DESCRIBE=off		# 
 	BEGINNER=off		#
 	VERBOSE=off		#
-	verbose=off		#
+	verbose=on 		#
 	COMPRESS=off		#
 	see_convergence=1