Changeset 6769


Ignore:
Timestamp:
Apr 8, 2009, 1:18:00 PM (15 years ago)
Author:
myall
Message:

running GBR models without LADS data

Location:
anuga_work/production/australia_ph2
Files:
9 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/australia_ph2/albany/project.py

    r6767 r6769  
    4141starttime=0             # start time for simulation
    4242finaltime=1000   #60000 # final time for simulation
    43 setup = 'trial'         # This can be one of three values
     43setup = 'final'         # This can be one of three values
    4444                        #    trial - coarsest mesh, fast
    4545                        #    basic - coarse mesh
  • anuga_work/production/australia_ph2/bamaga/project.py

    r6742 r6769  
    88from time import localtime, strftime, gmtime
    99from os.path import join, exists
    10 
     10from anuga.lib.add_csv_header.add_csv_header import add_csv_header
    1111
    1212#-------------------------------------------------------------------------------
     
    1818state = 'australia_ph2'
    1919scenario_name = 'bamaga'
    20 scenario_folder = scenario_name
    2120
    2221#-------------------------------------------------------------------------------
     
    2625# Model specific parameters.
    2726# One or all can be changed each time the run_model script is executed
    28 tide = 0                # difference between MSL and HAT
    29 ##zone = 53               # specify zone of model
    30 central_meridian = 143.5 # nonstandard projection
    31 
    32 
     27
     28central_meridian = 143.5 # Central meridian for projection (optional)
     29zone = None
    3330import sys
    3431if len(sys.argv) > 1:
     
    3633else:   
    3734    event_number = 31975    # the event number or the mux file name
    38    
     35
    3936event_number_list = [31975, 51390, 63735] # To piggy back multiple events
    40                                           # using run_multiple_events.py
    41 
     37
     38tide = 0                # difference between MSL and HAT
    4239alpha = 0.1             # smoothing parameter for mesh
    4340friction=0.01           # manning's friction coefficient
    4441starttime=0             # start time for simulation
    45 finaltime=80000          # final time for simulation
    46 
    47 # index is only used when wave = Tb
    48 index = 1889            # index from the PTHA
    49 wave = 'Tb'             # Bf (sts wave) Tb (index wave)
    50 
     42finaltime=1000   #60000 # final time for simulation
    5143setup = 'trial'         # This can be one of three values
    5244                        #    trial - coarsest mesh, fast
    5345                        #    basic - coarse mesh
    5446                        #    final - fine mesh, slowest
     47
     48# index is only used when wave = Tb
     49index = 2000            # index from the PTHA - Y2000 0.257m
     50wave = 'Tb'             # Bf (sts wave) Tb (index wave)
     51
    5552
    5653internal_polygon = False
     
    7067    internal_poly_comment = 'internal'
    7168else:
    72     internal_poly_comment = None
     69    internal_poly_comment = ''
    7370
    7471output_comment = [setup, tide, event_number, index, wave, internal_poly_comment]
     
    8380ascii_grid_filenames = ['grid250m_pro'] # 250m grid 2005
    8481
    85 # Format for point is x,y,elevation (with header)
    86 ##point_filenames = []
    87 point_filenames = ['ahs100_0.txt',
    88                    'ahs100_1.txt',
    89                    'ahs100_2.txt',
    90                    'ahs100_3.txt',
    91                    'ahs100_4.txt',
    92                    'ahs100_5.txt',
    93                    'ahs100_6.txt',
    94                    'ahs100_7.txt',
    95                    'ahs100_8.txt',
    96                    'ahs100_9.txt',
    97                    'ahs100_10.txt',
    98                    'ahs100_11.txt',
    99                    'ahs100_12.txt',
    100                    'ahs100_13.txt'] # LADS data
     82# Format for points is x,y,elevation (with header)
     83point_filenames = []
     84##point_filenames = ['ahs100_0.txt',
     85##                   'ahs100_1.txt',
     86##                   'ahs100_2.txt',
     87##                   'ahs100_3.txt',
     88##                   'ahs100_4.txt',
     89##                   'ahs100_5.txt',
     90##                   'ahs100_6.txt',
     91##                   'ahs100_7.txt',
     92##                   'ahs100_8.txt',
     93##                   'ahs100_9.txt',
     94##                   'ahs100_10.txt',
     95##                   'ahs100_11.txt',
     96##                   'ahs100_12.txt',
     97##                   'ahs100_13.txt'] # LADS data
     98
     99### Add csv header list to all files in point_filenames
     100##headerlist = ['x', 'y', 'elevation']
     101##for f in point_filenames:
     102##    add_csv_header(join(topographies_folder, f), headerlist)
    101103
    102104# BOUNDING POLYGON - for data clipping and estimate of triangles in mesh
    103105# Used in build_elevation.py
    104 # Format for points easting,northing (no header)
     106# Format for points: easting,northing (no header)
    105107bounding_polygon_filename = 'bounding_polygon.csv'
    106108bounding_polygon_maxarea = 125000
     
    131133# Used in run_building_inundation.py
    132134# Format latitude,longitude etc (geographic)
    133 building_exposure_filename = 'busselton_res_clip.csv' # from NEXIS
     135building_exposure_filename = '.csv' # from NEXIS
    134136
    135137# BOUNDING POLYGON - used in build_boundary.py and run_model.py respectively
     
    149151# If a meta-file from EventSelection is used, set 'multi-mux' to True.
    150152# If a single MUX stem filename (*.grd) is used, set 'multi-mux' to False.
    151 ##mux_input_filename = event_number # to be found in event_folder
    152                                     # (ie boundaries/event_number/)
     153##mux_input_filename = 'Java-0016-z.grd'
    153154##multi_mux = False
    154155mux_input_filename = 'event.list'
    155156multi_mux = True
     157
     158# Specify if share cache is to be used
     159# Whatever is specified here will be relative to INUNDATION_HOME/.cache
     160# If nothing is specified, local cache will be used.
     161cachedir = '.python_cache_phII'
    156162
    157163
     
    196202home = join(os.getenv(ENV_INUNDATIONHOME), 'data') # Absolute path for data folder
    197203muxhome = os.getenv(ENV_MUXHOME)
     204
     205# Create absolute pathname for cache directory
     206# and change caching to use it
     207if 'cachedir' in dir():
     208    cachedir = join(os.getenv(ENV_INUNDATIONHOME), '.cache', cachedir)
     209    from anuga.caching import caching
     210    caching.set_option('cachedir', cachedir)
    198211   
    199212# check various directories/files that must exist
    200 anuga_folder = join(home, state, scenario_folder, 'anuga')
     213anuga_folder = join(home, state, scenario_name, 'anuga')
    201214topographies_folder = join(anuga_folder, 'topographies')
    202215polygons_folder = join(anuga_folder, 'polygons')
     
    235248event_sts = join(event_folder, scenario_name)
    236249
    237 
    238 
    239250# The absolute pathname for the output folder names
    240251# Used for build_elevation.py
  • anuga_work/production/australia_ph2/bundaberg/project.py

    r6756 r6769  
    2626# One or all can be changed each time the run_model script is executed
    2727
    28 central_meridian = 153.0 # Central meridian for projection (optional)
    29 #zone = 56
     28central_meridian = None # Central meridian for projection (optional)
     29zone = 56
    3030import sys
    3131if len(sys.argv) > 1:
     
    8181
    8282# Format for points is x,y,elevation (with header)
    83 point_filenames = ['ahs100_0.txt',
    84                    'ahs100_1.txt',
    85                    'ahs100_2.txt',
    86                    'ahs100_3.txt',
    87                    'ahs100_4.txt',
    88                    'ahs100_5.txt',
    89                    'ahs100_6.txt',
    90                    'ahs100_7.txt',
    91                    'ahs100_8.txt',
    92                    'ahs100_9.txt',
    93                    'ahs100_10.txt',
    94                    'ahs100_11.txt',
    95                    'ahs100_12.txt',
    96                    'ahs100_13.txt'] # LADS data
     83point_filenames = []
     84##point_filenames = ['ahs100_0.txt',
     85##                   'ahs100_1.txt',
     86##                   'ahs100_2.txt',
     87##                   'ahs100_3.txt',
     88##                   'ahs100_4.txt',
     89##                   'ahs100_5.txt',
     90##                   'ahs100_6.txt',
     91##                   'ahs100_7.txt',
     92##                   'ahs100_8.txt',
     93##                   'ahs100_9.txt',
     94##                   'ahs100_10.txt',
     95##                   'ahs100_11.txt',
     96##                   'ahs100_12.txt',
     97##                   'ahs100_13.txt'] # LADS data
    9798
    9899### Add csv header list to all files in point_filenames
  • anuga_work/production/australia_ph2/cairns/project.py

    r6756 r6769  
    2626# One or all can be changed each time the run_model script is executed
    2727
    28 central_meridian = 147 # Central meridian for projection (optional)
    29 #zone = 55
     28central_meridian = None # Central meridian for projection (optional)
     29zone = 55
    3030import sys
    3131if len(sys.argv) > 1:
    3232    event_number = int(sys.argv[1])
    3333else:   
    34     event_number =     # the event number or the mux file name
    35 
    36 event_number_list = [] # To piggy back multiple events
     34    event_number = 31975    # the event number or the mux file name
     35
     36event_number_list = [31975, 51469, 63737] # To piggy back multiple events
    3737
    3838tide = 0                # difference between MSL and HAT
     
    4747
    4848# index is only used when wave = Tb
    49 index =             # index from the PTHA - Y2000 0.257m
     49index = 2172            # index from the PTHA - Y2000 0.257m
    5050wave = 'Tb'             # Bf (sts wave) Tb (index wave)
    5151
     
    8181
    8282# Format for points is x,y,elevation (with header)
    83 point_filenames = ['ahs100_0.txt',
    84                    'ahs100_1.txt',
    85                    'ahs100_2.txt',
    86                    'ahs100_3.txt',
    87                    'ahs100_4.txt',
    88                    'ahs100_5.txt',
    89                    'ahs100_6.txt',
    90                    'ahs100_7.txt',
    91                    'ahs100_8.txt',
    92                    'ahs100_9.txt',
    93                    'ahs100_10.txt',
    94                    'ahs100_11.txt',
    95                    'ahs100_12.txt',
    96                    'ahs100_13.txt'] # LADS data
     83point_filenames = []
     84##point_filenames = ['ahs100_0.txt',
     85##                   'ahs100_1.txt',
     86##                   'ahs100_2.txt',
     87##                   'ahs100_3.txt',
     88##                   'ahs100_4.txt',
     89##                   'ahs100_5.txt',
     90##                   'ahs100_6.txt',
     91##                   'ahs100_7.txt',
     92##                   'ahs100_8.txt',
     93##                   'ahs100_9.txt',
     94##                   'ahs100_10.txt',
     95##                   'ahs100_11.txt',
     96##                   'ahs100_12.txt',
     97##                   'ahs100_13.txt'] # LADS data
    9798
    9899### Add csv header list to all files in point_filenames
  • anuga_work/production/australia_ph2/cooktown/project.py

    r6756 r6769  
    2626# One or all can be changed each time the run_model script is executed
    2727
    28 central_meridian = 147.0 # Central meridian for projection (optional)
    29 #zone = 55
     28central_meridian = None # Central meridian for projection (optional)
     29zone = 55
    3030import sys
    3131if len(sys.argv) > 1:
     
    7878# Used in build_elevation.py
    7979# Format for ascii grids, as produced in ArcGIS + a projection file
    80 ascii_grid_filenames = []
     80ascii_grid_filenames = ['grid250m_pro']
    8181
    8282# Format for points is x,y,elevation (with header)
    83 point_filenames = ['ahs100_0.txt',
    84                    'ahs100_1.txt',
    85                    'ahs100_2.txt',
    86                    'ahs100_3.txt',
    87                    'ahs100_4.txt',
    88                    'ahs100_5.txt',
    89                    'ahs100_6.txt',
    90                    'ahs100_7.txt',
    91                    'ahs100_8.txt',
    92                    'ahs100_9.txt',
    93                    'ahs100_10.txt',
    94                    'ahs100_11.txt',
    95                    'ahs100_12.txt',
    96                    'ahs100_13.txt'] # LADS data
     83point_filenames = []
     84##point_filenames = ['ahs100_0.txt',
     85##                   'ahs100_1.txt',
     86##                   'ahs100_2.txt',
     87##                   'ahs100_3.txt',
     88##                   'ahs100_4.txt',
     89##                   'ahs100_5.txt',
     90##                   'ahs100_6.txt',
     91##                   'ahs100_7.txt',
     92##                   'ahs100_8.txt',
     93##                   'ahs100_9.txt',
     94##                   'ahs100_10.txt',
     95##                   'ahs100_11.txt',
     96##                   'ahs100_12.txt',
     97##                   'ahs100_13.txt'] # LADS data
    9798
    9899### Add csv header list to all files in point_filenames
  • anuga_work/production/australia_ph2/lockhart_river/project.py

    r6746 r6769  
    88from time import localtime, strftime, gmtime
    99from os.path import join, exists
    10 
     10from anuga.lib.add_csv_header.add_csv_header import add_csv_header
    1111
    1212#-------------------------------------------------------------------------------
     
    1818state = 'australia_ph2'
    1919scenario_name = 'lockhart_river'
    20 scenario_folder = scenario_name
    2120
    2221#-------------------------------------------------------------------------------
     
    2625# Model specific parameters.
    2726# One or all can be changed each time the run_model script is executed
    28 tide = 0                # difference between MSL and HAT
    29 ##zone = 55               # specify zone of model
    30 central_meridian = 143.5 # nonstandard projection
    31 
    32 
     27
     28central_meridian = 143.5 # Central meridian for projection (optional)
     29zone = None
    3330import sys
    3431if len(sys.argv) > 1:
     
    3633else:   
    3734    event_number = 51316    # the event number or the mux file name
    38    
     35
    3936event_number_list = [51316, 62786] # To piggy back multiple events
    40                                           # using run_multiple_events.py
    41 
     37
     38tide = 0                # difference between MSL and HAT
    4239alpha = 0.1             # smoothing parameter for mesh
    4340friction=0.01           # manning's friction coefficient
    4441starttime=0             # start time for simulation
    45 finaltime=80000          # final time for simulation
    46 
    47 # index is only used when wave = Tb
    48 index = 1889            # index from the PTHA
    49 wave = 'Tb'             # Bf (sts wave) Tb (index wave)
    50 
     42finaltime=1000   #60000 # final time for simulation
    5143setup = 'trial'         # This can be one of three values
    5244                        #    trial - coarsest mesh, fast
    5345                        #    basic - coarse mesh
    5446                        #    final - fine mesh, slowest
     47
     48# index is only used when wave = Tb
     49index = 1994            # index from the PTHA - Y2000 0.257m
     50wave = 'Tb'             # Bf (sts wave) Tb (index wave)
     51
    5552
    5653internal_polygon = False
     
    7067    internal_poly_comment = 'internal'
    7168else:
    72     internal_poly_comment = None
     69    internal_poly_comment = ''
    7370
    7471output_comment = [setup, tide, event_number, index, wave, internal_poly_comment]
     
    8380ascii_grid_filenames = ['grid250m_pro'] # 250m grid 2005
    8481
    85 # Format for point is x,y,elevation (with header)
    86 ##point_filenames = []
    87 point_filenames = ['ahs100_0.txt',
    88                    'ahs100_1.txt',
    89                    'ahs100_2.txt',
    90                    'ahs100_3.txt',
    91                    'ahs100_4.txt',
    92                    'ahs100_5.txt',
    93                    'ahs100_6.txt',
    94                    'ahs100_7.txt',
    95                    'ahs100_8.txt',
    96                    'ahs100_9.txt',
    97                    'ahs100_10.txt',
    98                    'ahs100_11.txt',
    99                    'ahs100_12.txt',
    100                    'ahs100_13.txt'] # LADS data
     82# Format for points is x,y,elevation (with header)
     83point_filenames = []
     84##point_filenames = ['ahs100_0.txt',
     85##                   'ahs100_1.txt',
     86##                   'ahs100_2.txt',
     87##                   'ahs100_3.txt',
     88##                   'ahs100_4.txt',
     89##                   'ahs100_5.txt',
     90##                   'ahs100_6.txt',
     91##                   'ahs100_7.txt',
     92##                   'ahs100_8.txt',
     93##                   'ahs100_9.txt',
     94##                   'ahs100_10.txt',
     95##                   'ahs100_11.txt',
     96##                   'ahs100_12.txt',
     97##                   'ahs100_13.txt'] # LADS data
     98
     99### Add csv header list to all files in point_filenames
     100##headerlist = ['x', 'y', 'elevation']
     101##for f in point_filenames:
     102##    add_csv_header(join(topographies_folder, f), headerlist)
    101103
    102104# BOUNDING POLYGON - for data clipping and estimate of triangles in mesh
    103105# Used in build_elevation.py
    104 # Format for points easting,northing (no header)
     106# Format for points: easting,northing (no header)
    105107bounding_polygon_filename = 'bounding_polygon.csv'
    106108bounding_polygon_maxarea = 125000
     
    131133# Used in run_building_inundation.py
    132134# Format latitude,longitude etc (geographic)
    133 building_exposure_filename = 'busselton_res_clip.csv' # from NEXIS
     135building_exposure_filename = '.csv' # from NEXIS
    134136
    135137# BOUNDING POLYGON - used in build_boundary.py and run_model.py respectively
     
    149151# If a meta-file from EventSelection is used, set 'multi-mux' to True.
    150152# If a single MUX stem filename (*.grd) is used, set 'multi-mux' to False.
    151 ##mux_input_filename = event_number # to be found in event_folder
    152                                     # (ie boundaries/event_number/)
     153##mux_input_filename = 'Java-0016-z.grd'
    153154##multi_mux = False
    154155mux_input_filename = 'event.list'
    155156multi_mux = True
     157
     158# Specify if share cache is to be used
     159# Whatever is specified here will be relative to INUNDATION_HOME/.cache
     160# If nothing is specified, local cache will be used.
     161cachedir = '.python_cache_phII'
    156162
    157163
     
    196202home = join(os.getenv(ENV_INUNDATIONHOME), 'data') # Absolute path for data folder
    197203muxhome = os.getenv(ENV_MUXHOME)
     204
     205# Create absolute pathname for cache directory
     206# and change caching to use it
     207if 'cachedir' in dir():
     208    cachedir = join(os.getenv(ENV_INUNDATIONHOME), '.cache', cachedir)
     209    from anuga.caching import caching
     210    caching.set_option('cachedir', cachedir)
    198211   
    199212# check various directories/files that must exist
    200 anuga_folder = join(home, state, scenario_folder, 'anuga')
     213anuga_folder = join(home, state, scenario_name, 'anuga')
    201214topographies_folder = join(anuga_folder, 'topographies')
    202215polygons_folder = join(anuga_folder, 'polygons')
     
    235248event_sts = join(event_folder, scenario_name)
    236249
    237 
    238 
    239250# The absolute pathname for the output folder names
    240251# Used for build_elevation.py
  • anuga_work/production/australia_ph2/lucinda/project.py

    r6756 r6769  
    2626# One or all can be changed each time the run_model script is executed
    2727
    28 central_meridian = 147.0 # Central meridian for projection (optional)
    29 #zone = 55
     28central_meridian = None # Central meridian for projection (optional)
     29zone = 55
    3030import sys
    3131if len(sys.argv) > 1:
     
    8181
    8282# Format for points is x,y,elevation (with header)
    83 ##point_filenames = []
    84 point_filenames = ['ahs100_0.txt',
    85                    'ahs100_1.txt',
    86                    'ahs100_2.txt',
    87                    'ahs100_3.txt',
    88                    'ahs100_4.txt',
    89                    'ahs100_5.txt',
    90                    'ahs100_6.txt',
    91                    'ahs100_7.txt',
    92                    'ahs100_8.txt',
    93                    'ahs100_9.txt',
    94                    'ahs100_10.txt',
    95                    'ahs100_11.txt',
    96                    'ahs100_12.txt',
    97                    'ahs100_13.txt'] # LADS data
     83point_filenames = []
     84##point_filenames = ['ahs100_0.txt',
     85##                   'ahs100_1.txt',
     86##                   'ahs100_2.txt',
     87##                   'ahs100_3.txt',
     88##                   'ahs100_4.txt',
     89##                   'ahs100_5.txt',
     90##                   'ahs100_6.txt',
     91##                   'ahs100_7.txt',
     92##                   'ahs100_8.txt',
     93##                   'ahs100_9.txt',
     94##                   'ahs100_10.txt',
     95##                   'ahs100_11.txt',
     96##                   'ahs100_12.txt',
     97##                   'ahs100_13.txt'] # LADS data
    9898
    9999### Add csv header list to all files in point_filenames
  • anuga_work/production/australia_ph2/rockhampton/project.py

    r6756 r6769  
    2626# One or all can be changed each time the run_model script is executed
    2727
    28 central_meridian = 153.0 # Central meridian for projection (optional)
    29 #zone = 56
     28central_meridian = None # Central meridian for projection (optional)
     29zone = 56
    3030import sys
    3131if len(sys.argv) > 1:
     
    8181
    8282# Format for points is x,y,elevation (with header)
    83 ##point_filenames = []
    84 point_filenames = ['ahs100_0.txt',
    85                    'ahs100_1.txt',
    86                    'ahs100_2.txt',
    87                    'ahs100_3.txt',
    88                    'ahs100_4.txt',
    89                    'ahs100_5.txt',
    90                    'ahs100_6.txt',
    91                    'ahs100_7.txt',
    92                    'ahs100_8.txt',
    93                    'ahs100_9.txt',
    94                    'ahs100_10.txt',
    95                    'ahs100_11.txt',
    96                    'ahs100_12.txt',
    97                    'ahs100_13.txt'] # LADS data
     83point_filenames = []
     84##point_filenames = ['ahs100_0.txt',
     85##                   'ahs100_1.txt',
     86##                   'ahs100_2.txt',
     87##                   'ahs100_3.txt',
     88##                   'ahs100_4.txt',
     89##                   'ahs100_5.txt',
     90##                   'ahs100_6.txt',
     91##                   'ahs100_7.txt',
     92##                   'ahs100_8.txt',
     93##                   'ahs100_9.txt',
     94##                   'ahs100_10.txt',
     95##                   'ahs100_11.txt',
     96##                   'ahs100_12.txt',
     97##                   'ahs100_13.txt'] # LADS data
    9898
    9999### Add csv header list to all files in point_filenames
  • anuga_work/production/australia_ph2/townsville/project.py

    r6749 r6769  
    88from time import localtime, strftime, gmtime
    99from os.path import join, exists
    10 
     10from anuga.lib.add_csv_header.add_csv_header import add_csv_header
    1111
    1212#-------------------------------------------------------------------------------
     
    1818state = 'australia_ph2'
    1919scenario_name = 'townsville'
    20 scenario_folder = scenario_name
    2120
    2221#-------------------------------------------------------------------------------
     
    2625# Model specific parameters.
    2726# One or all can be changed each time the run_model script is executed
    28 tide = 0                # difference between MSL and HAT
    29 ##zone = 55               # specify zone of model
    30 central_meridian = 149 # nonstandard projection
    31 
    32 
     27
     28central_meridian = 149.0 # Central meridian for projection (optional)
     29zone = None
    3330import sys
    3431if len(sys.argv) > 1:
     
    3633else:   
    3734    event_number = 31965    # the event number or the mux file name
    38    
     35
    3936event_number_list = [31965, 51469, 65840] # To piggy back multiple events
    40                                           # using run_multiple_events.py
    41 
     37
     38tide = 0                # difference between MSL and HAT
    4239alpha = 0.1             # smoothing parameter for mesh
    4340friction=0.01           # manning's friction coefficient
    4441starttime=0             # start time for simulation
    45 finaltime=80000          # final time for simulation
    46 
    47 # index is only used when wave = Tb
    48 index = 1889            # index from the PTHA
    49 wave = 'Tb'             # Bf (sts wave) Tb (index wave)
    50 
     42finaltime=1000   #60000 # final time for simulation
    5143setup = 'trial'         # This can be one of three values
    5244                        #    trial - coarsest mesh, fast
    5345                        #    basic - coarse mesh
    5446                        #    final - fine mesh, slowest
     47
     48# index is only used when wave = Tb
     49index = 2740            # index from the PTHA - Y2000 0.257m
     50wave = 'Tb'             # Bf (sts wave) Tb (index wave)
     51
    5552
    5653internal_polygon = False
     
    7067    internal_poly_comment = 'internal'
    7168else:
    72     internal_poly_comment = None
     69    internal_poly_comment = ''
    7370
    7471output_comment = [setup, tide, event_number, index, wave, internal_poly_comment]
     
    8380ascii_grid_filenames = ['grid250m_pro'] # 250m grid 2005
    8481
    85 # Format for point is x,y,elevation (with header)
    86 ##point_filenames = []
    87 point_filenames = ['ahs100_0.txt',
    88                    'ahs100_1.txt',
    89                    'ahs100_2.txt',
    90                    'ahs100_3.txt',
    91                    'ahs100_4.txt',
    92                    'ahs100_5.txt',
    93                    'ahs100_6.txt',
    94                    'ahs100_7.txt',
    95                    'ahs100_8.txt',
    96                    'ahs100_9.txt',
    97                    'ahs100_10.txt',
    98                    'ahs100_11.txt',
    99                    'ahs100_12.txt',
    100                    'ahs100_13.txt'] # LADS data
     82# Format for points is x,y,elevation (with header)
     83point_filenames = []
     84##point_filenames = ['ahs100_0.txt',
     85##                   'ahs100_1.txt',
     86##                   'ahs100_2.txt',
     87##                   'ahs100_3.txt',
     88##                   'ahs100_4.txt',
     89##                   'ahs100_5.txt',
     90##                   'ahs100_6.txt',
     91##                   'ahs100_7.txt',
     92##                   'ahs100_8.txt',
     93##                   'ahs100_9.txt',
     94##                   'ahs100_10.txt',
     95##                   'ahs100_11.txt',
     96##                   'ahs100_12.txt',
     97##                   'ahs100_13.txt'] # LADS data
     98
     99### Add csv header list to all files in point_filenames
     100##headerlist = ['x', 'y', 'elevation']
     101##for f in point_filenames:
     102##    add_csv_header(join(topographies_folder, f), headerlist)
    101103
    102104# BOUNDING POLYGON - for data clipping and estimate of triangles in mesh
    103105# Used in build_elevation.py
    104 # Format for points easting,northing (no header)
     106# Format for points: easting,northing (no header)
    105107bounding_polygon_filename = 'bounding_polygon.csv'
    106108bounding_polygon_maxarea = 125000
     
    131133# Used in run_building_inundation.py
    132134# Format latitude,longitude etc (geographic)
    133 building_exposure_filename = 'busselton_res_clip.csv' # from NEXIS
     135building_exposure_filename = '.csv' # from NEXIS
    134136
    135137# BOUNDING POLYGON - used in build_boundary.py and run_model.py respectively
     
    149151# If a meta-file from EventSelection is used, set 'multi-mux' to True.
    150152# If a single MUX stem filename (*.grd) is used, set 'multi-mux' to False.
    151 ##mux_input_filename = event_number # to be found in event_folder
    152                                     # (ie boundaries/event_number/)
     153##mux_input_filename = 'Java-0016-z.grd'
    153154##multi_mux = False
    154155mux_input_filename = 'event.list'
    155156multi_mux = True
     157
     158# Specify if share cache is to be used
     159# Whatever is specified here will be relative to INUNDATION_HOME/.cache
     160# If nothing is specified, local cache will be used.
     161cachedir = '.python_cache_phII'
    156162
    157163
     
    196202home = join(os.getenv(ENV_INUNDATIONHOME), 'data') # Absolute path for data folder
    197203muxhome = os.getenv(ENV_MUXHOME)
     204
     205# Create absolute pathname for cache directory
     206# and change caching to use it
     207if 'cachedir' in dir():
     208    cachedir = join(os.getenv(ENV_INUNDATIONHOME), '.cache', cachedir)
     209    from anuga.caching import caching
     210    caching.set_option('cachedir', cachedir)
    198211   
    199212# check various directories/files that must exist
    200 anuga_folder = join(home, state, scenario_folder, 'anuga')
     213anuga_folder = join(home, state, scenario_name, 'anuga')
    201214topographies_folder = join(anuga_folder, 'topographies')
    202215polygons_folder = join(anuga_folder, 'polygons')
     
    235248event_sts = join(event_folder, scenario_name)
    236249
    237 
    238 
    239250# The absolute pathname for the output folder names
    240251# Used for build_elevation.py
Note: See TracChangeset for help on using the changeset viewer.