예제 #1
0
import sys
from grid_3d import grid

if len(sys.argv) < 6:
        print("Usage:\t", sys.argv[0] , "<type> <output grid> <i_max> <j_max> <k_max>  [x_length] [y_length] [z_length]\n\t", "where <type> is one of: driven, step", file=sys.stderr)
        sys.exit(0)

i_max = int(sys.argv[3])
j_max = int(sys.argv[4])
k_max = int(sys.argv[5])

x_length = sys.argv[6] if len(sys.argv) > 6 else 1
y_length = sys.argv[7] if len(sys.argv) > 7 else 1
z_length = sys.argv[8] if len(sys.argv) > 8 else 1

type_string = sys.argv[1]
outfile_path = sys.argv[2]

b = grid(i_max, j_max, k_max, x_length, y_length, z_length)

if type_string.lower() == "driven":
        print("Creating 'Driven Cavity' grid at", outfile_path, ".")
elif type_string.lower() == "step":
        print("Creating 'Flow Over A Step' grid at", outfile_path, ".")
        b.insert_rectangle(0, 1 + i_max/4, 0, j_max + 1, 0, 1 + k_max/2)
else:
        print("Type not supported!", file=sys.stderr)
        sys.exit(0)
        
b.write_to(outfile_path)
예제 #2
0
in_grid = []

with open(infile_path, "rb") as csvfile:
    gridreader = csv.reader(csvfile, delimiter=",", quotechar='"')

    for row in gridreader:
        in_grid.append(row)


i_max = len(in_grid[0]) - 2
j_max = len(in_grid) - 2
k_max = 2 * distance + (num_repetitions - 1) * distance + height * num_repetitions

x_length = sys.argv[6] if len(sys.argv) > 6 else 1
y_length = sys.argv[7] if len(sys.argv) > 7 else 1
z_length = sys.argv[8] if len(sys.argv) > 8 else 1

b = grid(i_max, j_max, k_max, x_length, y_length, z_length, all_fluid=False)

k = 1 + distance

for _ in range(num_repetitions):
    for _ in range(height):
        b.add_xy_plane(k, in_grid)
        k += 1

    k += distance

print("Writing grid of size", i_max, j_max, k_max, "to", outfile_path, ".")
b.write_to(outfile_path)
예제 #3
0
in_grid = []

with open(infile_path, 'rb') as csvfile:
    gridreader = csv.reader(csvfile, delimiter=',', quotechar='"')

    for row in gridreader:
        in_grid.append(row)

i_max = len(in_grid[0]) - 2
j_max = len(in_grid) - 2
k_max = 2 * distance + (num_repetitions -
                        1) * distance + height * num_repetitions

x_length = sys.argv[6] if len(sys.argv) > 6 else 1
y_length = sys.argv[7] if len(sys.argv) > 7 else 1
z_length = sys.argv[8] if len(sys.argv) > 8 else 1

b = grid(i_max, j_max, k_max, x_length, y_length, z_length, all_fluid=False)

k = 1 + distance

for _ in range(num_repetitions):
    for _ in range(height):
        b.add_xy_plane(k, in_grid)
        k += 1

    k += distance

print("Writing grid of size", i_max, j_max, k_max, "to", outfile_path, ".")
b.write_to(outfile_path)
예제 #4
0
파일: make_3d.py 프로젝트: ltroska/nast_hpx
import sys
from grid_3d import grid

if len(sys.argv) < 6:
	print("Usage:\t", sys.argv[0] , "<type> <output grid> <i_max> <j_max> <k_max>  [x_length] [y_length] [z_length]\n\t", "where <type> is one of: driven, step", file=sys.stderr)
	sys.exit(0)

i_max = int(sys.argv[3])
j_max = int(sys.argv[4])
k_max = int(sys.argv[5])

x_length = sys.argv[6] if len(sys.argv) > 6 else 1
y_length = sys.argv[7] if len(sys.argv) > 7 else 1
z_length = sys.argv[8] if len(sys.argv) > 8 else 1

type_string = sys.argv[1]
outfile_path = sys.argv[2]

b = grid(i_max, j_max, k_max, x_length, y_length, z_length)

if type_string.lower() == "driven":
	print("Creating 'Driven Cavity' grid at", outfile_path, ".")
elif type_string.lower() == "step":
	print("Creating 'Flow Over A Step' grid at", outfile_path, ".")
	b.insert_rectangle(0, 1 + i_max/4, 0, j_max + 1, 0, 1 + k_max/2)
else:
	print("Type not supported!", file=sys.stderr)
	sys.exit(0)
								
b.write_to(outfile_path)