Changeset 3563 for anuga_core/source/anuga/shallow_water
- Timestamp:
- Sep 12, 2006, 4:35:11 PM (19 years ago)
- Location:
- anuga_core/source/anuga/shallow_water
- Files:
-
- 5 moved
Legend:
- Unmodified
- Added
- Removed
-
anuga_core/source/anuga/shallow_water/data_manager.py
r3560 r3563 3547 3547 from anuga.abstract_2d_finite_volumes.pmesh2domain import pmesh_to_domain_instance 3548 3548 import time, os 3549 from anuga.abstract_2d_finite_volumes.data_manager import get_dataobject3549 #from data_manager import get_dataobject 3550 3550 from os import sep, path 3551 3551 from anuga.utilities.numerical_tools import mean -
anuga_core/source/anuga/shallow_water/shallow_water_domain.py
r3560 r3563 55 55 56 56 57 from domain import Domain as Generic_Domain 58 from generic_boundary_conditions import Boundary 59 from generic_boundary_conditions import File_boundary 60 from generic_boundary_conditions import Dirichlet_boundary 61 from generic_boundary_conditions import Time_boundary 62 from generic_boundary_conditions import Transmissive_boundary 57 from anuga.abstract_2d_finite_volumes.domain import Domain as Generic_Domain 58 from anuga.abstract_2d_finite_volumes.generic_boundary_conditions\ 59 import Boundary 60 from anuga.abstract_2d_finite_volumes.generic_boundary_conditions\ 61 import File_boundary 62 from anuga.abstract_2d_finite_volumes.generic_boundary_conditions\ 63 import Dirichlet_boundary 64 from anuga.abstract_2d_finite_volumes.generic_boundary_conditions\ 65 import Time_boundary 66 from anuga.abstract_2d_finite_volumes.generic_boundary_conditions\ 67 import Transmissive_boundary 63 68 64 69 from anuga.utilities.numerical_tools import gradient, mean … … 334 339 """ 335 340 336 from anuga. abstract_2d_finite_volumes.data_manager import get_dataobject341 from anuga.shallow_water.data_manager import get_dataobject 337 342 338 343 #Initialise writer -
anuga_core/source/anuga/shallow_water/test_data_manager.py
r3560 r3563 10 10 from Scientific.IO.NetCDF import NetCDFFile 11 11 12 from anuga. abstract_2d_finite_volumes.data_manager import *13 from anuga. abstract_2d_finite_volumes.shallow_water import *12 from anuga.shallow_water.data_manager import * 13 from anuga.shallow_water import * 14 14 from anuga.config import epsilon 15 15 16 16 # This is needed to run the tests of local functions 17 import anuga.abstract_2d_finite_volumes.data_manager17 import data_manager 18 18 19 19 from anuga.coordinate_transforms.geo_reference import Geo_reference … … 2718 2718 ################################################ 2719 2719 from mesh_factory import rectangular 2720 from shallow_water import Domain, Reflective_boundary, Dirichlet_boundary,\2721 Constant_height, Time_boundary, Transmissive_boundary2722 2720 from Numeric import array 2723 2721 … … 2764 2762 #Import the example's file as a new domain 2765 2763 ########################################## 2766 from anuga.abstract_2d_finite_volumes.data_manager import sww2domain2764 from data_manager import sww2domain 2767 2765 from Numeric import allclose 2768 2766 import os … … 2864 2862 2865 2863 from mesh_factory import rectangular 2866 from shallow_water import Domain, Reflective_boundary, Dirichlet_boundary,\2867 Constant_height, Time_boundary, Transmissive_boundary2868 2864 from Numeric import array 2869 2865 … … 2908 2904 #Import the file as a new domain 2909 2905 ################################## 2910 from anuga.abstract_2d_finite_volumes.data_manager import sww2domain2906 from data_manager import sww2domain 2911 2907 from Numeric import allclose 2912 2908 import os … … 2946 2942 2947 2943 #def test_weed(self): 2948 from anuga.abstract_2d_finite_volumes.data_manager import weed2944 from data_manager import weed 2949 2945 2950 2946 coordinates1 = [[0.,0.],[1.,0.],[1.,1.],[1.,0.],[2.,0.],[1.,1.]] … … 2973 2969 ################################################ 2974 2970 from mesh_factory import rectangular 2975 from shallow_water import Domain, Reflective_boundary, Dirichlet_boundary,\2976 Constant_height, Time_boundary, Transmissive_boundary2977 2971 from Numeric import array 2978 2972 #Create basic mesh … … 3019 3013 #Import the example's file as a new domain 3020 3014 ########################################## 3021 from anuga.abstract_2d_finite_volumes.data_manager import sww2domain3015 from data_manager import sww2domain 3022 3016 from Numeric import allclose 3023 3017 import os … … 3320 3314 from Scientific.IO.NetCDF import NetCDFFile 3321 3315 3322 from anuga.abstract_2d_finite_volumes.data_manager import _read_asc3316 from data_manager import _read_asc 3323 3317 #Write test asc file 3324 3318 filename = tempfile.mktemp(".000") … … 4035 4029 # k - lat 4036 4030 # l - lon 4037 kmin, kmax, lmin, lmax = anuga.abstract_2d_finite_volumes.data_manager._get_min_max_indexes(4031 kmin, kmax, lmin, lmax = data_manager._get_min_max_indexes( 4038 4032 latitudes,longitudes, 4039 4033 -10,4,-10,31) … … 4050 4044 4051 4045 ## 2nd test 4052 kmin, kmax, lmin, lmax = anuga.abstract_2d_finite_volumes.data_manager._get_min_max_indexes(4046 kmin, kmax, lmin, lmax = data_manager._get_min_max_indexes( 4053 4047 latitudes,longitudes, 4054 4048 0.5,2.5,5,25) … … 4065 4059 4066 4060 ## 3rd test 4067 kmin, kmax, lmin, lmax = anuga.abstract_2d_finite_volumes.data_manager._get_min_max_indexes(\4061 kmin, kmax, lmin, lmax = data_manager._get_min_max_indexes(\ 4068 4062 latitudes, 4069 4063 longitudes, … … 4082 4076 4083 4077 ## 4th test 4084 kmin, kmax, lmin, lmax = anuga.abstract_2d_finite_volumes.data_manager._get_min_max_indexes(4078 kmin, kmax, lmin, lmax = data_manager._get_min_max_indexes( 4085 4079 latitudes,longitudes, 4086 4080 -0.1,1.9,-2,17) … … 4096 4090 'failed') 4097 4091 ## 5th test 4098 kmin, kmax, lmin, lmax = anuga.abstract_2d_finite_volumes.data_manager._get_min_max_indexes(4092 kmin, kmax, lmin, lmax = data_manager._get_min_max_indexes( 4099 4093 latitudes,longitudes, 4100 4094 0.1,1.9,2,17) … … 4112 4106 ## 6th test 4113 4107 4114 kmin, kmax, lmin, lmax = anuga.abstract_2d_finite_volumes.data_manager._get_min_max_indexes(4108 kmin, kmax, lmin, lmax = data_manager._get_min_max_indexes( 4115 4109 latitudes,longitudes, 4116 4110 1.5,4,18,32) … … 4129 4123 ## 7th test 4130 4124 m2d = array([[0,1,2,3],[4,5,6,7],[8,9,10,11],[12,13,14,15]]) 4131 kmin, kmax, lmin, lmax = anuga.abstract_2d_finite_volumes.data_manager._get_min_max_indexes(4125 kmin, kmax, lmin, lmax = data_manager._get_min_max_indexes( 4132 4126 latitudes,longitudes, 4133 4127 1.5,1.5,15,15) … … 4156 4150 # k - lat 4157 4151 # l - lon 4158 kmin, kmax, lmin, lmax = anuga.abstract_2d_finite_volumes.data_manager._get_min_max_indexes(4152 kmin, kmax, lmin, lmax = data_manager._get_min_max_indexes( 4159 4153 latitudes,longitudes, 4160 4154 -37,-27,147,149.5) … … 4185 4179 # k - lat 4186 4180 # l - lon 4187 kmin, kmax, lmin, lmax = anuga.abstract_2d_finite_volumes.data_manager._get_min_max_indexes(4181 kmin, kmax, lmin, lmax = data_manager._get_min_max_indexes( 4188 4182 latitudes,longitudes, 4189 4183 -43,-37,148.5,149.5) … … 4209 4203 # k - lat 4210 4204 # l - lon 4211 kmin, kmax, lmin, lmax = anuga.abstract_2d_finite_volumes.data_manager._get_min_max_indexes(4205 kmin, kmax, lmin, lmax = data_manager._get_min_max_indexes( 4212 4206 latitudes,longitudes) 4213 4207 -
anuga_core/source/anuga/shallow_water/test_shallow_water_domain.py
r3560 r3563 8 8 from anuga.utilities.numerical_tools import mean 9 9 10 from shallow_water import *10 from shallow_water_domain import * 11 11 12 12 … … 201 201 vertices = [ [1,0,2], [1,2,4], [4,2,5], [3,1,4]] 202 202 203 204 #from anuga.abstract_2d_finite_volumes.domain import Domain as Generic_domain 205 #msg = 'The class %s is not a subclass of the generic domain class %s'\ 206 # %(DomainClass, Domain) 207 #assert issubclass(DomainClass, Domain), msg 208 203 209 domain = Domain(points, vertices) 204 210 domain.check_integrity() … … 1029 1035 1030 1036 #Convert ASCII file to NetCDF (Which is what we really like!) 1031 from anuga.abstract_2d_finite_volumes.data_manager import timefile2netcdf1037 from data_manager import timefile2netcdf 1032 1038 timefile2netcdf(filename) 1033 1039 os.remove(filename + '.txt') … … 1584 1590 """ 1585 1591 from mesh_factory import rectangular 1586 from shallow_water import Domain, Reflective_boundary,\1587 Dirichlet_boundary, Constant_height1588 1592 from Numeric import array 1589 1593 … … 1638 1642 """ 1639 1643 from mesh_factory import rectangular 1640 from shallow_water import Domain, Reflective_boundary,\1641 Dirichlet_boundary, Constant_height1642 1644 from Numeric import array 1643 1645 … … 1690 1692 """ 1691 1693 from mesh_factory import rectangular 1692 from shallow_water import Domain, Reflective_boundary,\1693 Dirichlet_boundary, Constant_height1694 1694 from Numeric import array 1695 1695 … … 1767 1767 """ 1768 1768 from mesh_factory import rectangular 1769 from shallow_water import Domain, Reflective_boundary,\1770 Dirichlet_boundary, Constant_height1771 1769 from Numeric import array 1772 1770 … … 1855 1853 """ 1856 1854 from mesh_factory import rectangular 1857 from shallow_water import Domain, Reflective_boundary,\1858 Dirichlet_boundary, Constant_height1859 1855 from Numeric import array 1860 1856 … … 1929 1925 from os import sep; sys.path.append('..'+sep+'abstract_2d_finite_volumes') 1930 1926 from mesh_factory import rectangular 1931 from shallow_water import Domain, Reflective_boundary1932 1927 1933 1928 … … 1981 1976 1982 1977 from mesh_factory import rectangular 1983 from shallow_water import Domain, Reflective_boundary,\1984 Dirichlet_boundary, Constant_height1985 1978 from Numeric import array 1986 1979 … … 2051 2044 2052 2045 from mesh_factory import rectangular 2053 from shallow_water import Domain, Reflective_boundary,\2054 Dirichlet_boundary, Constant_height2055 2046 from Numeric import array 2056 2047 … … 2084 2075 def test_flatbed_first_order(self): 2085 2076 from mesh_factory import rectangular 2086 from shallow_water import Domain,\2087 Reflective_boundary, Dirichlet_boundary2088 2089 2077 from Numeric import array 2090 2078 … … 2131 2119 def test_flatbed_second_order(self): 2132 2120 from mesh_factory import rectangular 2133 from shallow_water import Domain,\2134 Reflective_boundary, Dirichlet_boundary2135 2136 2121 from Numeric import array 2137 2122 … … 2192 2177 def test_flatbed_second_order_vmax_0(self): 2193 2178 from mesh_factory import rectangular 2194 from shallow_water import Domain,\2195 Reflective_boundary, Dirichlet_boundary2196 2197 2179 from Numeric import array 2198 2180 … … 2241 2223 #painfully setup and extracted. 2242 2224 from mesh_factory import rectangular 2243 from shallow_water import Domain,\2244 Reflective_boundary, Dirichlet_boundary2245 2246 2225 from Numeric import array 2247 2226 … … 2391 2370 2392 2371 from mesh_factory import rectangular 2393 from shallow_water import Domain, Reflective_boundary, Constant_height2394 2372 from Numeric import array 2395 2373 … … 2429 2407 2430 2408 from mesh_factory import rectangular 2431 from shallow_water import Domain, Reflective_boundary, Constant_height2432 2409 from Numeric import array 2433 2410 … … 2493 2470 2494 2471 from mesh_factory import rectangular 2495 from shallow_water import Domain, Reflective_boundary, Constant_height2496 2472 from Numeric import array 2497 2473 … … 2580 2556 2581 2557 from mesh_factory import rectangular 2582 from shallow_water import Domain, Reflective_boundary, Constant_height2583 2558 from Numeric import array 2584 2559 … … 2673 2648 2674 2649 from mesh_factory import rectangular 2675 from shallow_water import Domain, Reflective_boundary, Constant_height2676 2650 from Numeric import array 2677 2651 … … 2762 2736 2763 2737 from mesh_factory import rectangular 2764 from shallow_water import Domain, Reflective_boundary, Constant_height2765 2738 from Numeric import array 2766 2739 … … 2879 2852 2880 2853 from mesh_factory import rectangular 2881 from shallow_water import Domain, Reflective_boundary, Constant_height2882 2854 from Numeric import array 2883 2855 … … 2921 2893 def test_complex_bed(self): 2922 2894 #No friction is tested here 2923 2924 from shallow_water import Domain, Reflective_boundary, Dirichlet_boundary,\2925 Transmissive_boundary, Time_boundary,\2926 Weir_simple as Weir, Constant_height2927 2895 2928 2896 from mesh_factory import rectangular
Note: See TracChangeset
for help on using the changeset viewer.