Changeset 6579


Ignore:
Timestamp:
Mar 23, 2009, 1:02:28 PM (15 years ago)
Author:
myall
Message:

putting events in for all models, and changing run_model.py to use Bt instead of Bd on the side boundaries

Location:
anuga_work/production/australia_ph2
Files:
15 edited

Legend:

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

    r6460 r6579  
    2626# One or all can be changed each time the run_model script is executed
    2727tide = 0                # difference between MSL and HAT
    28 event_number = 64481    # the event number
     28event_number = 64322    # the event number
     29# events: 7867 (altiplano), 27346 (java), 46591 (midAmerica),
     30# 64322 (sandwich), 65368 (shetland),
     31# other events: 51466 (newHebrides), 57537 (peru), 58334 (puysegur), 68791 (sumatra)
     32# 11100 (andaman),
    2933alpha = 0.1             # smoothing parameter for mesh
    3034friction=0.01           # manning's friction coefficient
     
    3337
    3438# index is only used when wave = Tb
    35 index = 1203            # index from the PTHA
     39index = 1937            # index from the PTHA - Y2000 0.257m
    3640wave = 'Tb'             # Bf (sts wave) Tb (index wave)
    3741
     
    217221mux_input = join(boundaries_folder, mux_input_filename)
    218222
    219 
    220 
     223# sts gauge with desired index number - used for wave 'Tb'
     224boundary_csv = join(event_folder, 'sts_gauge_', str(index), '.csv')
     225
     226
  • anuga_work/production/australia_ph2/adelaide/run_model.py

    r6511 r6579  
    156156                        verbose=True)
    157157    domain.set_boundary({'back': Br,
    158                      'side': Bd,
     158                     'side': Bt,
    159159                     'ocean': Bf})
    160160   
     
    163163
    164164    domain.set_boundary({'back': Br,
    165                          'side': Bd,
     165                         'side': Bt,
    166166                         'ocean': Tb})
    167167else:
  • anuga_work/production/australia_ph2/albany/export_results_max.py

    r6540 r6579  
    2222directory = project.output_folder
    2323
    24 time_dir1 = '20090310_171457_run_final_0_64344_1330_Tb_mhingee'
    25 time_dir2 = '20090310_171537_run_final_0_46697_1330_Tb_mhingee'
    26 time_dir3 = '20090310_175417_run_final_0_27319_1330_Tb_mhingee'
    27 ##time_dir4 = ''
    28 ##time_dir5 = ''
    29 ##time_dir6 = ''
    30 ##
    31 time_dirs = [time_dir1, time_dir2, time_dir3] # time_dir4, time_dir5, time_dir6]
     24time_dir1 = '20090318_171220_run_final_0_27319_1330_Tb_mhingee'
     25time_dir2 = '20090318_172003_run_final_0_64344_1330_Tb_mhingee'
     26time_dir3 = '20090318_172103_run_final_0_46697_1330_Tb_mhingee'
     27time_dir4 = '20090318_175318_run_final_1e-05_46697_1330_Tb_mhingee'
     28time_dir5 = '20090320_144313_run_final_1e-05_27319_1330_Tb_mhingee'
     29time_dir6 = '20090320_144404_run_final_1e-05_64344_1330_Tb_mhingee'
     30
     31time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
    3232 
    3333
  • anuga_work/production/australia_ph2/albany/project.py

    r6549 r6579  
    2828tide = 0.00001                # difference between MSL and HAT
    2929zone = 50               # specify zone of model
    30 event_number = 46697    # the event number or the mux file name
     30event_number = 64344    # the event number or the mux file name
    3131# selected events: 27319 (java), 64344 (sandwich), 46697 (midamerica),
    3232# other possible sources:
  • anuga_work/production/australia_ph2/albany/run_model.py

    r6511 r6579  
    156156                        verbose=True)
    157157    domain.set_boundary({'back': Br,
    158                      'side': Bd,
     158                     'side': Bt,
    159159                     'ocean': Bf})
    160160   
     
    163163
    164164    domain.set_boundary({'back': Br,
    165                          'side': Bd,
     165                         'side': Bt,
    166166                         'ocean': Tb})
    167167else:
  • anuga_work/production/australia_ph2/carnarvon/project.py

    r6460 r6579  
    3535
    3636# index is only used when wave = Tb
    37 index = 1203            # index from the PTHA
     37index = 109            # index from the PTHA - Y2000 1.6275m
    3838wave = 'Tb'             # Bf (sts wave) Tb (index wave)
    3939
  • anuga_work/production/australia_ph2/ceduna/project.py

    r6476 r6579  
    4141# index is only used when wave = Tb
    4242index = 1889            # index from the PTHA
    43 wave = 'Tb'             # Bf (sts wave) Tb (index wave)
     43wave = 'Bf'             # Bf (sts wave) Tb (index wave)
    4444
    4545setup = 'final'         # This can be one of three values
     
    5959#-------------------------------------------------------------------------------
    6060
    61 output_comment = [setup, tide, event_number, index, wave]
     61output_comment = [setup, tide, event_number, index, wave,]
    6262
    6363#-------------------------------------------------------------------------------
     
    225225
    226226# sts gauge with desired index number - used for wave 'Tb'
    227 boundary_csv = join(event_folder, 'sts_gauge_' + str(index) +'.csv')
    228 
     227boundary_csv = join(event_folder, 'sts_gauge_', str(index), '.csv')
     228
  • anuga_work/production/australia_ph2/ceduna/run_model.py

    r6511 r6579  
    156156                        verbose=True)
    157157    domain.set_boundary({'back': Br,
    158                      'side': Bd,
     158                     'side': Bt,
    159159                     'ocean': Bf})
    160160   
     
    163163
    164164    domain.set_boundary({'back': Br,
    165                          'side': Bd,
     165                         'side': Bt,
    166166                         'ocean': Tb})
    167167else:
  • anuga_work/production/australia_ph2/darwin/project.py

    r6549 r6579  
    8484# Format for points easting,northing (no header)
    8585bounding_polygon_filename = 'bounding_polygon_darwin2.csv'
    86 bounding_polygon_maxarea = 100000
     86bounding_polygon_maxarea = 125000 # 100000
    8787
    8888# INTERIOR REGIONS -  for designing the mesh
  • anuga_work/production/australia_ph2/darwin/run_model.py

    r6532 r6579  
    157157                        verbose=True)
    158158    domain.set_boundary({'back': Br,
    159                      'side': Bd,
     159                     'side': Bt,
    160160                     'ocean': Bf})
    161161   
     
    164164
    165165    domain.set_boundary({'back': Br,
    166                          'side': Bd,
     166                         'side': Bt,
    167167                         'ocean': Tb})
    168168else:
  • anuga_work/production/australia_ph2/esperance/run_model.py

    r6511 r6579  
    156156                        verbose=True)
    157157    domain.set_boundary({'back': Br,
    158                      'side': Bd,
     158                     'side': Bt,
    159159                     'ocean': Bf})
    160160   
     
    163163
    164164    domain.set_boundary({'back': Br,
    165                          'side': Bd,
     165                         'side': Bt,
    166166                         'ocean': Tb})
    167167else:
  • anuga_work/production/australia_ph2/eucla_motel/project.py

    r6476 r6579  
    3939# index is only used when wave = Tb
    4040index = 1837            # index from the PTHA
    41 wave = 'Tb'             # Bf (sts wave) Tb (index wave)
     41wave = 'Bf'             # Bf (sts wave) Tb (index wave)
    4242
    4343setup = 'final'         # This can be one of three values
  • anuga_work/production/australia_ph2/eucla_motel/run_model.py

    r6511 r6579  
    156156                        verbose=True)
    157157    domain.set_boundary({'back': Br,
    158                      'side': Bd,
     158                     'side': Bt,
    159159                     'ocean': Bf})
    160160   
     
    163163
    164164    domain.set_boundary({'back': Br,
    165                          'side': Bd,
     165                         'side': Bt,
    166166                         'ocean': Tb})
    167167else:
  • anuga_work/production/australia_ph2/perth/run_model.py

    r6511 r6579  
    156156                        verbose=True)
    157157    domain.set_boundary({'back': Br,
    158                      'side': Bd,
     158                     'side': Bt,
    159159                     'ocean': Bf})
    160160   
     
    163163
    164164    domain.set_boundary({'back': Br,
    165                          'side': Bd,
     165                         'side': Bt,
    166166                         'ocean': Tb})
    167167else:
  • anuga_work/production/australia_ph2/wyndham/run_model.py

    r6539 r6579  
    156156                        verbose=True)
    157157    domain.set_boundary({'back': Br,
    158                      'side': Bd,
     158                     'side': Bt,
    159159                     'ocean': Bf})
    160160   
     
    163163
    164164    domain.set_boundary({'back': Br,
    165                          'side': Bd,
     165                         'side': Bt,
    166166                         'ocean': Tb})
    167167else:
Note: See TracChangeset for help on using the changeset viewer.