extrapolate topography from interior to fictitious cells.

1 files changed, 8 insertions, 8 deletions

diff --git a/python/makeBoundary b/python/makeBoundary
index a0d3773..3684faf 100755
--- a/python/makeBoundary
+++ b/python/makeBoundary
@@ -183,20 +183,20 @@ print('dX =', dX)
 # blocks, with nrows lines per block.
 if args.location == 'top':
     xin = xtp[nrows-1:len(xtp):nrows]
-    yin = ytp[nrows-1:len(xtp):nrows]
-    zin = ztp[nrows-1:len(xtp):nrows]
+    yin = 2*ytp[nrows-1:len(xtp):nrows] - ytp[nrows-2:len(xtp):nrows]
+    zin = 2*ztp[nrows-1:len(xtp):nrows] - ztp[nrows-2:len(xtp):nrows]
 elif args.location == 'bottom':
     xin = xtp[0:len(xtp):nrows]
-    yin = ytp[0:len(xtp):nrows]
-    zin = ztp[0:len(xtp):nrows]
+    yin = 2*ytp[0:len(xtp):nrows] - ytp[1:len(xtp):nrows]
+    zin = 2*ztp[0:len(xtp):nrows] - ztp[1:len(xtp):nrows]
 elif args.location == 'left':
-    xin = xtp[:nrows]
+    xin = 2*xtp[:nrows] - xtp[nrows:2*nrows]
     yin = ytp[:nrows]
-    zin = ztp[:nrows]
+    zin = 2*ztp[:nrows] - ztp[nrows:2*nrows]
 elif args.location == 'right':
-    xin = xtp[nrows*(ncols-1):]
+    xin = 2*xtp[nrows*(ncols-1):] - xtp[nrows*(ncols-2):nrows*(ncols-1)]
     yin = ytp[nrows*(ncols-1):]
-    zin = ztp[nrows*(ncols-1):]
+    zin = 2*ztp[nrows*(ncols-1):] - ztp[nrows*(ncols-2):nrows*(ncols-1)]
     
 
 # print(xin)