Changeset 4372


Ignore:
Timestamp:
Apr 13, 2007, 4:59:58 PM (17 years ago)
Author:
nick
Message:

update broome

Location:
anuga_work/production/broome_2006
Files:
2 edited

Legend:

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

    r4368 r4372  
    2828print 'gtime: ', gtime
    2929
    30 tide = 0
     30#tide = -5.3
     31#tide = 0
     32tide = 4.9
    3133
    3234#Making assumptions about the location of scenario data
     
    5153combined_name ='broome_combined_elevation'
    5254combined_name1 ='broome_combined_elevation1'
     55combined_name_unclipped1 ='broome_combined_elevation_unclipped1'
    5356combined_small_name = 'broome_combined_elevation_small'
    5457
     
    7982#final topo files
    8083combined_dir_name = topographies_dir + combined_name
     84combined_dir_name_unclipped1 = topographies_dir + combined_name_unclipped1
    8185combined_dir_name1 = topographies_dir + combined_name1
    8286combined_small_name_dir = topographies_dir + combined_small_name
     
    120124from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon
    121125
    122 poly_all = read_polygon(polygons_dir+'extent.csv')
     126poly_all = read_polygon(polygons_dir+'extent_small.csv')
     127#poly_all = read_polygon(polygons_dir+'extent.csv')
    123128res_poly_all = 500000
     129#res_poly_all = 150000
    124130
    125131###############################
     
    129135poly_0 = read_polygon(polygons_dir+'neg20_coast_contour_pts.csv')
    130136res_0 = 100000
     137#res_0 = 20000
    131138
    132139poly_1 = read_polygon(polygons_dir+'broome_north_coast_inside_extent.csv')
    133140res_1 = 50000
     141#res_1 = 5000
    134142
    135143poly_2 = read_polygon(polygons_dir+'broome_south_coast_inside_extent.csv')
    136144res_2 = 50000
     145#res_2 = 5000
    137146
    138147poly_3 = read_polygon(polygons_dir+'Broome_town_pts.csv')
    139148res_3 = 20000
     149#res_3 = 2000
    140150
    141151poly_4 = read_polygon(polygons_dir+'Broome_inner_town_pts.csv')
    142152res_4 = 5000
     153#res_4 = 500
    143154#assert zone == refzone
    144155
  • anuga_work/production/broome_2006/run_broome.py

    r4367 r4372  
    3232from anuga_parallel.parallel_api import distribute, numprocs, myid, barrier
    3333from anuga.utilities.polygon import plot_polygons, polygon_area
     34from anuga_parallel.parallel_abstraction import get_processor_name
    3435
    3536# Application specific imports
     
    5657
    5758# creates copy of code in output dir
    58 
     59print "Processor Name:",get_processor_name()
    5960
    6061print 'USER: ', project.user
     
    7374    print 'start create mesh from regions'
    7475    create_mesh_from_regions(project.poly_all,
     76#                         boundary_tags={'back': [1, 2], 'side': [0,3],
     77#                                        'ocean': [4, 5, 6]},
    7578                         boundary_tags={'back': [1, 2], 'side': [0,3],
    76                                         'ocean': [4, 5, 6]},
     79                                        'ocean': [4]},
    7780                         maximum_triangle_area=project.res_poly_all,
    7881                         interior_regions=project.interior_regions,
     
    135138    #combined_time_dir_name = project.topographies_dir+build_time+project.combined_name
    136139    print 'Start Set quantity'
    137 
     140    print 'project.combined_dir_name_unclipped1',project.combined_dir_name_unclipped1+'.txt'
    138141    domain.set_quantity('elevation',
    139                     filename = project.combined_dir_name1+'.txt',
     142                    filename = project.combined_dir_name_unclipped1+'.txt',
    140143#                    filename = project.combined_smaller_name_dir+'.xya',
    141144                    use_cache = True,
Note: See TracChangeset for help on using the changeset viewer.