- Timestamp:
- Nov 12, 2008, 3:16:51 PM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_core/source_numpy_conversion/anuga/shallow_water/tsunami_okada.py
r5901 r5951 31 31 32 32 import numpy 33 from anuga.config import Float 34 33 35 34 36 def earthquake_tsunami(ns,NSMAX,length, width, strike, depth,\ … … 44 46 zrec=zrec0.get_vertex_values(xy=True) 45 47 46 x0= numpy.zeros(ns, numpy.float)47 y0= numpy.zeros(ns, numpy.float)48 x0= numpy.zeros(ns,Float) 49 y0= numpy.zeros(ns,Float) 48 50 if ns ==1: 49 51 x0[0]=xi … … 150 152 zrec=self.zrec 151 153 #initialization 152 disp0=numpy.zeros(3, numpy.float)153 strain0=numpy.zeros(6, numpy.float)154 tilt0 = numpy.zeros(2, numpy.float)155 dislocations=numpy.zeros(ns, numpy.float)156 depths=numpy.zeros(ns, numpy.float)157 strikes= numpy.zeros(ns, numpy.float)158 lengths= numpy.zeros(ns, numpy.float)159 slips= numpy.zeros(ns, numpy.float)160 rakes= numpy.zeros(ns, numpy.float)161 widths= numpy.zeros(ns, numpy.float)162 dips= numpy.zeros(ns, numpy.float)163 strikes= numpy.zeros(ns, numpy.float)164 strikes= numpy.zeros(ns, numpy.float)165 strain = numpy.zeros((N,6), numpy.float)166 disp = numpy.zeros((N,3), numpy.float)167 tilt = numpy.zeros((N,2), numpy.float)168 xs =numpy.zeros(ns, numpy.float)169 ys =numpy.zeros(ns, numpy.float)154 disp0=numpy.zeros(3,Float) 155 strain0=numpy.zeros(6,Float) 156 tilt0 = numpy.zeros(2,Float) 157 dislocations=numpy.zeros(ns,Float) 158 depths=numpy.zeros(ns,Float) 159 strikes= numpy.zeros(ns,Float) 160 lengths= numpy.zeros(ns,Float) 161 slips= numpy.zeros(ns,Float) 162 rakes= numpy.zeros(ns,Float) 163 widths= numpy.zeros(ns,Float) 164 dips= numpy.zeros(ns,Float) 165 strikes= numpy.zeros(ns,Float) 166 strikes= numpy.zeros(ns,Float) 167 strain = numpy.zeros((N,6),Float) 168 disp = numpy.zeros((N,3),Float) 169 tilt = numpy.zeros((N,2),Float) 170 xs =numpy.zeros(ns,Float) 171 ys =numpy.zeros(ns,Float) 170 172 z=[] 171 173 if ns==1: … … 213 215 xrec=y 214 216 yrec=x 217 ## print 'xrec=%s' % str(xrec) 218 ## print 'yrec=%s' % str(yrec) 219 ## print 'zrec=%s' % str(zrec) 215 220 for i in range(0,len(zrec[0])): 216 if zrec[0][i]==yrec and zrec[1][i]==xrec: 221 if numpy.alltrue(zrec[0][i]==yrec) and \ 222 numpy.alltrue(zrec[1][i]==xrec): 223 ## print 'Setting Z!' 217 224 Z=zrec[2][i] 218 225 Z=0.001*Z … … 377 384 378 385 F0=0.0 379 U=numpy.zeros((12,1), numpy.float)380 DUA=numpy.zeros((12,1), numpy.float)381 DUB=numpy.zeros((12,1), numpy.float)382 DUC=numpy.zeros((12,1), numpy.float)386 U=numpy.zeros((12,1),Float) 387 DUA=numpy.zeros((12,1),Float) 388 DUB=numpy.zeros((12,1),Float) 389 DUC=numpy.zeros((12,1),Float) 383 390 384 391 … … 523 530 WZ=self.WZ 524 531 525 DUA=numpy.zeros((12,1), numpy.float)526 DU=numpy.zeros((12,1), numpy.float)527 U=numpy.zeros((12,1), numpy.float)532 DUA=numpy.zeros((12,1),Float) 533 DU=numpy.zeros((12,1),Float) 534 U=numpy.zeros((12,1),Float) 528 535 #----- 529 536 for I in range(0,12): … … 633 640 # DATA PI2/6.283185307179586D0/ 634 641 635 DUB=numpy.zeros((12,1), numpy.float)636 DU=numpy.zeros((12,1), numpy.float)637 U=numpy.zeros((12,1), numpy.float)642 DUB=numpy.zeros((12,1),Float) 643 DU=numpy.zeros((12,1),Float) 644 U=numpy.zeros((12,1),Float) 638 645 639 646 F0=0.0 … … 807 814 # DATA PI2/6.283185307179586D0/ 808 815 809 DUC=numpy.zeros((12,1), numpy.float)810 DU=numpy.zeros((12,1), numpy.float)811 U=numpy.zeros((12,1), numpy.float)816 DUC=numpy.zeros((12,1),Float) 817 DU=numpy.zeros((12,1),Float) 818 U=numpy.zeros((12,1),Float) 812 819 813 820 F0=0.0 … … 1001 1008 EPS=0.000001 1002 1009 1003 XI=numpy.zeros(2, numpy.float)1004 ET=numpy.zeros(2, numpy.float)1005 KXI=numpy.zeros(2, numpy.float)1006 KET=numpy.zeros(2, numpy.float)1007 U=numpy.zeros(12, numpy.float)1008 DU=numpy.zeros(12, numpy.float)1009 DUA=numpy.zeros(12, numpy.float)1010 DUB=numpy.zeros(12, numpy.float)1011 DUC=numpy.zeros(12, numpy.float)1010 XI=numpy.zeros(2,Float) 1011 ET=numpy.zeros(2,Float) 1012 KXI=numpy.zeros(2,Float) 1013 KET=numpy.zeros(2,Float) 1014 U=numpy.zeros(12,Float) 1015 DU=numpy.zeros(12,Float) 1016 DUA=numpy.zeros(12,Float) 1017 DUB=numpy.zeros(12,Float) 1018 DUC=numpy.zeros(12,Float) 1012 1019 1013 1020 #----- … … 1206 1213 # DATA F0,F2,PI2/0.D0,2.D0,6.283185307179586D0/ 1207 1214 1208 U=numpy.zeros(12, numpy.float)1209 DU=numpy.zeros(12, numpy.float)1210 DUA=numpy.zeros(12, numpy.float)1215 U=numpy.zeros(12,Float) 1216 DU=numpy.zeros(12,Float) 1217 DUA=numpy.zeros(12,Float) 1211 1218 F0 =0.0 1212 1219 F2=2.0 … … 1337 1344 # DATA F0,F1,F2,PI2/0.D0,1.D0,2.D0,6.283185307179586D0/ 1338 1345 1339 DUB=numpy.zeros(12, numpy.float)1340 DU=numpy.zeros(12, numpy.float)1341 U=numpy.zeros(12, numpy.float)1346 DUB=numpy.zeros(12,Float) 1347 DU=numpy.zeros(12,Float) 1348 U=numpy.zeros(12,Float) 1342 1349 1343 1350 F0=0.0 … … 1502 1509 # DATA F0,F1,F2,F3,PI2/0.D0,1.D0,2.D0,3.D0,6.283185307179586D0/ 1503 1510 1504 DUC=numpy.zeros(12, numpy.float)1505 DU=numpy.zeros(12, numpy.float)1506 U=numpy.zeros(12, numpy.float)1511 DUC=numpy.zeros(12,Float) 1512 DU=numpy.zeros(12,Float) 1513 U=numpy.zeros(12,Float) 1507 1514 1508 1515 F0=0.0
Note: See TracChangeset
for help on using the changeset viewer.