Changeset 1014


Ignore:
Timestamp:
Mar 4, 2005, 7:57:38 PM (20 years ago)
Author:
ole
Message:
 
Location:
inundation/ga/storm_surge/pyvolution
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • inundation/ga/storm_surge/pyvolution/test_shallow_water.py

    r1013 r1014  
    15391539        domain = Domain(points, vertices, boundary)
    15401540        domain.smooth = False
    1541         domain.default_order=1
     1541        domain.default_order=2
    15421542
    15431543        #IC
     
    15921592        domain = Domain(points, vertices, boundary)
    15931593        domain.smooth = False
    1594         domain.default_order=1
     1594        domain.default_order=2
    15951595
    15961596        #IC
     
    16431643        domain = Domain(points, vertices, boundary)
    16441644        domain.smooth = False
    1645         domain.default_order=1
     1645        domain.default_order=1  #FIXME: 2 will fail
     1646        domain.beta_h = 0
    16461647
    16471648        #IC
     
    16551656                if 0.5 <= x[i] < 0.7: 
    16561657                    z[i] = 0.4       
     1658                    #z[i] = 0.7                     
    16571659                if 0.7 <= x[i]: 
    16581660                    z[i] = x[i]/3                   
     
    16891691
    16901692        #Evolution
    1691         for t in domain.evolve(yieldstep = 0.05, finaltime = 5.0):
     1693        for t in domain.evolve(yieldstep = 0.05, finaltime = 10.0):
    16921694            volume =  domain.quantities['stage'].get_integral()         
    16931695            #print t, volume, initial_volume
  • inundation/ga/storm_surge/pyvolution/twolevels.py

    r911 r1014  
    3636# beta_h == 0.0 reveals the problem
    3737# beta_h > 0.2 alleviates it
    38 domain.beta_h = 0.0   
     38domain.beta_h = 0.2   
    3939
    4040#IC
Note: See TracChangeset for help on using the changeset viewer.