Changeset 5704

Timestamp:

Aug 28, 2008, 9:42:21 AM (16 years ago)

Author:

kristy

Message:

 

File:

• anuga_work/development/boxingday08/export_results_GA.py (modified) (1 diff)

anuga_work/development/boxingday08/export_results_GA.py

@@ -6 +6 @@
 
 home = getenv('INUNDATIONHOME') + sep +'data'+sep #Sandpit's parent dir
+output_dir = home + 'anuga_validation/boxing_day_sceanario/anuga/outputs/20080813/good_polyline20080813_' 
 
-#name='good_simulation/good_simulation_polylinefriction20080529_115703'
-#name='good_simulation/good_simulation_ursfriction20080529_160110'
+name1= output_dir +'101348_0.45'
+name2= output_dir +'101521_0.55'
+name3= output_dir +'101610_0.35'
+name4= output_dir +'101741_0.25'
+name5= output_dir +'101819_0.15'
 
-name= home + 'anuga_validation/boxing_day_sceanario/anuga/outputs/20080723/good_polyline20080723_45_141453'
+names= [name1, name2, name3, name4, name5]
 
-print 'output dir:', name
+for the_name in names:
 
-bay_west = 417348.0
-bay_east = 425656.0
-bay_south = 870565.0
-bay_north = 877008.0
+    name = the_name[-11:]
+    print 'sww filename', name
 
-which_var = 2
+    bay_west = 417348.0
+    bay_east = 425656.0
+    bay_south = 870565.0
+    bay_north = 877008.0
 
-if which_var == 0:  # Stage
-    outname = name + '_stage-eq-thurs-test'
-    quantityname = 'stage'
+    which_var = 2
 
-if which_var == 1:  # Absolute Momentum
-    outname = name + '_momentum_i1'
-    quantityname = '(xmomentum**2 + ymomentum**2)**0.5'  
+    if which_var == 0:  # Stage
+        outname = name + '_stage-eq-thurs-test'
+        quantityname = 'stage'
 
-if which_var == 2:  # Depth
-    outname = name + '_depth'
-    quantityname = 'stage-elevation'  
+    if which_var == 1:  # Absolute Momentum
+        outname = name + '_momentum_i1'
+        quantityname = '(xmomentum**2 + ymomentum**2)**0.5'  
 
-if which_var == 3:  # Speed
-    outname = name + '_speed_i0'
-    quantityname = '(xmomentum**2 + ymomentum**2)**0.5/(stage-elevation+1.e-6/(stage-elevation))'  #Speed
+    if which_var == 2:  # Depth
+        outname = name + '_depth'
+        quantityname = 'stage-elevation'  
 
-if which_var == 4:  # Elevation
-    outname = name + '_elevation'
-    quantityname = 'elevation'  #Elevation
+    if which_var == 3:  # Speed
+        outname = name + '_speed_i0'
+        quantityname = '(xmomentum**2 + ymomentum**2)**0.5/(stage-elevation+1.e-6/(stage-elevation))'  #Speed
 
-print 'start sww2dem'
-sww2dem(name, basename_out = outname,
-        quantity = quantityname,
-        cellsize = 20,
-        easting_min = bay_west,
-        easting_max = bay_east,
-        northing_min = bay_south,
-        northing_max = bay_north,   
-        reduction = max, 
-        verbose = True,
-        format = 'asc')
+    if which_var == 4:  # Elevation
+        outname = name + '_elevation'
+        quantityname = 'elevation'  #Elevation
+
+    print 'start sww2dem'
+    sww2dem(the_name, basename_out = outname,
+            quantity = quantityname,
+            cellsize = 20,
+            easting_min = bay_west,
+            easting_max = bay_east,
+            northing_min = bay_south,
+            northing_max = bay_north,   
+            reduction = max, 
+            verbose = True,
+            format = 'asc')