Changeset 6476 for anuga_work/production/australia_ph2/eucla_motel

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/eucla_motel

Files:

• project.py (modified) (2 diffs)
• run_model.py (modified) (1 diff)

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,