Changeset 6511


Ignore:
Timestamp:
Mar 13, 2009, 1:25:37 PM (15 years ago)
Author:
jgriffin
Message:

Fized default boundary in Field_boundary

Location:
anuga_work/production/australia_ph2
Files:
14 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/australia_ph2/LHI/run_model.py

    r6502 r6511  
    8080# Reading the landward defined points, this incorporates the original clipping
    8181# polygon minus the 100m contour
    82 ##landward_boundary = read_polygon(project.landward_boundary)
     82landward_boundary = read_polygon(project.landward_boundary)
    8383
    8484# Combine sts polyline with landward points
    85 bounding_polygon_sts = event_sts #+ landward_boundary
     85bounding_polygon_sts = event_sts + landward_boundary
    8686
    8787# Number of boundary segments
    8888num_ocean_segments = len(event_sts) - 1
    8989# Number of landward_boundary points
    90 num_land_points = 0 #file_length(project.landward_boundary)
     90num_land_points = file_length(project.landward_boundary)
    9191
    9292# Boundary tags refer to project.landward_boundary
     
    151151                        domain, mean_stage=project.tide,
    152152                        time_thinning=1,
    153                         default_boundary=Bd,
     153                        default_boundary=Dirichlet_boundary([0, 0, 0]),
    154154                        boundary_polygon=bounding_polygon_sts,                   
    155155                        use_cache=True,
     
    160160   
    161161elif project.wave == 'Tb':
    162     Tb = Time_boundary(domain,f,default_boundary=Bd,
    163                        verbose=True)
     162    Tb = Time_boundary(domain,f,default_boundary=Bd )
    164163
    165164    domain.set_boundary({'back': Br,
  • anuga_work/production/australia_ph2/adelaide/run_model.py

    r6460 r6511  
    151151                        domain, mean_stage=project.tide,
    152152                        time_thinning=1,
    153                         default_boundary=Bd,
     153                        default_boundary=Dirichlet_boundary([0, 0, 0]),
    154154                        boundary_polygon=bounding_polygon_sts,                   
    155155                        use_cache=True,
  • anuga_work/production/australia_ph2/albany/run_model.py

    r6484 r6511  
    151151                        domain, mean_stage=project.tide,
    152152                        time_thinning=1,
    153                         default_boundary=Bd,
     153                        default_boundary=Dirichlet_boundary([0, 0, 0]),
    154154                        boundary_polygon=bounding_polygon_sts,                   
    155155                        use_cache=True,
     
    160160   
    161161elif project.wave == 'Tb':
    162     Tb = Time_boundary(domain,f,default_boundary=Bd,
    163                        verbose=True)
     162    Tb = Time_boundary(domain,f,default_boundary=Bd )
    164163
    165164    domain.set_boundary({'back': Br,
  • anuga_work/production/australia_ph2/broome/run_model.py

    r6460 r6511  
    151151                        domain, mean_stage=project.tide,
    152152                        time_thinning=1,
    153                         default_boundary=Bd,
     153                        default_boundary=Dirichlet_boundary([0, 0, 0]),
    154154                        boundary_polygon=bounding_polygon_sts,                   
    155155                        use_cache=True,
  • anuga_work/production/australia_ph2/carnarvon/run_model.py

    r6460 r6511  
    151151                        domain, mean_stage=project.tide,
    152152                        time_thinning=1,
    153                         default_boundary=Bd,
     153                        default_boundary=Dirichlet_boundary([0, 0, 0]),
    154154                        boundary_polygon=bounding_polygon_sts,                   
    155155                        use_cache=True,
  • anuga_work/production/australia_ph2/ceduna/run_model.py

    r6476 r6511  
    151151                        domain, mean_stage=project.tide,
    152152                        time_thinning=1,
    153                         default_boundary=Bd,
     153                        default_boundary=Dirichlet_boundary([0, 0, 0]),
    154154                        boundary_polygon=bounding_polygon_sts,                   
    155155                        use_cache=True,
     
    160160   
    161161elif project.wave == 'Tb':
    162     Tb = Time_boundary(domain,f,default_boundary=Bd,
    163                        verbose=True)
     162    Tb = Time_boundary(domain,f,default_boundary=Bd )
    164163
    165164    domain.set_boundary({'back': Br,
  • anuga_work/production/australia_ph2/darwin/run_model.py

    r6476 r6511  
    151151                        domain, mean_stage=project.tide,
    152152                        time_thinning=1,
    153                         default_boundary=Bd,
     153                        default_boundary=Dirichlet_boundary([0, 0, 0]),
    154154                        boundary_polygon=bounding_polygon_sts,                   
    155155                        use_cache=True,
     
    160160   
    161161elif project.wave == 'Tb':
    162     Tb = Time_boundary(domain,f,default_boundary=Bd,
    163                        verbose=True)
     162    Tb = Time_boundary(domain,f,default_boundary=Bd )
    164163
    165164    domain.set_boundary({'back': Br,
  • anuga_work/production/australia_ph2/esperance/run_model.py

    r6476 r6511  
    151151                        domain, mean_stage=project.tide,
    152152                        time_thinning=1,
    153                         default_boundary=Bd,
     153                        default_boundary=Dirichlet_boundary([0, 0, 0]),
    154154                        boundary_polygon=bounding_polygon_sts,                   
    155155                        use_cache=True,
     
    160160   
    161161elif project.wave == 'Tb':
    162     Tb = Time_boundary(domain,f,default_boundary=Bd,
    163                        verbose=True)
    164    
     162    Tb = Time_boundary(domain,f,default_boundary=Bd )
     163
    165164    domain.set_boundary({'back': Br,
    166165                         'side': Bd,
  • anuga_work/production/australia_ph2/eucla_motel/run_model.py

    r6476 r6511  
    151151                        domain, mean_stage=project.tide,
    152152                        time_thinning=1,
    153                         default_boundary=Bd,
     153                        default_boundary=Dirichlet_boundary([0, 0, 0]),
    154154                        boundary_polygon=bounding_polygon_sts,                   
    155155                        use_cache=True,
     
    160160   
    161161elif project.wave == 'Tb':
    162     Tb = Time_boundary(domain,f,default_boundary=Bd,
    163                        verbose=True)
     162    Tb = Time_boundary(domain,f,default_boundary=Bd )
    164163
    165164    domain.set_boundary({'back': Br,
  • anuga_work/production/australia_ph2/gold_coast/project.py

    r6496 r6511  
    2828tide = 0                # difference between MSL and HAT
    2929zone = 56               # specify zone of model
    30 event_number = 51453    # the event number or the mux file name
     30event_number = 50883    # the event number or the mux file name
    3131alpha = 0.1             # smoothing parameter for mesh
    3232friction=0.01           # manning's friction coefficient
     
    220220# full path to where MUX files (or meta-files) live
    221221mux_input = join(event_folder, mux_input_filename)
    222 
     222#full path to boundary polygon csv
     223boundary_csv = join(event_folder, 'sts_gauge_' + str(index) +'.csv')
     224
  • anuga_work/production/australia_ph2/hobart/run_model.py

    r6460 r6511  
    151151                        domain, mean_stage=project.tide,
    152152                        time_thinning=1,
    153                         default_boundary=Bd,
     153                        default_boundary=Dirichlet_boundary([0, 0, 0]),
    154154                        boundary_polygon=bounding_polygon_sts,                   
    155155                        use_cache=True,
  • anuga_work/production/australia_ph2/perth/run_model.py

    r6460 r6511  
    151151                        domain, mean_stage=project.tide,
    152152                        time_thinning=1,
    153                         default_boundary=Bd,
     153                        default_boundary=Dirichlet_boundary([0, 0, 0]),
    154154                        boundary_polygon=bounding_polygon_sts,                   
    155155                        use_cache=True,
  • anuga_work/production/australia_ph2/sydney/project.py

    r6498 r6511  
    2828# One or all can be changed each time the run_model script is executed
    2929tide = 0                # difference between MSL and HAT
    30 zone = 50               # specify zone of model
    31 event_number = 58326    # the event number or the mux file name
     30zone = 56              # specify zone of model
     31event_number = 58222    # the event number or the mux file name
    3232alpha = 0.1             # smoothing parameter for mesh
    3333friction=0.01           # manning's friction coefficient
     
    3939wave = 'Bf'             # Bf (sts wave) Tb (index wave)
    4040
    41 setup = 'final'         # This can be one of three values
     41setup = 'trial'         # This can be one of three values
    4242                        #    trial - coarsest mesh, fast
    4343                        #    basic - coarse mesh
     
    8888# LAND - used to set the initial stage/water to be offcoast only
    8989# Used in run_model.py.  Format for points easting,northing (no header)
    90 land_initial_conditions_filename = [['initial_condition_extend.csv', 0],
    91                                     ['initial_condition_comerong_island.csv', 0],
    92                                     ['initial_condition_gabo_island.csv', 0],
    93                                     ['initial_condition_montague_island.csv', 0]]
     90land_initial_conditions_filename = []
    9491
    9592# GAUGES - for creating timeseries at a specific point
     
    10198# Used in run_building_inundation.py
    10299# Format latitude,longitude etc (geographic)
    103 building_exposure_filename = 'busselton_res_clip.csv' # from NEXIS
     100building_exposure_filename = 'sydney_res_clip.csv' # from NEXIS
    104101
    105102# BOUNDING POLYGON - used in build_boundary.py and run_model.py respectively
  • anuga_work/production/australia_ph2/sydney/run_model.py

    r6460 r6511  
    151151                        domain, mean_stage=project.tide,
    152152                        time_thinning=1,
    153                         default_boundary=Bd,
     153                        default_boundary=Dirichlet_boundary([0, 0, 0]),
    154154                        boundary_polygon=bounding_polygon_sts,                   
    155155                        use_cache=True,
Note: See TracChangeset for help on using the changeset viewer.