Changeset 6065


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

general maintenance

Location:
anuga_work/production/geraldton
Files:
2 edited

Legend:

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

    r6041 r6065  
    2222gp.overwriteoutput = 1
    2323
    24 output_dir = "\\\\nas2\\gemd\\georisk_models\\inundation\\data\\western_australia\\"
    25 scenario_dir="geraldton_tsunami_scenario\\anuga\\outputs\\"
     24scenario_dir = "\\\\nas2\\gemd\\georisk_models\\inundation\\data\\western_australia\\geraldton_tsunami_scenario\\"
     25output_dir="anuga\\outputs\\"
    2626
     27##time_dir1 = '20081209_155330_run_final_0_27255_250m_none_lfountai'
     28##time_dir2 = '20081209_161229_run_final_0_68693_250m_none_lfountai'
     29##time_dirs = [time_dir1, time_dir2]
    2730
    28 time_dir1 = '20081117_141353_run_final_0_27283_alpha0.1_kvanputt'
    29 time_dir2 = '20081117_141443_run_final_0.6_27283_alpha0.1_kvanputt'
     31##time_dir1 = '20081117_141353_run_final_0_27283_alpha0.1_kvanputt'
     32##time_dir2 = '20081117_141443_run_final_0.6_27283_alpha0.1_kvanputt'
    3033time_dir3 = '20081117_141520_run_final_0_27255_alpha0.1_kvanputt'
    31 time_dir4 = '20081117_141558_run_final_0.6_27255_alpha0.1_kvanputt'
     34##time_dir4 = '20081117_141558_run_final_0.6_27255_alpha0.1_kvanputt'
    3235time_dir5 = '20081117_141634_run_final_0_68693_alpha0.1_kvanputt'
    33 time_dir6 = '20081117_141717_run_final_0.6_68693_alpha0.1_kvanputt'
    34 
    35 time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
     36##time_dir6 = '20081117_141717_run_final_0.6_68693_alpha0.1_kvanputt'
     37##
     38time_dirs = [time_dir3, time_dir5] #, time_dir3, time_dir4, time_dir5, time_dir6]
    3639
    3740for time_dir in time_dirs:
    3841
    3942    # Local variables...
    40     folder = output_dir + scenario_dir + time_dir +'\\'
     43    folder = scenario_dir + output_dir + time_dir +'\\'
    4144    raster_gbd = folder + 'raster.gdb'
    42     land = scenario_dir + "map_work\\\Geraldton.gdb\\outlines\\Clip_Topography"
     45    land = scenario_dir + "map_work\\\Geraldton.gdb\\outlines\\initial_condition"
     46    ocean = scenario_dir + "map_work\\\Geraldton.gdb\\outlines\\initial_conditions_ocean"
    4347
    44     # Process: Create File GDB...
    45     #gp.CreateFileGDB_management(folder, "raster")
     48##    print 'Process: Create File GDB'
     49##    gp.CreateFileGDB_management(folder, "raster")
    4650
    4751    gp.Workspace = raster_gbd
     
    5054   
    5155    #replication dictionary
    52     replicate = (('geraldton', ''),('_', ''),('max','M_'),
     56    replicate = (('geraldton', ''),('_', ''),('max','_M'),
    5357                 ('CBD', 'CDB'),
    5458                 ('depth','_depth'),('speed', '_speed'),
    55                  ('elevation', '_ele_'), ('stage','_sta_'))
     59                 ('elevation', '_ele_'), ('stage','_stage'))
    5660
    5761    generate_filename = []
    58     input_ascii = glob.glob(folder + '*max.asc')
     62    input_ascii = glob.glob(folder + '*stage_max.asc')
    5963
    6064    for infile in input_ascii:
     
    6973
    7074        print 'Output DEM ',output_DEM
    71         output_extract = output_DEM + 'E'
    72         print 'Output Extract ',output_extract
    7375       
    7476        print 'Process: ASCII to Raster'
     
    7981                                                   ",PROJECTION['Transverse_Mercator'],PARAMETER['False_Easting',500000.0],PARAMETER['False_Northing',10000000.0],PARAMETER['Central_Meridian',117.0],PARAMETER['Scale_Factor',0.9996]"
    8082                                                   ",PARAMETER['Latitude_Of_Origin',0.0],UNIT['Meter',1.0]]")
    81 
    82         print 'Process: Extract by Mask'
    83         gp.ExtractByMask_sa(output_DEM, land, output_extract)
     83##        output_extract = output_DEM + 'E'
     84##        print 'Output Extract ',output_extract
     85##        print 'Process: Extract by Mask'
     86##        gp.ExtractByMask_sa(output_DEM, land, output_extract)
    8487
    8588
  • anuga_work/production/geraldton/export_results_max.py

    r6040 r6065  
    2020
    2121directory = project.output_dir
     22##time_dir1 = '20081209_155330_run_final_0_27255_250m_none_lfountai'
     23##time_dir2 = '20081209_161229_run_final_0_68693_250m_none_lfountai'
     24##time_dirs = [time_dir1, time_dir2]
    2225
    23 time_dir1 = '20081117_141353_run_final_0_27283_alpha0.1_kvanputt'
    24 time_dir2 = '20081117_141443_run_final_0.6_27283_alpha0.1_kvanputt'
     26##time_dir1 = '20081117_141353_run_final_0_27283_alpha0.1_kvanputt'
     27##time_dir2 = '20081117_141443_run_final_0.6_27283_alpha0.1_kvanputt'
    2528time_dir3 = '20081117_141520_run_final_0_27255_alpha0.1_kvanputt'
    26 time_dir4 = '20081117_141558_run_final_0.6_27255_alpha0.1_kvanputt'
     29##time_dir4 = '20081117_141558_run_final_0.6_27255_alpha0.1_kvanputt'
    2730time_dir5 = '20081117_141634_run_final_0_68693_alpha0.1_kvanputt'
    28 time_dir6 = '20081117_141717_run_final_0.6_68693_alpha0.1_kvanputt'
    29 
    30 time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
     31##time_dir6 = '20081117_141717_run_final_0.6_68693_alpha0.1_kvanputt'
     32##
     33time_dirs = [time_dir3, time_dir5] #, time_dir3, time_dir4, time_dir5, time_dir6]
    3134
    3235 
    3336
    34 cellsize = 20
    35 #cellsize = 250
     37#cellsize = 10
     38cellsize = 250
    3639
    3740timestep = None    # None means no timestep!
     
    4245######
    4346
    44 area = ['CBD']  # strings must match keys in var_equations below
    45 #area = ['All']      # 'All' means no special areas - the whole thing
     47#area = ['CBD']  # strings must match keys in var_equations below
     48area = ['All']      # 'All' means no special areas - the whole thing
    4649
    4750######
     
    5659
    5760# one or more key strings from var_equations above
    58 var = ['depth', 'speed']
     61var = ['stage']
    5962
    6063######
Note: See TracChangeset for help on using the changeset viewer.