Changeset 7019


Ignore:
Timestamp:
May 12, 2009, 11:37:41 AM (11 years ago)
Author:
jgriffin
Message:
 
Location:
anuga_work/production/new_south_wales/batemans_bay
Files:
2 added
3 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/new_south_wales/batemans_bay/build_elevation.py

    r6427 r7019  
    6767if not project.point_filenames == []:
    6868    for filename in project.point_filenames:
    69         absolute_filename = join(project.topographies_folder, filename)
     69        print filename
     70        absolute_filename = project.topographies_folder+'/'+str(filename)
    7071        G_points = Geospatial_data(file_name=absolute_filename,
    7172                                                    verbose=True)
  • anuga_work/production/new_south_wales/batemans_bay/prepare_data.py

    r6941 r7019  
    1414import glob, fnmatch
    1515
    16 folder = '//nas/gemd/georisk_models/inundation/data/new_south_wales/gosford_tsunami_scenario_2009/original_reference_data/GosfordLidar/Ground/'
     16#folder = '//nas/gemd/georisk_models/inundation/data/new_south_wales/gosford_tsunami_scenario_2009/original_reference_data/GosfordLidar/Ground/'
     17folder = '//nas/gemd/georisk_models/inundation/data/new_south_wales/gosford_tsunami_scenario_2009/original_reference_data/topography/HCR-CMA_LiDAR/05.Raw_ALS_Data_Ground/ASCII_X,Y,Z,I/'
    1718
    1819
    19 
    20 input_grd = glob.glob(folder + '*.grd')
     20#input_grd = glob.glob(folder + '*.grd')
     21input_grd = glob.glob(folder + '*.XYZI')
    2122
    2223for file_in in input_grd:
    2324    #boundary_csv = file_in
    24     out_file = file_in[:-4] + '.txt'
     25    #out_file = file_in[:-4] + '.txt'
     26    out_file = file_in[:-5] + '.txt'
    2527    print 'Creating', out_file
    2628
     
    3032
    3133    # Skip first line
    32     line = fid.readline()
     34    #line = fid.readline()
    3335
    3436    # Read remaining lines
     
    5961
    6062filelist = []
    61 filelist_file = folder+'file_list.csv'
     63filelist_file = folder+'file_list2.csv'
    6264
    6365# Create csv file of filenames
  • anuga_work/production/new_south_wales/batemans_bay/project.py

    r6989 r7019  
    88from time import localtime, strftime, gmtime
    99from os.path import join, exists
     10import csv
    1011
    1112
     
    1718# Note, the user needs to set up the directory system accordingly
    1819state = 'new_south_wales'
    19 scenario_name = 'batemans_bay'
    20 scenario_folder = 'batemans_bay_tsunami_scenario_2009'
     20scenario_name = 'gosford'
     21scenario_folder = 'gosford_tsunami_scenario_2009'
    2122
    2223#-------------------------------------------------------------------------------
    2324# Initial Conditions
    2425#-------------------------------------------------------------------------------
     26
    2527
    2628# Model specific parameters.
    2729# One or all can be changed each time the run_model script is executed
    28 tide = 0                # difference between MSL and HAT (1.0)
    29 
    30 # the event number or the mux file name
    31 ##event_number = 58129  #1 in 200 yr Puyesgur
    32 ##event_number = 58115    #1 in 500 yr
    33 ##event_number = 58226    #1 in 1000 yr
    34 ##event_number = 58284    #1 in 2000 yr
    35 ##event_number = 58286    #1 in 5000 yr
    36 ##event_number = 58242    #1 in 10000 yr Puysegur
    37 event_number = 58368    #1 in 10000 yr Puysegur
    38 ##event_number = 51436    #1 in 10000 yr New Hebrides
    39 
     30tide = 0                # difference between MSL and HAT in metres
     31zone = 56               # specify zone of model
     32event_number = 31891    # the event number or the mux file name
    4033alpha = 0.1             # smoothing parameter for mesh
    4134friction=0.01           # manning's friction coefficient
    4235starttime=0             # start time for simulation
    43 finaltime=10000         # final time for simulation
    44 
    45 setup = 'final'         # This can be one of three values
     36finaltime=1000         # final time for simulation
     37
     38setup = 'trial'         # This can be one of three values
    4639                        #    trial - coarsest mesh, fast
    4740                        #    basic - coarse mesh
     
    5952#-------------------------------------------------------------------------------
    6053
    61 output_comment = [setup, tide, event_number, 'viz']
     54output_comment = [setup, tide, event_number, 'lidar']
    6255
    6356#-------------------------------------------------------------------------------
     
    6861# Used in build_elevation.py
    6962# Format for ascii grids, as produced in ArcGIS + a projection file
    70 ascii_grid_filenames = ['1a',# Topographic data
    71                         '1b',
    72                         '2b',
    73                         '2a_3',
    74                         '3b',
    75                         '3a',
    76                         '4a_2',
    77                         '4b',
    78                         'off1',
    79                         'off2',
    80                         'off3',
    81                         'bbhd',
    82                         'sd100031996_p',
    83                         'sd100031996_p2',
    84                         'sd100031996_p3',
    85                         'sd100031996_p4']
     63ascii_grid_filenames = ['e151_s34_clip1',
     64                        'aus197_topo2raster_umina',
     65                        'aus197_topo2raster_terrigal']   #GA data??
    8666                   
     67
     68
     69#Get long list of Lidar filenames.
     70ENV_INUNDATIONHOME = 'INUNDATIONHOME'
     71home = join(os.getenv(ENV_INUNDATIONHOME), 'data')
     72anuga_folder = join(home, state, scenario_folder, 'anuga')
     73topographies_folder = join(anuga_folder, 'topographies/original')
     74
     75file_list = join(topographies_folder, 'file_list.csv')
     76file_list2 = join(topographies_folder, 'file_list2.csv')
     77fid = csv.reader(open(file_list))
     78fid2 = csv.reader(open(file_list2))
     79lidar_filenames = []
     80for row in fid:
     81    lidar_filenames.append(','.join(row))
     82for row in fid2:
     83    lidar_filenames.append(','.join(row))
     84#print lidar_filenames
     85#lidar_filenames = list(lidar_filenames)
     86#print lidar_filenames
     87
    8788# Format for point is x,y,elevation (with header)
    88 point_filenames = ['SD100031996_jgriffin_clip.csv',
    89                    'tomaga_offshore_AHD_MGA_1997.csv',
    90                    'Batemans_BBHD_MGA_1995.csv',
    91                    'moruya_AHD_MGA_2000.csv']
     89point_filenames = ['estuaries.txt',
     90##                   'XYHI341_HSDB_SD100035103',
     91##                   'XYAHD',
     92##                   'XYSCHOOL_11_99_HSDB_SD100021002_MGA',
     93                   'ENT5C1S08_03_AHD_prepared.txt',
     94                   'ENT5C1S09_03_AHD_prepared.txt',
     95                   'ENT5C1S10_03_AHD_prepared.txt',
     96##                   'XYAHD',
     97##                   'XYHI341_HSDB_SD100035103_MGA',
     98##                   'XYHI339_F_HSDB_SD100031669_MGA',
     99                   'AUS197_MGA_AHD_1972.txt',
     100                   'hydro.txt'
     101                   ]
     102point_filenames.extend(lidar_filenames)
     103#print point_filenames
    92104         
    93105
     
    107119# Format for points easting,northing (no header)
    108120
    109 interior_regions_data = [['area_of_interest.csv', 500],
    110                         ['area_of_significance.csv', 2500],
    111                         ['shallow_water.csv', 10000]]
     121interior_regions_data =[['aoi_umina.csv', 500],
     122                        ['aoi_terrigal.csv', 500],
     123                        ['aos_umina.csv', 2500],
     124                        ['aos_terrigal.csv', 2500]]
    112125
    113126
     
    134147# Thinned ordering file from Hazard Map (geographic)
    135148# Format is index,latitude,longitude (with header)
    136 urs_order_filename = 'thinned_boundary_ordering_extend.csv'
     149urs_order_filename = 'urs_order.csv'
    137150
    138151# Landward bounding points
    139152# Format easting,northing (no header)
    140 landward_boundary_filename = 'landward_boundary_extend.csv'
     153landward_boundary_filename = 'landward_boundary.csv'
    141154
    142155# MUX input filename.
     
    210223# check various directories/files that must exist
    211224anuga_folder = join(home, state, scenario_folder, 'anuga')
    212 topographies_folder = join(anuga_folder, 'topographies')
     225topographies_folder = join(anuga_folder, 'topographies/original')
    213226polygons_folder = join(anuga_folder, 'polygons')
    214227boundaries_folder = join(anuga_folder, 'boundaries')
Note: See TracChangeset for help on using the changeset viewer.