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

tan2 at geodynamics.org tan2 at geodynamics.org
Thu Nov 8 15:35:33 PST 2007


Author: tan2
Date: 2007-11-08 15:35:32 -0800 (Thu, 08 Nov 2007)
New Revision: 8253

Modified:
   mc/3D/CitcomCU/trunk/src/Viscosity_structures.c
   mc/3D/CitcomCU/trunk/src/input1
Log:
Rollback accidental checkin

Modified: mc/3D/CitcomCU/trunk/src/Viscosity_structures.c
===================================================================
--- mc/3D/CitcomCU/trunk/src/Viscosity_structures.c	2007-11-08 23:33:04 UTC (rev 8252)
+++ mc/3D/CitcomCU/trunk/src/Viscosity_structures.c	2007-11-08 23:35:32 UTC (rev 8253)
@@ -345,28 +345,7 @@
 
 	}
 
-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:33:04 UTC (rev 8252)
+++ mc/3D/CitcomCU/trunk/src/input1	2007-11-08 23:35:32 UTC (rev 8253)
@@ -13,7 +13,7 @@
 
 	stokes_flow_only=0	# 1: only solve velocity once; 0: time-dependent problem.
 	maxstep=20000 		# max time steps
-	storage_spacing=10	# write data every ...
+	storage_spacing=50	# write data every ...
 
 # 2. Geometry, Ra numbers, Internal heating, Thermochemical/Purely thermal convection
 
@@ -34,16 +34,16 @@
 	
 # 3. Grid And Multiprocessor Information
 
-	nprocx=1		# number of processors in x or theta dir	
+	nprocx=2		# number of processors in x or theta dir	
 	nprocz=1		# number of processors in z or r dir
-	nprocy=1		# number of processors in y or fi dir
+	nprocy=2		# number of processors in y or fi dir
 
-	nodex=5  nodez=5  nodey=5 	# only relevant with conj-grad
+	nodex=33 nodez=33 nodey=33	# only relevant with conj-grad
 
-	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
+	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
 
 # 4. Coordinate Information
 
@@ -208,7 +208,7 @@
 	DESCRIBE=off		# 
 	BEGINNER=off		#
 	VERBOSE=off		#
-	verbose=on 		#
+	verbose=off		#
 	COMPRESS=off		#
 	see_convergence=1
 



More information about the cig-commits mailing list