Changeset 6391


Ignore:
Timestamp:
Feb 24, 2009, 9:18:23 AM (14 years ago)
Author:
kristy
Message:
 
Location:
anuga_work/production/perth
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/perth/Arc_asc2raster_GDA94z50.py

    r6067 r6391  
    2828##time_dir2 = '20081210_100528_run_final_0_68693_250m_none_kvanputt'
    2929##time_dirs = [time_dir1, time_dir2]
    30 
     30time_dir1 = '20090129_155526_run_final_0.6_27255_Rottnest_kvanputt'
    3131##time_dir1 = '20081031_133353_run_final_0.6_68693_alpha0.1_kvanputt'
    3232##time_dir2 = '20081031_133511_run_final_0_68693_alpha0.1_kvanputt'
    33 time_dir3 = '20081031_133624_run_final_0_27255_alpha0.1_kvanputt'
    34 time_dir4 = '20081031_133735_run_final_0.6_27255_alpha0.1_kvanputt'
     33##time_dir3 = '20081031_133624_run_final_0_27255_alpha0.1_kvanputt'
     34##time_dir4 = '20081031_133735_run_final_0.6_27255_alpha0.1_kvanputt'
    3535##time_dir5 = '20081031_133841_run_final_0_27283_alpha0.1_kvanputt'
    3636##time_dir6 = '20081031_133925_run_final_0.6_27283_alpha0.1_kvanputt'
    3737
    3838
    39 time_dirs = [time_dir3, time_dir4] #, time_dir3, time_dir4, time_dir5, time_dir6]
     39time_dirs = [time_dir1] #3, time_dir4] #, time_dir3, time_dir4, time_dir5, time_dir6]
    4040
    4141for time_dir in time_dirs:
     
    4848    ocean = scenario_dir + "map_work\\Perth.gdb\\Outlines\\initial_conditions_ocean1"
    4949   
    50 ##    print 'Process: Create File GDB'
    51 ##    gp.CreateFileGDB_management(folder, "raster")
     50    print 'Process: Create File GDB'
     51    gp.CreateFileGDB_management(folder, "raster")
    5252
    5353    gp.Workspace = raster_gbd
     
    6262
    6363    generate_filename = []
    64     input_ascii = glob.glob(folder + 'perth_Sorrento*speed_max.asc')
     64    input_ascii = glob.glob(folder + '*.asc')
    6565    print time_dir
    6666
  • anuga_work/production/perth/export_results_max.py

    r6067 r6391  
    2020directory = project.output_dir
    2121
     22time_dir1 = '20090129_155526_run_final_0.6_27255_Rottnest_kvanputt'
     23
     24
    2225##time_dir1 = '20081209_155431_run_final_0_27255_250m_none_lfountai'
    2326##time_dir2 = '20081210_100528_run_final_0_68693_250m_none_kvanputt'
     
    2932
    3033##time_dir1 = '20081031_133353_run_final_0.6_68693_alpha0.1_kvanputt'
    31 time_dir2 = '20081031_133511_run_final_0_68693_alpha0.1_kvanputt'
    32 time_dir3 = '20081031_133624_run_final_0_27255_alpha0.1_kvanputt'
     34##time_dir2 = '20081031_133511_run_final_0_68693_alpha0.1_kvanputt'
     35##time_dir3 = '20081031_133624_run_final_0_27255_alpha0.1_kvanputt'
    3336##time_dir4 = '20081031_133735_run_final_0.6_27255_alpha0.1_kvanputt'
    3437##time_dir5 = '20081031_133841_run_final_0_27283_alpha0.1_kvanputt'
    3538##time_dir6 = '20081031_133925_run_final_0.6_27283_alpha0.1_kvanputt'
    3639##
    37 time_dirs = [time_dir2, time_dir3] #, time_dir3, time_dir4, time_dir5, time_dir6]
     40time_dirs = [time_dir1] #2, time_dir3] #, time_dir3, time_dir4, time_dir5, time_dir6]
    3841
    3942 
    4043
    41 #cellsize = 5
    42 cellsize = 250
     44cellsize = 20
     45#cellsize = 250
    4346
    4447timestep = None    # None means no timestep!
     
    4952######
    5053
    51 #area = ['Sorrento', 'Geordie' ]  # strings must match keys in var_equations below
    52 area = ['All']      # 'All' means no special areas - the whole thing
     54area = ['Rottnest']  # strings must match keys in var_equations below
     55#area = ['All']      # 'All' means no special areas - the whole thing
    5356
    5457######
     
    6366
    6467# one or more key strings from var_equations above
    65 var = ['stage']
     68var = ['depth', 'speed']
    6669
    6770######
     
    9396
    9497            name1 = directory+time_dir+sep+project.scenario_name
    95             name2 = directory+time_dir+sep+project.scenario_name+'_time_39900_0'
    96             name3 = directory+time_dir+sep+project.scenario_name+'_time_79800_0'
     98##            name2 = directory+time_dir+sep+project.scenario_name+'_time_39900_0'
     99##            name3 = directory+time_dir+sep+project.scenario_name+'_time_79800_0'
    97100
    98             names = [name1, name2, name3]
     101            names = [name1]#, name2, name3]
    99102     
    100             asc_name = []   
     103       #     asc_name = []   
    101104               
    102105            for name in names:
     
    118121                            verbose = True,
    119122                            format = 'asc')
    120 
    121                 asc_name.append(outname + '.asc')
    122    
    123             maxasc_outname = directory+time_dir+sep+project.scenario_name+'_'+which_area+'_'+which_var+'_max.asc'
    124 
    125             print 'max asc outname ', maxasc_outname
    126             print 'asc_name ', str(asc_name)
    127            
    128             MaxAsc(maxasc_outname, asc_name)
     123##
     124##                asc_name.append(outname + '.asc')
     125##   
     126##            maxasc_outname = directory+time_dir+sep+project.scenario_name+'_'+which_area+'_'+which_var+'_max.asc'
     127##
     128##            print 'max asc outname ', maxasc_outname
     129##            print 'asc_name ', str(asc_name)
     130##           
     131##            MaxAsc(maxasc_outname, asc_name)
  • anuga_work/production/perth/project.py

    r6342 r6391  
    4444# run_scenario script is executed
    4545tide = 0.6               #0.6
    46 #event_number = 27255 # Java 9.3 worst case for Perth
     46event_number = 27255 # Java 9.3 worst case for Perth
    4747#event_number = 68693 # Sumatra 9.2
    48 event_number = 27283  # Java 9.3 original
     48#event_number = 27283  # Java 9.3 original
    4949alpha = 0.1             # smoothing parameter for mesh
    5050friction=0.01           # manning's friction coefficient
     
    7777# Important to distinguish each run - ensure str(user) is included!
    7878# Note, the user is free to include as many parameters as desired
    79 dir_comment='_'+setup+'_'+str(tide)+'_'+str(event_number)+'_'+ 'Rerun_'+str(user)
     79dir_comment='_'+setup+'_'+str(tide)+'_'+str(event_number)+'_'+ 'Rottnest_'+str(user)
    8080
    8181#------------------------------------------------------------------------------
     
    168168
    169169# where the mesh sits (this is created during the run_perth.py)
    170 meshes_dir_name = meshes_dir + scenario_name+'.msh'
     170meshes_dir_name = meshes_dir + scenario_name+'_rottnest.msh'
    171171
    172172# where the boundary ordering files sit (this is used within build_boundary.py)
     
    213213##res_aoi1 = 500*res_factor
    214214
    215 # Area of Interest 1 (Fremantle)
    216 poly_aoi1 = read_polygon(polygons_dir+'CBD_coastal.csv')
    217 res_aoi1 = 500*res_factor
    218 
    219 # Area of Interest 2 (Rockingham)
    220 poly_aoi2 = read_polygon(polygons_dir+'rockingham_penguin_increase.csv')
    221 res_aoi2 = 500*res_factor
    222 
    223 # Area of Interest 2 (Rockingham)
    224 poly_aoi2 = read_polygon(polygons_dir+'rockingham_penguin.csv')
    225 res_aoi2 = 500*res_factor
    226 
    227 # Area of Interest 2 (garden Island)
    228 poly_aoi2a = read_polygon(polygons_dir+'garden.csv')
    229 res_aoi2a= 500*res_factor
    230 
    231 # Area of Interest 3 (geordie bay - record of tsunami impact)
    232 poly_aoi3 = read_polygon(polygons_dir+'geordie_bay.csv')
    233 res_aoi3 = 500*res_factor
    234 
    235 # Area of Interest 4 (sorrento - record of tsunami impact)
    236 poly_aoi4 = read_polygon(polygons_dir+'sorrento_gauge.csv')
    237 res_aoi4 = 500*res_factor
    238 
    239 # Area of Significance 1 (Garden Island and sand bank infront of Rockingham)
    240 poly_aos1 = read_polygon(polygons_dir+'garden_rockingham_increase.csv')
    241 res_aos1 = 1000*res_factor
     215### Area of Interest 1 (Fremantle)
     216##poly_aoi1 = read_polygon(polygons_dir+'CBD_coastal.csv')
     217##res_aoi1 = 500*res_factor
     218##
     219### Area of Interest 2 (Rockingham)
     220##poly_aoi2 = read_polygon(polygons_dir+'rockingham_penguin_increase.csv')
     221##res_aoi2 = 500*res_factor
     222##
     223### Area of Interest 2 (Rockingham)
     224##poly_aoi2 = read_polygon(polygons_dir+'rockingham_penguin.csv')
     225##res_aoi2 = 500*res_factor
     226##
     227### Area of Interest 2 (garden Island)
     228##poly_aoi2a = read_polygon(polygons_dir+'garden.csv')
     229##res_aoi2a= 500*res_factor
     230##
     231### Area of Interest 3 (geordie bay - record of tsunami impact)
     232##poly_aoi3 = read_polygon(polygons_dir+'geordie_bay.csv')
     233##res_aoi3 = 500*res_factor
     234##
     235### Area of Interest 4 (sorrento - record of tsunami impact)
     236##poly_aoi4 = read_polygon(polygons_dir+'sorrento_gauge.csv')
     237##res_aoi4 = 500*res_factor
     238##
     239### Area of Significance 1 (Garden Island and sand bank infront of Rockingham)
     240##poly_aos1 = read_polygon(polygons_dir+'garden_rockingham_increase.csv')
     241##res_aos1 = 1000*res_factor
    242242
    243243# Area of Significance 2 (incorporate coastline of rottnest)
    244244poly_aos2 = read_polygon(polygons_dir+'rottnest_external.csv')
    245 res_aos2 = 1000*res_factor
    246 
    247 # Refined areas
    248 # Polygon designed to incorporate dredged area from Fremantle to
    249 # Rockingham as the steep incline was making the elevation go to 0
    250 poly_aos3 = read_polygon(polygons_dir+'DredgeArea.csv')
    251 res_aos3 = 1000*res_factor
    252 
    253 # Shallow water 1
    254 poly_sw1 = read_polygon(polygons_dir+'internal_h20mORd3km.csv')
    255 res_sw1 = 25000*res_factor
     245res_aos2 = 500*res_factor
     246
     247### Refined areas
     248### Polygon designed to incorporate dredged area from Fremantle to
     249### Rockingham as the steep incline was making the elevation go to 0
     250##poly_aos3 = read_polygon(polygons_dir+'DredgeArea.csv')
     251##res_aos3 = 1000*res_factor
     252##
     253### Shallow water 1
     254##poly_sw1 = read_polygon(polygons_dir+'internal_h20mORd3km.csv')
     255##res_sw1 = 25000*res_factor
     256
     257### Shallow water 1
     258poly_sw1 = read_polygon(polygons_dir+'rottnest_20m.csv')
     259res_sw1 = 1000*res_factor
    256260
    257261# Deep water (land of Rottnest, ANUGA does not do donuts!)
     
    264268##                     ,[poly_sw1,res_sw1], [poly_dw1,res_dw1]]
    265269
    266 interior_regions = [[poly_aoi1,res_aoi1],[poly_aoi2,res_aoi2]
    267                      ,[poly_aoi2a,res_aoi2a],[poly_aoi3,res_aoi3]
    268                      ,[poly_aoi4,res_aoi4],[poly_aos1,res_aos1]
    269                      ,[poly_aos2,res_aos2],[poly_aos3,res_aos3]
    270                      ,[poly_sw1,res_sw1], [poly_dw1,res_dw1]]
     270#For Rottnest zoom in
     271interior_regions = [[poly_aos2,res_aos2],[poly_sw1,res_sw1]
     272                     ,[poly_dw1,res_dw1]]
     273
     274##interior_regions = [[poly_aoi1,res_aoi1],[poly_aoi2,res_aoi2]
     275##                     ,[poly_aoi2a,res_aoi2a],[poly_aoi3,res_aoi3]
     276##                     ,[poly_aoi4,res_aoi4],[poly_aos1,res_aos1]
     277##                     ,[poly_aos2,res_aos2],[poly_aos3,res_aos3]
     278##                     ,[poly_sw1,res_sw1], [poly_dw1,res_dw1]]
    271279
    272280   
     
    309317ymaxPerth = 6467000
    310318
     319#Rottnest
     320xminRottnest = 352700
     321xmaxRottnest = 364800
     322yminRottnest = 6455000
     323ymaxRottnest = 6460500
     324
     325
Note: See TracChangeset for help on using the changeset viewer.