Changeset 5387

Timestamp:

Jun 2, 2008, 8:59:09 AM (17 years ago)

Author:

kristy

Message:

 

Location:

anuga_work/production

Files:

• busselton/build_busselton.py (modified) (2 diffs)
• busselton/project.py (modified) (3 diffs)
• perth/build_perth.py (modified) (1 diff)
• perth/project.py (modified) (1 diff)

anuga_work/production/busselton/build_busselton.py

@@ -90 +90 @@
 G2 = Geospatial_data(file_name = coast_in_dir_name + '.txt',verbose=True)
 G3 = Geospatial_data(file_name = coast_in_dir_name1 + '.txt',verbose=True)
-#G3 = Geospatial_data(file_name = island_dir_name + '.pts',verbose=True)
 G_off = Geospatial_data(file_name = offshore_in_dir_name + '.txt',verbose=True)
 
@@ -108 +107 @@
 
 print'project.combined_dir_name + .txt',project.combined_dir_name + '.txt'
-#G_all=Geospatial_data(file_name = project.combined_dir_name + '.txt')
-print'split'
-G_all_1 = G.split(.10)
-print'export 1'
-G_all_1.export_points_file(project.combined_dir_name+'_small' + '.txt')
-#print'export 2'
-#G_all_2.export_points_file(project.combined_dir_name+'_other1' + '.xya')
 
-'''
 #-------------------------------------------------------------------------
 # Convert URS to SWW file for boundary conditions 
 #-------------------------------------------------------------------------
 print 'starting to create boundary conditions'
+from anuga.shallow_water.data_manager import urs2sww, urs_ungridded2sww
+
 boundaries_in_dir_name = project.boundaries_in_dir_name
+print 'boundaries_in_dir_name',project.boundaries_in_dir_name
 
-from anuga.shallow_water.data_manager import urs2sww
-
-print 'minlat=project.north_boundary, maxlat=project.south_boundary',project.north_boundary, project.south_boundary
-print 'minlon= project.west_boundary, maxlon=project.east_boundary',project.west_boundary, project.east_boundary
 
 #import sys; sys.exit()
 
-#if access(project.boundaries_dir,F_OK) == 0:
-#    mkdir (project.boundaries_dir)
+urs_ungridded2sww(project.boundaries_in_dir_name, project.boundaries_in_dir_name,
+                  verbose=True, mint=4000, maxt=80000, zscale=1)
 
-from caching import cache
-cache(urs2sww,
-      (boundaries_in_dir_name,
-       project.boundaries_dir_name1), 
-      {'verbose': True,
-       'minlat': project.south_boundary,
-       'maxlat': project.north_boundary,
-       'minlon': project.west_boundary,
-       'maxlon': project.east_boundary,
-#       'minlat': project.south,
-#       'maxlat': project.north,
-#       'minlon': project.west,
-#       'maxlon': project.east,
-       'mint': 0, 'maxt': 40000,
-#       'origin': domain.geo_reference.get_origin(),
-       'mean_stage': project.tide,
-#       'zscale': 1,                 #Enhance tsunami
-       'fail_on_NaN': False},
-       verbose = True,
-       )
-#       dependencies = source_dir + project.boundary_basename + '.sww')
-
-'''
-
-
-
-
-
-
-

anuga_work/production/busselton/project.py

@@ -42 +42 @@
 finaltime=10000
 export_cellsize=50
-setup='trial'
-source='other'
+setup='final'
+source='exmouth'
 
 
@@ -108 +108 @@
 tide_dir = anuga_dir+'tide_data'+sep
 
-if source =='dampier':
-    boundaries_name = 'broome_3854_17042007' #Dampier gun
-    boundaries_in_dir = anuga_dir+'boundaries'+sep+'urs'+sep+'dampier'+sep+'1_10000'+sep
-
-if source=='onslow':
-    boundaries_name = 'broome_3859_16052007' #onslow_hedland_broome gun
-    boundaries_in_dir = anuga_dir+'boundaries'+sep+'urs'+sep+'onslow_hedland_broome'+sep+'1_10000'+sep
-    
 if source=='exmouth':
-    boundaries_name = 'broome_3103_18052007' #exmouth gun
+    boundaries_name = 'busselton_3103_30052008' #exmouth gun
     boundaries_in_dir = anuga_dir+'boundaries'+sep+'urs'+sep+'exmouth'+sep+'1_10000'+sep
 
-if source=='other':
-    boundaries_name = 'test' #exmouth gun
+if source=='test':
+    boundaries_name = 'other' #exmouth gun
     boundaries_in_dir = anuga_dir+'boundaries'+sep+'urs'+sep+'exmouth'+sep+'1_10000'+sep
 
@@ -190 +182 @@
 trigs_min = number_mesh_triangles(interior_regions, poly_all, res_poly_all)
 
-poly_mainland=read_polygon(topographies_in_dir+'initial_condition_south.csv')
+poly_mainland=read_polygon(polygons_dir+'initial_condition.csv')
 
 print 'min number triangles', trigs_min

anuga_work/production/perth/build_perth.py

@@ -141 +141 @@
 
 urs_ungridded2sww(project.boundaries_in_dir_name, project.boundaries_in_dir_name,
-                  verbose=True, mint=4000, maxt=10000, zscale=1)
+                  verbose=True, mint=4000, maxt=80000, zscale=1)
 
 

anuga_work/production/perth/project.py

@@ -130 +130 @@
 
 #boundaries locations
-boundaries_in_dir_name = boundaries_in_dir + scenario_name
+boundaries_in_dir_name = boundaries_in_dir + boundaries_name 
 boundaries_dir = anuga_dir+'boundaries'+sep
 boundaries_dir_name = boundaries_dir + scenario_name