Changeset 2441


Ignore:
Timestamp:
Feb 22, 2006, 6:08:30 PM (19 years ago)
Author:
nick
Message:
 
Location:
production/onslow_2006
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • production/onslow_2006/project.py

    r2436 r2441  
    6767polyAll = [d0, d1, d2, d3, d4, d5]
    6868
    69 #Interior region - Onslow - take 1
    70 i0 = [310000, 7610000]
    71 i1 = [300000, 7610000]
    72 i2 = [300000, 7600000]
    73 i3 = [310000, 7600000]
     69#Interior region - Onslow town
    7470
    75 polyInterior = [i0, i1, i2, i3]
     71i0 = [304000, 7608000]
     72i1 = [302000, 7605000]
     73i2 = [303000, 7602000]
     74i3 = [305000, 7601000]
     75i4 = [309000, 7603000]
     76i5 = [307000, 7606500]
     77
     78poly_onslow = [i0, i1, i2, i3, i4, i5]
     79
     80#Thevenard Island
     81j0 = [294000, 7629000]
     82j1 = [285000, 7625000]
     83j2 = [294000, 7621000]
     84j3 = [299000, 7625000]
     85
     86poly_thevenard = [j0, j1, j2, j3]
     87
     88# Direction Is
     89k0 = [309000, 7619000]
     90k1 = [304000, 7619000]
     91k2 = [304000, 7616500]
     92k3 = [309000, 7616500]
     93
     94poly_direction = [k0, k1, k2, k3]
    7695
    7796#!slump_origin = [385000.0, 6255000.0] #Absolute UTM
  • production/onslow_2006/run_onslow.py

    r2437 r2441  
    7676# original
    7777interior_res = 5000
    78 interior_regions = [[project.polyInterior, interior_res]]
     78interior_regions = [[project.poly_onslow, interior_res],
     79                    [project.poly_thevenard, interior_res],
     80                    [project.poly_direction, interior_res]]
    7981
    8082#FIXME: Fix caching of this one once the mesh_interface is ready
     
    150152                   f=lambda t: [(20<t<200)*6, 0, 0])
    151153
    152 domain.set_boundary( {'top': Bw, 'topleft': Br,
     154domain.set_boundary( {'top': Br, 'topleft': Br,
    153155                             'left': Br, 'bottom': Br,
    154156                             'bottomright': Br, 'topright': Br} )
     
    162164t0 = time.time()
    163165
    164 for t in domain.evolve(yieldstep = 10, finaltime = 300):
     166for t in domain.evolve(yieldstep = 10, finaltime = 20):
    165167    domain.write_time()
    166168    domain.write_boundary_statistics(tags = 'top')     
Note: See TracChangeset for help on using the changeset viewer.