From c1d741b8e1eb91dea84fdbe6e064ba92160dcc80 Mon Sep 17 00:00:00 2001 From: Paul Garlick Date: Mon, 25 Nov 2019 12:55:50 +0000 Subject: extrapolate topography from interior to fictitious cells. --- python/makeBoundary | 16 ++++++++-------- 1 file 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) -- cgit