Changeset 6982


Ignore:
Timestamp:
May 6, 2009, 12:06:04 PM (11 years ago)
Author:
kristy
Message:
 
Location:
anuga_work/production/australia_ph2
Files:
1 deleted
10 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/australia_ph2/carnarvon/build_elevation.py

    r6412 r6982  
    5555        convert_dem_from_ascii2netcdf(absolute_filename,
    5656                                      basename_out=absolute_filename,
    57                                       use_cache=True,
     57                                      use_cache=False,
    5858                                      verbose=True)
    59         dem2pts(absolute_filename, use_cache=True, verbose=True)
     59        dem2pts(absolute_filename, use_cache=False, verbose=True)
    6060
    6161        geospatial_data[filename] = Geospatial_data(file_name=absolute_filename+'.pts',
     
    8686##G_all = G_clip + geospatial_data[project.point_filenames[4]]
    8787#G_clipped = G_all.clip(project.poly_all)
    88 G_clip = G.clip(project.bounding_polygon)
     88#G_clip = G.clip(project.bounding_polygon)
    8989
    9090
    9191print 'Export combined DEM file'
    92 G_clip.export_points_file(project.combined_elevation + '.pts')
     92G.export_points_file(project.combined_elevation + '.pts')
    9393print 'Do txt version too'
    9494# Use for comparision in ARC
    95 G_clip.export_points_file(project.combined_elevation + '.txt')
     95G.export_points_file(project.combined_elevation + '.txt')
    9696
  • anuga_work/production/australia_ph2/carnarvon/project.py

    r6971 r6982  
    3232    event_number = int(sys.argv[1])
    3333else:   
    34     event_number = 27285    # the event number or the mux file name
     34    event_number = 68754    # the event number or the mux file name
    3535
    3636event_number_list = [27285, 68754, 70881] # To piggy back multiple events
     
    7878# Used in build_elevation.py
    7979# Format for ascii grids, as produced in ArcGIS + a projection file
    80 ascii_grid_filenames = ['grid250m_simple'] # 250m grid 2005
     80ascii_grid_filenames = ['grid_250m_simple'] # 250m grid 2005
    8181
    8282# Format for points is x,y,elevation (with header)
     
    132132# Landward bounding points
    133133# Format easting,northing (no header)
    134 landward_boundary_filename = 'landward_boundary_simple.csv'
     134landward_boundary_filename = 'landward_boundary_simple2.csv'
    135135
    136136# MUX input filename.
     
    145145# Whatever is specified here will be relative to INUNDATION_HOME/.cache
    146146# If nothing is specified, local cache will be used.
    147 cachedir = '.python_cache_phII'
     147#cachedir = '.python_cache_phII'
    148148
    149149
     
    173173# Output filename for elevation
    174174# this is a combination of all the data generated in build_elevation.py
    175 combined_elevation_basename = scenario_name + '_combined_elevation'
     175combined_elevation_basename = scenario_name + '_combined_elevation_simple'
    176176
    177177#-------------------------------------------------------------------------------
     
    191191# Create absolute pathname for cache directory
    192192# and change caching to use it
    193 if 'cachedir' in dir():
    194     cachedir = join(os.getenv(ENV_INUNDATIONHOME), '.cache', cachedir)
    195     from anuga.caching import caching
    196     caching.set_option('cachedir', cachedir)
     193##if 'cachedir' in dir():
     194##    cachedir = join(os.getenv(ENV_INUNDATIONHOME), '.cache', cachedir)
     195##    from anuga.caching import caching
     196##    caching.set_option('cachedir', cachedir)
    197197   
    198198# check various directories/files that must exist
     
    219219
    220220# The absolute pathname of the all elevation, generated in build_elevation.py
    221 combined_elevation = join(topographies_folder, combined_elevation_basename)+'_simple'
     221combined_elevation = join(topographies_folder, combined_elevation_basename)
    222222
    223223# The absolute pathname of the mesh, generated in run_model.py
  • anuga_work/production/australia_ph2/hobart/build_elevation.py

    r6407 r6982  
    8686##G_all = G_clip + geospatial_data[project.point_filenames[4]]
    8787#G_clipped = G_all.clip(project.poly_all)
    88 G_clip = G.clip(project.bounding_polygon)
     88#G_clip = G.clip(project.bounding_polygon)
    8989
    9090
    9191print 'Export combined DEM file'
    92 G_clip.export_points_file(project.combined_elevation + '.pts')
     92G.export_points_file(project.combined_elevation + '.pts')
    9393print 'Do txt version too'
    9494# Use for comparision in ARC
    95 G_clip.export_points_file(project.combined_elevation + '.txt')
     95G.export_points_file(project.combined_elevation + '.txt')
    9696
  • anuga_work/production/australia_ph2/hobart/export_results_max.py

    r6407 r6982  
    2222directory = project.output_folder
    2323
    24 time_dir1 = '20090130_165532_run_final_0.6_27255_extend_dt_kvanputt'
    25 ##time_dir2 = ''
    26 ##time_dir3 = ''
    27 ##time_dir4 = ''
    28 ##time_dir5 = ''
    29 ##time_dir6 = ''
    30 ##
    31 time_dirs = [time_dir1] #5, time_dir4]# , time_dir4, time_dir5, time_dir6]
     24time_dir1 = '20090424_121952_run_final_0_51204_2462_Tb__kvanputt'
     25time_dir2 = '20090424_123250_run_final_0_51204_2462_Tb_internal_kvanputt'
     26time_dir3 = '20090424_122951_run_final_0_58212_2462_Tb__kvanputt'
     27time_dir4 = '20090424_123104_run_final_0_64425_2462_Tb__kvanputt'
     28
     29time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4]
    3230 
    3331
     
    3533##cellsize = 5
    3634
    37 #timestep = None    # None means no timestep!
    38 timestep = 0
     35timestep = None    # None means no timestep!
     36#timestep = 0
    3937
    4038######
     
    4240######
    4341
    44 #area = ['Bunbury', 'Busselton']  # strings must match keys in var_equations below
    4542area = ['All']      # 'All' means no special areas - the whole thing
    4643
     
    5653
    5754# one or more key strings from var_equations above
    58 var = ['elevation']
     55var = ['stage']
    5956
    6057######
     
    8582        for time_dir in time_dirs:
    8683
    87             name1 = join(directory, time_dir, project.scenario_name)
    88             #name2 = join(directory, time_dir, project.scenario_name)+'_time_39600_0'
    89             #name3 = join(directory, time_dir, project.scenario_name)+'_time_79200_0'
     84            name = join(directory, time_dir, project.scenario_name)
     85                           
     86            outname = name + '_' + which_area + '_' + which_var
     87            quantityname = var_equations[which_var]
    9088
    91             names = [name1] #, name2]#, name3]
    92      
    93            # asc_name = []   
    94                
    95             for name in names:
    96                          
    97                 outname = name + '_' + which_area + '_' + which_var
    98                 quantityname = var_equations[which_var]
     89            print 'start sww2dem: time_dir=%s' % time_dir
     90           
     91            sww2dem(name, basename_out = outname,
     92                        quantity = quantityname,
     93                        timestep = timestep,
     94                        cellsize = cellsize,     
     95                        easting_min = easting_min,
     96                        easting_max = easting_max,
     97                        northing_min = northing_min,
     98                        northing_max = northing_max,       
     99                        reduction = max,
     100                        verbose = True,
     101                        format = 'asc')
    99102
    100                 print 'start sww2dem: time_dir=%s' % time_dir
    101                
    102                 sww2dem(name, basename_out = outname,
    103                             quantity = quantityname,
    104                             timestep = timestep,
    105                             cellsize = cellsize,     
    106                             easting_min = easting_min,
    107                             easting_max = easting_max,
    108                             northing_min = northing_min,
    109                             northing_max = northing_max,       
    110                             reduction = max,
    111                             verbose = True,
    112                             format = 'asc')
    113 
    114 ##                asc_name.append(outname + '.asc')
    115 ##   
    116 ##            maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc'
    117 ##
    118 ##            print 'max asc outname ', maxasc_outname
    119 ##            print 'asc_name ', str(asc_name)
    120 ##           
    121 ##            MaxAsc(maxasc_outname, asc_name)
  • anuga_work/production/australia_ph2/hobart/project.py

    r6772 r6982  
    3232    event_number = int(sys.argv[1])
    3333else:   
    34     event_number = 58082    # the event number or the mux file name
    35 
    36 event_number_list = [51268, 58082, 64309] # To piggy back multiple events
     34    event_number = 51204    # the event number or the mux file name
     35
     36event_number_list = [51204, 58212, 64425] # To piggy back multiple events
    3737
    3838tide = 0                # difference between MSL and HAT
     
    4040friction=0.01           # manning's friction coefficient
    4141starttime=0             # start time for simulation
    42 finaltime=60000 # final time for simulation
     42finaltime=60000         # final time for simulation
    4343setup = 'final'         # This can be one of three values
    4444                        #    trial - coarsest mesh, fast
     
    4747
    4848# index is only used when wave = Tb
    49 index = 2528            # index from the PTHA - Y2000 0.257m
     49index = 2462            # index from the PTHA - Y2000 0.257m
    5050wave = 'Tb'             # Bf (sts wave) Tb (index wave)
    5151
    5252
    53 internal_polygon = False
     53internal_polygon = True
    5454
    5555#-------------------------------------------------------------------------------
     
    7878# Used in build_elevation.py
    7979# Format for ascii grids, as produced in ArcGIS + a projection file
    80 ascii_grid_filenames = [] # 250m grid 2005
     80ascii_grid_filenames = ['grid_250m_simple'] # 250m grid 2005
    8181
    8282# Format for points is x,y,elevation (with header)
    83 point_filenames = ['grid_250m.txt']   # 250m grid 2005
     83point_filenames = []   # 250m grid 2005
    8484
    8585### Add csv header list to all files in point_filenames
     
    9191# Used in build_elevation.py
    9292# Format for points: easting,northing (no header)
    93 bounding_polygon_filename = 'bounding_polygon.csv'
     93bounding_polygon_filename = 'bounding_polygon_simple.csv'
    9494bounding_polygon_maxarea = 125000
    9595
     
    128128# Thinned ordering file from Hazard Map (geographic)
    129129# Format is index,latitude,longitude (with header)
    130 urs_order_filename = 'urs_order_simple.csv'
     130urs_order_filename = 'urs_order_simple2.csv'
    131131
    132132# Landward bounding points
     
    173173# Output filename for elevation
    174174# this is a combination of all the data generated in build_elevation.py
    175 combined_elevation_basename = scenario_name + '_combined_elevation'
     175combined_elevation_basename = scenario_name + '_combined_elevation_simple'
    176176
    177177#-------------------------------------------------------------------------------
  • anuga_work/production/australia_ph2/melbourne_east/project.py

    r6856 r6982  
    2929    event_number = int(sys.argv[1])
    3030else:   
    31     event_number = 58348    # the event number or the mux file name
     31    event_number = 31853    # the event number or the mux file name
    3232
    3333event_side = 'east' #'west'
    34 
    35 if event_side == 'west':
    36     index = 1979 # index from the PTHA - Y2000 0.257m
    37     event_number_list = [58348, 64288, 68793] # To piggy back multiple events
    38 
    39 elif event_side == 'east':
    40     index = 2649
    41     event_number_list = [31853, 51449, 58115] # To piggy back multiple events
     34index = 2649
     35event_number_list = [31853, 51449, 58115] # To piggy back multiple events
    4236
    4337central_meridian = 146.0 # Central meridian for projection (optional)
  • anuga_work/production/australia_ph2/sydney/export_results.py

    r6926 r6982  
    2929##cellsize = 5
    3030
    31 #timestep = None    # None means no timestep!
    32 timestep = 0
     31timestep = None    # None means no timestep!
     32#timestep = 0
    3333
    3434######
     
    5050
    5151# one or more key strings from var_equations above
    52 var = ['elevation']
     52var = ['stage']
    5353
    5454######
     
    7979        for time_dir in time_dirs:
    8080
    81             name1 = join(directory, time_dir, project.scenario_name)
    82 ##            name2 = join(directory, time_dir, project.scenario_name)+'_time_21060_0'
    83 ##            name3 = join(directory, time_dir, project.scenario_name)+'_time_42120_0'
    84 ##            name4 = join(directory, time_dir, project.scenario_name)+'_time_63180_0'
    85 ##            name5 = join(directory, time_dir, project.scenario_name)+'_time_53280_0'
    86 ##            name6 = join(directory, time_dir, project.scenario_name)+'_time_66600_0'
    87 ##            name7 = join(directory, time_dir, project.scenario_name)+'_time_79920_0'
     81            name = join(directory, time_dir, project.scenario_name)
    8882           
    89             names = [name1]#, name2, name3,name4]#,name5,name6,name7]
    90      
    91             asc_name = []   
    92                
    93             for name in names:
    94                          
    95                 outname = name + '_' + which_area + '_' + which_var
    96                 quantityname = var_equations[which_var]
     83            outname = name + '_' + which_area + '_' + which_var
     84            quantityname = var_equations[which_var]
     85           
     86            print 'start sww2dem: time_dir=%s' % time_dir
     87           
     88            sww2dem(name, basename_out = outname,
     89                        quantity = quantityname,
     90                        timestep = timestep,
     91                        cellsize = cellsize,     
     92                        easting_min = easting_min,
     93                        easting_max = easting_max,
     94                        northing_min = northing_min,
     95                        northing_max = northing_max,       
     96                        reduction = max,
     97                        verbose = True,
     98                        format = 'asc')
    9799
    98                 print 'start sww2dem: time_dir=%s' % time_dir
    99                
    100                 sww2dem(name, basename_out = outname,
    101                             quantity = quantityname,
    102                             timestep = timestep,
    103                             cellsize = cellsize,     
    104                             easting_min = easting_min,
    105                             easting_max = easting_max,
    106                             northing_min = northing_min,
    107                             northing_max = northing_max,       
    108                             reduction = max,
    109                             verbose = True,
    110                             format = 'asc')
    111 
    112 ##                asc_name.append(outname + '.asc')
    113 ##   
    114 ##            maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc'
    115 ##
    116 ##            print 'max asc outname ', maxasc_outname
    117 ##            print 'asc_name ', str(asc_name)
    118 ##           
    119 ##            MaxAsc(maxasc_outname, asc_name)
  • anuga_work/production/australia_ph2/sydney/export_results_max.py

    r6926 r6982  
    2222directory = project.output_folder
    2323
    24 time_dir1 = '20090421_150708_run_final_0_58152_2938_Tb_internal_mhingee'
    25 ##time_dir2 = '20090311_152929_run_final_0_70385_1203_Tb_jgriffin'
    26 ##time_dir3 = '20090311_153009_run_final_0_70454_1203_Bf_jgriffin'
    27 ##time_dir4 = '20090311_153054_run_final_0_70339_1203_Bf_jgriffin'
     24time_dir1 = '20090428_131539_run_final_30_0_58152_2938_Tb_internal_kvanputt'
     25time_dir2 = '20090311_152929_run_final_0_70385_1203_Tb_jgriffin'
     26time_dir3 = '20090311_153009_run_final_0_70454_1203_Bf_jgriffin'
     27time_dir4 = '20090311_153054_run_final_0_70339_1203_Bf_jgriffin'
    2828time_dirs = [time_dir1]#,time_dir2,time_dir3,time_dir4] 
    2929
     
    8686
    8787            name1 = join(directory, time_dir, project.scenario_name)
    88 ##            name2 = join(directory, time_dir, project.scenario_name)+'_time_66360_0'
    89 ##            name3 = join(directory, time_dir, project.scenario_name)+'_time_26640_0'
    90 ##            name4 = join(directory, time_dir, project.scenario_name)+'_time_39960_0'
    91 ##            name5 = join(directory, time_dir, project.scenario_name)+'_time_53280_0'
    92 ##            name6 = join(directory, time_dir, project.scenario_name)+'_time_66600_0'
    93 ##            name7 = join(directory, time_dir, project.scenario_name)+'_time_79920_0'
     88            name2 = join(directory, time_dir, project.scenario_name)+'_time_32460_0'
     89            name3 = join(directory, time_dir, project.scenario_name)+'_time_26640_0'
     90            name4 = join(directory, time_dir, project.scenario_name)+'_time_39960_0'
     91            name5 = join(directory, time_dir, project.scenario_name)+'_time_53280_0'
     92            name6 = join(directory, time_dir, project.scenario_name)+'_time_66600_0'
     93            name7 = join(directory, time_dir, project.scenario_name)+'_time_79920_0'
    9494           
    95             names = [name1]#, name2]#, name3,name4,name5,name6,name7]
     95            names = [name1, name2]#, name3,name4,name5,name6,name7]
    9696     
    9797            asc_name = []   
     
    116116                            format = 'asc')
    117117
    118 ##                asc_name.append(outname + '.asc')
    119 ##   
    120 ##            maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc'
    121 ##
    122 ##            print 'max asc outname ', maxasc_outname
    123 ##            print 'asc_name ', str(asc_name)
    124 ##           
    125 ##            MaxAsc(maxasc_outname, asc_name)
     118                asc_name.append(outname + '.asc')
     119   
     120            maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc'
     121
     122            print 'max asc outname ', maxasc_outname
     123            print 'asc_name ', str(asc_name)
     124           
     125
  • anuga_work/production/australia_ph2/sydney/project.py

    r6959 r6982  
    4545starttime=0             # start time for simulation
    4646finaltime=60000   #60000 # final time for simulation
    47 setup = 'final'         # This can be one of three values
     47setup = 'final_30'         # This can be one of three values
    4848                        #    trial - coarsest mesh, fast
    4949                        #    basic - coarse mesh
  • anuga_work/production/australia_ph2/sydney/setup_model.py

    r6461 r6982  
    109109    project.time_thinning = 12
    110110    project.yieldstep = 120
     111elif project.setup == 'final_120':
     112    print 'final_120'
     113    project.scale_factor = 1
     114    project.time_thinning = 4
     115    project.yieldstep = 120
     116elif project.setup == 'final_30':
     117    print 'final_120'
     118    project.scale_factor = 1
     119    project.time_thinning = 4
     120    project.yieldstep = 30
    111121elif project.setup == 'final':
    112122    print 'final'
Note: See TracChangeset for help on using the changeset viewer.