Changeset 1801

Timestamp:

Sep 7, 2005, 4:31:19 PM (20 years ago)

Author:

ole

Message:

 

Location:

production/karratha_2005

Files:

• prepare_boundary.py (modified) (1 diff)
• run_karratha.py (modified) (5 diffs)

production/karratha_2005/prepare_boundary.py

@@ -9 +9 @@
 """
 
+#FIXME: OBSOLETE
 
 import os

production/karratha_2005/run_karratha.py

@@ -7 +7 @@
 
 
-from pyvolution.shallow_water import Domain, Reflective_boundary
+from pyvolution.shallow_water import Domain, Reflective_boundary, File_boundary
 from pyvolution.data_manager import convert_dem_from_ascii2netcdf,\
      dem2pts, ferret2sww
@@ -27 +27 @@
       verbose = True)
 
+
 #Convert MOST boundary
-
-ferret2sww('test', verbose=False,
-           origin = (56, 0, 0))
+source_dir = project.boundarydir
 
 
-#Read mesh
+#Data
+from pyvolution.data_manager import ferret2sww
+
+
+south = project.south 
+north = project.north
+west = project.west
+east = project.east
+
+#Origin of existing dem (FIXME: Temporary measure)
+mesh_origin = (50, 421544.35127423, 7677669.5257159)
+
+ferret2sww(source_dir+project.boundary_basename,
+           project.boundary_basename, 
+           verbose=True,
+           minlat=south-1, maxlat=north+1,
+           minlon=west-1, maxlon=east+1,
+           origin = mesh_origin,
+           zscale = 1,
+           fail_on_NaN = False,
+           inverted_bathymetry = True),
+
+
+
+#Setup domain
 mesh = project.meshname + '.msh'
 domain = cache(pmesh_to_domain_instance, (mesh, Domain),
@@ -48 +71 @@
 
 
-#IC
+#Setup IC
 domain.set_quantity('stage', 0)
 domain.set_quantity('elevation',
@@ -55 +78 @@
                     verbose = True)
 
-#BC
+#Setup BC
+Bf = File_boundary(project.boundary_basename + '.sww', domain, verbose = True)
+
 Br = Reflective_boundary(domain)
 domain.set_boundary( {'wall': Br} )
@@ -61 +86 @@
 
 #Run
-for t in domain.evolve(yieldstep = 1, finaltime = 1):
+for t in domain.evolve(yieldstep = 1, finaltime = 100):
     domain.write_time()