Changeset 6695


Ignore:
Timestamp:
Apr 1, 2009, 8:55:37 PM (15 years ago)
Author:
sudi
Message:

Just modifying.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/development/anuga_1d/dam_h_sudi.py

    r5587 r6695  
    22from math import sqrt
    33#from shallow_water_h import *
    4 from shallow_water_domain import *
     4from shallow_water_domain_suggestion2 import *
    55from Numeric import zeros, Float
    66from analytic_dam_sudi import AnalyticDam
     
    3232
    3333L=2000.0
    34 N=400
     34N=200
    3535
    3636cell_len=L/N
     
    4242domain=Domain(points)
    4343
    44 domain.default_order = 2
    45 domain.default_time_order = 1
     44domain.order = 2
     45domain.set_timestepping_method('rk2')
    4646domain.cfl = 1.0
    4747domain.limiter = "vanleer"
     
    6060    return y
    6161
    62 domain.set_quantity('stage',height) #('height', height)
    63 domain.order=domain.default_order
    64 print "domain order", domain.order
     62domain.set_quantity('stage',height)
    6563
    6664domain.set_boundary({'exterior':Reflective_boundary(domain)})
     
    7472import time
    7573t0=time.time()
    76 yieldstep=30.0
     74yieldstep=20.0
    7775finaltime=20.0
    7876print "integral", domain.quantities['stage'].get_integral() #['height'].get_integral()
Note: See TracChangeset for help on using the changeset viewer.