Changeset 2084
- Timestamp:
- Nov 29, 2005, 10:13:34 AM (18 years ago)
- Location:
- production/gippsland_2005
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
production/gippsland_2005/create_mesh.py
r2076 r2084 10 10 from coordinate_transforms.geo_reference import Geo_reference 11 11 12 def create_mesh(max_area, mesh_file=None ):12 def create_mesh(max_area, mesh_file=None, triangles_in_name = False): 13 13 #Lower left corner 14 14 xllcorner = 570066 … … 52 52 m.addVertsSegs(dict) 53 53 54 m.generateMesh('q28.0z', maxArea = max_area) 54 m.generateMesh('q28.0z', maxArea = max_area) 55 triangle_count = len(m.getTriangulation()) 55 56 if mesh_file is None: 56 return m, len(m.getTriangulation())57 return m, triangle_count 57 58 else: 59 if triangles_in_name is True: 60 mesh_file = mesh_file[:-4] + '_' + str(triangle_count) \ 61 + mesh_file[-4:] 58 62 m.export_mesh_file(mesh_file) 59 return len(m.getTriangulation())63 return triangle_count 60 64 61 65 -
production/gippsland_2005/project.py
r2076 r2084 31 31 vcur_dir = boundarydir + 'vcur_expanded' 32 32 33 meshname = meshdir + basename + '. tsh'34 mesh_elevname = meshdir + basename + '_elev' + '. tsh'33 meshname = meshdir + basename + '.msh' 34 mesh_elevname = meshdir + basename + '_elev' + '.msh' 35 35 demname = datadir + basename + '.dem' 36 36 #demname = datadir + 'lakes_100' … … 70 70 refzone = 50 71 71 72 def update_names(new_info): 73 """ 74 Update the mesh_elevname and outputname, based on the string passed in. 75 Used to add the number of triangles to names. 76 77 """ 78 79 global mesh_elevname 80 global outputname 81 82 mesh_elevname = mesh_elevname[:-9] + '_' + str(new_info) + \ 83 mesh_elevname[-9:] 84 outputname = outputname[:-4] + '_' + str(new_info) + outputname[-4:] 85 86 #------------------------------------------------------------- 87 if __name__ == "__main__": 88 print "mesh_elevname",mesh_elevname 89 update_names(999) 90 print "mesh_elevname",mesh_elevname -
production/gippsland_2005/run_gippsland.py
r2082 r2084 51 51 52 52 #Create the mesh without elevation data 53 triagle_count = cache(create_mesh,(100000), {'mesh_file':project.meshname} 54 ,dependencies = ['create_mesh.py']) 53 # 100000 very course 54 triagle_count = cache(create_mesh,(5000), 55 {'mesh_file':project.meshname, 56 'triangles_in_name':True} 57 ,dependencies = ['create_mesh.py'] 58 #,evaluate = True 59 ) 60 print "project.meshname",project.meshname 61 project.update_names(triagle_count) 62 print "project.meshname",project.meshname 55 63 56 64 #Add elevation data to the mesh … … 65 73 'precrop': True} 66 74 ,dependencies = [project.meshname, project.pointname] 75 ,evaluate = True 67 76 ,verbose = False 68 77 ) … … 76 85 77 86 78 domain.set_name(project.basename ) # + '_%d' %resolution)87 domain.set_name(project.basename + '_%d' %triagle_count) 79 88 domain.set_datadir(project.outputdir) 80 89 domain.store = True
Note: See TracChangeset
for help on using the changeset viewer.