Changeset 1796

Timestamp:

Sep 7, 2005, 1:20:12 PM (20 years ago)

Author:

ole

Message:

Boundary from MOST

Location:

production/karratha_2005

Files:

• create_mesh.py (modified) (3 diffs)
• prepare_boundary.py (modified) (1 diff)
• project.py (modified) (2 diffs)
• run_karratha.py (modified) (3 diffs)

production/karratha_2005/create_mesh.py

@@ -8 +8 @@
 from pmesh.mesh import *
 from pyvolution.coordinate_transforms.geo_reference import Geo_reference
-from pyvolution.coordinate_transforms.redfearn import redfearn, degminsec2decimal_degrees
+from pyvolution.coordinate_transforms.redfearn import redfearn
+
+
+import project
 
 #-------------------------------------------------------------
@@ -14 +17 @@
 
 
-    refzone = 50
+    
     #Basic geometry
 
-    south = degminsec2decimal_degrees(-20,45,0)
-    north = degminsec2decimal_degrees(-20,15,0)
-    west = degminsec2decimal_degrees(116,30,0)
-    east = degminsec2decimal_degrees(117,0,0)        
+    south = project.south 
+    north = project.north
+    west = project.west
+    east = project.east
+
+    refzone = project.refzone
 
 
@@ -76 +81 @@
     m.generateMesh('pzq28.0za1000000a')
 
-    import project
     m.export_mesh_file(project.meshname + '.msh')
     

production/karratha_2005/prepare_boundary.py

@@ -21 +21 @@
 from pyvolution.data_manager import ferret2sww
 
-#Origin of existing tsh file (Temporary measure)
-mesh_origin = (56,308500.000000,6189000.000000) 
+
+south = project.south 
+north = project.north
+west = project.west
+east = project.east
+
+
+#Origin of existing dem (Temporary measure)
+mesh_origin = (50, 421544.35127423, 7677669.5257159)
+
+    
 
 ferret2sww(source_dir+project.boundary_basename,
            project.boundary_basename, 
            verbose=True,
-           minlat=-35, maxlat=-34,
-           minlon=150.5, maxlon=152,
+           minlat=south-1, maxlat=north+1,
+           minlon=west-1, maxlon=east+1,
            origin = mesh_origin,
            zscale = 1,
-           fail_on_NaN = False),
+           fail_on_NaN = False,
+           inverted_bathymetry = True),
 
 

production/karratha_2005/project.py

@@ -10 +10 @@
 
 meshdir = home+sep+'karratha_tsunami_scenario_2005'+sep+'meshes'+sep
-datadir = home+sep+'karratha_tsunami_scenario_2005'+sep+'data'+sep
+datadir = home+sep+'karratha_tsunami_scenario_2005'+sep+'topographies'+sep
 outputdir = home+sep+'karratha_tsunami_scenario_2005'+sep+'output'+sep
 boundarydir = home+sep+'karratha_tsunami_scenario_2005'+sep+'boundaries'+sep
@@ -18 +18 @@
 meshname = meshdir + basename
 demname = datadir + basename
+boundaryname = boundarydir + 'SU-AU'
 
 boundary_basename = 'SU-AU'
 
+#Georeferencing
+from pyvolution.coordinate_transforms.redfearn import degminsec2decimal_degrees
 
+south = degminsec2decimal_degrees(-20,45,0)
+north = degminsec2decimal_degrees(-20,15,0)
+west = degminsec2decimal_degrees(116,30,0)
+east = degminsec2decimal_degrees(117,0,0)
+
+
+refzone = 50

production/karratha_2005/run_karratha.py

@@ -8 +8 @@
 
 from pyvolution.shallow_water import Domain, Reflective_boundary
-from pyvolution.data_manager import convert_dem_from_ascii2netcdf, dem2pts
+from pyvolution.data_manager import convert_dem_from_ascii2netcdf,\
+     dem2pts, ferret2sww
 from pyvolution.pmesh2domain import pmesh_to_domain_instance 
 from caching import cache
 import project
 
+#Data preparation
 #Convert ASC 2 DEM 2 PTS using source data and store result in source data
-
 demname = project.demname
-
 
 cache(convert_dem_from_ascii2netcdf, demname, {'verbose': True},
       dependencies = [demname + '.asc'],
       verbose = True)
+      #evaluate = True)
 
 cache(dem2pts, demname, {'verbose': True},
@@ -26 +27 @@
       verbose = True)
 
+#Convert MOST boundary
 
+ferret2sww('test', verbose=False,
+           origin = (56, 0, 0))
+
+
+#Read mesh
 mesh = project.meshname + '.msh'
 domain = cache(pmesh_to_domain_instance, (mesh, Domain),
@@ -39 +46 @@
 print "Number of triangles = ", len(domain)
 print 'The extent is ', domain.get_extent()
+
 
 #IC