Changeset 5233


Ignore:
Timestamp:
Apr 23, 2008, 3:10:52 PM (15 years ago)
Author:
herve
Message:
 
File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_core/source/anuga/shallow_water/tsunami_okada.py

    r5206 r5233  
    4040    else:
    4141        x0=xi
    42         y10=yi
     42        y0=yi
    4343    if domain is not None:
    4444        xllcorner = domain.geo_reference.get_xllcorner()
     
    168168            widths[0]=width
    169169            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]
    172176        else:
    173177            dislocations=dislocation
     
    189193        AL1 = 0.0
    190194        AW2 = 0.0
    191         rad=1.745329252e-2
     195        #rad=1.745329252e-2
    192196        #zrec = domain.get_quantity('elevation').get_values(interpolation_points=[[x, y]])   #get Zrec... has to be constant!!!
    193197        #zrec=zreci.get_values(interpolation_points=[[x, y]])
     
    10791083             self.UA(XI[J],ET[K],Q,DD1,DD2,DD3)
    10801084             DUA=self.DUA
    1081              for I in range(0,10,3):                                                 
     1085             for I in range(0,10,3):
    10821086                  DU[I]  =-DUA[I]                                               
    10831087                  DU[I+1]=-DUA[I+1]*CD+DUA[I+2]*SD                             
     
    11981202      from Numeric import zeros, Float
    11991203       
    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)
    12031207      F0 =0.0
    12041208      F2=2.0
     
    13041308      #for I in range (0,12):
    13051309          #DUA[I]=U[I]     
    1306 
    13071310      self.DUA=U
    13081311
     
    13321335      from Numeric import zeros, Float
    13331336       
    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)
    13371340     
    13381341      F0=0.0
     
    14991502      from Numeric import zeros, Float
    15001503       
    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)
    15041507     
    15051508      F0=0.0
Note: See TracChangeset for help on using the changeset viewer.