Changeset 6476

Timestamp:

Mar 10, 2009, 9:29:00 AM (16 years ago)

Author:

myall

Message:

running final models for various events

Location:

anuga_work/production/australia_ph2

Files:

• ceduna/project.py (modified) (3 diffs)
• ceduna/run_model.py (modified) (1 diff)
• darwin/project.py (modified) (3 diffs)
• darwin/run_model.py (modified) (1 diff)
• esperance/project.py (modified) (2 diffs)
• esperance/run_model.py (modified) (1 diff)
• eucla_motel/project.py (modified) (2 diffs)
• eucla_motel/run_model.py (modified) (1 diff)
• hobart/project.py (modified) (2 diffs)

anuga_work/production/australia_ph2/ceduna/project.py

@@ -28 +28 @@
 tide = 0                # difference between MSL and HAT
 zone = 53               # specify zone of model
-event_number = 64448    # the event number or the mux file name
+event_number = 27347    # the event number or the mux file name
 # ceduna: 27347 (Java), 64448 (sandwich), 58331 (Puysegur)
 # other possible sources:
@@ -37 +37 @@
 friction=0.01           # manning's friction coefficient
 starttime=0             # start time for simulation
-finaltime=1000          # final time for simulation
+finaltime=80000          # final time for simulation
 
 # index is only used when wave = Tb
@@ -224 +224 @@
 mux_input = join(event_folder, mux_input_filename)
 
+# sts gauge with desired index number - used for wave 'Tb'
 boundary_csv = join(event_folder, 'sts_gauge_' + str(index) +'.csv')
 

anuga_work/production/australia_ph2/ceduna/run_model.py

@@ -160 +160 @@
     
 elif project.wave == 'Tb':
-    Tb = Time_boundary(domain,f,default_boundary=Bd )
+    Tb = Time_boundary(domain,f,default_boundary=Bd,
+                       verbose=True)
 
     domain.set_boundary({'back': Br,

anuga_work/production/australia_ph2/darwin/project.py

@@ -28 +28 @@
 tide = 0                # difference between MSL and HAT
 zone = 52               # specify zone of model
-event_number = 17890    # the event number or the mux file name
+event_number = 71261    # the event number or the mux file name
+# contributing events: 17898 (Flores), 64977 (Seram), 70906 (Sumba_normal),
+# 71261 (Timor_OTF)
 alpha = 0.1             # smoothing parameter for mesh
 friction=0.01           # manning's friction coefficient
@@ -35 +37 @@
 
 # index is only used when wave = Tb
-index = 1203            # index from the PTHA
+index = 1833            # index from the PTHA - 1 in 2000 wave height = 0.64m
 wave = 'Tb'             # Bf (sts wave) Tb (index wave)
 
@@ -220 +222 @@
 mux_input = join(event_folder, mux_input_filename)
 
+# sts gauge with desired index number - used for wave 'Tb'
+boundary_csv = join(event_folder, 'sts_gauge_' + str(index) +'.csv')

anuga_work/production/australia_ph2/darwin/run_model.py

@@ -160 +160 @@
     
 elif project.wave == 'Tb':
-    Tb = Time_boundary(domain,f,default_boundary=Bd )
+    Tb = Time_boundary(domain,f,default_boundary=Bd,
+                       verbose=True)
 
     domain.set_boundary({'back': Br,

anuga_work/production/australia_ph2/esperance/project.py

@@ -28 +28 @@
 tide = 0                # difference between MSL and HAT
 zone = 51               # specify zone of model
-event_number = 27255    # the event number or the mux file name
+event_number = 27329    # the event number or the mux file name
+# 7858 (altiplano), 27329 (java), 31951 (kermandec), 46483 (midamerica),
+# 53960 (Puna), 57616 (peru), 58367 (puysegur), 64449 (sandwich), 65356 (shetland)
 alpha = 0.1             # smoothing parameter for mesh
 friction=0.01           # manning's friction coefficient
 starttime=0             # start time for simulation
-finaltime=1000          # final time for simulation
+finaltime=80000          # final time for simulation
 
 # index is only used when wave = Tb
-index = 1203            # index from the PTHA
+index = 1699            # index from the PTHA
 wave = 'Tb'             # Bf (sts wave) Tb (index wave)
 
@@ -220 +222 @@
 mux_input = join(event_folder, mux_input_filename)
 
+# sts gauge with desired index number - used for wave 'Tb'
+boundary_csv = join(event_folder, 'sts_gauge_' + str(index) +'.csv')

anuga_work/production/australia_ph2/esperance/run_model.py

@@ -160 +160 @@
     
 elif project.wave == 'Tb':
-    Tb = Time_boundary(domain,f,default_boundary=Bd )
-
+    Tb = Time_boundary(domain,f,default_boundary=Bd,
+                       verbose=True)
+    
     domain.set_boundary({'back': Br,
                          'side': Bd,

anuga_work/production/australia_ph2/eucla_motel/project.py

@@ -28 +28 @@
 tide = 0                # difference between MSL and HAT
 zone = 52               # specify zone of model
-event_number = 65371    # the event number or the mux file name
+event_number = 64449    # the event number or the mux file name
+# sources: 64449 (sandwich), 68755 (sumatra), 58332 (Puysegur)
+# others: 7859 (altiplano), 11016 (andaman), 17017 (columbia), 27252 (java), 46699 (midamerica)
+# 53961 (puna), 57582 (peru), 58332 (Puysegur), 64449 (sandwich), 68755 (sumatra)
 alpha = 0.1             # smoothing parameter for mesh
 friction=0.01           # manning's friction coefficient
 starttime=0             # start time for simulation
-finaltime=1000          # final time for simulation
+finaltime=80000          # final time for simulation
 
 # index is only used when wave = Tb
-index = 1203            # index from the PTHA
+index = 1837            # index from the PTHA
 wave = 'Tb'             # Bf (sts wave) Tb (index wave)
 
@@ -219 +222 @@
 mux_input = join(event_folder, mux_input_filename)
 
+# sts gauge with desired index number - used for wave 'Tb'
+boundary_csv = join(event_folder, 'sts_gauge_' + str(index) +'.csv')

anuga_work/production/australia_ph2/eucla_motel/run_model.py

@@ -160 +160 @@
     
 elif project.wave == 'Tb':
-    Tb = Time_boundary(domain,f,default_boundary=Bd )
+    Tb = Time_boundary(domain,f,default_boundary=Bd,
+                       verbose=True)
 
     domain.set_boundary({'back': Br,

anuga_work/production/australia_ph2/hobart/project.py

@@ -28 +28 @@
 tide = 0                # difference between MSL and HAT
 zone = 55               # specify zone of model
-event_number = 58364    # the event number or the mux file name
+
 alpha = 0.1             # smoothing parameter for mesh
 friction=0.01           # manning's friction coefficient
 starttime=0             # start time for simulation
-finaltime=1000          # final time for simulation
+finaltime=80000          # final time for simulation
 
 # index is only used when wave = Tb
-index = 1203            # index from the PTHA
 wave = 'Tb'             # Bf (sts wave) Tb (index wave)
+# for melbourne (west) side:
+event_number = 58349    # the event number or the mux file name
+# contributing sources:
+# 51328 (newhebrides), 58349 (puysegur), 64535 (sandwich)
+# other possible: 31966 (kermandec), 65350 (shetland)
+index = 1983            # index from the PTHA - west side, close to melbourne; 1983 has 1 in 2000 wave height 0.534m
+
+# for east, Hobart, side - point opposite hobart
+##event_number = 51202    # the event number or the mux file name
+# contributing events: 64434 (sandwich), 58117 (puysegur), 51202 (new hebrides)
+# Others: 31899 (kermandec), 46657 (midamerica), 57586 (peru),
+##index = 2473            # index from PTHA - east, opposite Hobart; 1.112m for 1 in 2000
 
 setup = 'final'         # This can be one of three values
@@ -220 +231 @@
 mux_input = join(event_folder, mux_input_filename)
 
+# sts gauge with desired index number - used for wave 'Tb'
+boundary_csv = join(event_folder, 'sts_gauge_' + str(index) +'.csv')