Changeset 4525 for anuga_work/production/test_dampier_2006/test.py

Timestamp:

Jun 1, 2007, 4:06:46 PM (17 years ago)

Author:

nick

Message:

update run

File:

• anuga_work/production/test_dampier_2006/test.py (modified) (3 diffs)

anuga_work/production/test_dampier_2006/test.py

@@ -3 +3 @@
 import os
 
-import project
+#import project
+#from Numeric import allclose
 from anuga.pmesh.mesh_interface import create_mesh_from_regions
 from anuga.shallow_water import Domain
 from anuga.shallow_water import File_boundary
+from time import localtime, strftime, gmtime
 
 home = getenv('INUNDATIONHOME') #Sandpit's parent dir   
@@ -13 +15 @@
 state = 'western_australia'
 scenario_name = 'dampier'
-scenario = 'dampier_tsunami_scenario_2006'
+scenario = 'test_dampier'
 an = 'anuga'
 bo = 'boundaries'
+topo = 'topographies'
+out = 'outputs'
 urs = 'urs'
 gridded = '1_10000'
-best_boundary_sww_file = 'o_test_8500_12000.sww'
+#best_boundary_sww_file = 'o_test_8500_12000.sww'
+best_boundary_sww_file = 'o_8500_12000_no_zone.sww'
+
+time = strftime('%Y%m%d_%H%M%S',localtime()) #gets time for new dir
+gtime = strftime('%Y%m%d_%H%M%S',gmtime()) #gets time for new dir
+run_time = time+'_run'
+print 'gtime: ', gtime
+print 'time: ', time
+
 
 mesh_filename = 'small.tsh' # this will be local
-best_boundary_sww = os.path.join(data,state,scenario_name,
-                                 scenario,an,bo,best_boundary_sww_file)
+best_boundary_sww = os.path.join(home,data,state,scenario,an,bo,best_boundary_sww_file)
+#J:\inundation\data\western_australia\dampier_tsunami_scenario_2006\anuga\topographies
+                                 
+elevation_file = 'dampier_combined_elevation_smallest.txt'
+smallest_elevation= os.path.join(home,data,state,scenario,an,topo,elevation_file)
+
+output= os.path.join(home,data,state,scenario,an,out,run_time)
 
 #pro_instance = project.Project([data,state,scenario_name],
 #                                       outputdir_name=outputdir_name)
 
-poly_all = [[250000,7660000],[280000,7660000],
-                             [280000,7630000]]
+#poly_all = [[250000,7660000],[280000,7660000], [280000,7630000]]
+poly_all = [[469000,7760000],[470000,7758000],[468000,7758000]]
+
+print 'create mesh'
 create_mesh_from_regions(poly_all,
                              boundary_tags={'back': [0,1,2]},
@@ -38 +57 @@
 print domain.statistics()
 
-domain.set_quantity('stage', 0.0)
+domain.set_quantity('stage', 2.4)
 domain.set_quantity('friction', 0.01)
+domain.set_name(scenario_name)
+domain.set_datadir(output)
+
+print 'set elevation',smallest_elevation
+
 domain.set_quantity('elevation', 
-#                    filename = project.combined_dir_name + '.pts',
-# MUST USE TXT FILES FOR CACHING TO WORK!
-                    filename = project.combined_smaller_dir_name + '.txt',
+                    filename = smallest_elevation,
                     use_cache = True,
                     verbose = True,
                     alpha = 0.1)
-    
+                    
+
+print 'best_boundary_sww',best_boundary_sww
+
+print 'get file boundary'
 Bf = File_boundary(best_boundary_sww,
-                  domain,  use_cache=False, verbose=True)
+                  domain,  time_thinning=24, use_cache=False, verbose=True)
+                  
+domain.set_boundary({'back': Bf})
+                     
+                     
+#t0 = time.time()
+
+print 'evolve'
+for t in domain.evolve(yieldstep = 120, finaltime = 3400):
+    domain.write_time()
+    domain.write_boundary_statistics(tags = 'back')
+
+
+
+