Changeset 2270

Timestamp:

Jan 24, 2006, 2:08:25 PM (19 years ago)

Author:

jakeman

Message:

Cairns_2006 now contains create_mesh file hich creates mesh with multiple refined regions. Still not working correctly though.
Also contains sww files from two basic runs.

Location:

development/cairns_2006

Files:

• cairns_fault.xya (modified) (1 diff)
• run_cairns.py (modified) (7 diffs)

development/cairns_2006/cairns_fault.xya

@@ -1 @@
-5481413.1609155666, 1558472.3824136395
-4933271.8448240105, 1558472.3824136395
-5444342.9345278544, 1335833.4706402624
-4899908.64107507, 1335833.4706402624
+5481413.1609155666, 1001875.1029801967
+4933271.8448240105, 1001875.1029801967
+5481413.1609155666, 779236.19120681973
+4933271.8448240105, 779236.19120681973

development/cairns_2006/run_cairns.py

@@ -14 +14 @@
 from pyvolution.pmesh2domain import pmesh_to_domain_instance
 from Numeric import array, zeros, Float, allclose
-import cairns_project
+#import cairns_project
+import project
 from caching import cache
 from utilities.polygon import read_polygon, Polygon_function
@@ -22 +23 @@
 #-------------------------------
 print 'Creating domain'
-M = 72 #float(len1)/m
-N = 20 #float(len2)/n
+"""
+M = 36 #float(len1)/m
+N = 10 #float(len2)/n
 
 points, elements, boundary = rectangular_cross(M, N, len1=6393693.9504678631,
@@ -31 +33 @@
 #points, elements, boundary = rectangular_cross(N, M, len1=58.0, len2=1.0, origin=(145.0,-24.0))
                                                
-domain = Domain(points, elements, boundary)
+domain = cache(Domain, (points, elements, boundary))
+
+"""
+domain = cache(pmesh_to_domain_instance,
+               (project.mesh_filename, Domain),
+               dependencies = [project.mesh_filename])
 
 domain.check_integrity()
@@ -49 +56 @@
 #p0 = read_polygon('cairns_fault_degrees.xya')
 
-domain.set_quantity('stage',Polygon_function([(p0,1.0)]))
+domain.set_quantity('stage',Polygon_function([(p0,1000.0)]))
 
 #print domain.quantities['stage'].vertex_values
 
 domain.set_quantity('elevation',
-                    filename = cairns_project.bathymetry_filename[:-4] + '.xya',
+                    filename = project.bathymetry_filename[:-4] + '.xya',
                     alpha = 10.0,
                     verbose = True,
@@ -70 +77 @@
 print 'Boundaries'
 
-#   Tidal cycle recorded at Eden as open
-#print 'Open sea boundary condition from ',cairns_project.boundary_filename
-#from pyvolution.util import file_function
-#tide_function = file_function(cairns_project.boundary_filename[:-4] + '.tms',
-#                               domain,
-#                               verbose = True)
+
 #Bts = Transmissive_Momentum_Set_Stage_boundary(domain, tide_function)
 
-#   All other boundaries are reflective
-#Br = Reflective_boundary(domain)
-#domain.set_boundary({'left': Br, 'right': Br, 'top': Br, 'bottom': Br})
+Br = Reflective_boundary(domain)
+domain.set_boundary({'border': Br, 'border': Br, 'border': Br, 'border': Br})
 
-Bt = Transmissive_boundary(domain)
-domain.set_boundary({'left': Bt, 'right': Bt, 'top': Bt, 'bottom': Bt})
+#Bt = Transmissive_boundary(domain)
+#domain.set_boundary({'left': Bt, 'right': Bt, 'top': Bt, 'bottom': Bt})
 
 #-------------------------------
 # Setup domain runtime parameters
 #-------------------------------
+
 domain.visualise = True
 domain.visualise_color_stage = True
+
+#domain.visualise = False
+#domain.visualise_color_stage = False
 
 base = os.path.basename(sys.argv[0])
@@ -96 +101 @@
 domain.store = True    #Store for visualisation purposes
 
-
-
 #-------------------------------
 # Evolve
@@ -103 +106 @@
 import time
 t0 = time.time()
-yieldstep = 10
-finaltime = 120
-
+yieldstep = 60.0
+finaltime = 480.0
 
 for t in domain.evolve(yieldstep = yieldstep, finaltime = finaltime):