Changeset 5233
- Timestamp:
- Apr 23, 2008, 3:10:52 PM (17 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_core/source/anuga/shallow_water/tsunami_okada.py
r5206 r5233 40 40 else: 41 41 x0=xi 42 y 10=yi42 y0=yi 43 43 if domain is not None: 44 44 xllcorner = domain.geo_reference.get_xllcorner() … … 168 168 widths[0]=width 169 169 dips[0]=dip 170 xs[0]=x0 171 ys[0]=y0 170 try: 171 xs[0]=x0 172 ys[0]=y0 173 except: 174 xs[0]=x0[0] 175 ys[0]=y0[0] 172 176 else: 173 177 dislocations=dislocation … … 189 193 AL1 = 0.0 190 194 AW2 = 0.0 191 rad=1.745329252e-2195 #rad=1.745329252e-2 192 196 #zrec = domain.get_quantity('elevation').get_values(interpolation_points=[[x, y]]) #get Zrec... has to be constant!!! 193 197 #zrec=zreci.get_values(interpolation_points=[[x, y]]) … … 1079 1083 self.UA(XI[J],ET[K],Q,DD1,DD2,DD3) 1080 1084 DUA=self.DUA 1081 for I in range(0,10,3): 1085 for I in range(0,10,3): 1082 1086 DU[I] =-DUA[I] 1083 1087 DU[I+1]=-DUA[I+1]*CD+DUA[I+2]*SD … … 1198 1202 from Numeric import zeros, Float 1199 1203 1200 U=zeros( (12,1),Float)1201 DU=zeros( (12,1),Float)1202 DUA=zeros( (12,1),Float)1204 U=zeros(12,Float) 1205 DU=zeros(12,Float) 1206 DUA=zeros(12,Float) 1203 1207 F0 =0.0 1204 1208 F2=2.0 … … 1304 1308 #for I in range (0,12): 1305 1309 #DUA[I]=U[I] 1306 1307 1310 self.DUA=U 1308 1311 … … 1332 1335 from Numeric import zeros, Float 1333 1336 1334 DUB=zeros( (12,1),Float)1335 DU=zeros( (12,1),Float)1336 U=zeros( (12,1),Float)1337 DUB=zeros(12,Float) 1338 DU=zeros(12,Float) 1339 U=zeros(12,Float) 1337 1340 1338 1341 F0=0.0 … … 1499 1502 from Numeric import zeros, Float 1500 1503 1501 DUC=zeros( (12,1),Float)1502 DU=zeros( (12,1),Float)1503 U=zeros( (12,1),Float)1504 DUC=zeros(12,Float) 1505 DU=zeros(12,Float) 1506 U=zeros(12,Float) 1504 1507 1505 1508 F0=0.0
Note: See TracChangeset
for help on using the changeset viewer.