<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>