diff options
| -rwxr-xr-x | makeBoundary | 5 |
1 files changed, 1 insertions, 4 deletions
diff --git a/makeBoundary b/makeBoundary index b16512e..6695404 100755 --- a/makeBoundary +++ b/makeBoundary @@ -158,6 +158,7 @@ numH = definition_dict["numH"] # number of height intervals # m, c = np.polyfit(ych, zch, 1) # print('gradient =', m, 'intercept =', c) +# read topography: with open(height_data, "r") as topo: xtp, ytp, ztp = np.loadtxt(topo, delimiter=' ', unpack=True) @@ -177,22 +178,18 @@ 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] - outputFilename = "BCTop.txt" elif args.location == 'bottom': xin = xtp[0:len(xtp):nrows] yin = ytp[0:len(xtp):nrows] zin = ztp[0:len(xtp):nrows] - outputFilename = "BCBottom.txt" elif args.location == 'left': xin = xtp[:nrows] yin = ytp[:nrows] zin = ztp[:nrows] - outputFilename = "BCLeft.txt" elif args.location == 'right': xin = xtp[nrows*(ncols-1):] yin = ytp[nrows*(ncols-1):] zin = ztp[nrows*(ncols-1):] - outputFilename = "BCRight.txt" # print(xin) |