Changeset 449
- Timestamp:
- Oct 26, 2004, 6:01:37 PM (20 years ago)
- Location:
- inundation/ga/storm_surge
- Files:
-
- 6 edited
Legend:
- Unmodified
- Added
- Removed
-
inundation/ga/storm_surge/examples/run_tsh.py
r443 r449 50 50 51 51 domain.checkpoint = False #True 52 domain.default_order = 252 domain.default_order = 1 53 53 domain.visualise = visualise 54 54 … … 75 75 76 76 #Set friction 77 #manning = 0.07 78 manning = 0.3 77 manning = 0.07 79 78 inflow_stage = 10.0 80 79 domain.set_quantity('friction', manning) 80 81 #FIXME: Make flat 82 #domain.set_quantity('elevation', 0.0) 83 #domain.set_quantity('level', 0.0) 81 84 82 85 ###################### -
inundation/ga/storm_surge/pyvolution/config.py
r444 r449 2 2 """ 3 3 4 epsilon = 1.0e-1 34 epsilon = 1.0e-12 5 5 6 6 default_boundary_tag = 'exterior' -
inundation/ga/storm_surge/pyvolution/issues.txt
r433 r449 1 1 CURRENT PYVOLUTION ISSUES 2 2 3 3 4 OPEN - issues that have not beeen resolved … … 18 19 to == the bed elevation, if level < bed elevation. This will make 19 20 setting things up easier. 21 22 23 Issue: Transmissive boundary condition can cause unphysical momentums if 24 material is not moving out of the domain. 25 Importance: Med 26 Suggested Action: Add some sort of warning in (subclass of) Transmissive_boundary 27 alerting inflow of material. 20 28 21 29 Issue: Run-time error; … … 106 114 _________________________________________________________ 107 115 116 108 117 **************** 109 118 CLOSED - issues that have beeen resolved 110 119 ************ 120 121 122 New example 26/10 16:20 123 124 V:\7\cit\workspaces\ws_nielsen\projects\storm_surge\examples>python run_ 125 cns_ground_pts_triangulation_wave.tsh 10 500 126 Creating domain from cns_ground_pts_triangulation_wave.tsh 127 Number of triangles = 4962 128 Output being written to C:\grohm_output\cns_ground_pts_triangulation_wav 129 _ys10.0_ft500.0.sww 130 Boundaries 131 Available boundary tags are ['external', 'wave'] 132 Time = 0.0000, steps=0 (0) 133 Time = 10.0000, delta t in [0.00018562, 0.47635714], steps=296 (0) 134 Time = 20.0000, delta t in [0.00058530, 0.37116719], steps=267 (0) 135 Time = 30.0000, delta t in [0.00014949, 0.22819244], steps=282 (0) 136 Time = 40.0000, delta t in [0.00000000, 0.15756478], steps=478 (0) 137 Time = 50.0000, delta t in [0.00000000, 0.32346737], steps=648 (1) 138 Time = 60.0000, delta t in [0.00000000, 0.24878640], steps=664 (0) 139 Time = 70.0000, delta t in [0.00000000, 0.08966319], steps=819 (0) 140 Time = 80.0000, delta t in [0.00000000, 0.05706256], steps=1518 (1) 141 Time = 90.0000, delta t in [0.00000000, 0.06957202], steps=2129 (3) 142 Time = 100.0000, delta t in [0.00000000, 0.08955953], steps=3644 (18) 143 Time = 110.0000, delta t in [0.00000000, 0.05102101], steps=2120 (0) 144 Time = 120.0000, delta t in [0.00000000, 0.07261881], steps=2423 (0) 145 Time = 130.0000, delta t in [0.00000000, 0.07346308], steps=2366 (1) 146 Time = 140.0000, delta t in [0.00000000, 0.17995247], steps=2825 (5) 147 Time = 150.0000, delta t in [0.00000000, 0.08268458], steps=1477 (0) 148 Time = 160.0000, delta t in [0.00000000, 0.07423284], steps=3631 (0) 149 Time = 170.0000, delta t in [0.00000000, 0.07305694], steps=4369 (31) 150 Time = 180.0000, delta t in [0.00000000, 0.06342947], steps=3235 (0) 151 Time = 190.0000, delta t in [0.00000000, 0.10134596], steps=2235 (0) 152 Time = 200.0000, delta t in [0.00000000, 0.06387375], steps=4239 (10) -
inundation/ga/storm_surge/pyvolution/netherlands.py
r443 r449 93 93 N = 600 #Size = 720000 94 94 N = 20 95 #N = 150 95 96 N = 130 96 97 … … 104 105 105 106 domain.check_integrity() 106 domain.default_order = 1 107 domain.smooth = True 107 domain.default_order = 2 108 109 #Output params 110 domain.smooth = False 108 111 domain.reduction = min #Looks a lot better on top of steep slopes 109 112 … … 150 153 151 154 152 #Time dependent inflow153 from math import sin, pi, exp154 Bw = Time_boundary(domain=domain,155 f=lambda t: [0.7*(1 + sin(t*pi/2.5))*\156 (inflow_stage*(sin(t*pi)+0.7)),0,0])157 158 155 #Set boundary conditions 159 #domain.set_boundary({'left': Bw, 'right': Br, 'bottom': Br, 'top': Br})160 domain.set_boundary({'left': Bd, 'right': Bt, 'bottom': Bt, 'top': Bt})156 domain.set_boundary({'left': Bd, 'right': Br, 'bottom': Br, 'top': Br}) 157 #domain.set_boundary({'left': Bd, 'right': Bt, 'bottom': Bt, 'top': Bt}) 161 158 162 159 … … 172 169 t0 = time.time() 173 170 174 for t in domain.evolve(yieldstep = 0.1, finaltime = 2.0):171 for t in domain.evolve(yieldstep = 0.1, finaltime = 5.0): 175 172 domain.write_time() 176 173 -
inundation/ga/storm_surge/pyvolution/shallow_water.py
r443 r449 450 450 balance_deep_and_shallow(domain) 451 451 452 #Protect against infinitesimal depths at vertices453 ####dry(domain)454 455 452 #Compute edge values by interpolation 456 453 for name in domain.conserved_quantities: -
inundation/ga/storm_surge/pyvolution/shallow_water_ext.c
r443 r449 226 226 alpha = 1.0; //Flat bed 227 227 228 228 229 229 //Weighted balance between stage parallel to bed elevation 230 230 //(wvi = zvi + hc) and stage as computed by 1st or 2nd
Note: See TracChangeset
for help on using the changeset viewer.