Changeset 912


Ignore:
Timestamp:
Feb 17, 2005, 11:27:49 PM (19 years ago)
Author:
ole
Message:

Played more with 'lapping' issue.
I think it was largely due to missing
distribute_to_vertices_and_edges in shallow_water evolve

Location:
inundation/ga/storm_surge/pyvolution
Files:
2 edited

Legend:

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

    r910 r912  
    1818
    1919#Create basic mesh
    20 N = 10
     20N = 20
    2121points, vertices, boundary = rectangular(N, N, 100, 100)
    2222
     
    2525domain.store = True
    2626domain.smooth = False
    27 domain.visualise = True
     27domain.visualise = False
    2828domain.set_name('island')
    2929print 'Output being written to ' + domain.get_datadir() + sep + \
     
    3636#I tried to introduce this parameter top control the h-limiter,
    3737#but it doesn't remove the 'lapping water'
     38#NEW (ON): I think it has been fixed (or at least reduced significantly) now
    3839#
    3940# beta_h == 1.0 means that the largest gradients (on h) are allowed
    4041# beta_h == 0.0 means that constant (1st order) gradients are introduced
    4142# on h. This is equivalent to the constant depth used previously.
    42 domain.beta_h = 0.0
     43domain.beta_h = 0.2
    4344
    4445
     
    6768######################
    6869#Evolution
    69 for t in domain.evolve(yieldstep = 0.05, finaltime = 0.5):
     70import time
     71for t in domain.evolve(yieldstep = 1, finaltime = 20):
    7072    domain.write_time()
    7173
  • inundation/ga/storm_surge/pyvolution/shallow_water.py

    r908 r912  
    172172        #self.check_integrity()
    173173       
     174       
     175        #Initial update of vertex and edge values before any storage
     176        #and or visualisation
     177        self.distribute_to_vertices_and_edges()
     178       
     179               
    174180        #Initialise real time viz if requested
    175181        if self.visualise is True and self.time == 0.0:
     
    782788        # Momentum is balanced between constant and limited
    783789
     790       
     791        #for i in range(3):
     792        #    wv[k,i] = zv[k,i] + hvbar[k,i]
     793       
     794        #return
     795       
    784796        if alpha < 1:         
    785797       
Note: See TracChangeset for help on using the changeset viewer.