[cig-commits] r5645 - in long/3D/Gale/trunk/src/Gale: . plugins/SurfaceProcess

[walter at geodynamics.org]

Author: walter
Date: 2007-01-05 11:19:07 -0800 (Fri, 05 Jan 2007)
New Revision: 5645

Modified:
   long/3D/Gale/trunk/src/Gale/
   long/3D/Gale/trunk/src/Gale/plugins/SurfaceProcess/SurfaceProcess.c
Log:
 r96 at earth:  boo | 2007-01-05 11:15:05 -0800
  r88 at earth (orig r73):  LukeHodkinson | 2007-01-04 20:17:30 -0800
  Found an index mapping bug.
  
  
 



Property changes on: long/3D/Gale/trunk/src/Gale
___________________________________________________________________
Name: svk:merge
   - 4e4aea6e-fd02-0410-981f-80ab108a5659:/trunk:72
8f887497-cf10-0410-afc2-df8ae19c7fb0:/cig:95
   + 4e4aea6e-fd02-0410-981f-80ab108a5659:/trunk:73
8f887497-cf10-0410-afc2-df8ae19c7fb0:/cig:96

Modified: long/3D/Gale/trunk/src/Gale/plugins/SurfaceProcess/SurfaceProcess.c
===================================================================
--- long/3D/Gale/trunk/src/Gale/plugins/SurfaceProcess/SurfaceProcess.c	2007-01-05 19:19:06 UTC (rev 5644)
+++ long/3D/Gale/trunk/src/Gale/plugins/SurfaceProcess/SurfaceProcess.c	2007-01-05 19:19:07 UTC (rev 5645)
@@ -244,8 +244,10 @@
 		if( vDomain >= Mesh_GetLocalSize( spCtx->mesh, MT_VERTEX ) )
 			continue;
 
-		vGlobal = Mesh_DomainToGlobal( spCtx->mesh, MT_VERTEX, v_i );
+		vGlobal = Mesh_DomainToGlobal( spCtx->mesh, MT_VERTEX, vDomain );
 		RegularMeshUtils_Node_1DTo3D( spCtx->mesh, vGlobal, ijk );
+		if( nSurfDims == 2 )
+			ijk[1] = ijk[2];
 		surfaceInd = Grid_Project( surfaceGrid, ijk );
 		UIntMap_Insert( map, topVerts[v_i], surfaceInd );
 	}