Changeset 3535
- Timestamp:
- Aug 29, 2006, 11:30:14 AM (19 years ago)
- Files:
-
- 52 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_core/source/anuga/examples/run_gong.py
r951 r3535 20 20 Transmissive_boundary, Time_boundary, Add_value_to_region, File_boundary 21 21 from mesh_factory import rectangular 22 from pmesh2domain import pmesh_to_domain, pmesh_to_domain_instance22 from anuga.pyvolution.pmesh2domain import pmesh_to_domain, pmesh_to_domain_instance 23 23 24 24 from Numeric import array -
anuga_core/source/anuga/examples/run_tsh.py
r3514 r3535 21 21 from region import Add_value_to_region, Set_region 22 22 from mesh_factory import rectangular 23 from pmesh2domain import pmesh_to_domain_instance23 from anuga.pyvolution.pmesh2domain import pmesh_to_domain_instance 24 24 25 25 from Numeric import array -
anuga_core/source/anuga/examples/runup_convergence.py
r3530 r3535 14 14 import sys 15 15 16 from pmesh.mesh_interface import create_mesh_from_regions16 from anuga.pmesh.mesh_interface import create_mesh_from_regions 17 17 from pyvolution.mesh_factory import rectangular_cross 18 18 from pyvolution.shallow_water import Domain -
anuga_core/source/anuga/pmesh/timing.py
r3437 r3535 3 3 """ 4 4 from mesh import * 5 from pmesh import *5 from anuga.pmesh import * 6 6 import time 7 7 -
anuga_core/source/anuga/pyvolution/check_sww_tsh.py
r3514 r3535 8 8 from load_mesh.loadASCII import import_mesh_file 9 9 from shallow_water import Domain 10 from pmesh2domain import pmesh_to_domain_instance10 from anuga.pyvolution.pmesh2domain import pmesh_to_domain_instance 11 11 12 12 def check_sww_tsh(sww_file, tsh_file, verbose = False): -
anuga_core/source/anuga/pyvolution/data_manager.py
r3532 r3535 31 31 Manually created files: 32 32 ASC, PRJ: Digital elevation models (gridded) 33 TSH: Triangular meshes (e.g. created from pmesh)33 TSH: Triangular meshes (e.g. created from anuga.pmesh) 34 34 NC Model outputs for use as boundary conditions (e.g from MOST) 35 35 … … 3545 3545 3546 3546 from shallow_water import Domain 3547 from pmesh2domain import pmesh_to_domain_instance3547 from anuga.pyvolution.pmesh2domain import pmesh_to_domain_instance 3548 3548 import time, os 3549 3549 from anuga.pyvolution.data_manager import get_dataobject -
anuga_core/source/anuga/pyvolution/run_tsh_weir.py
r773 r3535 16 16 17 17 from mesh_factory import rectangular 18 from pmesh2domain import pmesh_to_domain18 from anuga.pyvolution.pmesh2domain import pmesh_to_domain 19 19 20 20 from Numeric import array -
anuga_core/source/anuga/pyvolution/test_pmesh2domain.py
r3514 r3535 8 8 from shallow_water import Domain 9 9 10 #from anuga.pyvolution.pmesh2domain import * 11 10 12 from pmesh2domain import * 11 13 … … 14 16 15 17 from anuga.coordinate_transforms.geo_reference import Geo_reference 16 17 #This is making pyvolution dependent on pmesh.18 # not good. this should be in a seperate package.(directory)19 18 from anuga.pmesh.mesh import importMeshFromFile 20 19 -
anuga_core/source/anuga/pyvolution/test_shallow_water.py
r3514 r3535 3341 3341 tags["3"] = b3 3342 3342 3343 #from pmesh2domain import pmesh_to_domain_instance3343 #from anuga.pyvolution.pmesh2domain import pmesh_to_domain_instance 3344 3344 #domain = pmesh_to_domain_instance(fileName, Domain) 3345 3345 -
anuga_core/source/anuga/pyvolution/tsh2sww.py
r3514 r3535 11 11 12 12 from shallow_water import Domain 13 from pmesh2domain import pmesh_to_domain_instance13 from anuga.pyvolution.pmesh2domain import pmesh_to_domain_instance 14 14 import time, os 15 15 from anuga.pyvolution.data_manager import get_dataobject -
anuga_core/source/anuga/pyvolution/view_tsh.py
r3514 r3535 8 8 9 9 from shallow_water import Domain 10 from pmesh2domain import pmesh_to_domain_instance10 from anuga.pyvolution.pmesh2domain import pmesh_to_domain_instance 11 11 from anuga.pyvolution.util import file_function 12 12 from anuga.utilities.polygon import Polygon_function, read_polygon -
development/cairns_2006/create_mesh.py
r3434 r3535 6 6 #!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 7 7 8 from pmesh.mesh import *9 from coordinate_transforms.geo_reference import Geo_reference8 from anuga.pmesh.mesh import * 9 from anuga.coordinate_transforms.geo_reference import Geo_reference 10 10 from conversion import convert_latlon_to_xycoords 11 11 -
development/dam_2006/create_mesh.py
r3446 r3535 2 2 """ 3 3 4 #!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!5 # Assume that the root AnuGA dir is included in your pythonpath6 #!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!7 4 8 from pmesh.mesh import *9 from coordinate_transforms.geo_reference import Geo_reference5 from anuga.pmesh.mesh import * 6 from anuga.coordinate_transforms.geo_reference import Geo_reference 10 7 11 8 -
development/dam_2006/run_dam.py
r3533 r3535 26 26 # Application specific imports 27 27 import project # Definition of file names and polygons 28 28 29 29 30 def main(): … … 71 72 domain.set_quantities_to_be_stored(['stage', 'xmomentum', 'ymomentum']) 72 73 domain.set_minimum_sww_depth(0.01) 74 domain.set_store_vertices_uniquely(True) # for writting to sww 73 75 74 76 #------------------------------------------------------------------------- … … 88 90 domain.set_region('dam','stage',0.2, 89 91 location = 'unique vertices') 90 domain.set_region(Set_region('dam','stage',0.2, 91 location = 'unique vertices')) 92 domain.set_store_vertices_uniquely(True) # for writting to sww 92 #domain.set_region(Set_region('dam','stage',0.2, 93 # location = 'unique vertices')) 93 94 94 95 Br = Reflective_boundary(domain) 95 Bd = Dirichlet_boundary([0,0,0]) 96 Bd = Dirichlet_boundary([0,0,0]) # to drain the water out. 96 97 domain.set_boundary( {'wall': Br, 'edge': Bd} ) 97 98 … … 116 117 117 118 #------------------------------------------------------------------------- 118 # Calc gauge info119 # Calculate gauge info 119 120 #------------------------------------------------------------------------- 120 121 interpolate_sww2csv(project.outputtimedir + project.basename +".sww", -
development/momentum_sink/scripts/building_generator.py
r2927 r3535 3 3 4 4 # add inundation dir to your pythonpath 5 from pmesh.mesh import Mesh6 from coordinate_transforms.geo_reference import Geo_reference5 from anuga.pmesh.mesh import Mesh 6 from anuga.coordinate_transforms.geo_reference import Geo_reference 7 7 8 8 -
development/momentum_sink/scripts/create_buildings.py
r2927 r3535 6 6 import array 7 7 # add inundation dir to your pythonpath 8 from pmesh.mesh import Mesh9 from coordinate_transforms.geo_reference import Geo_reference8 from anuga.pmesh.mesh import Mesh 9 from anuga.coordinate_transforms.geo_reference import Geo_reference 10 10 11 11 WidtH = 200 # width of boudary in metres -
development/momentum_sink/scripts/create_two_buildings.py
r2927 r3535 4 4 import random 5 5 # add inundation dir to your pythonpath 6 from pmesh.mesh import Mesh7 from coordinate_transforms.geo_reference import Geo_reference6 from anuga.pmesh.mesh import Mesh 7 from anuga.coordinate_transforms.geo_reference import Geo_reference 8 8 9 9 WidtH = 200 # width of boudary in metres -
development/momentum_sink/scripts/friction_block.py
r2927 r3535 3 3 4 4 # add inundation dir to your pythonpath 5 from pmesh.mesh import Mesh6 from coordinate_transforms.geo_reference import Geo_reference5 from anuga.pmesh.mesh import Mesh 6 from anuga.coordinate_transforms.geo_reference import Geo_reference 7 7 8 8 -
development/momentum_sink/scripts/loop_buildings.py
r3514 r3535 15 15 import project 16 16 from create_buildings import create_mesh 17 from pmesh.mesh import importMeshFromFile17 from anuga.pmesh.mesh import importMeshFromFile 18 18 import Numeric 19 19 -
development/momentum_sink/scripts/loop_friction.py
r3514 r3535 18 18 from anuga.pyvolution.least_squares import fit_to_mesh_file, DEFAULT_ALPHA 19 19 from friction_block import create_mesh 20 from pmesh.mesh import importMeshFromFile20 from anuga.pmesh.mesh import importMeshFromFile 21 21 22 22 # generation of friction value list -
development/momentum_sink/scripts/plot_cross_section.py
r3514 r3535 7 7 from anuga.pyvolution.util import file_function 8 8 #from anuga.pyvolution.coordinate_transforms.redfearn import degminsec2decimal_degrees, redfearn 9 #from coordinate_transforms.redfearn import degminsec2decimal_degrees, redfearn9 #from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees, redfearn 10 10 from pylab import * 11 11 #from compare_sww import gauge_locations -
development/momentum_sink/scripts/run_buildings.py
r3514 r3535 18 18 #from building_generator import create_mesh 19 19 20 from pmesh.mesh import importMeshFromFile20 from anuga.pmesh.mesh import importMeshFromFile 21 21 22 22 depth = 3.5 # depth of building side to oncoming wave -
development/momentum_sink/scripts/run_friction.py
r3514 r3535 25 25 #from building_generator import create_mesh 26 26 27 from pmesh.mesh import importMeshFromFile27 from anuga.pmesh.mesh import importMeshFromFile 28 28 29 29 meshname = project_friction.meshname -
development/momentum_sink/scripts/two_friction_block.py
r2927 r3535 3 3 4 4 # add inundation dir to your pythonpath 5 from pmesh.mesh import Mesh6 from coordinate_transforms.geo_reference import Geo_reference5 from anuga.pmesh.mesh import Mesh 6 from anuga.coordinate_transforms.geo_reference import Geo_reference 7 7 8 8 -
development/okushiri_2005/create_mesh.py
r2427 r3535 6 6 #!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 7 7 8 from pmesh.mesh import *9 from coordinate_transforms.geo_reference import Geo_reference8 from anuga.pmesh.mesh import * 9 from anuga.coordinate_transforms.geo_reference import Geo_reference 10 10 11 11 #------------------------------------------------------------- -
development/pyvolution-1d/data_manager.py
r3514 r3535 54 54 55 55 56 from coordinate_transforms.geo_reference import Geo_reference56 from anuga.coordinate_transforms.geo_reference import Geo_reference 57 57 58 58 def make_filename(s): … … 1053 1053 1054 1054 #Georeferencing 1055 from coordinate_transforms.geo_reference import Geo_reference1055 from anuga.coordinate_transforms.geo_reference import Geo_reference 1056 1056 if zone is None: 1057 1057 assert xllcorner is None, 'xllcorner must be None' … … 2309 2309 2310 2310 #Store 2311 from coordinate_transforms.redfearn import redfearn2311 from anuga.coordinate_transforms.redfearn import redfearn 2312 2312 x = zeros(number_of_points, Float) #Easting 2313 2313 y = zeros(number_of_points, Float) #Northing … … 3439 3439 from Scientific.IO.NetCDF import NetCDFFile 3440 3440 3441 from coordinate_transforms.redfearn import redfearn3441 from anuga.coordinate_transforms.redfearn import redfearn 3442 3442 3443 3443 precision = Float # So if we want to change the precision its done here … … 3591 3591 3592 3592 #Store 3593 from coordinate_transforms.redfearn import redfearn3593 from anuga.coordinate_transforms.redfearn import redfearn 3594 3594 x = zeros(number_of_points, Float) #Easting 3595 3595 y = zeros(number_of_points, Float) #Northing -
development/pyvolution-1d/domain.py
r3510 r3535 8 8 """ 9 9 from generic_boundary_conditions import * 10 from coordinate_transforms.geo_reference import Geo_reference10 from anuga.coordinate_transforms.geo_reference import Geo_reference 11 11 12 12 class Domain: -
development/stochastic_study/create_mesh.py
r2979 r3535 6 6 #!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 7 7 8 from pmesh.mesh import *9 from coordinate_transforms.geo_reference import Geo_reference8 from anuga.pmesh.mesh import * 9 from anuga.coordinate_transforms.geo_reference import Geo_reference 10 10 11 11 #------------------------------------------------------------- -
production/broome_2006/project.py
r3514 r3535 6 6 #from anuga.utilities.polygon import read_polygon 7 7 import sys 8 from pmesh.create_mesh import convert_points_from_latlon_to_utm9 from coordinate_transforms.redfearn import degminsec2decimal_degrees8 from anuga.pmesh.create_mesh import convert_points_from_latlon_to_utm 9 from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees 10 10 from time import localtime, strftime 11 11 from anuga.geospatial_data.geospatial_data import * … … 87 87 refzone = 50 88 88 89 from coordinate_transforms.redfearn import redfearn89 from anuga.coordinate_transforms.redfearn import redfearn 90 90 # boundary up to 50 m contour 91 91 lat1_50 = degminsec2decimal_degrees(-19,20,0) -
production/broome_2006/run_broome.py
r3514 r3535 97 97 #------------------------------------------------------------------------------- 98 98 99 from pmesh.mesh_interface import create_mesh_from_regions99 from anuga.pmesh.mesh_interface import create_mesh_from_regions 100 100 101 101 region_res = 500000 -
production/gippsland_2005/create_mesh2.py
r2301 r3535 4 4 5 5 # add inundation dir to your pythonpath 6 from pmesh.mesh import *7 from coordinate_transforms.geo_reference import Geo_reference6 from anuga.pmesh.mesh import * 7 from anuga.coordinate_transforms.geo_reference import Geo_reference 8 8 9 9 -
production/gippsland_2005/get_timeseries.py
r3514 r3535 6 6 from anuga.pyvolution.util import file_function 7 7 #from anuga.pyvolution.coordinate_transforms.redfearn import degminsec2decimal_degrees, redfearn 8 from coordinate_transforms.redfearn import degminsec2decimal_degrees, redfearn8 from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees, redfearn 9 9 from pylab import * 10 10 -
production/karratha_2005/project.py
r2291 r3535 38 38 39 39 #Georeferencing 40 from coordinate_transforms.redfearn import degminsec2decimal_degrees40 from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees 41 41 42 42 #Origin of existing dem (FIXME: Temporary measure) -
production/karratha_2005/run_karratha.py
r3514 r3535 84 84 85 85 #Create Triangular Mesh 86 from pmesh.create_mesh import create_mesh_from_regions86 from anuga.pmesh.create_mesh import create_mesh_from_regions 87 87 88 88 resolution = 4000 -
production/onslow_2006/get_timeseries.py
r3514 r3535 5 5 import project 6 6 from anuga.pyvolution.util import file_function 7 from coordinate_transforms.redfearn import degminsec2decimal_degrees, redfearn7 from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees, redfearn 8 8 from pylab import * 9 9 from os import sep -
production/onslow_2006/make_gauges.py
r3514 r3535 7 7 from anuga.utilities.polygon import poly_xy, is_inside_polygon 8 8 from Numeric import arange 9 from pmesh.create_mesh import convert_points_from_latlon_to_utm9 from anuga.pmesh.create_mesh import convert_points_from_latlon_to_utm 10 10 11 11 x_bound, y_bound = poly_xy(project.polyAll) -
production/onslow_2006/plot_data_extent.py
r3514 r3535 1 1 from pylab import * 2 2 from anuga.geospatial_data.geospatial_data import * 3 from pmesh.create_mesh import convert_points_from_latlon_to_utm4 from coordinate_transforms.redfearn import degminsec2decimal_degrees3 from anuga.pmesh.create_mesh import convert_points_from_latlon_to_utm 4 from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees 5 5 from os import sep 6 6 from anuga.utilities.polygon import inside_polygon … … 17 17 18 18 def read_community_file(filename1, filename2, polygon): 19 from coordinate_transforms.redfearn import redfearn19 from anuga.coordinate_transforms.redfearn import redfearn 20 20 from anuga.utilities.polygon import inside_polygon 21 21 fid1 = open(filename1) -
production/onslow_2006/project.py
r3514 r3535 8 8 import sys 9 9 10 from pmesh.create_mesh import convert_points_from_latlon_to_utm11 12 from coordinate_transforms.redfearn import degminsec2decimal_degrees10 from anuga.pmesh.create_mesh import convert_points_from_latlon_to_utm 11 12 from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees 13 13 14 14 from time import localtime, strftime, gmtime … … 135 135 136 136 #Georeferencing 137 from coordinate_transforms.redfearn import degminsec2decimal_degrees137 from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees 138 138 139 139 refzone = 50 -
production/onslow_2006/project_new.py
r3514 r3535 8 8 import sys 9 9 10 from pmesh.create_mesh import convert_points_from_latlon_to_utm10 from anuga.pmesh.create_mesh import convert_points_from_latlon_to_utm 11 11 12 from coordinate_transforms.redfearn import degminsec2decimal_degrees12 from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees 13 13 14 14 from time import localtime, strftime … … 85 85 86 86 #Georeferencing 87 from coordinate_transforms.redfearn import degminsec2decimal_degrees87 from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees 88 88 89 89 refzone = 50 -
production/onslow_2006/project_old.py
r3514 r3535 8 8 import sys 9 9 10 from pmesh.create_mesh import convert_points_from_latlon_to_utm10 from anuga.pmesh.create_mesh import convert_points_from_latlon_to_utm 11 11 12 from coordinate_transforms.redfearn import degminsec2decimal_degrees12 from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees 13 13 14 14 from time import localtime, strftime … … 87 87 88 88 #Georeferencing 89 from coordinate_transforms.redfearn import degminsec2decimal_degrees89 from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees 90 90 91 91 refzone = 50 -
production/onslow_2006/run_onslow.py
r3514 r3535 115 115 #------------------------------------------------------------------------------- 116 116 117 from pmesh.mesh_interface import create_mesh_from_regions117 from anuga.pmesh.mesh_interface import create_mesh_from_regions 118 118 119 119 #new -
production/onslow_2006/run_onslow_new.py
r3514 r3535 103 103 #------------------------------------------------------------------------------- 104 104 105 from pmesh.mesh_interface import create_mesh_from_regions105 from anuga.pmesh.mesh_interface import create_mesh_from_regions 106 106 107 107 # original -
production/onslow_2006/run_onslow_old.py
r3514 r3535 128 128 # resolutions (maximal area of per triangle) for each polygon 129 129 #------------------------------------------------------------------------------ 130 from pmesh.mesh_interface import create_mesh_from_regions130 from anuga.pmesh.mesh_interface import create_mesh_from_regions 131 131 132 132 # original -
production/pt_hedland_2006/project.py
r3514 r3535 7 7 #from anuga.utilities.polygon import read_polygon 8 8 import sys 9 from pmesh.create_mesh import convert_points_from_latlon_to_utm10 from coordinate_transforms.redfearn import degminsec2decimal_degrees9 from anuga.pmesh.create_mesh import convert_points_from_latlon_to_utm 10 from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees 11 11 from time import localtime, strftime 12 12 from anuga.geospatial_data.geospatial_data import * … … 111 111 refzone = 50 112 112 113 from coordinate_transforms.redfearn import redfearn113 from anuga.coordinate_transforms.redfearn import redfearn 114 114 # boundary up to 50 m contour 115 115 lat1_50 = degminsec2decimal_degrees(-19,20,0) -
production/pt_hedland_2006/run_pt_hedland.py
r3514 r3535 99 99 #------------------------------------------------------------------------------- 100 100 101 from pmesh.mesh_interface import create_mesh_from_regions101 from anuga.pmesh.mesh_interface import create_mesh_from_regions 102 102 103 103 region_res = 500000 -
production/swan_river_tin_2005/thin_generate_alpha.py
r3514 r3535 4 4 5 5 from anuga.pyvolution.load_mesh.loadASCII import reduce_pts 6 from pmesh.mesh import importMeshFromFile6 from anuga.pmesh.mesh import importMeshFromFile 7 7 8 8 thin_points_file = 'thinned_cross_sections.pts' -
production/sydney_2006/get_timeseries.py
r3514 r3535 5 5 import project 6 6 from anuga.pyvolution.util import file_function 7 #from coordinate_transforms.redfearn import degminsec2decimal_degrees, redfearn7 #from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees, redfearn 8 8 from pylab import * 9 9 from matplotlib.ticker import MultipleLocator, FormatStrFormatter -
production/sydney_2006/project.py
r3514 r3535 8 8 import sys 9 9 10 from pmesh.create_mesh import convert_points_from_latlon_to_utm10 from anuga.pmesh.create_mesh import convert_points_from_latlon_to_utm 11 11 12 12 … … 77 77 78 78 #Georeferencing 79 from coordinate_transforms.redfearn import degminsec2decimal_degrees79 from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees 80 80 81 81 # define clipping polygon -
production/sydney_2006/run_sydney.py
r3514 r3535 45 45 46 46 #Create Triangular Mesh 47 from pmesh.create_mesh import create_mesh_from_regions47 from anuga.pmesh.create_mesh import create_mesh_from_regions 48 48 49 49 if mytest == 0: -
production/sydney_2006/run_sydney_smf.py
r3514 r3535 66 66 project.combineddemname + '.pts') 67 67 68 #from pmesh.create_mesh import create_mesh_from_regions68 #from anuga.pmesh.create_mesh import create_mesh_from_regions 69 69 #new interface 70 from pmesh.mesh_interface import create_mesh_from_regions70 from anuga.pmesh.mesh_interface import create_mesh_from_regions 71 71 72 72 # original issue to Benfield -
production/sydney_2006/run_sydney_smf_polyingo.py
r3514 r3535 78 78 # resolutions (maximal area of per triangle) for each polygon 79 79 #------------------------------------------------------------------------------- 80 #from pmesh.create_mesh import create_mesh_from_regions81 from pmesh.mesh_interface import create_mesh_from_regions80 #from anuga.pmesh.create_mesh import create_mesh_from_regions 81 from anuga.pmesh.mesh_interface import create_mesh_from_regions 82 82 83 83 meshname = project.meshname4+'.msh' -
production/wollongong_2006/run_flagstaff.py
r3519 r3535 29 29 from anuga.pyvolution.shallow_water import Reflective_boundary 30 30 from anuga.pyvolution.data_manager import convert_dem_from_ascii2netcdf, dem2pts 31 from pmesh.mesh_interface import create_mesh_from_regions32 from pmesh.mesh import importUngenerateFile, Segment31 from anuga.pmesh.mesh_interface import create_mesh_from_regions 32 from anuga.pmesh.mesh import importUngenerateFile, Segment 33 33 34 34 # Parallelism
Note: See TracChangeset
for help on using the changeset viewer.