Changeset 6064

Timestamp:

Dec 11, 2008, 4:10:54 PM (16 years ago)

Author:

kristy

Message:

general maintenance

Location:

anuga_work/production/carnarvon

Files:

• build_carnarvon.py (modified) (1 diff)
• export_results_max.py (modified) (5 diffs)
• project.py (modified) (5 diffs)

anuga_work/production/carnarvon/build_carnarvon.py

@@ -103 +103 @@
 G_off = G_off1 + G_off2
 
-print 'clipping data to coast'
-G_ocean = G_off.clip(project.poly_ocean, verbose=True) #mainland
-G_ocean_b = G_ocean.clip_outside(project.poly_island1, verbose=True) #bernier
-G_ocean_bd = G_ocean_b.clip_outside(project.poly_island2, verbose=True) #dorne
-
-G1_clip_m = G1.clip(project.poly_mainland, verbose=True)
-G1_clip_b = G1.clip(project.poly_island1, verbose=True)
-G1_clip_d = G1.clip(project.poly_island2, verbose=True)
+##print 'clipping data to coast'
+##G_ocean = G_off.clip(project.poly_ocean, verbose=True) #mainland
+##G_ocean_b = G_ocean.clip_outside(project.poly_island1, verbose=True) #bernier
+##G_ocean_bd = G_ocean_b.clip_outside(project.poly_island2, verbose=True) #dorne
+##
+##G1_clip_m = G1.clip(project.poly_mainland, verbose=True)
+##G1_clip_b = G1.clip(project.poly_island1, verbose=True)
+##G1_clip_d = G1.clip(project.poly_island2, verbose=True)
+##G1_clip = G1_clip_m + G1_clip_b + G1_clip_d
 
 print'add all geospatial objects'
-G_aos = G1_clip_m+ G1_clip_b+G1_clip_d+ G_coast + G_ocean_bd
+G_aos = G1 + G_off
 G_aosclip = G_aos.clip_outside(project.poly_aos_extract, verbose=True)
 G_aoi = G_aosclip.clip_outside(project.poly_aoi1, verbose=True)
-G = G_aoi +Goffaoi+Goffaos
+G = G_aoi +G_offaoi+G_offaos
 
 print'export combined DEM file'

anuga_work/production/carnarvon/export_results_max.py

@@ -13 +13 @@
 import project, os
 import sys
+from anuga.lib.maxasc.maxasc import MaxAsc
 from anuga.shallow_water.data_manager import convert_dem_from_ascii2netcdf, dem2pts
 from anuga.shallow_water.data_manager import sww2dem
 from os import sep
-import maxasc
 
 directory = project.output_dir
 
-time_dir1 = '20081110_105300_run_final_0_27283_alpha0.1_kvanputt'
-time_dir2 = '20081002_111207_run_final_0_27283_250m_all_kvanputt'
-time_dir3 = '20081002_111053_run_final_0_27283_250m_none_kvanputt'
+##time_dir1 = '20081209_155247_run_final_0_27255_250m_none_lfountai'
+##time_dir2 = '20081209_161626_run_final_0_68693_250m_none_lfountai'
+##time_dirs = [time_dir1, time_dir2]
 
-time_dirs = [time_dir1, time_dir2, time_dir3]
+##time_dir1 = '20081202_084220_run_final_1_27283_alpha0.1_kvanputt'
+##time_dir2 = '20081202_084202_run_final_1_27255_alpha0.1_kvanputt'
+##time_dir3 = '20081202_084132_run_final_1_68693_alpha0.1_kvanputt'
+time_dir4 = '20081202_084025_run_final_0_68693_alpha0.1_kvanputt'
+time_dir5 = '20081202_083932_run_final_0_27255_alpha0.1_kvanputt'
+##time_dir6 = '20081201_103449_run_final_0_27283_alpha0.1_kvanputt'
+
+time_dirs = [time_dir4, time_dir5] #, time_dir3, time_dir4, time_dir5, time_dir6]
+ 
 
 #cellsize = 20
@@ -36 +44 @@
 ######
 
-#area = ['Geordie', 'Sorrento', 'Fremantle', 'Rockingham']  # strings must match keys in var_equations below
+#area = ['CBD']  # strings must match keys in var_equations below
 area = ['All']      # 'All' means no special areas - the whole thing
 
@@ -50 +58 @@
 
 # one or more key strings from var_equations above
-var = ['stage', 'elevation']
+var = ['stage']
 
 ######
 # Start running the various conversions we require.
 ######
+
 for which_var in var:
     if which_var not in var_equations:
@@ -60 +69 @@
         break
 
-    asc_name = []
+    for which_area in area:
+        if which_area == 'All':
+            easting_min = None
+            easting_max = None
+            northing_min = None
+            northing_max = None
+        else:
+            try:
+                easting_min = eval('project.xmin%s' % which_area)
+                easting_max = eval('project.xmax%s' % which_area)
+                northing_min = eval('project.ymin%s' % which_area)
+                northing_max = eval('project.ymax%s' % which_area)
+            except AttributeError:
+                print 'Unrecognized area name: %s' % which_area
+                break
 
-    for time_dir in time_dirs:
-        if time_dir == '20081002_111053_run_final_0_27283_250m_none_kvanputt':
+        for time_dir in time_dirs:
+
             name1 = directory+time_dir+sep+project.scenario_name
-            name2 = directory+time_dir+sep+project.scenario_name+'_time_77160_0'
-            names = [name1, name2]
-        else:
-            name1 = directory+time_dir+sep+project.scenario_name
-            name2 = directory+time_dir+sep+project.scenario_name+'_time_39540_0'
-            name3 = directory+time_dir+sep+project.scenario_name+'_time_79080_0'
-            names = [name1, name2, name3]     
-       
-        for name in names:
-            for which_area in area:
-                if which_area == 'All':
-                    easting_min = None
-                    easting_max = None
-                    northing_min = None
-                    northing_max = None
-                else:
-                    try:
-                        easting_min = eval('project.xmin%s' % which_area)
-                        easting_max = eval('project.xmax%s' % which_area)
-                        northing_min = eval('project.ymin%s' % which_area)
-                        northing_max = eval('project.ymax%s' % which_area)
-                    except AttributeError:
-                        print 'Unrecognized area name: %s' % which_area
-                        break
-                  
+##            name2 = directory+time_dir+sep+project.scenario_name+'_time_77160_0'
+
+            name2 = directory+time_dir+sep+project.scenario_name+'_time_39300_0'
+            name3 = directory+time_dir+sep+project.scenario_name+'_time_78600_0'
+
+            names = [name1, name2, name3]
+      
+            asc_name = []   
+               
+            for name in names:
+                         
                 outname = name + '_' + which_area + '_' + which_var
                 quantityname = var_equations[which_var]
@@ -109 +118 @@
                 asc_name.append(outname + '.asc')
     
-        maxasc_outname = directory+time_dir+sep+project.scenario_name+'_'+which_area+'_'+which_var+'_max.asc'
+            maxasc_outname = directory+time_dir+sep+project.scenario_name+'_'+which_area+'_'+which_var+'_max.asc'
 
-        maxasc.MaxAsc(maxasc_outname, asc_name)
-
-
+            print 'max asc outname ', maxasc_outname
+            print 'asc_name ', str(asc_name)
+            
+            MaxAsc(maxasc_outname, asc_name)

anuga_work/production/carnarvon/project.py

@@ -43 +43 @@
 # Model specific parameters. One or all can be changed each time the
 # run_scenario script is executed
-tide = 0                #0.6
+tide = 1                #1.0
 #event_number = 27255 # Java 9.3 worst case for Perth
 #event_number = 68693 # Sumatra 9.2
@@ -96 +96 @@
 gauge_name = 'carnarvon.csv'
 gauge_name2 = 'thinned_MGA50.csv'
+#buildings - used in run_building_inundation.py
+building = 'carnarvon_res'
 
 # BOUNDING POLYGON - used in build_boundary.py and run_carnarvon.py respectively
@@ -110 +112 @@
 # Output filename for elevation
 # this is a combination of all the data (utilisied in build_boundary)
-combined_name ='carnarvon_combined_elevation'
+combined_name ='carnarvon_combined_elevation_quick'
 combined_smaller_name = 'carnarvon_combined_elevation_smaller'
 
@@ -168 +170 @@
 gauges_dir_name = gauges_dir + gauge_name       #used for get_timeseries.py
 gauges_dir_name2 = gauges_dir + gauge_name2     #used for get_timeseries.py
-
+building_in_dir_name = gauges_dir + building + '.csv' 
 #------------------------------------------------------------------------------
 # Interior region definitions
@@ -215 +217 @@
 
 # carnarvon CBD extract ascii grid
-##xminCBD = 
-##xmaxCBD = 
-##yminCBD = 
-##ymaxCBD = 
+xminCBD = 763000
+xmaxCBD = 775000
+yminCBD = 7239000
+ymaxCBD = 7251000