Changeset 6617


Ignore:
Timestamp:
Mar 25, 2009, 4:26:41 PM (15 years ago)
Author:
myall
Message:

putting in the option for a internal polygon with same triangle area (to test discretisation issues) in project scripts.
adding required scripts for broome

Location:
anuga_work/production/australia_ph2
Files:
2 added
10 edited

Legend:

Unmodified
Added
Removed

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

    r6595 r6617  
    2626# One or all can be changed each time the run_model script is executed
    2727
    28 central_meridian = 137.5 # This is the nonstandard projection needed
     28central_meridian = 137.5 # Central meridian for projection (optional)
    2929tide = 0                # difference between MSL and HAT
    3030event_number = 64322    # the event number
    … …  
    4747                        #    final - fine mesh, slowest
    4848
     49internal_polygon = False
     50
    4951#-------------------------------------------------------------------------------
    5052# Output filename
    … …  
    5860#-------------------------------------------------------------------------------
    5961
    60 output_comment = [setup, tide, event_number, index, wave]
     62if internal_polygon:
     63    internal_poly_comment = 'internal'
     64else:
     65    internal_poly_comment = None
     66
     67output_comment = [setup, tide, event_number, index, wave, internal_poly_comment]
    6168
    6269#-------------------------------------------------------------------------------
    … …  
    8996# Format for points easting,northing (no header)                   
    9097interior_regions_data = []
     98
     99# add an internal polygon to force different mesh generation
     100# used to test for discretisation error when building elevation
     101# make sure file is in same folder as interior regions and bouding polygon;
     102# format is same (2 column .csv; easting, northing; no header)
     103if internal_polygon:
     104    interior_regions_data.append(['internal_polygon.csv',
     105                                  bounding_polygon_maxarea])
    91106
    92107# LAND - used to set the initial stage/water to be offcoast only

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

    r6579 r6617  
    2626# Model specific parameters.
    2727# One or all can be changed each time the run_model script is executed
    28 tide = 0.00001                # difference between MSL and HAT
    29 zone = 50               # specify zone of model
     28tide = 0                # difference between MSL and HAT
    3029event_number = 64344    # the event number or the mux file name
    3130# selected events: 27319 (java), 64344 (sandwich), 46697 (midamerica),
    … …  
    4746                        #    basic - coarse mesh
    4847                        #    final - fine mesh, slowest
     48
     49internal_polygon = True
    4950
    5051#-------------------------------------------------------------------------------
    … …  
    5960#-------------------------------------------------------------------------------
    6061
    61 output_comment = [setup, tide, event_number, index, wave]
     62if internal_polygon:
     63    internal_poly_comment = 'internal'
     64else:
     65    internal_poly_comment = None
     66
     67output_comment = [setup, tide, event_number, index, wave, internal_poly_comment]
     68
    6269
    6370#-------------------------------------------------------------------------------
    … …  
    8895# Format for points easting,northing (no header)                   
    8996##interior_regions_data = []
    90 interior_regions_data = [['coast_to_20km_polygon.csv',
    91                             bounding_polygon_maxarea]]
     97interior_regions_data = []
     98
     99# add an internal polygon to force different mesh generation
     100# used to test for discretisation error when building elevation
     101# make sure file is in same folder as interior regions and bouding polygon;
     102# format is same (2 column .csv; easting, northing; no header)
     103if internal_polygon:
     104    interior_regions_data.append(['coast_to_20km_polygon.csv',
     105                                  bounding_polygon_maxarea])
    92106
    93107# LAND - used to set the initial stage/water to be offcoast only

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

    r6601 r6617  
    2727# One or all can be changed each time the run_model script is executed
    2828tide = 0                # difference between MSL and HAT
    29 zone = 51               # specify zone of model
     29central_meridian = 123.0 # central meridian for projection (optional)
    3030event_number = 70468    # (sumba) the event number or the mux file name
    3131alpha = 0.1             # smoothing parameter for mesh
    … …  
    4343                        #    basic - coarse mesh
    4444                        #    final - fine mesh, slowest
     45
     46internal_polygon = False
    4547
    4648#-------------------------------------------------------------------------------
    … …  
    5557#-------------------------------------------------------------------------------
    5658
    57 output_comment = [setup, tide, event_number, index, wave]
     59if internal_polygon:
     60    internal_poly_comment = 'internal'
     61else:
     62    internal_poly_comment = None
     63
     64output_comment = [setup, tide, event_number, index, wave, internal_poly_comment]
    5865
    5966#-------------------------------------------------------------------------------
    … …  
    7986# Format for points easting,northing (no header)                   
    8087interior_regions_data = []
     88
     89# add an internal polygon to force different mesh generation
     90# used to test for discretisation error when building elevation
     91# make sure file is in same folder as interior regions and bouding polygon;
     92# format is same (2 column .csv; easting, northing; no header)
     93if internal_polygon:
     94    interior_regions_data.append(['internal_polygon.csv',
     95                                  bounding_polygon_maxarea])
    8196
    8297# LAND - used to set the initial stage/water to be offcoast only

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

    r6601 r6617  
    2727# One or all can be changed each time the run_model script is executed
    2828tide = 0                # difference between MSL and HAT
    29 zone = 49               # specify zone of model
     29central_meridian = 111.0 # Central meridian for projection (optional)
    3030event_number = 27285    # the event number or the mux file name
    3131# events: 27285 (java), 68794 (sumatra), 70375 (sumba)
    … …  
    4343                        #    basic - coarse mesh
    4444                        #    final - fine mesh, slowest
     45
     46internal_polygon = False
    4547
    4648#-------------------------------------------------------------------------------
    … …  
    5557#-------------------------------------------------------------------------------
    5658
    57 output_comment = [setup, tide, event_number, index, wave]
     59if internal_polygon:
     60    internal_poly_comment = 'internal'
     61else:
     62    internal_poly_comment = None
     63
     64output_comment = [setup, tide, event_number, index, wave, internal_poly_comment]
    5865
    5966#-------------------------------------------------------------------------------
    … …  
    8592# Format for points easting,northing (no header)                   
    8693interior_regions_data = []
     94
     95# add an internal polygon to force different mesh generation
     96# used to test for discretisation error when building elevation
     97# make sure file is in same folder as interior regions and bouding polygon;
     98# format is same (2 column .csv; easting, northing; no header)
     99if internal_polygon:
     100    interior_regions_data.append(['internal_polygon.csv',
     101                                  bounding_polygon_maxarea])
    87102
    88103# LAND - used to set the initial stage/water to be offcoast only

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

    r6610 r6617  
    4949                        #    final - fine mesh, slowest
    5050
     51internal_polygon = False
     52
    5153#-------------------------------------------------------------------------------
    5254# Output filename
    … …  
    6062#-------------------------------------------------------------------------------
    6163
    62 output_comment = [setup, tide, event_number, index, wave,]
     64if internal_polygon:
     65    internal_poly_comment = 'internal'
     66else:
     67    internal_poly_comment = None
     68
     69output_comment = [setup, tide, event_number, index, wave, internal_poly_comment]
    6370
    6471#-------------------------------------------------------------------------------
    … …  
    8996# Format for points easting,northing (no header)                   
    9097interior_regions_data = []
     98
     99# add an internal polygon to force different mesh generation
     100# used to test for discretisation error when building elevation
     101# make sure file is in same folder as interior regions and bouding polygon;
     102# format is same (2 column .csv; easting, northing; no header)
     103if internal_polygon:
     104    interior_regions_data.append(['internal_polygon.csv',
     105                                  bounding_polygon_maxarea])
    91106
    92107# LAND - used to set the initial stage/water to be offcoast only

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

    r6579 r6617  
    2727# One or all can be changed each time the run_model script is executed
    2828tide = 0                # difference between MSL and HAT
    29 zone = 52               # specify zone of model
     29central_meridian = 129.0 # for projection (necessary for nonstandard)
    3030event_number = 70906    # the event number or the mux file name
    3131# contributing events: 17898 (Flores), 64977 (Seram), 70906 (Sumba_normal),
    … …  
    4848# for 1860 - point to north of Darwin, large wave from event off Timor
    4949# event_number 71065    #(timor - mag 7.4, prob 1e-4)
     50
     51internal_polygon = False
    5052
    5153#-------------------------------------------------------------------------------
    … …  
    6062#-------------------------------------------------------------------------------
    6163
    62 output_comment = [setup, tide, event_number, index, wave]
     64if internal_polygon:
     65    internal_poly_comment = 'internal'
     66else:
     67    internal_poly_comment = None
     68
     69output_comment = [setup, tide, event_number, index, wave, internal_poly_comment]
    6370
    6471#-------------------------------------------------------------------------------
    … …  
    9198#interior_regions_data = []
    9299interior_regions_data = [['poly_island.csv',200000000]]  # large triangles on island
     100
     101# add an internal polygon to force different mesh generation
     102# used to test for discretisation error when building elevation
     103# make sure file is in same folder as interior regions and bouding polygon;
     104# format is same (2 column .csv; easting, northing; no header)
     105if internal_polygon:
     106    interior_regions_data.append(['internal_polygon.csv',
     107                                  bounding_polygon_maxarea])
    93108
    94109# LAND - used to set the initial stage/water to be offcoast only

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

    r6484 r6617  
    2727# One or all can be changed each time the run_model script is executed
    2828tide = 0                # difference between MSL and HAT
    29 zone = 51               # specify zone of model
     29central_meridian = 123.0 # Central meridian for projection (optional)
    3030event_number = 58367    # the event number or the mux file name
    3131# main events: 27329 (java), 64449 (sandwich), 58367 (puysegur)
    … …  
    4646                        #    final - fine mesh, slowest
    4747
     48internal_polygon = False
     49
    4850#-------------------------------------------------------------------------------
    4951# Output filename
    … …  
    5759#-------------------------------------------------------------------------------
    5860
    59 output_comment = [setup, tide, event_number, index, wave]
     61if internal_polygon:
     62    internal_poly_comment = 'internal'
     63else:
     64    internal_poly_comment = None
     65
     66output_comment = [setup, tide, event_number, index, wave, internal_poly_comment]
    6067
    6168#-------------------------------------------------------------------------------
    … …  
    8794# Format for points easting,northing (no header)                   
    8895interior_regions_data = []
     96
     97# add an internal polygon to force different mesh generation
     98# used to test for discretisation error when building elevation
     99# make sure file is in same folder as interior regions and bouding polygon;
     100# format is same (2 column .csv; easting, northing; no header)
     101if internal_polygon:
     102    interior_regions_data.append(['internal_polygon.csv',
     103                                  bounding_polygon_maxarea])
    89104
    90105# LAND - used to set the initial stage/water to be offcoast only

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

    r6581 r6617  
    2727# One or all can be changed each time the run_model script is executed
    2828tide = 0                # difference between MSL and HAT
    29 zone = 52               # specify zone of model
     29central_meridian = 129.0 # Central meridian for projection (optional)
    3030event_number = 64449    # the event number or the mux file name
    3131# sources: 64449 (sandwich), 68755 (sumatra), 58332 (Puysegur)
    … …  
    4646                        #    final - fine mesh, slowest
    4747
     48internal_polygon = False
     49
    4850#-------------------------------------------------------------------------------
    4951# Output filename
    … …  
    5759#-------------------------------------------------------------------------------
    5860
    59 output_comment = [setup, tide, event_number, index, wave]
     61if internal_polygon:
     62    internal_poly_comment = 'internal'
     63else:
     64    internal_poly_comment = None
     65
     66output_comment = [setup, tide, event_number, index, wave, internal_poly_comment]
    6067
    6168#-------------------------------------------------------------------------------
    … …  
    8693# Format for points easting,northing (no header)                   
    8794interior_regions_data = []
     95
     96# add an internal polygon to force different mesh generation
     97# used to test for discretisation error when building elevation
     98# make sure file is in same folder as interior regions and bouding polygon;
     99# format is same (2 column .csv; easting, northing; no header)
     100if internal_polygon:
     101    interior_regions_data.append(['internal_polygon.csv',
     102                                  bounding_polygon_maxarea])
    88103
    89104# LAND - used to set the initial stage/water to be offcoast only

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

    r6549 r6617  
    2727# One or all can be changed each time the run_model script is executed
    2828tide = 0                # difference between MSL and HAT
    29 zone = 50               # specify zone of model
     29central_meridian = 117.0 # Central meridian for projection (optional)
    3030event_number = 27240    # the event number or the mux file name
    3131# events: 11126 (andaman), 27240 (java), 64454 (sandwich)
    … …  
    4646                        #    final - fine mesh, slowest
    4747
     48internal_polygon = False
     49
    4850#-------------------------------------------------------------------------------
    4951# Output filename
    … …  
    5759#-------------------------------------------------------------------------------
    5860
    59 output_comment = [setup, tide, event_number, index, wave]
     61if internal_polygon:
     62    internal_poly_comment = 'internal'
     63else:
     64    internal_poly_comment = None
     65
     66output_comment = [setup, tide, event_number, index, wave, internal_poly_comment]
    6067
    6168#-------------------------------------------------------------------------------
    … …  
    8693# Format for points easting,northing (no header)                   
    8794interior_regions_data = []
     95
     96# add an internal polygon to force different mesh generation
     97# used to test for discretisation error when building elevation
     98# make sure file is in same folder as interior regions and bouding polygon;
     99# format is same (2 column .csv; easting, northing; no header)
     100if internal_polygon:
     101    interior_regions_data.append(['internal_polygon.csv',
     102                                  bounding_polygon_maxarea])
    88103
    89104# LAND - used to set the initial stage/water to be offcoast only

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

    r6601 r6617  
    2727# One or all can be changed each time the run_model script is executed
    2828tide = 0                # difference between MSL and HAT
    29 zone = 52               # specify zone of model
     29central_meridian = 129.0 # Central meridian for projection (optional)
    3030event_number = 64962    # the event number or the mux file name
    3131# contributing events: 17918 (Flores), 64962 (Seram), 70920 (Sumba_normal),
    … …  
    4444                        #    basic - coarse mesh
    4545                        #    final - fine mesh, slowest
     46
     47internal_polygon = False
    4648
    4749#-------------------------------------------------------------------------------
    … …  
    5658#-------------------------------------------------------------------------------
    5759
    58 output_comment = [setup, tide, event_number, index, wave]
     60if internal_polygon:
     61    internal_poly_comment = 'internal'
     62else:
     63    internal_poly_comment = None
     64
     65output_comment = [setup, tide, event_number, index, wave, internal_poly_comment]
    5966
    6067#-------------------------------------------------------------------------------
    … …  
    8693# Format for points easting,northing (no header)                   
    8794interior_regions_data = []
     95
     96# add an internal polygon to force different mesh generation
     97# used to test for discretisation error when building elevation
     98# make sure file is in same folder as interior regions and bouding polygon;
     99# format is same (2 column .csv; easting, northing; no header)
     100if internal_polygon:
     101    interior_regions_data.append(['internal_polygon.csv',
     102                                  bounding_polygon_maxarea])
    88103
    89104# LAND - used to set the initial stage/water to be offcoast only
Note: See TracChangeset for help on using the changeset viewer.