Changeset 6321


Ignore:
Timestamp:
Feb 12, 2009, 3:21:39 PM (15 years ago)
Author:
kristy
Message:

Adjusted for larger area, not running yet

Location:
anuga_work/production/patong
Files:
1 deleted
2 edited

Legend:

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

    r6271 r6321  
    9292
    9393# gauges - used in get_timeseries.py
    94 gauge_name = 'patong.csv'
     94gauge_name = 'gauges.csv'
    9595
    9696# BOUNDING POLYGON - used in build_boundary.py and run_patong.py respectively
     
    9898# Check the run_patong.py for boundary_tags
    9999# thinned ordering file from Hazard Map: format is index,latitude,longitude (with title)
    100 order_filename = 'thinned_boundary_ordering.csv'
     100order_filename = 'urs_order.csv'
    101101#landward bounding points
    102 landward = 'landward_bounding_polygon.csv'
     102landward = 'landward_bounding_polygon2.csv'
    103103
    104104#------------------------------------------------------------------------------
     
    107107# Output filename for elevation
    108108# this is a combination of all the data (utilisied in build_boundary)
    109 combined_name ='patong_combined_elevation'
     109combined_name ='patong_combined_elevation_large'
    110110combined_smaller_name = 'patong_combined_elevation_smaller'
    111111
     
    174174
    175175# Initial bounding polygon for data clipping
    176 poly_all = read_polygon(polygons_dir+'poly_all.csv')
     176poly_all = read_polygon(polygons_dir+'bounding_polygon.csv')
    177177res_poly_all = 150000*res_factor
    178178
  • anuga_work/production/patong/run_patong.py

    r6264 r6321  
    3333from anuga.interface import create_sts_boundary
    3434from anuga.interface import csv2building_polygons
     35from file_length import file_length
    3536
    3637from anuga.caching import cache
     
    8081# Counting segments
    8182N = len(urs_bounding_polygon)-1
    82 
    83 # boundary tags refer to project.landward 4 points equals 5 segments start at N
    84 boundary_tags={'back': [N+1,N+2,N+3], 'side': [N,N+4], 'ocean': range(N)}
     83# Number of landward_boundary points
     84M = file_length(project.landward_dir)
     85
     86# Boundary tags refer to project.landward_boundary
     87# 4 points equals 5 segments start at N
     88boundary_tags={'back': range(N+1,N+M),
     89               'side': [N,N+M],
     90               'ocean': range(N)}
    8591
    8692#--------------------------------------------------------------------------
     
    176182
    177183Br = Reflective_boundary(domain)
    178 #Bd = Dirichlet_boundary([project.tide,0,0])
    179 Bs = Transmissive_stage_zero_momentum_boundary(domain)
     184Bd = Dirichlet_boundary([project.tide,0,0])
     185#Bs = Transmissive_stage_zero_momentum_boundary(domain)
    180186Bf = Field_boundary(urs_boundary_name+'.sts',
    181187                    domain,
    182188                    mean_stage= project.tide,
    183189                    time_thinning=project.time_thinning,
    184                     default_boundary=Bs,
     190                    default_boundary=Bd,
    185191                    use_cache=True,
    186192                    verbose=True,
     
    188194
    189195domain.set_boundary({'back': Br,
    190                      'side': Bs,
     196                     'side': Bd,
    191197                     'ocean': Bf})
    192198
Note: See TracChangeset for help on using the changeset viewer.