Changeset 7082


Ignore:
Timestamp:
May 26, 2009, 10:46:11 AM (11 years ago)
Author:
kristy
Message:
 
Location:
anuga_work/production/hobart_2009
Files:
5 added
3 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/hobart_2009/export_results_max.py

    r6984 r7082  
    2121
    2222directory = project.output_folder
    23 time_dir1 = '20090420_153539_run_final_0_58274_kvanputt'
    24 time_dir2 = '20090420_153654_run_final_0_58279_kvanputt'
    25 time_dir3 = '20090420_153739_run_final_0_58281_kvanputt'
    26 time_dir4 = '20090420_153826_run_final_0_58078_kvanputt'
    27 time_dir5 = '20090330_093828_run_final_0_58292_kvanputt'
    28 time_dir6 = '20090504_123724_run_final_0_58282_small_kvanputt'
    29 time_dir7 = '20090504_124158_run_final_0_58068_small_kvanputt'
     23time_dir1 = '20090505_150430_run_final_0.8_58292_None_kvanputt'
     24time_dir2 = '20090505_150517_run_final_0_58292_None_kvanputt'
     25time_dir3 = '20090505_150711_run_final_0_58280_None_kvanputt'
     26time_dir4 = '20090505_150805_run_final_0.8_58280_None_kvanputt'
     27time_dir5 = '20090505_151322_run_final_0.8_64477_None_kvanputt'
     28time_dir6 = '20090505_151447_run_final_0_64477_None_kvanputt'
    3029
    31 time_dirs = [time_dir6, time_dir7]
    32 #time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5]
     30time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
    3331 
    3432
     
    4341######
    4442
    45 area = ['SArms', 'Bruny']  # strings must match keys in var_equations below
     43area = ['Hobart', 'NW', 'South']  # strings must match keys in var_equations below
    4644#area = ['All']      # 'All' means no special areas - the whole thing
    4745
     
    5755
    5856# one or more key strings from var_equations above
    59 var = ['depth', 'speed']
     57var = ['speed']
    6058
    6159######
     
    8684        for time_dir in time_dirs:
    8785            name1 = join(directory, time_dir, project.scenario_name)
    88             name2 = join(directory, time_dir, project.scenario_name)+'_time_22680_0'
    89             name3 = join(directory, time_dir, project.scenario_name)+'_time_45360_0'
    90             name4 = join(directory, time_dir, project.scenario_name)+'_time_68040_0'
     86            name2 = join(directory, time_dir, project.scenario_name)+'_time_16980_0'
     87            name3 = join(directory, time_dir, project.scenario_name)+'_time_33960_0'
     88            name4 = join(directory, time_dir, project.scenario_name)+'_time_50940_0'
    9189           
    92             names = [name1] #, name2, name3, name4]
     90            names = [name1, name2, name3, name4]
    9391     
    94 ##            asc_name = []   
     92            asc_name = []   
    9593               
    9694            for name in names:
     
    113111                            format = 'asc')
    114112
    115 ##                asc_name.append(outname + '.asc')
    116 ##   
    117 ##            maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc'
    118 ##
    119 ##            print 'max asc outname ', maxasc_outname
    120 ##            print 'asc_name ', str(asc_name)
    121 ##           
    122 ##            MaxAsc(maxasc_outname, asc_name)
     113                asc_name.append(outname + '.asc')
     114   
     115            maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc'
     116
     117            print 'max asc outname ', maxasc_outname
     118            print 'asc_name ', str(asc_name)
     119           
     120            MaxAsc(maxasc_outname, asc_name)
  • anuga_work/production/hobart_2009/project.py

    r6984 r7082  
    128128# Used in get_timeseries.py. 
    129129# Format easting,northing,name,elevation (with header)
    130 gauges_filename = 'tsunamipointsMGA.csv'
     130gauges_filename = 'time_of_arrival.csv' #'TideGaugesPoints.csv' #'tsunamipointsMGA.csv'
    131131
    132132# BUILDINGS EXPOSURE - for identifying inundated houses
     
    134134# Format latitude,longitude etc (geographic)
    135135building_exposure_filename = 'busselton_res_clip.csv' # from NEXIS
     136
     137# AREA OF IMAGES - Extent of each image to find out highest runup
     138# Used in get_runup.py
     139images_filename = 'images.csv'
    136140
    137141# BOUNDING POLYGON - used in build_boundary.py and run_model.py respectively
     
    173177yminSArms = 5237028
    174178ymaxSArms = 5238465
     179
     180# ASCII export grid for Hobart Large
     181xminHobart = 520000
     182xmaxHobart = 545000
     183yminHobart = 5230000
     184ymaxHobart = 5260000
     185
     186# ASCII export grid for North West Communities
     187xminNW = 548000
     188xmaxNW = 561000
     189yminNW = 5250000
     190ymaxNW = 5258000
     191
     192# ASCII export grid for South Communities
     193xminSouth = 523000
     194xmaxSouth = 533300
     195yminSouth = 5197000
     196ymaxSouth = 5214000
    175197
    176198
     
    273295    building_exposure = join(gauges_folder, building_exposure_filename)
    274296
     297# The absolute pathname for the image file
     298# Used for get_runup.py
     299if images_filename:
     300    images = join(polygons_folder, images_filename)
     301
     302
    275303# full path to where MUX files (or meta-files) live
    276304mux_input = join(event_folder, mux_input_filename)
  • anuga_work/production/hobart_2009/setup_model.py

    r6984 r7082  
    154154project.interior_regions = []
    155155
    156 if not project.PriorityArea_filename == None:
     156if project.PriorityArea_filename is not None:
    157157    polygons, maxareas = csv2polygons(project.PriorityAreas)
    158158    print 'Creating %d internal polygons' % len(polygons)
Note: See TracChangeset for help on using the changeset viewer.