Ignore:
Timestamp:
Mar 2, 2010, 2:31:30 PM (14 years ago)
Author:
fountain
Message:

updates to mandurah and bunbury storm surge models

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/bunbury_storm_surge_2009/build_elevation.py

    r7613 r7647  
    8484keylist = project.ascii_grid_filenames + project.point_filenames
    8585
    86 G1 = geospatial_data['swwa_10m'].clip_outside(project.elevation_clip_box)
    87 G2 = geospatial_data['bunbury_5m'].clip(project.elevation_clip_box)
     86G1 = geospatial_data['swwa_10m_IC'].clip_outside(project.elevation_clip_box)
     87G2 = geospatial_data['bunbury_5m_IC'].clip(project.elevation_clip_box)
    8888#~ G2 = geospatial_data['bunbury_nth_a'].clip(project.elevation_clip_box)
    8989#~ G3 = geospatial_data['bunbury_nth_b'].clip(project.elevation_clip_box)
     
    9494G5 = geospatial_data['Busselton_NavyFinal_Clip_ss.txt']
    9595G6 = geospatial_data['Leschenault_TIN.txt']
    96 G7 = geospatial_data['DPI5U1A02_01a_edited.txt']
    97 G8 = geospatial_data['DPI5U1A02_01b_edited.txt']
    98 G9 = geospatial_data['DPI5U1A02_01c_edited.txt']
    99 G10 = geospatial_data['DPI5U1A02_01d_edited.txt']
    100 G11 = geospatial_data['DPI5U1A02_01e_edited.txt']
     96G7 = geospatial_data['Leschenault_Inlet_TIN.txt']
     97G8 = geospatial_data['DPI5U1A02_01a_edited.txt']
     98G9 = geospatial_data['DPI5U1A02_01b_edited.txt']
     99G10 = geospatial_data['DPI5U1A02_01c_edited.txt']
     100G11 = geospatial_data['DPI5U1A02_01d_edited.txt']
     101G12 = geospatial_data['DPI5U1A02_01e_edited.txt']
    101102
    102103#~ G_list = [G1, G2, G3, G4, G5, G6, G7, G8, G9, G10]
     
    113114#~ print 'G10', G10.attributes
    114115
    115 G = G1 + G2 + G3 + G4 + G5 + G6 + G7 + G8 + G9 + G10 + G11
     116G = G1 + G2 + G3 + G4 + G5 + G6 + G7 + G8 + G9 + G10 + G11 + G12
    116117
    117118# G = None
     
    145146G10.export_points_file(join(project.topographies_folder, keylist[9] +'_export.txt'))
    146147G11.export_points_file(join(project.topographies_folder, keylist[10] +'_export.txt'))
     148G12.export_points_file(join(project.topographies_folder, keylist[11] +'_export.txt'))
Note: See TracChangeset for help on using the changeset viewer.