Changeset 449


Ignore:
Timestamp:
Oct 26, 2004, 6:01:37 PM (20 years ago)
Author:
ole
Message:

Bug huntin'

Location:
inundation/ga/storm_surge
Files:
6 edited

Legend:

Unmodified
Added
Removed
  • inundation/ga/storm_surge/examples/run_tsh.py

    r443 r449  
    5050
    5151    domain.checkpoint = False #True
    52     domain.default_order = 2
     52    domain.default_order = 1
    5353    domain.visualise = visualise
    5454
     
    7575
    7676    #Set friction
    77     #manning = 0.07
    78     manning = 0.3   
     77    manning = 0.07
    7978    inflow_stage = 10.0
    8079    domain.set_quantity('friction', manning)
     80
     81    #FIXME: Make flat
     82    #domain.set_quantity('elevation', 0.0)
     83    #domain.set_quantity('level', 0.0)   
    8184
    8285    ######################
  • inundation/ga/storm_surge/pyvolution/config.py

    r444 r449  
    22"""
    33
    4 epsilon = 1.0e-13
     4epsilon = 1.0e-12
    55
    66default_boundary_tag = 'exterior'
  • inundation/ga/storm_surge/pyvolution/issues.txt

    r433 r449  
    11CURRENT PYVOLUTION ISSUES
     2
    23
    34OPEN - issues that have not beeen resolved
     
    1819to == the bed elevation, if level < bed elevation.  This will make
    1920setting things up easier.
     21
     22
     23Issue: Transmissive boundary condition can cause unphysical momentums if
     24material is not moving out of the domain.
     25Importance: Med
     26Suggested Action: Add some sort of warning in (subclass of) Transmissive_boundary
     27alerting inflow of material.
    2028
    2129Issue:  Run-time error;
     
    106114_________________________________________________________
    107115
     116
    108117****************
    109118CLOSED - issues that have  beeen resolved
    110119************
     120
     121
     122New example 26/10 16:20
     123
     124V:\7\cit\workspaces\ws_nielsen\projects\storm_surge\examples>python run_
     125 cns_ground_pts_triangulation_wave.tsh 10 500
     126Creating domain from cns_ground_pts_triangulation_wave.tsh
     127Number of triangles =  4962
     128Output being written to C:\grohm_output\cns_ground_pts_triangulation_wav
     129_ys10.0_ft500.0.sww
     130Boundaries
     131Available boundary tags are ['external', 'wave']
     132Time = 0.0000, steps=0 (0)
     133Time = 10.0000, delta t in [0.00018562, 0.47635714], steps=296 (0)
     134Time = 20.0000, delta t in [0.00058530, 0.37116719], steps=267 (0)
     135Time = 30.0000, delta t in [0.00014949, 0.22819244], steps=282 (0)
     136Time = 40.0000, delta t in [0.00000000, 0.15756478], steps=478 (0)
     137Time = 50.0000, delta t in [0.00000000, 0.32346737], steps=648 (1)
     138Time = 60.0000, delta t in [0.00000000, 0.24878640], steps=664 (0)
     139Time = 70.0000, delta t in [0.00000000, 0.08966319], steps=819 (0)
     140Time = 80.0000, delta t in [0.00000000, 0.05706256], steps=1518 (1)
     141Time = 90.0000, delta t in [0.00000000, 0.06957202], steps=2129 (3)
     142Time = 100.0000, delta t in [0.00000000, 0.08955953], steps=3644 (18)
     143Time = 110.0000, delta t in [0.00000000, 0.05102101], steps=2120 (0)
     144Time = 120.0000, delta t in [0.00000000, 0.07261881], steps=2423 (0)
     145Time = 130.0000, delta t in [0.00000000, 0.07346308], steps=2366 (1)
     146Time = 140.0000, delta t in [0.00000000, 0.17995247], steps=2825 (5)
     147Time = 150.0000, delta t in [0.00000000, 0.08268458], steps=1477 (0)
     148Time = 160.0000, delta t in [0.00000000, 0.07423284], steps=3631 (0)
     149Time = 170.0000, delta t in [0.00000000, 0.07305694], steps=4369 (31)
     150Time = 180.0000, delta t in [0.00000000, 0.06342947], steps=3235 (0)
     151Time = 190.0000, delta t in [0.00000000, 0.10134596], steps=2235 (0)
     152Time = 200.0000, delta t in [0.00000000, 0.06387375], steps=4239 (10)
  • inundation/ga/storm_surge/pyvolution/netherlands.py

    r443 r449  
    9393N = 600 #Size = 720000
    9494N = 20
     95#N = 150
    9596N = 130
    9697
     
    104105   
    105106domain.check_integrity()
    106 domain.default_order = 1
    107 domain.smooth = True
     107domain.default_order = 2
     108
     109#Output params
     110domain.smooth = False
    108111domain.reduction = min  #Looks a lot better on top of steep slopes
    109112
     
    150153
    151154
    152 #Time dependent inflow
    153 from math import sin, pi, exp
    154 Bw = Time_boundary(domain=domain,
    155                    f=lambda t: [0.7*(1 + sin(t*pi/2.5))*\
    156                                 (inflow_stage*(sin(t*pi)+0.7)),0,0])
    157 
    158155#Set boundary conditions
    159 #domain.set_boundary({'left': Bw, 'right': Br, 'bottom': Br, 'top': Br})
    160 domain.set_boundary({'left': Bd, 'right': Bt, 'bottom': Bt, 'top': Bt})
     156domain.set_boundary({'left': Bd, 'right': Br, 'bottom': Br, 'top': Br})
     157#domain.set_boundary({'left': Bd, 'right': Bt, 'bottom': Bt, 'top': Bt})
    161158
    162159                 
     
    172169t0 = time.time()
    173170
    174 for t in domain.evolve(yieldstep = 0.1, finaltime = 2.0):
     171for t in domain.evolve(yieldstep = 0.1, finaltime = 5.0):
    175172    domain.write_time()
    176173   
  • inundation/ga/storm_surge/pyvolution/shallow_water.py

    r443 r449  
    450450    balance_deep_and_shallow(domain)
    451451
    452     #Protect against infinitesimal depths at vertices   
    453     ####dry(domain)
    454    
    455452    #Compute edge values by interpolation
    456453    for name in domain.conserved_quantities:
  • inundation/ga/storm_surge/pyvolution/shallow_water_ext.c

    r443 r449  
    226226      alpha = 1.0;  //Flat bed
    227227
    228 
     228     
    229229    //Weighted balance between stage parallel to bed elevation
    230230    //(wvi = zvi + hc) and stage as computed by 1st or 2nd
Note: See TracChangeset for help on using the changeset viewer.