Changeset 1784

Timestamp:

Aug 30, 2005, 4:44:56 PM (20 years ago)

Author:

ole

Message:

Basic scripts for Karratha scenario

Location:

production/karratha_2005

Files:

• create_mesh.py (added)
• project.py (modified) (1 diff)
• run_karratha.py (modified) (1 diff)

production/karratha_2005/project.py

@@ -11 +11 @@
 meshdir = home+sep+'karratha_tsunami_scenario_2005'+sep+'meshes'+sep
 datadir = home+sep+'karratha_tsunami_scenario_2005'+sep+'data'+sep
-basename = 'karratha_250m'
-
-meshfile = meshdir + basename + '.msh'
-demfile = datadir + basename + '.dem'
+outputdir = home+sep+'karratha_tsunami_scenario_2005'+sep+'output'+sep
 
 
-print demfile
+basename = 'karratha_250m'
+meshname = meshdir + basename
+demname = datadir + basename
 
-#boundary_filename = 
-#bathymetry_filename = 
-#mesh_filename = 

production/karratha_2005/run_karratha.py

@@ -6 +6 @@
 import time
 
+
+from pyvolution.shallow_water import Domain, Reflective_boundary
 from pyvolution.data_manager import convert_dem_from_ascii2netcdf, dem2pts
+from pyvolution.pmesh2domain import pmesh_to_domain_instance 
+from caching import cache
 import project
 
 #Convert ASC 2 DEM 2 PTS using source data and store result in source data
-convert_dem_from_ascii2netcdf(project.demfile[:-4], verbose=True)
-#dem2pts(source_dir2+rootfilename2, verbose = True)
 
+demname = project.demname
+
+
+cache(convert_dem_from_ascii2netcdf, demname, {'verbose': True},
+      dependencies = [demname + '.asc'],
+      verbose = True)
+
+cache(dem2pts, demname, {'verbose': True},
+      dependencies = [demname + '.dem'],      
+      verbose = True)
+
+
+mesh = project.meshname + '.msh'
+domain = cache(pmesh_to_domain_instance, (mesh, Domain),
+               dependencies = [mesh],                     
+               verbose = True)               
+
+
+domain.set_name(project.basename)
+
+
+#IC
+domain.set_quantity('stage', 0)
+
+
+#BC
+Br = Reflective_boundary(domain)
+domain.set_boundary( {'wall': Br} )
+
+
+#Run
+for t in domain.evolve(yieldstep = 1, finaltime = 1):
+    domain.write_time()