Changeset 6067


Ignore:
Timestamp:
Dec 11, 2008, 4:12:14 PM (15 years ago)
Author:
kristy
Message:

general maintenance

Location:
anuga_work/production/perth
Files:
5 edited

Legend:

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

    r6033 r6067  
    2525output_dir = "anuga\\outputs\\"
    2626
     27##time_dir1 = '20081209_155431_run_final_0_27255_250m_none_lfountai'
     28##time_dir2 = '20081210_100528_run_final_0_68693_250m_none_kvanputt'
     29##time_dirs = [time_dir1, time_dir2]
    2730
    28 time_dir1 = '20081031_133353_run_final_0.6_68693_alpha0.1_kvanputt'
     31##time_dir1 = '20081031_133353_run_final_0.6_68693_alpha0.1_kvanputt'
    2932##time_dir2 = '20081031_133511_run_final_0_68693_alpha0.1_kvanputt'
    30 ##time_dir3 = '20081031_133624_run_final_0_27255_alpha0.1_kvanputt'
    31 ##time_dir4 = '20081031_133735_run_final_0.6_27255_alpha0.1_kvanputt'
     33time_dir3 = '20081031_133624_run_final_0_27255_alpha0.1_kvanputt'
     34time_dir4 = '20081031_133735_run_final_0.6_27255_alpha0.1_kvanputt'
    3235##time_dir5 = '20081031_133841_run_final_0_27283_alpha0.1_kvanputt'
    3336##time_dir6 = '20081031_133925_run_final_0.6_27283_alpha0.1_kvanputt'
    3437
    3538
    36 time_dirs = [time_dir1]#, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
     39time_dirs = [time_dir3, time_dir4] #, time_dir3, time_dir4, time_dir5, time_dir6]
    3740
    3841for time_dir in time_dirs:
     
    4144    folder =  scenario_dir + output_dir + time_dir + '\\'
    4245    raster_gbd = folder + 'raster.gdb'
    43     land = scenario_dir + "map_work\\Perth.gdb\\poly_all_cut"
     46    contour = raster_gbd + '\\contour_dep'
     47    land = scenario_dir + "map_work\\Perth.gdb\\Outlines\\initial_conditions_rottnest"
     48    ocean = scenario_dir + "map_work\\Perth.gdb\\Outlines\\initial_conditions_ocean1"
    4449   
    45     # Process: Create File GDB...
    46     #gp.CreateFileGDB_management(folder, "raster")
     50##    print 'Process: Create File GDB'
     51##    gp.CreateFileGDB_management(folder, "raster")
    4752
    4853    gp.Workspace = raster_gbd
     
    5459                 ('_', ''),('Geordie', 'Geo'),('Sorrento', 'Sor'), ('max','M_'),
    5560                 ('Fremantle', 'Fre'),('Rockingham', 'Roc'),('depth','_dep_'),
    56                  ('speed', '_spe_'), ('elevation', '_ele_'), ('stage','_sta_'))
     61                 ('speed', '_spe_'), ('elevation', '_ele_'), ('stage','_stage'))
    5762
    5863    generate_filename = []
    59     input_ascii = glob.glob(folder + '*elevation.asc')
     64    input_ascii = glob.glob(folder + 'perth_Sorrento*speed_max.asc')
     65    print time_dir
    6066
    6167    for infile in input_ascii:
     
    6975        generate_filename.append(output_DEM)
    7076        print 'Output DEM ',output_DEM
    71         output_extract = output_DEM + 'E'
    72         print 'Output Extract ',output_extract
    7377       
    7478        print 'Process: ASCII to Raster'
     
    8084                                                   ",PARAMETER['Latitude_Of_Origin',0.0],UNIT['Meter',1.0]]")
    8185
     86   
     87##        output_extract = output_DEM + 'E'
     88##        print 'Output Extract ',output_extract
     89##   
    8290##        print 'Process: Extract by Mask'
    8391##        gp.ExtractByMask_sa(output_DEM, land, output_extract)
    84 
    85 
     92##
     93##       
  • anuga_work/production/perth/export_results_max.py

    r6032 r6067  
    2020directory = project.output_dir
    2121
    22 #time_dir1 = '20081031_133841_run_final_0_27283_alpha0.1_kvanputt'
    23 #time_dir1 = '20080924_123626_run_final_0_27283_250m_all_kvanputt' # This uses the 250m bathymetry with all interior polygons
    24 #time_dir2 = '20080912_154716_run_final_0_27283_alpha0.1_kvanputt' # This uses original bathyemetry data
    25 #time_dir3 = '20080924_123601_run_final_0_27283_250m_none_kvanputt' # This uses the 250m bathymetry without any interior polygons
    26 #time_dirs = [time_dir1] #, time_dir2]#, time_dir3]
     22##time_dir1 = '20081209_155431_run_final_0_27255_250m_none_lfountai'
     23##time_dir2 = '20081210_100528_run_final_0_68693_250m_none_kvanputt'
    2724
    28 time_dir1 = '20081031_133353_run_final_0.6_68693_alpha0.1_kvanputt'
    29 ##time_dir2 = '20081031_133511_run_final_0_68693_alpha0.1_kvanputt'
    30 ##time_dir3 = '20081031_133624_run_final_0_27255_alpha0.1_kvanputt'
     25###time_dir1 = '20080924_123626_run_final_0_27283_250m_all_kvanputt' # This uses the 250m bathymetry with all interior polygons
     26###time_dir2 = '20080912_154716_run_final_0_27283_alpha0.1_kvanputt' # This uses original bathyemetry data
     27###time_dir3 = '20080924_123601_run_final_0_27283_250m_none_kvanputt' # This uses the 250m bathymetry without any interior polygons
     28##time_dirs = [time_dir1, time_dir2]#, time_dir3]
     29
     30##time_dir1 = '20081031_133353_run_final_0.6_68693_alpha0.1_kvanputt'
     31time_dir2 = '20081031_133511_run_final_0_68693_alpha0.1_kvanputt'
     32time_dir3 = '20081031_133624_run_final_0_27255_alpha0.1_kvanputt'
    3133##time_dir4 = '20081031_133735_run_final_0.6_27255_alpha0.1_kvanputt'
    3234##time_dir5 = '20081031_133841_run_final_0_27283_alpha0.1_kvanputt'
    3335##time_dir6 = '20081031_133925_run_final_0.6_27283_alpha0.1_kvanputt'
    34 
    35 time_dirs = [time_dir1]#, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
     36##
     37time_dirs = [time_dir2, time_dir3] #, time_dir3, time_dir4, time_dir5, time_dir6]
    3638
    3739 
    3840
    39 cellsize = 20
    40 #cellsize = 250
     41#cellsize = 5
     42cellsize = 250
    4143
    42 #timestep = None    # None means no timestep!
    43 timestep = 0
     44timestep = None    # None means no timestep!
     45#timestep = 0
    4446
    4547######
     
    4749######
    4850
    49 #area = ['Perth'] #Rockingham', 'Sorrento', 'Fremantle','Geordie' ]  # strings must match keys in var_equations below
     51#area = ['Sorrento', 'Geordie' ]  # strings must match keys in var_equations below
    5052area = ['All']      # 'All' means no special areas - the whole thing
    5153
     
    6163
    6264# one or more key strings from var_equations above
    63 var = ['elevation']
     65var = ['stage']
    6466
    6567######
     
    9193
    9294            name1 = directory+time_dir+sep+project.scenario_name
    93             #name2 = directory+time_dir+sep+project.scenario_name+'_time_39180_0'
    94             #name3 = directory+time_dir+sep+project.scenario_name+'_time_78360_0'
     95            name2 = directory+time_dir+sep+project.scenario_name+'_time_39900_0'
     96            name3 = directory+time_dir+sep+project.scenario_name+'_time_79800_0'
    9597
    96             names = [name1] #, name2, name3]
     98            names = [name1, name2, name3]
    9799     
    98100            asc_name = []   
     
    124126            print 'asc_name ', str(asc_name)
    125127           
    126             maxasc.MaxAsc(maxasc_outname, asc_name)
     128            MaxAsc(maxasc_outname, asc_name)
  • anuga_work/production/perth/get_timeseries.py

    r6032 r6067  
    1818directory = project.output_dir
    1919
    20 ##time_dir1 = '20081031_133353_run_final_0.6_68693_alpha0.1_kvanputt'
    21 ##time_dir2 = '20081031_133511_run_final_0_68693_alpha0.1_kvanputt'
    22 ##time_dir3 = '20081031_133624_run_final_0_27255_alpha0.1_kvanputt'
    23 ##time_dir4 = '20081031_133735_run_final_0.6_27255_alpha0.1_kvanputt'
    24 ##time_dir5 = '20081031_133841_run_final_0_27283_alpha0.1_kvanputt'
    25 time_dir6 = '20081117_141914_run_final_0.6_27283_Rockingham_kvanputt'
     20time_dir1 = '20081031_133353_run_final_0.6_68693_alpha0.1_kvanputt'+sep+'gauge'
     21time_dir2 = '20081031_133511_run_final_0_68693_alpha0.1_kvanputt'+sep+'gauge'
     22time_dir3 = '20081031_133624_run_final_0_27255_alpha0.1_kvanputt'+sep+'gauge'
     23time_dir4 = '20081031_133735_run_final_0.6_27255_alpha0.1_kvanputt'+sep+'gauge'
     24time_dir5 = '20081031_133841_run_final_0_27283_alpha0.1_kvanputt'+sep+'gauge'
     25time_dir6 = '20081031_133925_run_final_0.6_27283_alpha0.1_kvanputt'+sep+'gauge'
    2626
    27 time_dirs = [time_dir6] #1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
     27time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
    2828
    2929for time_dir in time_dirs:
     30    name = directory+time_dir+sep+project.scenario_name
     31    gauge = project.gauges_dir_name
     32    sww2csv_gauges(name+'.sww',gauge,
     33                   quantities = ['stage','speed','depth','elevation'],
     34                   verbose=True)
     35
    3036   
    31     name1 = directory+time_dir+sep+project.scenario_name
    32     name2 = directory+time_dir+sep+project.scenario_name+'_time_43380_0'
    33     #name3 = directory+time_dir+sep+project.scenario_name+'_time_78360_0'
    34 
    35     names = [name1, name2] #, name3]
    36 
    37     for i, name in enumerate(names):
    38         gauge = project.gauges_dir_name
    39         out_name = 'gauge_%d_' % i
    40         sww2csv_gauges(name+'.sww',gauge,out_name = out_name,
    41                         quantities = ['stage','speed','depth','elevation'],
    42                         verbose=True)
    43 ##        for out_name in listdir(directory+time_dir):
    44 ##            #total_filename = directory+time_dir+sep+'All_' + out_name
    45 ##            if i == 0:
    46 ##               
    47 ##            else:
    48 ##                system('cat %s %s > %s' %(,out_name,))
    49 ##   
     37##    name1 = directory+time_dir+sep+project.scenario_name
     38####    name2 = directory+time_dir+sep+project.scenario_name+'_time_43380_0'
     39####    name3 = directory+time_dir+sep+project.scenario_name+'_time_78360_0'
     40##
     41##    names = [name1]#, name2, name3]
     42##
     43##    for i, name in enumerate(names):
     44##        gauge = project.gauges_dir_name
     45####        out_name = 'gauge_%d_' % i
     46##        print name
     47##        sww2csv_gauges(name+'.sww',gauge,#out_name = out_name,
     48##                        quantities = ['stage','speed','depth','elevation'],
     49##                        verbose=True)
     50##
  • anuga_work/production/perth/project.py

    r6032 r6067  
    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)+'_'+ 'Rockingham_'+str(user)
     79dir_comment='_'+setup+'_'+str(tide)+'_'+str(event_number)+'_'+ 'Rerun_'+str(user)
    8080
    8181#------------------------------------------------------------------------------
  • anuga_work/production/perth/run_building_inundation.py

    r6031 r6067  
    5656       newline = record[-8:]
    5757       newline.insert(0,record[0])
     58       newline.insert(1,record[1])
     59       newline.insert(2,record[2])
    5860       writer.writerow(newline)
    5961       
Note: See TracChangeset for help on using the changeset viewer.