Changeset 4066


Ignore:
Timestamp:
Dec 7, 2006, 9:49:27 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

    r4063 r4066  
    160160#west_boundary = degminsec2decimal_degrees(116,0,0)
    161161#east_boundary = degminsec2decimal_degrees(118,00,0)
    162 west_boundary = degminsec2decimal_degrees(116,10,0)
    163 east_boundary = degminsec2decimal_degrees(117,15,0)
     162west_boundary = degminsec2decimal_degrees(116,00,0)
     163east_boundary = degminsec2decimal_degrees(117,20,0)
    164164
    165165
     
    178178#bounding_polygon, zone =\
    179179#                  convert_from_latlon_to_utm([p1, p2, p3, p4, p5, p6, p7])
     180print bounding_polygon
    180181refzone = zone
    181182from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon
  • anuga_work/production/dampier_2006/run_dampier.py

    r4049 r4066  
    7070                 dirname(project.__file__)+sep+ project.__name__+'.py' )
    7171barrier()
    72 start_screen_catcher(project.output_run_time_dir, myid, numprocs)
     72#start_screen_catcher(project.output_run_time_dir, myid, numprocs)
    7373
    7474print 'USER: ', project.user
     
    8989                    [project.poly_facility, 500]]   
    9090#                    [project.poly_interior, 1000]]       
    91    
     91    from anuga.utilities.polygon import plot_polygons
     92    figname = project.output_run_time_dir + 'poly_pic'
     93    plot_polygons([project.poly_coast,project.poly_pipeline,project.poly_facility,
     94               project.bounding_polygon],
     95               figname,
     96               verbose = True)
     97#    if access(project.meshes_time_dir,F_OK) == 0:
     98#        mkdir(project.meshes_time_dir)
     99#    if access(project.meshes_dir,F_OK) == 0:
     100#        mkdir(project.meshes_dir)
     101    import sys; sys.exit()
     102    print 'start create mesh from regions'
     103    interior_regions = [#[project.karratha_polygon, 25000],
     104#                    [project.cipma_polygon, 1000],
     105#                    [project.poly_pipeline, 5000],
     106#                    [project.poly_facility, 500]]   
     107                    [project.poly_interior, 1000]]   
     108#    meshes_dir_name = project.meshes_dir_name + '.msh'
    92109
    93110    create_mesh_from_regions(project.bounding_polygon,
     
    194211domain.set_store_vertices_uniquely(False)
    195212domain.set_quantities_to_be_stored(['stage', 'xmomentum', 'ymomentum'])
    196 domain.set_maximum_allowed_speed(0.1) # Allow a little runoff (0.1 is OK)
     213domain.set_maximum_allowed_speed(0.0) # Allow a little runoff (0.1 is OK)
    197214
    198215#-------------------------------------------------------------------------
Note: See TracChangeset for help on using the changeset viewer.