Changeset 2077

Timestamp:

Nov 25, 2005, 3:39:56 PM (19 years ago)

Author:

duncan

Message:

update

File:

• production/gippsland_2005/run_gippsland.py (modified) (4 diffs)

production/gippsland_2005/run_gippsland.py

@@ -10 +10 @@
 import project
 from create_mesh import create_mesh
-
+from pyvolution.shallow_water import Domain, Reflective_boundary,\
+     File_boundary, Dirichlet_boundary, Time_boundary, Transmissive_boundary
+from pyvolution.least_squares import fit_to_mesh_file, DEFAULT_ALPHA 
 
 #Data preparation
@@ -20 +22 @@
 #      #evaluate = True)
 
-cache(dem2pts, project.demname, {'verbose': True},
-      dependencies = [project.demname + '.dem']      
+cache(dem2pts, project.demname[:-4], {'verbose': True},
+      dependencies = [project.demname]      
       ,verbose = False
       )
@@ -34 +36 @@
        project.ucur_dir,
        project.vcur_dir,
-       project.boundaryname + '.sww'),
+       project.boundaryname),
        {'verbose': True,
        'minlat': -38.25,
@@ -54 +56 @@
 #fit_to_mesh_file(mesh_file, point_file, mesh_output_file, DEFAULT_ALPHA,
 #                 verbose=True, expand_search=True, precrop=True)      
-cache(fit_to_mesh_file,(project.mesh_file,
-                 point_file,
-                 mesh_output_file,
+cache(fit_to_mesh_file,(project.meshname,
+                 project.pointname,
+                 project.mesh_elevname,
                  DEFAULT_ALPHA),
       {'verbose': True,
        'expand_search': True,
        'precrop': True}
-      ,dependencies = [project.mesh_file, point_file]
-      #,verbose = True
+      ,dependencies = [project.meshname, project.pointname]
+      ,verbose = False
       )
     
+
+#Setup domain
+domain = cache(pmesh_to_domain_instance, (project.mesh_elevname, Domain),
+               dependencies = [project.mesh_elevname]                    
+               ,verbose = False
+               )               
+
+
+domain.set_name(project.basename) # + '_%d' %resolution)
+domain.set_datadir(project.outputdir)
+domain.store = True
+domain.quantities_to_be_stored = ['stage']
+
+print 'Number of triangles = ', len(domain)
+print 'The extent is ', domain.get_extent()
+
+#Setup Initial Conditions
+domain.set_quantity('friction', 0)
+domain.set_quantity('stage', project.tide)
+
+
+#Setup Boundary Conditions
+print domain.get_boundary_tags()
+
+print "****** run  ****"
+print "project.boundaryname",project.boundaryname 
+print "**********"
+
+Bf = File_boundary(project.boundaryname, domain, verbose = True)
+#domain.starttime = 3000  #Obtained from MOST
+domain.starttime = 0  #Obtained from MOST
+
+Br = Reflective_boundary(domain)
+Bt = Transmissive_boundary(domain)
+Bd = Dirichlet_boundary([project.tide,0,0])
+Bw = Time_boundary(domain=domain,
+                   f=lambda t: [(60<t<660)*4, 0, 0])
+
+domain.set_boundary( {'back': Br,'side': Bd, 'ocean': Bf} ) #CSIRO storm surge
+
+#Evolve
+import time
+t0 = time.time()
+
+for t in domain.evolve(yieldstep = 60, finaltime = 1200):
+#                       skip_initial_step = True):
+    domain.write_time()
+    domain.write_boundary_statistics(tags = 'ocean') #quantities = 'stage')       
+
+print 'That took %.2f seconds' %(time.time()-t0)