Changeset 1839


Ignore:
Timestamp:
Sep 16, 2005, 12:11:41 PM (19 years ago)
Author:
steve
Message:

Investigating Circular Hydraulic jump.

Location:
inundation
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • inundation/analytical solutions/Analytical_solution_circular_hydraulic_jump.py

    r1834 r1839  
    113113# Order of accuracy
    114114domain.default_order = 2
     115domain.CFL = 1.0
     116#domain.beta_w = 0.5
     117#domain.beta_h = 0.2
    115118domain.smooth = True
     119
     120
    116121
    117122
     
    119124    #domain.visualiser.coloring['stage'] = True
    120125domain.visualiser.scale_z['stage'] = 2.0
    121 #domain.visualiser.scale_z['elevation'] = 0.05
     126domain.visualiser.scale_z['elevation'] = 0.05
     127
     128
     129from realtime_visualisation_new import Visualiser
     130vxmom = Visualiser(domain,title='xmomentum',scale_z=10.0)
     131vymom = Visualiser(domain,title='ymomentum',scale_z=10.0)
     132
     133
    122134#----------
    123135# Evolution
    124136import time
     137
     138print 'PROBLEM WITH INSTABILITY AFTER t'
    125139t0 = time.time()
    126140for t in domain.evolve(yieldstep = .01, finaltime = 10):
    127141    domain.write_time()
     142    vxmom.update_quantity('xmomentum')
     143    vymom.update_quantity('ymomentum')
    128144
    129145print 'That took %.2f seconds' %(time.time()-t0)
  • inundation/pyvolution/config.py

    r1697 r1839  
    7878beta_h = 0.2
    7979CFL = 1.0  #FIXME (ole): Is this in use yet??
     80           #(Steve) yes, change domain.CFL to
     81           #make changes
    8082
    8183
  • inundation/pyvolution/realtime_visualisation_new.py

    r1697 r1839  
    2222        self.frame  = frame()
    2323
    24         self.domain   = domain
    25         self.scale_z = scale_z
    26         self.vertices = domain.vertex_coordinates
     24        self.domain          = domain
     25        self.default_scale_z = scale_z
     26        self.vertices        = domain.vertex_coordinates
    2727
    2828
     
    176176
    177177        #print 'update '+qname+' arrays'
    178 
    179         qcolor   = self.qcolor[qname]
    180         scale_z  = self.scale_z[qname]
    181         coloring = self.coloring[qname]
    182         updating = self.updating[qname]
     178        try:
     179            qcolor   = self.qcolor[qname]
     180            scale_z  = self.scale_z[qname]
     181            coloring = self.coloring[qname]
     182            updating = self.updating[qname]
     183        except:
     184            qcolor   = None
     185            scale_z  = None
     186            coloring = False
     187            updating = True
    183188
    184189
     
    196201
    197202        if scale_z is None:
    198             scale_z = self.scale_z
     203            scale_z = self.default_scale_z
    199204
    200205        #try:
Note: See TracChangeset for help on using the changeset viewer.