Changeset 7385


Ignore:
Timestamp:
Aug 20, 2009, 9:20:31 AM (16 years ago)
Author:
kristy
Message:
 
Location:
anuga_work/production
Files:
3 deleted
15 edited

Legend:

Unmodified
Added
Removed

  • anuga_work/production/busselton/busselton_rerun/Arc_asc2raster_GDA94z50.py

    r6984 r7385  
    3232time_dir6 = '20090330_120458_run_final_0_27283_Bt_kvanputt'
    3333
    34 time_dirs = [time_dir6] #1, time_dir2, time_dir3, time_dir4, time_dir5]
     34time_dirs = [time_dir1] #1, time_dir2, time_dir3, time_dir4, time_dir5]
    3535   
    3636for time_dir in time_dirs:
    … …  
    5757
    5858    generate_filename = []
    59     input_ascii = glob.glob(folder + '*depth_max.asc')
     59    input_ascii = glob.glob(folder + '*elevation.asc')
    6060
    6161    for infile in input_ascii:
    … …  
    7878                                                   ",PROJECTION['Transverse_Mercator'],PARAMETER['False_Easting',500000.0],PARAMETER['False_Northing',10000000.0],PARAMETER['Central_Meridian',117.0],PARAMETER['Scale_Factor',0.9996]"
    7979                                                   ",PARAMETER['Latitude_Of_Origin',0.0],UNIT['Meter',1.0]]")
    80         output_extract = output_DEM + '_E'
    81         print 'Output Extract ',output_extract
    82         print 'Process: Extract by Mask'
    83         gp.ExtractByMask_sa(output_DEM, land, output_extract)
     80##        output_extract = output_DEM + '_E'
     81##        print 'Output Extract ',output_extract
     82##        print 'Process: Extract by Mask'
     83##        gp.ExtractByMask_sa(output_DEM, land, output_extract)
    8484
    8585

  • anuga_work/production/carnarvon/carnarvon_rerun/Arc_asc2raster_GDA94z49.py

    r6984 r7385  
    3232time_dir6 = '20090327_160753_run_final_0_27283_Bt_kvanputt'
    3333
    34 time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]   
     34time_dirs = [time_dir1] #, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]   
    3535
    3636for time_dir in time_dirs:
    … …  
    5555
    5656    generate_filename = []
    57     input_ascii = glob.glob(folder + '*depth_max.asc')
     57    input_ascii = glob.glob(folder + '*elevation.asc')
    5858
    5959    for infile in input_ascii:
    … …  
    8080
    8181 
    82         print 'Process: Extract by Mask'
    83         output_extract = output_DEM + '_E'
    84         print 'Output Extract ',output_extract
    85         gp.ExtractByMask_sa(output_DEM, land, output_extract)
     82##        print 'Process: Extract by Mask'
     83##        output_extract = output_DEM + '_E'
     84##        print 'Output Extract ',output_extract
     85##        gp.ExtractByMask_sa(output_DEM, land, output_extract)
    8686
    8787

  • anuga_work/production/dampier_2009/Arc_asc2raster_GDA94z50.py

    r7272 r7385  
    2222gp.overwriteoutput = 1
    2323
    24 scenario_dir = "\\\\nas2\\gemd\\georisk_models\\inundation\\data\\western_australia\\onslow_tsunami_scenario_2009\\"
     24scenario_dir = "\\\\nas2\\gemd\\georisk_models\\inundation\\data\\western_australia\\dampier_tsunami_scenario_2009\\"
    2525output_dir="anuga\\outputs\\"
    2626
    27 time_dir1 = '20090616_154908_run_final_0.0_27283_mhingee'
    28 ##time_dir2 = '20090416_072538_run_final_0_27338_1307_Tb__kvanputt'
    29 ##time_dir3 = '20090417_021634_run_final_0_70337_1307_Tb__kvanputt'
     27time_dir1 = '20090817_150959_run_final_0.0_70408_kvanputt'
    3028
    31 ##time_dir2 = '20090422_102457_run_final_0_70844_1307_Tb_internal_mhingee'
    32 
    33 time_dirs = [time_dir1]#, time_dir2]#, time_dir3]
     29time_dirs = [time_dir1]
    3430
    3531for time_dir in time_dirs:
    … …  
    4844   
    4945    #replication dictionary
    50     replicate = (('onslow', ''),('_', ''),('max','_M'),
    51                  ('CBD', 'CDB'),('All',''),
     46    replicate = (('dampier', ''),('_', ''),('max','_M'),
     47                 ('CBD', 'CDB'),('Dampier', 'Damp'),('Karratha', 'Karr')
    5248                 ('depth','_depth'),('speed', '_speed'),
    53                  ('elevation', '_ele_'), ('stage','_stage'))
     49                 ('elevation', '_ele'), ('stage','_stage'))
    5450
    5551    generate_filename = []
    56     input_ascii = glob.glob(folder + '*elevation.asc')
     52    input_ascii = glob.glob(folder + '*.asc')
    5753
    5854    for infile in input_ascii:
    … …  
    7268
    7369        print 'Process: Define Projection'
    74         #GDA_1994_MGA_Zone_54
    7570        gp.DefineProjection_management(output_DEM, "PROJCS['GDA_1994_MGA_Zone_50',GEOGCS['GCS_GDA_1994',DATUM['D_GDA_1994',SPHEROID['GRS_1980',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]]"
    7671                                                   ",PROJECTION['Transverse_Mercator'],PARAMETER['False_Easting',500000.0],PARAMETER['False_Northing',10000000.0],PARAMETER['Central_Meridian',117.0],PARAMETER['Scale_Factor',0.9996]"

  • anuga_work/production/dampier_2009/project.py

    r7272 r7385  
    2929zone = 50               # specify UTM zone of model
    3030central_meridian = None
    31 event_number = 27283    # 27255,27283 the event number or the mux file name
     31event_number = 70408    # 27255,27283 the event number or the mux file name
    3232alpha = 0.1             # smoothing parameter for mesh
    3333friction=0.01           # manning's friction coefficient
    … …  
    6060# Used in build_elevation.py
    6161# Format for ascii grids, as produced in ArcGIS + a projection file
    62 ascii_grid_filenames = []
     62ascii_grid_filenames = ['topo_to_raster_0_30', 'land_final_240107']
    6363
    6464# Format for point is x,y,elevation (with header)
    65 point_filenames = []
     65point_filenames = ['Hydro_fair_chart.txt', 'coast_line.txt']
    6666
    6767### Add csv header list to all files in point_filenames
    … …  
    7979# Used in run_model.py
    8080# Format for points easting,northing (no header)                   
    81 interior_regions_data = [['dampier_region.csv', 3000],
    82                          ['dampier_aos.csv', 800],
    83                          ['dampier_aoi.csv', 500],
    84                          ['dampier_north_aos.csv', 800],
    85                          ['dampier_north_aoi.csv', 500],
    86                          ['karratha_aos.csv', 800],
    87                          ['karratha_aoi.csv', 500],
    88                          ['ref_ne_islands1.csv', 1500],  ## islands to the NE and headland between dampier and karratha
    89                          ['ref_ne_islands2.csv', 1500],  ## island in entrance to bay for Karratha
    90 ##                         ['nwislands.csv', 1500],   ## all islands to the NW:
    91                          ['ref_nw4.csv', 1500],
    92                          ['ref_nw5.csv', 1500],
    93                          ['ref_nw6.csv', 1500],
    94                          ['ref_nw7.csv', 1500],
    95                          ['ref_nw8.csv', 1500],
    96                          ['island0.csv', 10000],   ## make large triangles for land areas of islands to reduce no. of triangles
    97                          ['island1.csv', 10000],
    98                          ['island2.csv', 10000],
    99                          ['island3.csv', 10000],
    100                          ['island4.csv', 10000],
    101                          ['island5.csv', 10000],
    102                          ['island6.csv', 10000],
    103                          ['island7.csv', 10000],
    104                          ['island8.csv', 10000]]
    105 
    106 PriorityArea_filename = None
     81interior_regions_data = [['dampier_aoi.csv', 500],
     82                         ['karratha_aoi.csv', 500]]
     83
     84PriorityArea_filename = 'aos.csv'
    10785   
    10886# LAND - used to set the initial stage/water to be offcoast only
    10987# Used in run_model.py.  Format for points easting,northing,id,value
    110 land_initial_conditions_filename = [['mainland_only.csv',0],  ## these are all placeholders until proper initial conditions made
    111                                     ['island0.csv',0]]
     88land_initial_conditions_filename = 'initial_conditions_simple.csv'
    11289
    11390# GAUGES - for creating timeseries at a specific point
    … …  
    153130
    154131# ASCII export grid for Dampier
    155 ##xmin =
    156 ##xmax =
    157 ##ymin =
    158 ##ymax =
     132xminDampier = 465000
     133xmaxDampier = 471000
     134yminDampier = 7713000
     135ymaxDampier = 7719000
    159136
    160137# ASCII export grid for Karratha
    161 ##xmin =
    162 ##xmax =
    163 ##ymin =
    164 ##ymax =
    165 
    166 # ASCII export grid for Dampier northern area
    167 ##xmin =
    168 ##xmax =
    169 ##ymin =
    170 ##ymax =
    171 
     138xminKarratha = 475700
     139xmaxKarratha = 489600
     140yminKarratha = 7707000
     141ymaxKarratha = 7718000
    172142
    173143################################################################################

  • anuga_work/production/dampier_2009/setup_model.py

    r7272 r7385  
    2727    sanity_error = True
    2828
    29 if os.getenv(project.ENV_MUXHOME) is None:
    30     print "Environment variable '%s' is not set" % project.ENV_MUXHOME
    31     sanity_error = True
    32 
    3329#-------------------------------------------------------------------------------
    3430# Check filename strings that MUST be set
    … …  
    5248if not exists(project.home):
    5349    print "Sorry, data directory '%s' doesn't exist" % project.home
    54     sanity_error = True
    55    
    56 if not exists(project.muxhome):
    57     print "Sorry, MUX directory '%s' doesn't exist" % project.muxhome
    5850    sanity_error = True
    5951   
    … …  
    8476    print "Sorry, gauges directory '%s' doesn't exist" % project.gauges_folder
    8577    sanity_error = True
    86 
    87 if not exists(project.mux_data_folder):
    88     print "Sorry, mux data directory '%s' doesn't exist" % project.mux_data_folder
    89     sanity_error = True
    90 
    91 # if multi_mux is True, check if multi-mux file exists
    92 if project.multi_mux:
    93     if not exists(project.mux_input):
    94         print ("Sorry, MUX input file '%s' doesn't exist"
    95                % project.mux_input)
    96         sanity_error = True
    9778
    9879#-----
    … …  
    166147           % type(project.land_initial_conditions_filename))
    167148    raise Exception, msg
     149
    168150# Create list of interior polygons with scaling factor
    169151project.interior_regions = []

  • anuga_work/production/exmouth_2009/project.py

    r7274 r7385  
    88from time import localtime, strftime, gmtime
    99from os.path import join, exists
    10 from anuga.lib.add_csv_header.add_csv_header import add_csv_header
    1110
    1211#-------------------------------------------------------------------------------
    … …  
    2625# Model specific parameters.
    2726# One or all can be changed each time the run_model script is executed
    28 tide = 0.0              # 3.6 difference between MSL and HAT in metres
    29 zone = None               # specify UTM zone of model
    30 central_meridian = 114.0
    31 event_number = 27283    # 27255,27283 the event number or the mux file name
     27tide = 0                # difference between MSL and HAT in metres
     28zone = None
     29central_meridian = 114   # specify the central merdian for the data
     30event_number = 27283    # the event number or the mux file name
    3231alpha = 0.1             # smoothing parameter for mesh
    3332friction=0.01           # manning's friction coefficient
    … …  
    6059# Used in build_elevation.py
    6160# Format for ascii grids, as produced in ArcGIS + a projection file
    62 ##ascii_grid_filenames = ['extract_dli',
    63 ##                        'extract_dted']
    64 ascii_grid_filenames = []
    65 
    66 # Format for point is x,y,elevation (with header)
     61ascii_grid_filenames = ['extract_dli', 'extract_dted'] # onshore elevation data
     62
     63
     64    # Format for point is x,y,elevation (with header)
    6765point_filenames = ['coastline.txt',  # coastline points with 0m elevation
    6866                   'Exmouth.txt',    # combined bathymetry for most of model
    69                    'extract_dgapcs100.txt',  # interpolated bathy. for gap in data just NE of Exmouth itself
    70                    'extract_rivoli.txt',     # interpolated bathy. for gap in data near Rivoli islands
    71                    'extract_dli.txt',        # onshore elevation data
    72                    'extract_dted.txt',       # more onshore elevation data
    73                    'exmouth_aoi_aos_gridded.txt']  # fine grid of data interpolated from Exmouth.txt for use in refined grids (aoi and aos)
     67                   'extract_dgap_cs100.txt',  # interpolated bathy. for gap in data just NE of Exmouth itself
     68                   'extract_rivoli100.txt']     # interpolated bathy. for gap in data near Rivoli islands
     69                   #'exmouth_aoi_aos_gridded.txt']  # fine grid of data interpolated from Exmouth.txt for use in refined grids (aoi and aos)
    7470
    7571### Add csv header list to all files in point_filenames
    … …  
    8177# Used in build_elevation.py
    8278# Format for points easting,northing (no header)
    83 bounding_polygon_filename = 'poly_all_MH.csv'
    84 bounding_polygon_maxarea = 100000
     79bounding_polygon_filename = 'bounding_polygon.csv'
     80bounding_polygon_maxarea = 1000000
    8581
    8682# INTERIOR REGIONS -  for designing the mesh
    8783# Used in run_model.py
    8884# Format for points easting,northing (no header)                   
    89 interior_regions_data = [['sw_MH.csv', 3000],
    90                          ['aos_MH.csv', 800],
    91                          ['aoi_small_MH.csv', 500],
    92                          ['island1_MH.csv', 3000],
    93                          ['island2_MH.csv', 3000],
    94                          ['island3_MH.csv', 3000],
    95                          ['island4_MH.csv', 3000],
    96                          ['island5_MH.csv', 3000]]#,
    97 ##                         ['island6_MH.csv', 3000]]
    98 
    99 PriorityArea_filename = None
     85
     86interior_regions_data = [['aoi.csv', 500]]
     87PriorityArea_filename = 'aos.csv'
    10088   
    10189# LAND - used to set the initial stage/water to be offcoast only
    102 # Used in run_model.py.  Format for points easting,northing,id,value
    103 land_initial_conditions_filename = 'initial_condition_comb.csv'
     90# Used in run_model.py.  Format for points easting,northing (no header)
     91land_initial_conditions_filename = 'initial_conditions_simple.csv'
    10492
    10593# GAUGES - for creating timeseries at a specific point
    10694# Used in get_timeseries.py. 
    10795# Format easting,northing,name,elevation (with header)
    108 gauges_filename = ''
     96gauges_filename = '' #'TideGaugesPoints.csv' #'tsunamipointsMGA.csv'
     97
    10998# BUILDINGS EXPOSURE - for identifying inundated houses
    11099# Used in run_building_inundation.py
    111100# Format latitude,longitude etc (geographic)
    112 building_exposure_filename = 'busselton_res_clip.csv' # from NEXIS
     101building_exposure_filename = '' # from NEXIS
    113102
    114103# AREA OF IMAGES - Extent of each image to find out highest runup
    115104# Header - easting,northing,id,value
    116105# Used in get_runup.py
    117 images_filename = ''
     106images_filename = 'images.csv'
    118107
    119108# BOUNDING POLYGON - used in build_boundary.py and run_model.py respectively
    … …  
    145134
    146135# ASCII export grid for Exmouth
    147 xminExmouth = 511616
    148 xmaxExmouth = 516764
    149 yminExmouth = 7568731
    150 ymaxExmouth = 7577540
    151 
     136xminExmouth = 511000
     137xmaxExmouth = 516000
     138yminExmouth = 7567700
     139ymaxExmouth = 7576500
    152140
    153141################################################################################
    … …  
    201189event_folder = join(boundaries_folder, str(event_number))
    202190
    203 ### Add csv header list to all files in point_filenames
    204 ##headerlist = ['x', 'y', 'elevation']
    205 ##for f in point_filenames:
    206 ##    add_csv_header(join(topographies_folder, f), headerlist)
    207 
    208191# MUX data files
    209192# Directory containing the MUX data files to be used with EventSelection.

  • anuga_work/production/exmouth_2009/run_model.py

    r7209 r7385  
    137137Br = Reflective_boundary(domain)
    138138Bt = Transmissive_stage_zero_momentum_boundary(domain)
    139 Bd = Dirichlet_boundary([project.tide, 0, 0])
    140139Bf = Field_boundary(project.event_sts+'.sts',
    141140                    domain, mean_stage=project.tide,

  • anuga_work/production/geraldton/geraldton_rerun/Arc_asc2raster_GDA94z50.py

    r6682 r7385  
    3232time_dir6 = '20090327_160043_run_final_0_27255_Bt_kvanputt'
    3333
    34 time_dirs = [time_dir1, time_dir2, time_dir3]
     34time_dirs = [time_dir1] #, time_dir2, time_dir3]
    3535
    3636for time_dir in time_dirs:
    … …  
    5656
    5757    generate_filename = []
    58     input_ascii = glob.glob(folder + '*speed_max.asc')
     58    input_ascii = glob.glob(folder + '*elevation.asc')
    5959
    6060    for infile in input_ascii:

  • anuga_work/production/hobart_2009/export_results_max.py

    r7355 r7385  
    2121
    2222directory = project.output_folder
    23 time_dir1 = '20090714_162136_run_final_0_58260_None_kvanputt'
     23time_dir1 = '20090504_124158_run_final_0_58068_small_kvanputt'
    2424time_dir2 = '20090714_162041_run_final_0.8_58292_None_kvanputt'
    2525time_dir3 = '20090714_161944_run_final_0.8_58280_None_kvanputt'
    2626time_dir4 = '20090714_161740_run_final_0.8_58260_None_kvanputt'
    2727
    28 time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4] 
     28time_dirs = [time_dir1] 
    2929
    3030# sww filename extensions ie. hobart_time_17640_0.sww, input into list 17640
    3131# make sure numbers are in sequential order
    32 times = [16980, 33960, 50940]
     32times = [22680, 45360, 68040]
    3333           
    3434#Modify the cellsize value to set the size of the raster you require
    3535#Take into account mesh size when aplying this paramater
    36 cellsize = 20 #250
     36cellsize = 250 #250
    3737
    3838#Now set the timestep at which you want the raster generated.
    … …  
    120120            print 'which var=%s' % which_var
    121121            #print 'asc_name ', str(asc_name)
    122             if times == []:     
    123                 print 'Only one sww file, no Max file required'
     122            if len(names) >1:
     123                maxasc_outname = join(directory, time_dir, 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)
    124129            else:
    125                 print 'Multiple sww files'
    126                 maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc'
    127                 print 'maxasc_outname = ', maxasc_outname
    128                 MaxAsc(maxasc_outname, asc_name)
     130                print 'Only one sww file, no Max file required'
     131

  • anuga_work/production/onslow_2009_june/Arc_asc2raster_GDA94z50.py

    r7261 r7385  
    2525output_dir="anuga\\outputs\\"
    2626
    27 time_dir1 = '20090617_110255_run_final_0.0_27283_mhingee'
    28 ##time_dir2 = '20090416_072538_run_final_0_27338_1307_Tb__kvanputt'
    29 ##time_dir3 = '20090417_021634_run_final_0_70337_1307_Tb__kvanputt'
     27time_dir1 = '20090817_171542_run_final_0.0_27283_kvanputt'
    3028
    31 ##time_dir2 = '20090422_102457_run_final_0_70844_1307_Tb_internal_mhingee'
    32 
    33 time_dirs = [time_dir1]#, time_dir2]#, time_dir3]
     29time_dirs = [time_dir1]
    3430
    3531for time_dir in time_dirs:
    … …  
    4036##    ocean = scenario_dir + "map_work\\\port_hedland.gdb\\outlines\\initial_conditions_ocean"
    4137
    42 ##    print 'Process: Create File GDB'
    43 ##    gp.CreateFileGDB_management(folder, "raster")
     38    print 'Process: Create File GDB'
     39    gp.CreateFileGDB_management(folder, "raster")
    4440
    4541    gp.Workspace = raster_gbd
    … …  
    4844   
    4945    #replication dictionary
    50     replicate = (('onslow', 'o'),('_', ''),('max','_M'),
    51                  ('CBD', 'CDB'),('All',''),
     46    replicate = (('onslow', ''),('_', ''),('max','_M'),
     47                 ('CBD', 'CDB'),
    5248                 ('depth','_depth'),('speed', '_speed'),
    53                  ('elevation', '_ele_'), ('stage','_stage'))
     49                 ('elevation', '_ele'), ('stage','_stage'))
    5450
    5551    generate_filename = []
    56     input_ascii = glob.glob(folder + '*w_O*.asc')
     52    input_ascii = glob.glob(folder + '*.asc')
    5753
    5854    for infile in input_ascii:
    … …  
    7268
    7369        print 'Process: Define Projection'
    74         #GDA_1994_MGA_Zone_54
    7570        gp.DefineProjection_management(output_DEM, "PROJCS['GDA_1994_MGA_Zone_50',GEOGCS['GCS_GDA_1994',DATUM['D_GDA_1994',SPHEROID['GRS_1980',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]]"
    7671                                                   ",PROJECTION['Transverse_Mercator'],PARAMETER['False_Easting',500000.0],PARAMETER['False_Northing',10000000.0],PARAMETER['Central_Meridian',117.0],PARAMETER['Scale_Factor',0.9996]"

  • anuga_work/production/onslow_2009_june/project.py

    r7261 r7385  
    2626# One or all can be changed each time the run_model script is executed
    2727tide = 0.0              # 3.6 difference between MSL and HAT in metres
    28 zone = 50               # specify UTM zone of model
     28zone= 50               # specify UTM zone of model
    2929central_meridian = None
    3030event_number = 27283    # 27255,27283 the event number or the mux file name
    … …  
    3838                        #    basic - coarse mesh
    3939                        #    final - fine mesh, slowest
     40
     41# MUX input filename.
     42# If a meta-file from EventSelection is used, set 'multi-mux' to True.
     43# If a single MUX stem filename (*.grd) is used, set 'multi-mux' to False.
     44##mux_input_filename = event_number # to be found in event_folder
     45                                    # (ie boundaries/event_number/)
     46##multi_mux = False
     47mux_input_filename = 'event.list'
     48multi_mux = True
    4049
    4150#-------------------------------------------------------------------------------
    … …  
    6473# Format for point is x,y,elevation (with header)
    6574point_filenames = ['onslow_offshore_points.txt',
    66                    'revised_coastline.txt',
    6775                   'beach_survey.txt',
    6876                   'ONS16FINAL.txt',
    … …  
    7583##    add_csv_header(join(topographies_folder, f), headerlist)
    7684
    77 # BOUNDING POLYGON - for data clipping and estimate of triangles in mesh
     85# SIMPLE BOUNDING POLYGON - for data clipping and estimate of triangles in mesh
    7886# Used in build_elevation.py
    7987# Format for points easting,northing (no header)
    … …  
    92100                         ['island4_MH.csv', 3000]]
    93101
    94 PriorityArea_filename = None
     102#Format for points easting,northing,id,value
     103interior_regions_combined_filename = None
    95104   
    96105# LAND - used to set the initial stage/water to be offcoast only
    97 # Used in run_model.py.  Format for points easting,northing,id,value
    98 land_initial_conditions_filename = 'initial_conditions_mainland.csv'
     106# Used in run_model.py.
     107# If a list format for points easting,northing (no header)
     108
     109##land_initial_conditions_filename = []
     110
     111#Format for points easting,northing,id,value
     112land_initial_conditions_filename = 'initial_conditions_simple.csv'
     113
     114# MODEL BOUNDING POLYGON - used in build_boundary.py and run_model.py respectively
     115# NOTE: when files are put together the points must be in sequence
     116# For ease go clockwise!
     117# Check the run_model.py for boundary_tags
     118
     119# Thinned ordering file from Hazard Map (geographic)
     120# Format is index,latitude,longitude (with header)
     121urs_order_filename = 'urs_order.csv'
     122
     123# Landward bounding points
     124# Format easting,northing (no header)
     125landward_boundary_filename = 'landward_boundary.csv'
     126
     127#-------------------------------------------------------------------------------
     128# For output sww file
     129#-------------------------------------------------------------------------------
    99130
    100131# GAUGES - for creating timeseries at a specific point
    … …  
    105136# Used in run_building_inundation.py
    106137# Format latitude,longitude etc (geographic)
    107 building_exposure_filename = 'busselton_res_clip.csv' # from NEXIS
     138building_exposure_filename = '' # from NEXIS
    108139
    109140# AREA OF IMAGES - Extent of each image to find out highest runup
    … …  
    112143images_filename = ''
    113144
    114 # BOUNDING POLYGON - used in build_boundary.py and run_model.py respectively
    115 # NOTE: when files are put together the points must be in sequence
    116 # For ease go clockwise!
    117 # Check the run_model.py for boundary_tags
    118 
    119 # Thinned ordering file from Hazard Map (geographic)
    120 # Format is index,latitude,longitude (with header)
    121 urs_order_filename = 'urs_order.csv'
    122 
    123 # Landward bounding points
    124 # Format easting,northing (no header)
    125 landward_boundary_filename = 'landward_boundary.csv'
    126 
    127 # MUX input filename.
    128 # If a meta-file from EventSelection is used, set 'multi-mux' to True.
    129 # If a single MUX stem filename (*.grd) is used, set 'multi-mux' to False.
    130 ##mux_input_filename = event_number # to be found in event_folder
    131                                     # (ie boundaries/event_number/)
    132 ##multi_mux = False
    133 mux_input_filename = 'event.list'
    134 multi_mux = True
    135 
    136 #-------------------------------------------------------------------------------
     145
    137146# Clipping regions for export to asc and regions for clipping data
    138147# Final inundation maps should only be created in regions of the finest mesh
    139 #-------------------------------------------------------------------------------
    140148
    141149# ASCII export grid for Onslow
    142 xminOnslow = 301199
    143 xmaxOnslow = 311353
    144 yminOnslow = 7600921
    145 ymaxOnslow = 7609833
     150xminOnslow = 301200
     151xmaxOnslow = 311500
     152yminOnslow = 7601000
     153ymaxOnslow = 7609200
    146154
    147155
    … …  
    254262
    255263#Multiple polygons in one CSV file to make internal polygons
    256 if not PriorityArea_filename == None:
    257     PriorityAreas = join(polygons_folder, PriorityArea_filename)
    258 
    259 
     264if not interior_regions_combined_filename == None:
     265    interior_regions_combined = join(polygons_folder, interior_regions_combined_filename)
     266
     267

  • anuga_work/production/onslow_2009_june/setup_model.py

    r7209 r7385  
    169169project.interior_regions = []
    170170
    171 if project.PriorityArea_filename is not None:
     171if project.interior_regions_combined_filename is not None:
    172172    polygons, maxareas = csv2polygons(project.PriorityAreas)
    173173    print 'Creating %d internal polygons' % len(polygons)

  • anuga_work/production/perth/perth_rerun/Arc_asc2raster_GDA94z50.py

    r7113 r7385  
    3232time_dir6 = '20090327_155803_run_final_0_68693_Bt_kvanputt'
    3333
    34 time_dirs = [time_dir4, time_dir5, time_dir6] #time_dir1, time_dir2, time_dir3,
     34time_dirs = [time_dir1] #4, time_dir5, time_dir6] #time_dir1, time_dir2, time_dir3,
    3535 
    3636for time_dir in time_dirs:
    … …  
    5757
    5858    generate_filename = []
    59     input_ascii = glob.glob(folder + '*speed_max.asc')
     59    input_ascii = glob.glob(folder + '*elevation.asc')
    6060    print time_dir
    6161

  • anuga_work/production/pt_hedland_2009/Arc_asc2raster_GDA94z50.py

    r7261 r7385  
    2525output_dir="anuga\\outputs\\"
    2626
    27 time_dir1 = '20090615_151314_run_final_0.0_27283_mhingee'
    28 ##time_dir2 = '20090416_072538_run_final_0_27338_1307_Tb__kvanputt'
    29 ##time_dir3 = '20090417_021634_run_final_0_70337_1307_Tb__kvanputt'
     27time_dir1 = '20090818_133107_run_final_0.0_27283_kvanputt'
    3028
    31 ##time_dir2 = '20090422_102457_run_final_0_70844_1307_Tb_internal_mhingee'
    32 
    33 time_dirs = [time_dir1]#, time_dir2]#, time_dir3]
     29time_dirs = [time_dir1]
    3430
    3531for time_dir in time_dirs:
    … …  
    4036##    ocean = scenario_dir + "map_work\\\port_hedland.gdb\\outlines\\initial_conditions_ocean"
    4137
    42 ##    print 'Process: Create File GDB'
    43 ##    gp.CreateFileGDB_management(folder, "raster")
     38    print 'Process: Create File GDB'
     39    gp.CreateFileGDB_management(folder, "raster")
    4440
    4541    gp.Workspace = raster_gbd
    … …  
    4945    #replication dictionary
    5046    replicate = (('pt_hedland', ''),('_', ''),('max','_M'),
    51                  ('CBD', 'CDB'),('All',''),
     47                 ('CBD', 'CDB'), ('PtHedland', 'PH'),
    5248                 ('depth','_depth'),('speed', '_speed'),
    53                  ('elevation', '_ele_'), ('stage','_stage'))
     49                 ('elevation', '_ele'), ('stage','_stage'))
    5450
    5551    generate_filename = []
    56     input_ascii = glob.glob(folder + '*Town*.asc')
     52    input_ascii = glob.glob(folder + '*.asc')
    5753
    5854    for infile in input_ascii:
    … …  
    7268
    7369        print 'Process: Define Projection'
    74         #GDA_1994_MGA_Zone_54
    7570        gp.DefineProjection_management(output_DEM, "PROJCS['GDA_1994_MGA_Zone_50',GEOGCS['GCS_GDA_1994',DATUM['D_GDA_1994',SPHEROID['GRS_1980',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]]"
    7671                                                   ",PROJECTION['Transverse_Mercator'],PARAMETER['False_Easting',500000.0],PARAMETER['False_Northing',10000000.0],PARAMETER['Central_Meridian',117.0],PARAMETER['Scale_Factor',0.9996]"

  • anuga_work/production/pt_hedland_2009/project.py

    r7261 r7385  
    3232friction=0.01           # manning's friction coefficient
    3333starttime=0             # start time for simulation
    34 finaltime=80000         # final time for simulation
     34finaltime=1000         # final time for simulation
    3535
    3636setup = 'final'         # This can be one of three values
    … …  
    8383# LAND - used to set the initial stage/water to be offcoast only
    8484# Used in run_model.py.  Format for points easting,northing,id,value
    85 land_initial_conditions_filename = 'initial_conditions_comb.csv'
     85land_initial_conditions_filename = 'initial_conditions_simple.csv'
    8686
    8787# GAUGES - for creating timeseries at a specific point
    … …  
    9292# Used in run_building_inundation.py
    9393# Format latitude,longitude etc (geographic)
    94 building_exposure_filename = 'busselton_res_clip.csv' # from NEXIS
     94building_exposure_filename = '' # from NEXIS
    9595
    9696# AREA OF IMAGES - Extent of each image to find out highest runup
    … …  
    127127
    128128# ASCII export grid for Pt Hedland
    129 xminTown = 659643
    130 xmaxTown = 676236
    131 yminTown = 7744216
    132 ymaxTown = 7757347
     129xminPtHedland = 659400
     130xmaxPtHedland = 676600
     131yminPtHedland = 7744000
     132ymaxPtHedland = 7757300
    133133
    134134
Note: See TracChangeset for help on using the changeset viewer.