Ignore:
Timestamp:
Oct 3, 2006, 5:28:47 PM (18 years ago)
Author:
sexton
Message:

updates (clipping 25m data into 50m grid, putting MOST boundary on relevant boundary segments, change saving unique vertices to False)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/hobart_2006/run_hobart.py

    r3679 r3683  
    6262# filenames
    6363onshore_offshore_dem_name = project.onshore_offshore_dem_name
     64onshore_offshore_dem_name_25 = project.onshore_offshore_dem_name_25
    6465meshname = project.meshname+'.msh'
    6566source_dir = project.boundarydir
     
    6768copied_files = False
    6869
    69 # create DEM from asc data - 25m data
     70# create DEM from 50m asc data
    7071convert_dem_from_ascii2netcdf(onshore_offshore_dem_name, use_cache=True, verbose=True)
    7172
    72 #creates pts file for combined DEM - 25
     73# creates pts file for combined 50m DEM
    7374dem2pts(onshore_offshore_dem_name, use_cache=True, verbose=True)
    7475
     76# 25m data (clipping the around the Hobart area)
     77convert_dem_from_ascii2netcdf(onshore_offshore_dem_name_25, use_cache=True, verbose=True)
     78# creates pts file for 25m data around Hobart
     79dem2pts(onshore_offshore_dem_name_25, project.hobart_dem_name_25,
     80        easting_min=project.eastingmin25,
     81        easting_max=project.eastingmax25,
     82        northing_min=project.northingmin25,
     83        northing_max= project.northingmax25,
     84        use_cache=True,
     85        verbose=True)
     86
     87# combining the 50m and Hobart 25m data
     88combine_rectangular_points_files(project.hobart_dem_name_25 + '.pts',
     89                                 project.onshore_offshore_dem_name + '.pts',
     90                                 project.combined_dem_name + '.pts')
     91
     92# 25m data (clipping the around site 24 on Bruny Island)
     93convert_dem_from_ascii2netcdf(onshore_offshore_dem_name_25, use_cache=True, verbose=True)
     94# creates pts file for 25m data around site 24 at Bruny Island
     95dem2pts(onshore_offshore_dem_name_25, project.bruny_dem_name_25,
     96        easting_min=project.eastingmin25_2,
     97        easting_max=project.eastingmax25_2,
     98        northing_min=project.northingmin25_2,
     99        northing_max= project.northingmax25_2,
     100        use_cache=True,
     101        verbose=True)
     102
     103# combining the 50m and Hobart 25m data with Bruny Island 25m data
     104combine_rectangular_points_files(project.bruny_dem_name_25 + '.pts',
     105                                 project.combined_dem_name + '.pts',
     106                                 project.combined_dem_name_2 + '.pts')
     107
    75108# create geospatial data set and export
    76 G = Geospatial_data(file_name = project.onshore_offshore_dem_name + '.pts')
    77 G.export_points_file(project.combined_dem_name + '.pts')
     109#G = Geospatial_data(file_name = project.onshore_offshore_dem_name + '.pts')
     110#G.export_points_file(project.combined_dem_name + '.pts')
    78111
    79112#----------------------------------------------------------------------------
     
    104137                              'e12': [12], 'e13': [13], 'e14': [14],
    105138                              'e15': [15]},
    106            'maximum_triangle_area': 250000,
    107            'filename': meshname,           
     139           'maximum_triangle_area': 2000000,
     140           'filename': meshname,
    108141           'interior_regions': interior_regions},
    109142          verbose = True, evaluate=False)
     
    123156domain.set_quantities_to_be_stored(['stage', 'xmomentum', 'ymomentum'])
    124157domain.set_minimum_storable_height(0.01)
    125 domain.set_store_vertices_uniquely(True)  # for writting to sww
     158domain.set_store_vertices_uniquely(False)  # for writting to sww
    126159
    127160#-------------------------------------------------------------------------------                                 
     
    136169domain.set_quantity('elevation',
    137170#                    filename = project.onshore_dem_name + '.pts',
    138                     filename = project.combined_dem_name + '.pts',
     171                    filename = project.combined_dem_name_2 + '.pts',
    139172#                    filename = project.offshore_dem_name + '.pts',
    140173                    use_cache = True,
     
    195228domain.set_boundary( {'e0': Bd, 'e1': Bd, 'e2': Bd, 'e3': Bd, 'e4': Bd,
    196229                        'e5': Bd, 'e6': Bd, 'e7': Bd, 'e8': Bd, 'e9': Bd,
    197                         'e10': Bd, 'e11': Bd, 'e12': Bd, 'e13': Bd, 'e14': Bf,
    198                         'e15': Bd} )
     230                        'e10': Bd, 'e11': Bd, 'e12': Bf, 'e13': Bf, 'e14': Bf,
     231                        'e15': Bf} )
    199232
    200233#-------------------------------------------------------------------------------                                 
Note: See TracChangeset for help on using the changeset viewer.