<p><b>mpetersen@lanl.gov</b> 2011-07-13 07:58:10 -0600 (Wed, 13 Jul 2011)</p><p>Add lines to have solid boundaries in both x and y in the basin code.  This is needed for the double gyre simulation.<br>
</p><hr noshade><pre><font color="gray">Modified: branches/ocean_projects/basin/src/basin.F
===================================================================
--- branches/ocean_projects/basin/src/basin.F        2011-07-12 16:24:08 UTC (rev 917)
+++ branches/ocean_projects/basin/src/basin.F        2011-07-13 13:58:10 UTC (rev 918)
@@ -1072,7 +1072,7 @@
 real (kind=4), allocatable, dimension(:) :: x,y, work_kmt
 real (kind=4), allocatable, dimension(:,:) :: ztopo
 integer :: nx, ny, inx, iny, ix, iy
-real :: pi, dtr, zdata, rlon, rlat, r, ymin, ymax
+real :: pi, dtr, zdata, rlon, rlat, r, ymin, ymax, xmin, xmax
 real :: latmin, latmax, lonmin, lonmax
 logical :: flag, kmt_flag
 pi = 4.0*atan(1.0)
@@ -1096,12 +1096,21 @@
         where(lonCell.lt.lonmin) kmt = 0
         where(lonCell.gt.lonmax) kmt = 0
     else
+        ! solid boundary in y
         ymin = minval(yCell)
         write(6,*) ' minimum yCell ', ymin
         ymax = maxval(yCell)
         write(6,*) ' maximum yCell ', ymax
         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.1.001*xmin) kmt = 0
+        where(xCell.gt.0.999*xmax) kmt = 0
     endif
 
     

</font>
</pre>