Changeset 6476


Ignore:
Timestamp:
Mar 10, 2009, 9:29:00 AM (14 years ago)
Author:
myall
Message:

running final models for various events

Location:
anuga_work/production/australia_ph2
Files:
9 edited

Legend:

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

    r6468 r6476  
    2828tide = 0                # difference between MSL and HAT
    2929zone = 53               # specify zone of model
    30 event_number = 64448    # the event number or the mux file name
     30event_number = 27347    # the event number or the mux file name
    3131# ceduna: 27347 (Java), 64448 (sandwich), 58331 (Puysegur)
    3232# other possible sources:
     
    3737friction=0.01           # manning's friction coefficient
    3838starttime=0             # start time for simulation
    39 finaltime=1000          # final time for simulation
     39finaltime=80000          # final time for simulation
    4040
    4141# index is only used when wave = Tb
     
    224224mux_input = join(event_folder, mux_input_filename)
    225225
     226# sts gauge with desired index number - used for wave 'Tb'
    226227boundary_csv = join(event_folder, 'sts_gauge_' + str(index) +'.csv')
    227228
  • anuga_work/production/australia_ph2/ceduna/run_model.py

    r6460 r6476  
    160160   
    161161elif project.wave == 'Tb':
    162     Tb = Time_boundary(domain,f,default_boundary=Bd )
     162    Tb = Time_boundary(domain,f,default_boundary=Bd,
     163                       verbose=True)
    163164
    164165    domain.set_boundary({'back': Br,
  • anuga_work/production/australia_ph2/darwin/project.py

    r6461 r6476  
    2828tide = 0                # difference between MSL and HAT
    2929zone = 52               # specify zone of model
    30 event_number = 17890    # the event number or the mux file name
     30event_number = 71261    # the event number or the mux file name
     31# contributing events: 17898 (Flores), 64977 (Seram), 70906 (Sumba_normal),
     32# 71261 (Timor_OTF)
    3133alpha = 0.1             # smoothing parameter for mesh
    3234friction=0.01           # manning's friction coefficient
     
    3537
    3638# index is only used when wave = Tb
    37 index = 1203            # index from the PTHA
     39index = 1833            # index from the PTHA - 1 in 2000 wave height = 0.64m
    3840wave = 'Tb'             # Bf (sts wave) Tb (index wave)
    3941
     
    220222mux_input = join(event_folder, mux_input_filename)
    221223
     224# sts gauge with desired index number - used for wave 'Tb'
     225boundary_csv = join(event_folder, 'sts_gauge_' + str(index) +'.csv')
  • anuga_work/production/australia_ph2/darwin/run_model.py

    r6461 r6476  
    160160   
    161161elif project.wave == 'Tb':
    162     Tb = Time_boundary(domain,f,default_boundary=Bd )
     162    Tb = Time_boundary(domain,f,default_boundary=Bd,
     163                       verbose=True)
    163164
    164165    domain.set_boundary({'back': Br,
  • anuga_work/production/australia_ph2/esperance/project.py

    r6460 r6476  
    2828tide = 0                # difference between MSL and HAT
    2929zone = 51               # specify zone of model
    30 event_number = 27255    # the event number or the mux file name
     30event_number = 27329    # the event number or the mux file name
     31# 7858 (altiplano), 27329 (java), 31951 (kermandec), 46483 (midamerica),
     32# 53960 (Puna), 57616 (peru), 58367 (puysegur), 64449 (sandwich), 65356 (shetland)
    3133alpha = 0.1             # smoothing parameter for mesh
    3234friction=0.01           # manning's friction coefficient
    3335starttime=0             # start time for simulation
    34 finaltime=1000          # final time for simulation
     36finaltime=80000          # final time for simulation
    3537
    3638# index is only used when wave = Tb
    37 index = 1203            # index from the PTHA
     39index = 1699            # index from the PTHA
    3840wave = 'Tb'             # Bf (sts wave) Tb (index wave)
    3941
     
    220222mux_input = join(event_folder, mux_input_filename)
    221223
     224# sts gauge with desired index number - used for wave 'Tb'
     225boundary_csv = join(event_folder, 'sts_gauge_' + str(index) +'.csv')
  • anuga_work/production/australia_ph2/esperance/run_model.py

    r6460 r6476  
    160160   
    161161elif project.wave == 'Tb':
    162     Tb = Time_boundary(domain,f,default_boundary=Bd )
    163 
     162    Tb = Time_boundary(domain,f,default_boundary=Bd,
     163                       verbose=True)
     164   
    164165    domain.set_boundary({'back': Br,
    165166                         'side': Bd,
  • anuga_work/production/australia_ph2/eucla_motel/project.py

    r6460 r6476  
    2828tide = 0                # difference between MSL and HAT
    2929zone = 52               # specify zone of model
    30 event_number = 65371    # the event number or the mux file name
     30event_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)
    3134alpha = 0.1             # smoothing parameter for mesh
    3235friction=0.01           # manning's friction coefficient
    3336starttime=0             # start time for simulation
    34 finaltime=1000          # final time for simulation
     37finaltime=80000          # final time for simulation
    3538
    3639# index is only used when wave = Tb
    37 index = 1203            # index from the PTHA
     40index = 1837            # index from the PTHA
    3841wave = 'Tb'             # Bf (sts wave) Tb (index wave)
    3942
     
    219222mux_input = join(event_folder, mux_input_filename)
    220223
     224# sts gauge with desired index number - used for wave 'Tb'
     225boundary_csv = join(event_folder, 'sts_gauge_' + str(index) +'.csv')
  • anuga_work/production/australia_ph2/eucla_motel/run_model.py

    r6460 r6476  
    160160   
    161161elif project.wave == 'Tb':
    162     Tb = Time_boundary(domain,f,default_boundary=Bd )
     162    Tb = Time_boundary(domain,f,default_boundary=Bd,
     163                       verbose=True)
    163164
    164165    domain.set_boundary({'back': Br,
  • anuga_work/production/australia_ph2/hobart/project.py

    r6460 r6476  
    2828tide = 0                # difference between MSL and HAT
    2929zone = 55               # specify zone of model
    30 event_number = 58364    # the event number or the mux file name
     30
    3131alpha = 0.1             # smoothing parameter for mesh
    3232friction=0.01           # manning's friction coefficient
    3333starttime=0             # start time for simulation
    34 finaltime=1000          # final time for simulation
     34finaltime=80000          # final time for simulation
    3535
    3636# index is only used when wave = Tb
    37 index = 1203            # index from the PTHA
    3837wave = 'Tb'             # Bf (sts wave) Tb (index wave)
     38# for melbourne (west) side:
     39event_number = 58349    # the event number or the mux file name
     40# contributing sources:
     41# 51328 (newhebrides), 58349 (puysegur), 64535 (sandwich)
     42# other possible: 31966 (kermandec), 65350 (shetland)
     43index = 1983            # index from the PTHA - west side, close to melbourne; 1983 has 1 in 2000 wave height 0.534m
     44
     45# for east, Hobart, side - point opposite hobart
     46##event_number = 51202    # the event number or the mux file name
     47# contributing events: 64434 (sandwich), 58117 (puysegur), 51202 (new hebrides)
     48# Others: 31899 (kermandec), 46657 (midamerica), 57586 (peru),
     49##index = 2473            # index from PTHA - east, opposite Hobart; 1.112m for 1 in 2000
    3950
    4051setup = 'final'         # This can be one of three values
     
    220231mux_input = join(event_folder, mux_input_filename)
    221232
     233# sts gauge with desired index number - used for wave 'Tb'
     234boundary_csv = join(event_folder, 'sts_gauge_' + str(index) +'.csv')
Note: See TracChangeset for help on using the changeset viewer.