[cig-commits] commit: Make HRS_Erosion not skip the edges of the grid

[Mercurial]

changeset:   829:3bb9aae4d6a1
tag:         tip
user:        Walter Landry <wlandry at caltech.edu>
date:        Tue Mar 16 16:03:48 2010 -0700
files:       plugins/HRS_Erosion/HRS_Erosion.c
description:
Make HRS_Erosion not skip the edges of the grid


diff -r bd95c8514b39 -r 3bb9aae4d6a1 plugins/HRS_Erosion/HRS_Erosion.c
--- a/plugins/HRS_Erosion/HRS_Erosion.c	Fri Mar 12 15:48:39 2010 -0800
+++ b/plugins/HRS_Erosion/HRS_Erosion.c	Tue Mar 16 16:03:48 2010 -0700
@@ -157,7 +157,7 @@ void Underworld_HRS_Erosion_Execute( Tim
 
       if(!Mesh_GlobalToDomain
          (mesh,MT_VERTEX,
-          RegularMeshUtils_Node_3DTo1D(mesh,ijk_left),&n_left)) 
+          RegularMeshUtils_Node_3DTo1D(mesh,ijk_left),&n_left))
         {
           printf("Can not map to local domain %d %d\n",
                  ijk_left[0],ijk_left[1]);
@@ -170,7 +170,7 @@ void Underworld_HRS_Erosion_Execute( Tim
           ijk_left[0]=i;
           if(!Mesh_GlobalToDomain
              (mesh,MT_VERTEX,
-              RegularMeshUtils_Node_3DTo1D(mesh,ijk_left),&n_left)) 
+              RegularMeshUtils_Node_3DTo1D(mesh,ijk_left),&n_left))
             {
               printf("Can not map to local domain %d %d\n",
                      ijk_left[0],ijk_left[1]);
@@ -238,8 +238,7 @@ void Underworld_HRS_Erosion_Execute( Tim
           RegularMeshUtils_Node_1DTo3D
             ( mesh, Mesh_DomainToGlobal( mesh, MT_VERTEX, n_i ), ijk );
 					
-          if(ijk[1]==grid->sizes[1]-1
-             && ijk[0]!=0 && ijk[0]!=grid->sizes[0]-1) 
+          if(ijk[1]==grid->sizes[1]-1)
             {
               double y_old,delta_v, v[3];
               double x;