Changeset 4310


Ignore:
Timestamp:
Mar 21, 2007, 6:25:37 PM (18 years ago)
Author:
sexton
Message:

update damage scripts for FESA areas

Location:
anuga_work/production
Files:
11 edited

Legend:

Unmodified
Added
Removed

  • anuga_work/production/broome_2006/run_building_inundation.py

    r4228 r4310  
    3333
    3434add_depth_and_momentum2csv(swwfile, buildings_filename, buildings_filename_out)
     35inundation_damage(swwfile, buildings_filename, buildings_filename_out)
    3536
    3637

  • anuga_work/production/busselton_2006/run_building_inundation.py

    r4228 r4310  
    3333
    3434add_depth_and_momentum2csv(swwfile, buildings_filename, buildings_filename_out)
    35 
     35inundation_damage(swwfile, buildings_filename, buildings_filename_out)
    3636
    3737print '\n Augmented building file written to %s \n' \

  • anuga_work/production/dampier_2006/run_building_inundation.py

    r4228 r4310  
    3333
    3434add_depth_and_momentum2csv(swwfile, buildings_filename, buildings_filename_out)
    35 
     35inundation_damage(swwfile, buildings_filename, buildings_filename_out)
    3636
    3737print '\n Augmented building file written to %s \n' \

  • anuga_work/production/exmouth_2006/run_building_inundation.py

    r4228 r4310  
    3333
    3434add_depth_and_momentum2csv(swwfile, buildings_filename, buildings_filename_out)
    35 
     35inundation_damage(swwfile, buildings_filename, buildings_filename_out)
    3636
    3737print '\n Augmented building file written to %s \n' \

  • anuga_work/production/gold_coast_2007/export_results.py

    r4294 r4310  
    55from os import sep
    66
    7 time_dir = '20070228_071238'
     7#time_dir = '20070228_071238' # 6.0 m
     8#time_dir = '20070305_051718' # 2.5m
     9#time_dir = '20070305_051728' # 3.5m
     10#time_dir = '20070305_051817' # 4.5m
     11time_dir = '20070314_231144' # 4.5 tide
     12#time_dir = '20070314_231156' # 5m tide
     13
    814
    915directory = project.outputdir
    … …  
    2430
    2531if which_var == 2:  # Depth
    26     outname = name + '_depth-eq'
     32    outname = name + '_depth'
    2733    quantityname = 'stage-elevation' 
    2834
    … …  
    5763                quantity = quantityname,
    5864                #timestep = 0,
    59                 cellsize = 50,     
     65                cellsize = 20,     
    6066                easting_min = project.eastingmin,
    6167                easting_max = project.eastingmax,

  • anuga_work/production/gold_coast_2007/project.py

    r4294 r4310  
    2828gc250 = 'dem_250'
    2929
     30# better DEM
     31gc5 = 'dem_5m_ascii'
     32offshore = 'xy_AHD'
     33coastline = 'coastline'
     34
    3035#swollen/ all data output
    3136basename = 'source'
    … …  
    4247outputtimedir = outputdir + local_time + sep
    4348polygondir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'polygons'+sep
     49alex_datadir = home+sep+state+sep+scenario_dir_name+sep+'elevation_final'+sep+'points'+sep
    4450
    45 #gauge_filename = gaugedir + 'sydney_slide_gauges.csv'
    4651codedir = getcwd()+sep                           
    4752codedirname = codedir + 'project.py'
    … …  
    4954
    5055# Necessary if using point datasets, rather than grid
    51 gc_dem_name = datadir + gc250
     56gc_dem_name = alex_datadir + gc5
     57offshore_name = datadir + offshore
     58coast_name = datadir + coastline
    5259combined_dem_name   = datadir + 'gc_elevation'
    5360
     61gauge_name = gaugedir + 'gc_gauges.csv'
     62buildings_filename = gaugedir + 'gc_nexis.csv'
     63buildings_filename_out = 'gc_nexis_modified.csv'
    5464###############################
    5565# Domain definitions
    … …  
    5767
    5868# bounding polygon for study area
    59 polyAll = read_polygon(polygondir+'extent.csv')
     69polyAll = read_polygon(polygondir+'extentv2.csv')
    6070
    6171print 'Area of bounding polygon', polygon_area(polyAll)/1000000.0
    … …  
    7080print 'Area of local polygon', polygon_area(poly_int)/1000000.0
    7181
     82poly_bathy = read_polygon(polygondir+'mainland_only.csv')
     83
    7284###################################################################
    7385# Clipping regions for export to asc and regions for clipping data
    … …  
    7587
    7688# exporting asc grid
    77 eastingmin = 538450
     89eastingmin = 534300
    7890eastingmax = 546775
    79 northingmin = 6890540
    80 northingmax = 6902925
     91northingmin = 6890000
     92northingmax = 6912925

  • anuga_work/production/gold_coast_2007/run_goldcoast.py

    r4294 r4310  
    5555gc_dem_name = project.gc_dem_name
    5656meshname = project.meshname+'.msh'
    57 
     57'''
    5858# creates DEM from asc data
    5959convert_dem_from_ascii2netcdf(gc_dem_name, use_cache=True, verbose=True)
    … …  
    6262dem2pts(gc_dem_name, use_cache=True, verbose=True)
    6363
    64 G = Geospatial_data(file_name = project.gc_dem_name + '.pts')
     64G = Geospatial_data(file_name = project.gc_dem_name + '.pts') + \
     65    Geospatial_data(file_name = project.offshore_name + '.xya') +\
     66    Geospatial_data(file_name = project.coast_name + '.xya')
    6567
    6668G.export_points_file(project.combined_dem_name + '.pts')
    67 
    68 
     69'''
    6970#----------------------------------------------------------------------------
    7071# Create the triangular mesh based on overall clipping polygon with a tagged
    … …  
    7475
    7576from anuga.pmesh.mesh_interface import create_mesh_from_regions
    76 remainder_res = 150000.
    77 int_res = 1000.
     77remainder_res = 1000000.
     78int_res = 10000.
    7879interior_regions = [[project.poly_int, int_res]]
    7980
    … …  
    105106#-------------------------------------------------------------------------------
    106107
    107 tide = 0.0
    108 domain.set_quantity('stage', tide)
     108#tide = 3.0
     109tide = 1.0
     110from polygon import *
     111#IC = Polygon_function( [(project.poly_bathy, 0.)], default = tide)
     112#IC = Polygon_function( [(poly, 0.)], default = tide)
     113from anuga.utilities.polygon import read_polygon
     114poly = read_polygon('test.csv')
     115IC = Polygon_function( [(poly, 1000.)], default = 0.0)
     116domain.set_quantity('stage', IC)
    109117domain.set_quantity('friction', 0.0)
    110118domain.set_quantity('elevation',
    … …  
    122130Br = Reflective_boundary(domain)
    123131Bd = Dirichlet_boundary([tide,0,0])
     132# 10 min period = 1/600 frequency
     133from math import sin
    124134Bw = Time_boundary(domain=domain, # Time dependent boundary
    125                    f=lambda t: [(60<t<660)*6., 0, 0])
     135                   #f=lambda t: [1.0*sin(t/600.), 0, 0])
     136                   f=lambda t: [(60<t<3660)*0.1, 0, 0])
    126137
    127138domain.set_boundary( {'e0': Bw,  'e1': Bd, 'e2': Bd, 'e3': Bd} )
    … …  
    135146t0 = time.time()
    136147
    137 for t in domain.evolve(yieldstep = 10, finaltime = 5000):
     148for t in domain.evolve(yieldstep = 10, finaltime = 15000):
    138149    domain.write_time()
    139150    domain.write_boundary_statistics(tags = 'e0')

  • anuga_work/production/onslow_2006/run_building_inundation.py

    r3316 r4310  
    3333
    3434add_depth_and_momentum2csv(swwfile, buildings_filename, buildings_filename_out)
    35 
     35inundation_damage(swwfile, buildings_filename, buildings_filename_out)
    3636
    3737print '\n Augmented building file written to %s \n' \

  • anuga_work/production/perth_2006/run_building_inundation.py

    r4228 r4310  
    3333
    3434add_depth_and_momentum2csv(swwfile, buildings_filename, buildings_filename_out)
    35 
     35inundation_damage(swwfile, buildings_filename, buildings_filename_out)
    3636
    3737print '\n Augmented building file written to %s \n' \

  • anuga_work/production/pt_hedland_2006/run_building_inundation.py

    r4228 r4310  
    3333
    3434add_depth_and_momentum2csv(swwfile, buildings_filename, buildings_filename_out)
    35 
     35inundation_damage(swwfile, buildings_filename, buildings_filename_out)
    3636
    3737print '\n Augmented building file written to %s \n' \
Note: See TracChangeset for help on using the changeset viewer.