from src.utilities import DataInput, MeshGenerator from data.data_factory import DataFactory # create meshgrid for contour : vara = DataFactory.get_bamber() H = vara['H']['map_data'] S = vara['h']['map_data'] B = vara['b']['map_data'] H[S <= 0] = 0.0 vara['H']['map_data'] = H # process the data : dd = DataInput(None, vara, flip=False, gen_space=False) #dd.set_data_max('mask', 2, 0) m = MeshGenerator(dd, 'mesh', './') # 20 works well for the last arg. below m.create_contour('H', 200.0, 3) m.plot_contour() m.eliminate_intersections() m.plot_contour() # 10000 works well on the following line m.write_gmsh_contour(1000) m.finish(None) m.create_2D_mesh('mesh') #FIXME: fails m.convert_msh_to_xml('mesh', 'mesh')
from src.utilities import DataInput, MeshGenerator from data.data_factory import DataFactory # create meshgrid for contour : vara = DataFactory.get_bamber() H = vara['H']['map_data'] S = vara['h']['map_data'] B = vara['b']['map_data'] H[S<=0] = 0.0 vara['H']['map_data'] = H # process the data : dd = DataInput(None, vara, flip=False, gen_space=False) #dd.set_data_max('mask', 2, 0) m = MeshGenerator(dd, 'mesh', './') # 20 works well for the last arg. below m.create_contour('H', 200.0,3) m.plot_contour() m.eliminate_intersections() m.plot_contour() # 10000 works well on the following line m.write_gmsh_contour(1000) m.finish(None) m.create_2D_mesh('mesh') #FIXME: fails m.convert_msh_to_xml('mesh', 'mesh')
import sys import os src_directory = '../../../' sys.path.append(src_directory) from src.utilities import DataInput, MeshGenerator from data.data_factory import DataFactory # create meshgrid for contour : vara = DataFactory.get_bedmap2() # process the data : dd = DataInput(None, vara, flip=True, gen_space=False) #dd.set_data_max('mask', 2, 0) m = MeshGenerator(dd, 'mesh', '../meshes/') m.create_contour('mask', 0.999, 100) m.plot_contour() m.eliminate_intersections() m.write_gmsh_contour(10000.) m.finish()