Changeset 3997


Ignore:
Timestamp:
Nov 16, 2006, 9:27:20 AM (17 years ago)
Author:
nick
Message:

updates to dampier

Location:
anuga_work/production/dampier_2006
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/dampier_2006/project.py

    r3996 r3997  
    204204poly_facility = read_polygon(polygons_dir+'facility.csv')
    205205
    206 poly_pipeline = read_polygon(polygons_dir+'pipeline.csv')
     206poly_pipeline = read_polygon(polygons_dir+'pipeline2.csv')
     207
     208poly_interior = read_polygon(polygons_dir+'interior.csv')
     209
     210poly_coast = read_polygon(polygons_dir+'coast_final.csv')
    207211
    208212clip_poly_e = read_polygon(polygons_dir+'gap_e.csv')
  • anuga_work/production/dampier_2006/run_dampier.py

    r3940 r3997  
    8181
    8282if myid == 0:
     83   
     84    print 'start create mesh from regions'
     85    interior_regions = [#[project.karratha_polygon, 25000],
     86#                    [project.poly_coast, 10000],
     87#                    [project.poly_pipeline, 2000],
     88#                    [project.poly_facility, 500]]   
     89                    [project.poly_interior, 1000]]       
     90   
     91    from anuga.utilities.polygon import plot_polygons
     92    figname = project.output_run_time_dir + 'poly_pic'
     93    plot_polygons([interior_regions,
     94               project.bounding_polygon],
     95               figname,
     96               verbose = True)
    8397#    if access(project.meshes_time_dir,F_OK) == 0:
    8498#        mkdir(project.meshes_time_dir)
    85     if access(project.meshes_dir,F_OK) == 0:
    86         mkdir(project.meshes_dir)
    87     print 'start create mesh from regions'
    88     interior_regions = [#[project.karratha_polygon, 25000],
    89 #                    [project.cipma_polygon, 1000],
    90                     [project.poly_pipeline, 5000],
    91                     [project.poly_facility, 500]]   
     99#    if access(project.meshes_dir,F_OK) == 0:
     100#        mkdir(project.meshes_dir)
     101 
    92102#    meshes_dir_name = project.meshes_dir_name + '.msh'
    93103
     
    112122domain = Domain(meshes_dir_name, use_cache=True, verbose=True)
    113123print domain.statistics()
    114 
    115 
    116124
    117125print 'starting to create boundary conditions'
     
    238246    domain.write_boundary_statistics(tags = 'ocean')
    239247    print 'time: ',time.time()
     248 
     249    if allclose(t, 6000):
     250        domain.set_quantity('xmomentum', 0)
     251        domain.set_quantity('ymomentum', 0)
     252#import sys; sys.exit()
    240253     
    241254for t in domain.evolve(yieldstep = 60, finaltime = 28800
Note: See TracChangeset for help on using the changeset viewer.