Changeset 6755


Ignore:
Timestamp:
Apr 8, 2009, 7:48:21 AM (16 years ago)
Author:
myall
Message:
 
Location:
anuga_work/production/australia_ph2
Files:
8 added
10 edited

Legend:

Unmodified
Added
Removed

  • anuga_work/production/australia_ph2/adelaide/project.py

    r6753 r6755  
    3232    event_number = int(sys.argv[1])
    3333else:   
    34     event_number = 46591    # the event number or the mux file name
     34    event_number = 64322    # the event number or the mux file name
    3535
    3636event_number_list = [27346, 58334, 64322] # To piggy back multiple events
    … …  
    8181# Used in build_elevation.py
    8282# Format for ascii grids, as produced in ArcGIS + a projection file
    83 ##ascii_grid_filenames = ['busselton_v2',   # Topo
    84 ##                        'grid_250m_clip'] # Busselton Topo
     83ascii_grid_filenames = []
    8584
    8685# Format for points is x,y,elevation (with header)

  • anuga_work/production/australia_ph2/albany/project.py

    r6724 r6755  
    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 = 'albany'
    20 scenario_folder = scenario_name
    21 
    2220
    2321#-------------------------------------------------------------------------------
    … …  
    2725# Model specific parameters.
    2826# One or all can be changed each time the run_model script is executed
     27
     28central_meridian = 117.0 # Central meridian for projection (optional)
     29#zone =
     30import sys
     31if len(sys.argv) > 1:
     32    event_number = int(sys.argv[1])
     33else:   
     34    event_number = 64344    # the event number or the mux file name
     35
     36event_number_list = [27319, 64344, 46697] # To piggy back multiple events
     37
    2938tide = 0                # difference between MSL and HAT
    30 central_meridian = 117.0 # central meridian for projection (optional)
    31 event_number = 64344    # the event number or the mux file name
    32 # selected events: 27319 (java), 64344 (sandwich), 46697 (midamerica),
    33 # other possible sources:
    34 # 7873 (altiplano), 11061 (andaman),
    35 # 53986 (puna), 57486 (peru), 58353 (puysegur), 68821 (sumatra)
    36 
    3739alpha = 0.1             # smoothing parameter for mesh
    3840friction=0.01           # manning's friction coefficient
    3941starttime=0             # start time for simulation
    40 finaltime=1000          # final time for simulation
    41 
    42 # index is only used when wave = Tb
    43 index = 1330            # index from the PTHA - 1 in 2000 wave 0.389m
    44 wave = 'Tb'             # Bf (sts wave) Tb (index wave)
    45 
    46 setup = 'basic'         # This can be one of three values
     42finaltime=1000   #60000 # final time for simulation
     43setup = 'trial'         # This can be one of three values
    4744                        #    trial - coarsest mesh, fast
    4845                        #    basic - coarse mesh
    4946                        #    final - fine mesh, slowest
     47
     48# index is only used when wave = Tb
     49index = 1330            # index from the PTHA
     50wave = 'Tb'             # Bf (sts wave) Tb (index wave)
     51
    5052
    5153internal_polygon = False
    … …  
    6567    internal_poly_comment = 'internal'
    6668else:
    67     internal_poly_comment = None
     69    internal_poly_comment = ''
    6870
    6971output_comment = [setup, tide, event_number, index, wave, internal_poly_comment]
    70 
    7172
    7273#-------------------------------------------------------------------------------
    … …  
    7980ascii_grid_filenames = [] # 250m grid 2005
    8081
    81 # Format for point is x,y,elevation (with header)
    82 point_filenames = ['grid_250m.txt'] # 250m grid 2005
     82# Format for points is x,y,elevation (with header)
     83point_filenames = ['grid_250m.txt']   # 250m grid 2005
    8384
    8485### Add csv header list to all files in point_filenames
    … …  
    8788##    add_csv_header(join(topographies_folder, f), headerlist)
    8889
    89     # BOUNDING POLYGON - for data clipping and estimate of triangles in mesh
     90# BOUNDING POLYGON - for data clipping and estimate of triangles in mesh
    9091# Used in build_elevation.py
    91 # Format for points easting,northing (no header)
     92# Format for points: easting,northing (no header)
    9293bounding_polygon_filename = 'bounding_polygon.csv'
    9394bounding_polygon_maxarea = 100000
    … …  
    9697# Used in run_model.py
    9798# Format for points easting,northing (no header)                   
    98 ##interior_regions_data = []
    9999interior_regions_data = []
    100100
    … …  
    104104# format is same (2 column .csv; easting, northing; no header)
    105105if internal_polygon:
    106     interior_regions_data.append(['coast_to_20km_polygon.csv',
     106    interior_regions_data.append(['internal_polygon.csv',
    107107                                  bounding_polygon_maxarea])
    108108
    … …  
    119119# Used in run_building_inundation.py
    120120# Format latitude,longitude etc (geographic)
    121 building_exposure_filename = 'busselton_res_clip.csv' # from NEXIS
     121building_exposure_filename = '.csv' # from NEXIS
    122122
    123123# BOUNDING POLYGON - used in build_boundary.py and run_model.py respectively
    … …  
    137137# If a meta-file from EventSelection is used, set 'multi-mux' to True.
    138138# If a single MUX stem filename (*.grd) is used, set 'multi-mux' to False.
    139 ##mux_input_filename = event_number # to be found in event_folder
    140                                     # (ie boundaries/event_number/)
     139##mux_input_filename = 'Java-0016-z.grd'
    141140##multi_mux = False
    142141mux_input_filename = 'event.list'
    143142multi_mux = True
    144 
    145143
    146144# Specify if share cache is to be used
    … …  
    168166# Path to MUX data
    169167ENV_MUXHOME = 'MUXHOME'
    170 
    171168
    172169#-------------------------------------------------------------------------------
    … …  
    198195    from anuga.caching import caching
    199196    caching.set_option('cachedir', cachedir)
    200 
    201197   
    202198# check various directories/files that must exist
    203 anuga_folder = join(home, state, scenario_folder, 'anuga')
     199anuga_folder = join(home, state, scenario_name, 'anuga')
    204200topographies_folder = join(anuga_folder, 'topographies')
    205201polygons_folder = join(anuga_folder, 'polygons')
    … …  
    258254
    259255# sts gauge with desired index number - used for wave 'Tb'
    260 boundary_csv = join(event_folder, 'sts_gauge_' + str(index) +'.csv')
    261 
     256boundary_csv = join(event_folder, 'sts_gauge_' + str(index) + '.csv')
     257
     258

  • anuga_work/production/australia_ph2/broome/project.py

    r6731 r6755  
    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 = 'broome'
    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
     27
     28central_meridian = 123.0 # Central meridian for projection (optional)
     29#zone =
     30import sys
     31if len(sys.argv) > 1:
     32    event_number = int(sys.argv[1])
     33else:   
     34    event_number = 70869    # the event number or the mux file name
     35
     36event_number_list = [27337, 70336, 70869] # To piggy back multiple events
     37
    2838tide = 0                # difference between MSL and HAT
    29 central_meridian = 123.0 # central meridian for projection (optional)
    30 event_number = 70468    # (sumba) the event number or the mux file name
    3139alpha = 0.1             # smoothing parameter for mesh
    3240friction=0.01           # manning's friction coefficient
    3341starttime=0             # start time for simulation
    34 finaltime=80000          # final time for simulation
    35 
    36 # index is only used when wave = Tb
    37 ##index = 1693            # index from the PTHA - Y2000 5.4259m
    38 index = 1661            # index from the PTHA - TO DO
    39 
    40 wave = 'Tb'             # Bf (sts wave) Tb (index wave)
    41 
    42 
    43 setup = 'final'         # This can be one of three values
     42finaltime=1000   #60000 # final time for simulation
     43setup = 'trial'         # This can be one of three values
    4444                        #    trial - coarsest mesh, fast
    4545                        #    basic - coarse mesh
    4646                        #    final - fine mesh, slowest
     47
     48# index is only used when wave = Tb
     49index = 1661            # index from the PTHA - Y2000 0.257m
     50wave = 'Tb'             # Bf (sts wave) Tb (index wave)
     51
    4752
    4853internal_polygon = False
    … …  
    6267    internal_poly_comment = 'internal'
    6368else:
    64     internal_poly_comment = None
     69    internal_poly_comment = ''
    6570
    6671output_comment = [setup, tide, event_number, index, wave, internal_poly_comment]
    … …  
    7378# Used in build_elevation.py
    7479# Format for ascii grids, as produced in ArcGIS + a projection file
    75 ascii_grid_filenames = [] # Grid data
    76 
    77 # Format for point is x,y,elevation (with header)
    78 point_filenames = ['grid.txt']            # 250m grid 2005
     80ascii_grid_filenames = ['grid250m_simple']
     81
     82# Format for points is x,y,elevation (with header)
     83point_filenames = []   # 250m grid 2005
     84
     85### Add csv header list to all files in point_filenames
     86##headerlist = ['x', 'y', 'elevation']
     87##for f in point_filenames:
     88##    add_csv_header(join(topographies_folder, f), headerlist)
    7989
    8090# BOUNDING POLYGON - for data clipping and estimate of triangles in mesh
    8191# Used in build_elevation.py
    82 # Format for points easting,northing (no header)
     92# Format for points: easting,northing (no header)
    8393bounding_polygon_filename = 'bounding_polygon.csv'
    8494bounding_polygon_maxarea = 100000
    … …  
    109119# Used in run_building_inundation.py
    110120# Format latitude,longitude etc (geographic)
    111 building_exposure_filename = 'busselton_res_clip.csv' # from NEXIS
     121building_exposure_filename = '.csv' # from NEXIS
    112122
    113123# BOUNDING POLYGON - used in build_boundary.py and run_model.py respectively
    … …  
    118128# Thinned ordering file from Hazard Map (geographic)
    119129# Format is index,latitude,longitude (with header)
    120 urs_order_filename = 'urs_order.csv'
     130urs_order_filename = 'urs_order_simple.csv'
    121131
    122132# Landward bounding points
    … …  
    127137# If a meta-file from EventSelection is used, set 'multi-mux' to True.
    128138# If a single MUX stem filename (*.grd) is used, set 'multi-mux' to False.
    129 ##mux_input_filename = event_number # to be found in event_folder
    130                                     # (ie boundaries/event_number/)
     139##mux_input_filename = 'Java-0016-z.grd'
    131140##multi_mux = False
    132141mux_input_filename = 'event.list'
    133142multi_mux = True
     143
     144# Specify if share cache is to be used
     145# Whatever is specified here will be relative to INUNDATION_HOME/.cache
     146# If nothing is specified, local cache will be used.
     147cachedir = '.python_cache_phII'
    134148
    135149
    … …  
    174188home = join(os.getenv(ENV_INUNDATIONHOME), 'data') # Absolute path for data folder
    175189muxhome = os.getenv(ENV_MUXHOME)
     190
     191# Create absolute pathname for cache directory
     192# and change caching to use it
     193if 'cachedir' in dir():
     194    cachedir = join(os.getenv(ENV_INUNDATIONHOME), '.cache', cachedir)
     195    from anuga.caching import caching
     196    caching.set_option('cachedir', cachedir)
    176197   
    177198# check various directories/files that must exist
    178 anuga_folder = join(home, state, scenario_folder, 'anuga')
     199anuga_folder = join(home, state, scenario_name, 'anuga')
    179200topographies_folder = join(anuga_folder, 'topographies')
    180201polygons_folder = join(anuga_folder, 'polygons')
    … …  
    234255# sts gauge with desired index number - used for wave 'Tb'
    235256boundary_csv = join(event_folder, 'sts_gauge_' + str(index) + '.csv')
     257
     258

  • anuga_work/production/australia_ph2/darwin/project.py

    r6741 r6755  
    128128# Thinned ordering file from Hazard Map (geographic)
    129129# Format is index,latitude,longitude (with header)
    130 urs_order_filename = 'urs_order.csv'
     130urs_order_filename = 'urs_order_simple.csv'
    131131
    132132# Landward bounding points

  • anuga_work/production/australia_ph2/esperance/project.py

    r6617 r6755  
    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 = 'esperance'
    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
     27
     28central_meridian = 123.0 # Central meridian for projection (optional)
     29#zone =
     30import sys
     31if len(sys.argv) > 1:
     32    event_number = int(sys.argv[1])
     33else:   
     34    event_number = 64449    # the event number or the mux file name
     35
     36event_number_list = [27329, 58367, 64449] # To piggy back multiple events
     37
    2838tide = 0                # difference between MSL and HAT
    29 central_meridian = 123.0 # Central meridian for projection (optional)
    30 event_number = 58367    # the event number or the mux file name
    31 # main events: 27329 (java), 64449 (sandwich), 58367 (puysegur)
    32 # others: 7858 (altiplano), 31951 (kermandec), 46483 (midamerica),
    33 # 53960 (Puna), 57616 (peru), 65356 (shetland)
    3439alpha = 0.1             # smoothing parameter for mesh
    3540friction=0.01           # manning's friction coefficient
    3641starttime=0             # start time for simulation
    37 finaltime=80000          # final time for simulation
    38 
    39 # index is only used when wave = Tb
    40 index = 1699            # index from the PTHA
    41 wave = 'Tb'             # Bf (sts wave) Tb (index wave)
    42 
    43 setup = 'final'         # This can be one of three values
     42finaltime=1000   #60000 # final time for simulation
     43setup = 'trial'         # This can be one of three values
    4444                        #    trial - coarsest mesh, fast
    4545                        #    basic - coarse mesh
    4646                        #    final - fine mesh, slowest
     47
     48# index is only used when wave = Tb
     49index = 1699            # index from the PTHA - Y2000 0.257m
     50wave = 'Tb'             # Bf (sts wave) Tb (index wave)
     51
    4752
    4853internal_polygon = False
    … …  
    6267    internal_poly_comment = 'internal'
    6368else:
    64     internal_poly_comment = None
     69    internal_poly_comment = ''
    6570
    6671output_comment = [setup, tide, event_number, index, wave, internal_poly_comment]
    … …  
    7580ascii_grid_filenames = [] # 250m grid 2005
    7681
    77 # Format for point is x,y,elevation (with header)
     82# Format for points is x,y,elevation (with header)
    7883point_filenames = ['grid250m_pts.txt',
    7984                   'coast_250.txt'] # 250m grid 2005
    … …  
    8489##    add_csv_header(join(topographies_folder, f), headerlist)
    8590
    86     # BOUNDING POLYGON - for data clipping and estimate of triangles in mesh
     91# BOUNDING POLYGON - for data clipping and estimate of triangles in mesh
    8792# Used in build_elevation.py
    88 # Format for points easting,northing (no header)
     93# Format for points: easting,northing (no header)
    8994bounding_polygon_filename = 'bounding_polygon.csv'
    9095bounding_polygon_maxarea = 100000
    … …  
    115120# Used in run_building_inundation.py
    116121# Format latitude,longitude etc (geographic)
    117 building_exposure_filename = 'busselton_res_clip.csv' # from NEXIS
     122building_exposure_filename = '.csv' # from NEXIS
    118123
    119124# BOUNDING POLYGON - used in build_boundary.py and run_model.py respectively
    … …  
    133138# If a meta-file from EventSelection is used, set 'multi-mux' to True.
    134139# If a single MUX stem filename (*.grd) is used, set 'multi-mux' to False.
    135 ##mux_input_filename = event_number # to be found in event_folder
    136                                     # (ie boundaries/event_number/)
     140##mux_input_filename = 'Java-0016-z.grd'
    137141##multi_mux = False
    138142mux_input_filename = 'event.list'
    139143multi_mux = True
     144
     145# Specify if share cache is to be used
     146# Whatever is specified here will be relative to INUNDATION_HOME/.cache
     147# If nothing is specified, local cache will be used.
     148cachedir = '.python_cache_phII'
    140149
    141150
    … …  
    180189home = join(os.getenv(ENV_INUNDATIONHOME), 'data') # Absolute path for data folder
    181190muxhome = os.getenv(ENV_MUXHOME)
     191
     192# Create absolute pathname for cache directory
     193# and change caching to use it
     194if 'cachedir' in dir():
     195    cachedir = join(os.getenv(ENV_INUNDATIONHOME), '.cache', cachedir)
     196    from anuga.caching import caching
     197    caching.set_option('cachedir', cachedir)
    182198   
    183199# check various directories/files that must exist
    184 anuga_folder = join(home, state, scenario_folder, 'anuga')
     200anuga_folder = join(home, state, scenario_name, 'anuga')
    185201topographies_folder = join(anuga_folder, 'topographies')
    186202polygons_folder = join(anuga_folder, 'polygons')
    … …  
    239255
    240256# sts gauge with desired index number - used for wave 'Tb'
    241 boundary_csv = join(event_folder, 'sts_gauge_' + str(index) +'.csv')
     257boundary_csv = join(event_folder, 'sts_gauge_' + str(index) + '.csv')
     258
     259

  • anuga_work/production/australia_ph2/eucla_motel/project.py

    r6617 r6755  
    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 = 'eucla_motel'
    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
     27
     28central_meridian = 129.0 # Central meridian for projection (optional)
     29#zone =
     30import sys
     31if len(sys.argv) > 1:
     32    event_number = int(sys.argv[1])
     33else:   
     34    event_number = 64449    # the event number or the mux file name
     35
     36event_number_list = [58332, 64449, 68755] # To piggy back multiple events
     37
    2838tide = 0                # difference between MSL and HAT
    29 central_meridian = 129.0 # Central meridian for projection (optional)
    30 event_number = 64449    # the event number or the mux file name
    31 # sources: 64449 (sandwich), 68755 (sumatra), 58332 (Puysegur)
    32 # others: 7859 (altiplano), 11016 (andaman), 17017 (columbia), 27252 (java), 46699 (midamerica)
    33 # 53961 (puna), 57582 (peru), 58332 (Puysegur), 64449 (sandwich), 68755 (sumatra)
    3439alpha = 0.1             # smoothing parameter for mesh
    3540friction=0.01           # manning's friction coefficient
    3641starttime=0             # start time for simulation
    37 finaltime=80000          # final time for simulation
    38 
    39 # index is only used when wave = Tb
    40 index = 1837            # index from the PTHA
    41 wave = 'Bf'             # Bf (sts wave) Tb (index wave)
    42 
    43 setup = 'final'         # This can be one of three values
     42finaltime=1000   #60000 # final time for simulation
     43setup = 'trial'         # This can be one of three values
    4444                        #    trial - coarsest mesh, fast
    4545                        #    basic - coarse mesh
    4646                        #    final - fine mesh, slowest
     47
     48# index is only used when wave = Tb
     49index = 1837            # index from the PTHA - Y2000 0.257m
     50wave = 'Tb'             # Bf (sts wave) Tb (index wave)
     51
    4752
    4853internal_polygon = False
    … …  
    6267    internal_poly_comment = 'internal'
    6368else:
    64     internal_poly_comment = None
     69    internal_poly_comment = ''
    6570
    6671output_comment = [setup, tide, event_number, index, wave, internal_poly_comment]
    … …  
    7378# Used in build_elevation.py
    7479# Format for ascii grids, as produced in ArcGIS + a projection file
    75 ascii_grid_filenames = [] # 250m grid 2005
    76 
    77 # Format for point is x,y,elevation (with header)
    78 point_filenames = ['grid250m_pts.txt'] # 250m grid 2005
     80ascii_grid_filenames = []
     81
     82# Format for points is x,y,elevation (with header)
     83point_filenames = ['grid250m_pts.txt']   # 250m grid 2005
    7984
    8085### Add csv header list to all files in point_filenames
    … …  
    8388##    add_csv_header(join(topographies_folder, f), headerlist)
    8489
    85     # BOUNDING POLYGON - for data clipping and estimate of triangles in mesh
     90# BOUNDING POLYGON - for data clipping and estimate of triangles in mesh
    8691# Used in build_elevation.py
    87 # Format for points easting,northing (no header)
     92# Format for points: easting,northing (no header)
    8893bounding_polygon_filename = 'bounding_polygon.csv'
    8994bounding_polygon_maxarea = 100000
    … …  
    114119# Used in run_building_inundation.py
    115120# Format latitude,longitude etc (geographic)
    116 building_exposure_filename = 'busselton_res_clip.csv' # from NEXIS
     121building_exposure_filename = '.csv' # from NEXIS
    117122
    118123# BOUNDING POLYGON - used in build_boundary.py and run_model.py respectively
    … …  
    132137# If a meta-file from EventSelection is used, set 'multi-mux' to True.
    133138# If a single MUX stem filename (*.grd) is used, set 'multi-mux' to False.
    134 ##mux_input_filename = event_number # to be found in event_folder
    135                                     # (ie boundaries/event_number/)
     139##mux_input_filename = 'Java-0016-z.grd'
    136140##multi_mux = False
    137141mux_input_filename = 'event.list'
    138142multi_mux = True
     143
     144# Specify if share cache is to be used
     145# Whatever is specified here will be relative to INUNDATION_HOME/.cache
     146# If nothing is specified, local cache will be used.
     147cachedir = '.python_cache_phII'
    139148
    140149
    … …  
    179188home = join(os.getenv(ENV_INUNDATIONHOME), 'data') # Absolute path for data folder
    180189muxhome = os.getenv(ENV_MUXHOME)
     190
     191# Create absolute pathname for cache directory
     192# and change caching to use it
     193if 'cachedir' in dir():
     194    cachedir = join(os.getenv(ENV_INUNDATIONHOME), '.cache', cachedir)
     195    from anuga.caching import caching
     196    caching.set_option('cachedir', cachedir)
    181197   
    182198# check various directories/files that must exist
    183 anuga_folder = join(home, state, scenario_folder, 'anuga')
     199anuga_folder = join(home, state, scenario_name, 'anuga')
    184200topographies_folder = join(anuga_folder, 'topographies')
    185201polygons_folder = join(anuga_folder, 'polygons')
    … …  
    239255# sts gauge with desired index number - used for wave 'Tb'
    240256boundary_csv = join(event_folder, 'sts_gauge_' + str(index) + '.csv')
     257
     258

  • anuga_work/production/australia_ph2/hobart/project.py

    r6716 r6755  
    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 = 'hobart'
    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
     27
     28central_meridian = 147.0 # Central meridian for projection (optional)
     29#zone =
     30import sys
     31if len(sys.argv) > 1:
     32    event_number = int(sys.argv[1])
     33else:   
     34    event_number = 58082    # the event number or the mux file name
     35
     36event_number_list = [51268, 58082, 64309] # To piggy back multiple events
     37
    2838tide = 0                # difference between MSL and HAT
    29 central_meridian = 147.0 # Central meridian for projection (optional)
    30 
    3139alpha = 0.1             # smoothing parameter for mesh
    3240friction=0.01           # manning's friction coefficient
    3341starttime=0             # start time for simulation
    34 finaltime=80000          # final time for simulation
    35 
    36 # index is only used when wave = Tb
    37 wave = 'Tb'             # Bf (sts wave) Tb (index wave)
    38 
    39 # for melbourne (west) side:
    40 ##event_number = 58349    # the event number or the mux file name
    41 # contributing sources:
    42 # 51328 (newhebrides), 58349 (puysegur), 64535 (sandwich)
    43 # other possible: 31966 (kermandec), 65350 (shetland)
    44 ##index = 1983            # index from the PTHA - west side, close to melbourne; 1983 has 1 in 2000 wave height 0.534m
    45 
    46 # for east, Hobart, side - point opposite hobart
    47 event_number = 64461    # the event number or the mux file name
    48 # contributing events:  51450 (newHebrides 7.3e-6), 58078 (puysegur 5.2e-6), 64461 (sandwich 8.2e-6)
    49 # Others: 31816 (kermadec 4.7e-7), 46651 (midamerica 3.5e-7),
    50 #
    51 index = 2479            # index from PTHA - Y2000 1.0977m
    52 
    53 # index = 2460   # for Bichen
    54 
    55 setup = 'final'         # This can be one of three values
     42finaltime=1000   #60000 # final time for simulation
     43setup = 'trial'         # This can be one of three values
    5644                        #    trial - coarsest mesh, fast
    5745                        #    basic - coarse mesh
    5846                        #    final - fine mesh, slowest
     47
     48# index is only used when wave = Tb
     49index = 2528            # index from the PTHA - Y2000 0.257m
     50wave = 'Tb'             # Bf (sts wave) Tb (index wave)
     51
    5952
    6053internal_polygon = False
    … …  
    7467    internal_poly_comment = 'internal'
    7568else:
    76     internal_poly_comment = None
     69    internal_poly_comment = ''
    7770
    7871output_comment = [setup, tide, event_number, index, wave, internal_poly_comment]
    … …  
    8780ascii_grid_filenames = [] # 250m grid 2005
    8881
    89 # Format for point is x,y,elevation (with header)
    90 point_filenames = ['grid_250m.txt'] # 250m grid 2005
     82# Format for points is x,y,elevation (with header)
     83point_filenames = ['grid_250m.txt']   # 250m grid 2005
    9184
    9285### Add csv header list to all files in point_filenames
    … …  
    9588##    add_csv_header(join(topographies_folder, f), headerlist)
    9689
    97     # BOUNDING POLYGON - for data clipping and estimate of triangles in mesh
     90# BOUNDING POLYGON - for data clipping and estimate of triangles in mesh
    9891# Used in build_elevation.py
    99 # Format for points easting,northing (no header)
     92# Format for points: easting,northing (no header)
    10093bounding_polygon_filename = 'bounding_polygon.csv'
    10194bounding_polygon_maxarea = 100000
    … …  
    121114# Used in get_timeseries.py. 
    122115# Format easting,northing,name,elevation (with header)
    123 gauges_filename = 'MH_gauges_10m.csv'
     116gauges_filename = 'gauges.csv'
    124117
    125118# BUILDINGS EXPOSURE - for identifying inundated houses
    126119# Used in run_building_inundation.py
    127120# Format latitude,longitude etc (geographic)
    128 building_exposure_filename = 'busselton_res_clip.csv' # from NEXIS
     121building_exposure_filename = '.csv' # from NEXIS
    129122
    130123# BOUNDING POLYGON - used in build_boundary.py and run_model.py respectively
    … …  
    135128# Thinned ordering file from Hazard Map (geographic)
    136129# Format is index,latitude,longitude (with header)
    137 urs_order_filename = 'urs_order.csv'
     130urs_order_filename = 'urs_order_simple.csv'
    138131
    139132# Landward bounding points
    … …  
    144137# If a meta-file from EventSelection is used, set 'multi-mux' to True.
    145138# If a single MUX stem filename (*.grd) is used, set 'multi-mux' to False.
    146 ##mux_input_filename = event_number # to be found in event_folder
    147                                     # (ie boundaries/event_number/)
     139##mux_input_filename = 'Java-0016-z.grd'
    148140##multi_mux = False
    149141mux_input_filename = 'event.list'
    150142multi_mux = True
     143
     144# Specify if share cache is to be used
     145# Whatever is specified here will be relative to INUNDATION_HOME/.cache
     146# If nothing is specified, local cache will be used.
     147cachedir = '.python_cache_phII'
    151148
    152149
    … …  
    191188home = join(os.getenv(ENV_INUNDATIONHOME), 'data') # Absolute path for data folder
    192189muxhome = os.getenv(ENV_MUXHOME)
     190
     191# Create absolute pathname for cache directory
     192# and change caching to use it
     193if 'cachedir' in dir():
     194    cachedir = join(os.getenv(ENV_INUNDATIONHOME), '.cache', cachedir)
     195    from anuga.caching import caching
     196    caching.set_option('cachedir', cachedir)
    193197   
    194198# check various directories/files that must exist
    195 anuga_folder = join(home, state, scenario_folder, 'anuga')
     199anuga_folder = join(home, state, scenario_name, 'anuga')
    196200topographies_folder = join(anuga_folder, 'topographies')
    197201polygons_folder = join(anuga_folder, 'polygons')
    … …  
    250254
    251255# sts gauge with desired index number - used for wave 'Tb'
    252 boundary_csv = join(event_folder, 'sts_gauge_' + str(index) +'.csv')
     256boundary_csv = join(event_folder, 'sts_gauge_' + str(index) + '.csv')
     257
     258

  • anuga_work/production/australia_ph2/perth/project.py

    r6617 r6755  
    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 = 'perth'
    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
     27
     28central_meridian = 117.0 # Central meridian for projection (optional)
     29#zone =
     30import sys
     31if len(sys.argv) > 1:
     32    event_number = int(sys.argv[1])
     33else:   
     34    event_number = 70921    # the event number or the mux file name
     35
     36event_number_list = [64479, 70921] # To piggy back multiple events
     37
    2838tide = 0                # difference between MSL and HAT
    29 central_meridian = 117.0 # Central meridian for projection (optional)
    30 event_number = 27240    # the event number or the mux file name
    31 # events: 11126 (andaman), 27240 (java), 64454 (sandwich)
    32 # 68840 (sumatra), 70903 (sumba_normal)
    33 # other events: 13215 (arakan),
    3439alpha = 0.1             # smoothing parameter for mesh
    3540friction=0.01           # manning's friction coefficient
    3641starttime=0             # start time for simulation
    37 finaltime=80000          # final time for simulation
    38 
    39 # index is only used when wave = Tb
    40 index = 1039            # index from the PTHA - 1 in 2000 yr 0.5357m
    41 wave = 'Bf'             # Bf (sts wave) Tb (index wave)
    42 
    43 setup = 'final'         # This can be one of three values
     42finaltime=1000   #60000 # final time for simulation
     43setup = 'trial'         # This can be one of three values
    4444                        #    trial - coarsest mesh, fast
    4545                        #    basic - coarse mesh
    4646                        #    final - fine mesh, slowest
     47
     48# index is only used when wave = Tb
     49index = 803            # index from the PTHA - Y2000 0.257m
     50wave = 'Tb'             # Bf (sts wave) Tb (index wave)
     51
    4752
    4853internal_polygon = False
    … …  
    6267    internal_poly_comment = 'internal'
    6368else:
    64     internal_poly_comment = None
     69    internal_poly_comment = ''
    6570
    6671output_comment = [setup, tide, event_number, index, wave, internal_poly_comment]
    … …  
    7378# Used in build_elevation.py
    7479# Format for ascii grids, as produced in ArcGIS + a projection file
    75 ascii_grid_filenames = [] # 250m grid 2005
    76 
    77 # Format for point is x,y,elevation (with header)
    78 point_filenames = ['grid_250m_small.txt'] # 250m grid 2005
     80ascii_grid_filenames = ['grid250m_simple'] # 250m grid 2005
     81
     82# Format for points is x,y,elevation (with header)
     83point_filenames = []  # 250m grid 2005
    7984
    8085### Add csv header list to all files in point_filenames
    … …  
    8388##    add_csv_header(join(topographies_folder, f), headerlist)
    8489
    85     # BOUNDING POLYGON - for data clipping and estimate of triangles in mesh
     90# BOUNDING POLYGON - for data clipping and estimate of triangles in mesh
    8691# Used in build_elevation.py
    87 # Format for points easting,northing (no header)
     92# Format for points: easting,northing (no header)
    8893bounding_polygon_filename = 'bounding_polygon.csv'
    8994bounding_polygon_maxarea = 100000
    … …  
    114119# Used in run_building_inundation.py
    115120# Format latitude,longitude etc (geographic)
    116 building_exposure_filename = 'busselton_res_clip.csv' # from NEXIS
     121building_exposure_filename = '.csv' # from NEXIS
    117122
    118123# BOUNDING POLYGON - used in build_boundary.py and run_model.py respectively
    … …  
    123128# Thinned ordering file from Hazard Map (geographic)
    124129# Format is index,latitude,longitude (with header)
    125 urs_order_filename = 'urs_order.csv'
     130urs_order_filename = 'urs_order_simple.csv'
    126131
    127132# Landward bounding points
    … …  
    132137# If a meta-file from EventSelection is used, set 'multi-mux' to True.
    133138# If a single MUX stem filename (*.grd) is used, set 'multi-mux' to False.
    134 ##mux_input_filename = event_number # to be found in event_folder
    135                                     # (ie boundaries/event_number/)
     139##mux_input_filename = 'Java-0016-z.grd'
    136140##multi_mux = False
    137141mux_input_filename = 'event.list'
    138142multi_mux = True
     143
     144# Specify if share cache is to be used
     145# Whatever is specified here will be relative to INUNDATION_HOME/.cache
     146# If nothing is specified, local cache will be used.
     147cachedir = '.python_cache_phII'
    139148
    140149
    … …  
    179188home = join(os.getenv(ENV_INUNDATIONHOME), 'data') # Absolute path for data folder
    180189muxhome = os.getenv(ENV_MUXHOME)
     190
     191# Create absolute pathname for cache directory
     192# and change caching to use it
     193if 'cachedir' in dir():
     194    cachedir = join(os.getenv(ENV_INUNDATIONHOME), '.cache', cachedir)
     195    from anuga.caching import caching
     196    caching.set_option('cachedir', cachedir)
    181197   
    182198# check various directories/files that must exist
    183 anuga_folder = join(home, state, scenario_folder, 'anuga')
     199anuga_folder = join(home, state, scenario_name, 'anuga')
    184200topographies_folder = join(anuga_folder, 'topographies')
    185201polygons_folder = join(anuga_folder, 'polygons')
    … …  
    238254
    239255# sts gauge with desired index number - used for wave 'Tb'
    240 boundary_csv = join(event_folder, 'sts_gauge_' + str(index) +'.csv')
     256boundary_csv = join(event_folder, 'sts_gauge_' + str(index) + '.csv')
     257
     258

  • anuga_work/production/australia_ph2/portland/project.py

    r6716 r6755  
    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 = 'portland'
    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
     27
     28central_meridian = 141.0 # Central meridian for projection (optional)
     29#zone =
     30import sys
     31if len(sys.argv) > 1:
     32    event_number = int(sys.argv[1])
     33else:   
     34    event_number = 64322    # the event number or the mux file name
     35
     36event_number_list = [27289, 58358, 64322] # To piggy back multiple events
     37
    2838tide = 0                # difference between MSL and HAT
    29 central_meridian = 141.0 # Central meridian for projection (optional)
    30 event_number = 64409    # the event number or the mux file name
    31 # main events: 64409 (sandwich 4.9e-6)
    32 # others:
    33 
    3439alpha = 0.1             # smoothing parameter for mesh
    3540friction=0.01           # manning's friction coefficient
    3641starttime=0             # start time for simulation
    37 finaltime=80000          # final time for simulation
    38 
    39 # index is only used when wave = Tb
    40 index = 1970            # index from the PTHA - Y2000 is 0.43074m
    41 wave = 'Bf'             # Bf (sts wave) Tb (index wave)
    42 
    43 setup = 'final'         # This can be one of three values
     42finaltime=1000   #60000 # final time for simulation
     43setup = 'trial'         # This can be one of three values
    4444                        #    trial - coarsest mesh, fast
    4545                        #    basic - coarse mesh
    4646                        #    final - fine mesh, slowest
     47
     48# index is only used when wave = Tb
     49index = 1969            # index from the PTHA
     50wave = 'Tb'             # Bf (sts wave) Tb (index wave)
     51
    4752
    4853internal_polygon = False
    … …  
    6267    internal_poly_comment = 'internal'
    6368else:
    64     internal_poly_comment = None
     69    internal_poly_comment = ''
    6570
    6671output_comment = [setup, tide, event_number, index, wave, internal_poly_comment]
    … …  
    7580ascii_grid_filenames = [] # 250m grid 2005
    7681
    77 # Format for point is x,y,elevation (with header)
     82# Format for points is x,y,elevation (with header)
    7883point_filenames = ['bight_grid250m.txt',
    7984                   'melbourne_grid250m.txt'] # 250m grid 2005
    … …  
    8489##    add_csv_header(join(topographies_folder, f), headerlist)
    8590
    86     # BOUNDING POLYGON - for data clipping and estimate of triangles in mesh
     91# BOUNDING POLYGON - for data clipping and estimate of triangles in mesh
    8792# Used in build_elevation.py
    88 # Format for points easting,northing (no header)
     93# Format for points: easting,northing (no header)
    8994bounding_polygon_filename = 'bounding_polygon.csv'
    9095bounding_polygon_maxarea = 100000
    … …  
    115120# Used in run_building_inundation.py
    116121# Format latitude,longitude etc (geographic)
    117 building_exposure_filename = 'busselton_res_clip.csv' # from NEXIS
     122building_exposure_filename = '.csv' # from NEXIS
    118123
    119124# BOUNDING POLYGON - used in build_boundary.py and run_model.py respectively
    … …  
    124129# Thinned ordering file from Hazard Map (geographic)
    125130# Format is index,latitude,longitude (with header)
    126 urs_order_filename = 'urs_order.csv'
     131urs_order_filename = 'urs_order_simple.csv'
    127132
    128133# Landward bounding points
    … …  
    133138# If a meta-file from EventSelection is used, set 'multi-mux' to True.
    134139# If a single MUX stem filename (*.grd) is used, set 'multi-mux' to False.
    135 ##mux_input_filename = event_number # to be found in event_folder
    136                                     # (ie boundaries/event_number/)
     140##mux_input_filename = 'Java-0016-z.grd'
    137141##multi_mux = False
    138142mux_input_filename = 'event.list'
    139143multi_mux = True
     144
     145# Specify if share cache is to be used
     146# Whatever is specified here will be relative to INUNDATION_HOME/.cache
     147# If nothing is specified, local cache will be used.
     148cachedir = '.python_cache_phII'
    140149
    141150
    … …  
    180189home = join(os.getenv(ENV_INUNDATIONHOME), 'data') # Absolute path for data folder
    181190muxhome = os.getenv(ENV_MUXHOME)
     191
     192# Create absolute pathname for cache directory
     193# and change caching to use it
     194if 'cachedir' in dir():
     195    cachedir = join(os.getenv(ENV_INUNDATIONHOME), '.cache', cachedir)
     196    from anuga.caching import caching
     197    caching.set_option('cachedir', cachedir)
    182198   
    183199# check various directories/files that must exist
    184 anuga_folder = join(home, state, scenario_folder, 'anuga')
     200anuga_folder = join(home, state, scenario_name, 'anuga')
    185201topographies_folder = join(anuga_folder, 'topographies')
    186202polygons_folder = join(anuga_folder, 'polygons')
    … …  
    239255
    240256# sts gauge with desired index number - used for wave 'Tb'
    241 boundary_csv = join(event_folder, 'sts_gauge_' + str(index) +'.csv')
     257boundary_csv = join(event_folder, 'sts_gauge_' + str(index) + '.csv')
     258
     259

  • anuga_work/production/australia_ph2/strahan/project.py

    r6672 r6755  
    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 = 'strahan'
    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
     27
     28central_meridian = 147.0 # Central meridian for projection (optional)
     29#zone =
     30import sys
     31if len(sys.argv) > 1:
     32    event_number = int(sys.argv[1])
     33else:   
     34    event_number = 58337    # the event number or the mux file name
     35
     36event_number_list = [58337, 64214, 68779] # To piggy back multiple events
     37
    2838tide = 0                # difference between MSL and HAT
    29 central_meridian = 147.0 # Central meridian for projection (optional)
    30 event_number = 68779    # the event number or the mux file name
    31 # main events: 64214 (sandwich 1.4e-5), 58337 (puysegur 3e-6), 68779 (sumatra 7.53-7)
    32 # others: 27295 (java 3e-7), 31817 (kermadec 5e-7), 46559 (midAmerica 4.2e-7)
    33 # 51460 (newHebrides 7e-7), 57477 (peru 5.4e-7), 65322 (shetland 3.6e-7),
    34 
    3539alpha = 0.1             # smoothing parameter for mesh
    3640friction=0.01           # manning's friction coefficient
    3741starttime=0             # start time for simulation
    38 finaltime=80000          # final time for simulation
    39 
    40 # index is only used when wave = Tb
    41 index = 2044            # index from the PTHA - Y2000 is 0.60374m
    42 wave = 'Bf'             # Bf (sts wave) Tb (index wave)
    43 
    44 setup = 'final'         # This can be one of three values
     42finaltime=1000   #60000 # final time for simulation
     43setup = 'trial'         # This can be one of three values
    4544                        #    trial - coarsest mesh, fast
    4645                        #    basic - coarse mesh
    4746                        #    final - fine mesh, slowest
     47
     48# index is only used when wave = Tb
     49index = 2044            # index from the PTHA - Y2000 0.257m
     50wave = 'Tb'             # Bf (sts wave) Tb (index wave)
     51
    4852
    4953internal_polygon = False
    … …  
    6367    internal_poly_comment = 'internal'
    6468else:
    65     internal_poly_comment = None
     69    internal_poly_comment = ''
    6670
    6771output_comment = [setup, tide, event_number, index, wave, internal_poly_comment]
    … …  
    7680ascii_grid_filenames = [] # 250m grid 2005
    7781
    78 # Format for point is x,y,elevation (with header)
    79 point_filenames = ['grid_250m.txt'] # 250m grid 2005
     82# Format for points is x,y,elevation (with header)
     83point_filenames = ['grid_250m.txt']   # 250m grid 2005
    8084
    8185### Add csv header list to all files in point_filenames
    … …  
    8488##    add_csv_header(join(topographies_folder, f), headerlist)
    8589
    86     # BOUNDING POLYGON - for data clipping and estimate of triangles in mesh
     90# BOUNDING POLYGON - for data clipping and estimate of triangles in mesh
    8791# Used in build_elevation.py
    88 # Format for points easting,northing (no header)
     92# Format for points: easting,northing (no header)
    8993bounding_polygon_filename = 'bounding_polygon.csv'
    9094bounding_polygon_maxarea = 100000
    … …  
    115119# Used in run_building_inundation.py
    116120# Format latitude,longitude etc (geographic)
    117 building_exposure_filename = 'busselton_res_clip.csv' # from NEXIS
     121building_exposure_filename = '.csv' # from NEXIS
    118122
    119123# BOUNDING POLYGON - used in build_boundary.py and run_model.py respectively
    … …  
    124128# Thinned ordering file from Hazard Map (geographic)
    125129# Format is index,latitude,longitude (with header)
    126 urs_order_filename = 'urs_order.csv'
     130urs_order_filename = 'urs_order_simple.csv'
    127131
    128132# Landward bounding points
    … …  
    133137# If a meta-file from EventSelection is used, set 'multi-mux' to True.
    134138# If a single MUX stem filename (*.grd) is used, set 'multi-mux' to False.
    135 ##mux_input_filename = event_number # to be found in event_folder
    136                                     # (ie boundaries/event_number/)
     139##mux_input_filename = 'Java-0016-z.grd'
    137140##multi_mux = False
    138141mux_input_filename = 'event.list'
    139142multi_mux = True
     143
     144# Specify if share cache is to be used
     145# Whatever is specified here will be relative to INUNDATION_HOME/.cache
     146# If nothing is specified, local cache will be used.
     147cachedir = '.python_cache_phII'
    140148
    141149
    … …  
    180188home = join(os.getenv(ENV_INUNDATIONHOME), 'data') # Absolute path for data folder
    181189muxhome = os.getenv(ENV_MUXHOME)
     190
     191# Create absolute pathname for cache directory
     192# and change caching to use it
     193if 'cachedir' in dir():
     194    cachedir = join(os.getenv(ENV_INUNDATIONHOME), '.cache', cachedir)
     195    from anuga.caching import caching
     196    caching.set_option('cachedir', cachedir)
    182197   
    183198# check various directories/files that must exist
    184 anuga_folder = join(home, state, scenario_folder, 'anuga')
     199anuga_folder = join(home, state, scenario_name, 'anuga')
    185200topographies_folder = join(anuga_folder, 'topographies')
    186201polygons_folder = join(anuga_folder, 'polygons')
    … …  
    239254
    240255# sts gauge with desired index number - used for wave 'Tb'
    241 boundary_csv = join(event_folder, 'sts_gauge_' + str(index) +'.csv')
     256boundary_csv = join(event_folder, 'sts_gauge_' + str(index) + '.csv')
     257
     258
Note: See TracChangeset for help on using the changeset viewer.