- Timestamp:
- Aug 4, 2008, 12:08:33 PM (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_core/source/anuga/shallow_water/test_tsunami_okada.py
r5421 r5601 24 24 rectangular sources 25 25 """ 26 # get path where this test is run26 # Get path where this test is run 27 27 path = get_pathname_from_package('anuga.shallow_water') 28 #choose what test to proceed 28 29 # Choose what test to proceed 29 30 T = 1 30 31 31 32 32 33 33 if T==0: 34 # fotran output file 35 34 # Fortran output file 36 35 filename = path+sep+'fullokada_SP.txt' 37 # initial condition of earthquake for multiple source 36 37 # Initial condition of earthquake for multiple source 38 38 x0 = 7000.0 39 39 y0 = 10000.0 … … 48 48 NSMAX=1 49 49 elif T==1: 50 # fotran output file50 # Fortran output file 51 51 filename = path+sep+'fullokada_SS.txt' 52 # initial condition of earthquake for multiple source 52 53 # Initial condition of earthquake for multiple source 53 54 x0 = 7000.0 54 55 y0 = 10000.0 … … 65 66 elif T==2: 66 67 67 # fotran output file68 # Fortran output file 68 69 filename = path+sep+'fullokada_MS.txt' 69 # initial condition of earthquake for multiple source 70 71 # Initial condition of earthquake for multiple source 70 72 x0 = [7000.0,10000.0] 71 73 y0 = [10000.0,7000.0] … … 82 84 83 85 84 # get output file from original okada fortran script.85 # vertical displacement is listed under tmp.86 # Get output file from original okada fortran script. 87 # Vertical displacement is listed under tmp. 86 88 polyline_file=open(filename,'r') 87 89 lines=polyline_file.readlines() … … 160 162 """ 161 163 162 # get path where this test is run164 # Get path where this test is run 163 165 path= get_pathname_from_package('anuga.shallow_water') 164 166 165 # choose what test to proceed167 # Choose what test to proceed 166 168 T=1 167 169 168 170 if T==0: 169 # fotran output file 170 171 # Fortran output file 171 172 filename = path+sep+'fullokada_SP.txt' 172 # initial condition of earthquake for multiple source 173 174 # Initial condition of earthquake for multiple source 173 175 x0 = 7000.0 174 176 y0 = 10000.0 … … 183 185 NSMAX=1 184 186 elif T==1: 185 # fotran output file187 # Fortran output file 186 188 filename = path+sep+'fullokada_SS.txt' 187 # initial condition of earthquake for multiple source 189 190 # Initial condition of earthquake for multiple source 188 191 x0 = 7000.0 189 192 y0 = 10000.0 … … 200 203 elif T==2: 201 204 202 # fotran output file205 # Fortran output file 203 206 filename = path+sep+'fullokada_MS.txt' 204 # initial condition of earthquake for multiple source 207 208 # Initial condition of earthquake for multiple source 205 209 x0 = [7000.0,10000.0] 206 210 y0 = [10000.0,7000.0] … … 217 221 218 222 219 # get output file from original okada fortran script.220 # vertical displacement is listed under tmp.223 # Get output file from original okada fortran script. 224 # Vertical displacement is listed under tmp. 221 225 polyline_file=open(filename,'r') 222 226 lines=polyline_file.readlines()
Note: See TracChangeset
for help on using the changeset viewer.