aboutsummaryrefslogtreecommitdiff
path: root/python/makeBoundary
diff options
context:
space:
mode:
authorPaul Garlick <pgarlick@tourbillion-technology.com>2019-11-25 12:55:50 +0000
committerPaul Garlick <pgarlick@tourbillion-technology.com>2019-11-25 12:55:50 +0000
commitc1d741b8e1eb91dea84fdbe6e064ba92160dcc80 (patch)
tree07f28b508297370251bd7bd3d574493099a6e87d /python/makeBoundary
parentce8f415da2cc27348a6deb83d090d974d9162a06 (diff)
downloadfullSWOF-utils-c1d741b8e1eb91dea84fdbe6e064ba92160dcc80.tar.gz
extrapolate topography from interior to fictitious cells.
Diffstat (limited to 'python/makeBoundary')
-rwxr-xr-xpython/makeBoundary16
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)