<p><b>ringler@lanl.gov</b> 2011-10-14 09:20:08 -0600 (Fri, 14 Oct 2011)</p><p><br>
adding connections needed to eliminate loops, currently turned off<br>
</p><hr noshade><pre><font color="gray">Modified: branches/ocean_projects/basin/src/Makefile
===================================================================
--- branches/ocean_projects/basin/src/Makefile        2011-10-14 15:12:55 UTC (rev 1074)
+++ branches/ocean_projects/basin/src/Makefile        2011-10-14 15:20:08 UTC (rev 1075)
@@ -46,11 +46,12 @@
        module_read_netcdf.o \
        module_read_topo.o \
        module_read_TS.o \
+       module_cullLoops.o \
        module_write_netcdf.o
 
 all: map
 
-basin.o: utilities.o module_write_netcdf.o module_read_netcdf.o module_read_topo.o module_read_TS.o
+basin.o: utilities.o module_write_netcdf.o module_read_netcdf.o module_read_topo.o module_read_TS.o module_cullLoops.o
 
 map: $(OBJS)
         $(FC) $(LDFLAGS) -o $@ $(OBJS) $(LIBS)

Modified: branches/ocean_projects/basin/src/basin.F
===================================================================
--- branches/ocean_projects/basin/src/basin.F        2011-10-14 15:12:55 UTC (rev 1074)
+++ branches/ocean_projects/basin/src/basin.F        2011-10-14 15:20:08 UTC (rev 1075)
@@ -5,6 +5,7 @@
 use read_TS
 use write_netcdf
 use utilities
+!use cullLoops
 
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 ! Program: basin.F
@@ -1107,13 +1108,13 @@
         where(yCell.lt.1.001*ymin) kmt = 0
         where(yCell.gt.0.999*ymax) kmt = 0
 
-        ! solid boundary in x
-        xmin = minval(xCell)
-        write(6,*) ' minimum xCell ', xmin
-        xmax = maxval(xCell)
-        write(6,*) ' maximum xCell ', xmax
-        where(xCell.lt.xmin+dc/1.5) kmt = 0
-        where(xCell.gt.xmax-dc/1.5) kmt = 0
+     !  ! solid boundary in x
+     !  xmin = minval(xCell)
+     !  write(6,*) ' minimum xCell ', xmin
+     !  xmax = maxval(xCell)
+     !  write(6,*) ' maximum xCell ', xmax
+     !  where(xCell.lt.xmin+dc/1.5) kmt = 0
+     !  where(xCell.gt.xmax-dc/1.5) kmt = 0
     endif
 
     
@@ -1192,8 +1193,14 @@
     if(flag) kmt(iCell)=0
 enddo
 
+!call eliminateLoops(nCells,nEdges,nVertices,maxEdges,vertexDegree, &amp;
+                    !nEdgesOnCell, cellsOnCell, verticesOnEdge, cellsOnVertex, edgesOnCell, lonCell, latCell, &amp;
+                    !xCell, yCell, zCell, xEdge, yEdge, zEdge, xVertex, yVertex, zVertex, &amp;
+                    !KMT)
+
 if(real_bathymetry) then
-    where(kmt.eq.1) kmt=2
+    where(kmt.eq.1) kmt=3
+    where(kmt.eq.2) kmt=3
 endif
 
 end subroutine define_kmt

Modified: branches/ocean_projects/basin/src/utilities.F
===================================================================
--- branches/ocean_projects/basin/src/utilities.F        2011-10-14 15:12:55 UTC (rev 1074)
+++ branches/ocean_projects/basin/src/utilities.F        2011-10-14 15:20:08 UTC (rev 1075)
@@ -232,7 +232,7 @@
  20   format(e20.10)
  21   format(i20)
  22   format(3e20.10)
- 23   format(3i6, 3e20.10)
+ 23   format(3i8, 3e20.10)
 
       close(10)
       close(11)

</font>
</pre>