Changeset 6062


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

updated scripts for new data in bunbury and marina in busselton

Location:
anuga_work/production/busselton
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/busselton/build_busselton.py

    r5829 r6062  
    7676print "creates DEMs from asc data"
    7777convert_dem_from_ascii2netcdf(onshore_in_dir_name, basename_out=onshore_dir_name, use_cache=True, verbose=True)
     78convert_dem_from_ascii2netcdf(offshore_in_dir_name5, basename_out=offshore_dir_name5, use_cache=True, verbose=True)
    7879
    7980# creates pts file for onshore DEM
    8081print "creates pts file for onshore DEM"
    8182dem2pts(onshore_dir_name ,use_cache=True,verbose=True)
     83dem2pts(offshore_dir_name5 ,use_cache=True,verbose=True)
    8284
    8385# create onshore pts files
     
    101103G_off4 = Geospatial_data(file_name = offshore_in_dir_name4)
    102104print'create Geospatial data9 objects from topographies'
    103 G_off5 = Geospatial_data(file_name = offshore_in_dir_name5)
     105G_off5 = Geospatial_data(file_name = offshore_dir_name5 + '.pts')
    104106
    105107
     
    109111
    110112print'add all geospatial objects'
    111 G = G1 + G2 + G3 + G_off + G_off1 + G_off2 + G_off3 + G_off4 + G_off5
     113G = G1 + G2 + G3 + G_off + G_off1 + G_off2 + G_off3 + G_off4
    112114
    113115print'clip combined geospatial object by bounding polygon'
    114 G_clipped = G.clip(project.poly_all)
     116G_clip = G.clip_outside(project.poly_aoi1)
     117G_all = G_clip + G_off5
     118G_clipped = G_all.clip(project.poly_all)
    115119
    116120print'export combined DEM file'
     
    118122    mkdir (project.topographies_dir)
    119123G_clipped.export_points_file(project.combined_dir_name + '.pts')
    120 #G_clipped.export_points_file(project.combined_dir_name + '.txt') #Use for comparision in ARC
     124G_clipped.export_points_file(project.combined_dir_name + '.txt') #Use for comparision in ARC
    121125
  • anuga_work/production/busselton/project.py

    r6018 r6062  
    4343# Model specific parameters. One or all can be changed each time the
    4444# run_scenario script is executed
    45 tide = 0.6                #0.6
    46 event_number = 27255 # Java 9.3 worst case for Perth
     45tide = 0              #0.6
     46#event_number = 27255 # Java 9.3 worst case for Perth
    4747#event_number = 68693 # Sumatra 9.2
    48 #event_number = 27283  # Java 9.3 original
     48event_number = 27283  # Java 9.3 original
    4949alpha = 0.1             # smoothing parameter for mesh
    5050friction=0.01           # manning's friction coefficient
     
    113113offshore_name2 = 'Busselton_Digitised.txt'
    114114offshore_name3 = 'Busselton_250m.txt' # for areas that were heading to zero - 2005 Bathymetry grid
    115 offshore_name4 = 'Bunbury_TIN.txt' # for area within Bunbury 500 mesh less than zero generated from TIN
    116 offshore_name5 = 'Busselton_TIN.txt' # for area within Busselton 500 mesh less than zero generated from TIN
     115offshore_name4 = 'DPI.txt' # for area within Bunbury 500 mesh less than zero generated from TIN
     116offshore_name5 = 'topo_20m_buss_1km' # for area within Busselton 500 mesh less than zero generated from TIN
    117117
    118118# gauges - used in get_timeseries.py
     
    162162offshore_in_dir_name2 = topographies_in_dir + offshore_name2 #Digitised Fairsheet
    163163offshore_in_dir_name3 = topographies_in_dir + offshore_name3 #250m
    164 offshore_in_dir_name4 = topographies_in_dir + offshore_name4 #Bunbury TIN
    165 offshore_in_dir_name5 = topographies_in_dir + offshore_name5 #Busselton TIN
     164offshore_in_dir_name4 = topographies_in_dir + offshore_name4 #Bunbury DPI
     165offshore_in_dir_name5 = topographies_in_dir + offshore_name5 #Busselton Topo
    166166
    167167# where the output data sits
     
    211211poly_mainland = read_polygon(polygons_dir+'initial_condition.csv')
    212212
     213#Land, to set the initial stage/water to be offcoast only
     214poly_marina = read_polygon(polygons_dir+'initial_condition_marina.csv')
     215
    213216# Initial bounding polygon for data clipping
    214217poly_all = read_polygon(polygons_dir+'poly_all_extend.csv')
  • anuga_work/production/busselton/run_busselton.py

    r6023 r6062  
    132132
    133133    # sets the initial stage in the offcoast region only
    134     IC = Polygon_function( [(project.poly_mainland, 0)], default = kwargs['tide'],
     134    IC = Polygon_function( [(project.poly_mainland, 0),(project.poly_marina, 0)], default = kwargs['tide'],
    135135                             geo_reference = domain.geo_reference)
    136136    domain.set_quantity('stage', IC)
Note: See TracChangeset for help on using the changeset viewer.