Changeset 2539


Ignore:
Timestamp:
Mar 13, 2006, 4:47:16 PM (18 years ago)
Author:
nick
Message:
 
Location:
production/onslow_2006
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • production/onslow_2006/project.py

    r2519 r2539  
    3838
    3939#swollen/ all data output
    40 basename = 'source'
     40basename = 'source_fine_clipped'
    4141
    4242if sys.platform == 'win32':
  • production/onslow_2006/run_onslow.py

    r2522 r2539  
    2121# Related major packages
    2222from pyvolution.shallow_water import Domain, Reflective_boundary, \
    23                             Dirichlet_boundary, Time_boundary
     23                            Dirichlet_boundary, Time_boundary, File_boundary
    2424from pyvolution.data_manager import convert_dem_from_ascii2netcdf, dem2pts
    2525from pyvolution.combine_pts import combine_rectangular_points_files
     
    5555convert_dem_from_ascii2netcdf(onshore_dem_name, use_cache=True, verbose=True)
    5656dem2pts(onshore_dem_name,
    57 #        easting_min=project.eastingmin,
    58 #        easting_max=project.eastingmax,
    59 #        northing_min=project.northingmin,
    60 #        northing_max= project.northingmax,
     57        easting_min=project.eastingmin,
     58        easting_max=project.eastingmax,
     59        northing_min=project.northingmin,
     60        northing_max= project.northingmax,
    6161        use_cache=True,
    6262        verbose=True)
     
    6464
    6565# combining the coarse and fine data
    66 combine_rectangular_points_files(project.onshore_dem_name + '.pts',
     66combine_rectangular_points_files(
    6767                                 project.coarsedemname + '.pts',
     68                                 project.onshore_dem_name + '.pts',
    6869                                 project.combineddemname + '.pts')
    6970
     
    140141domain.set_quantity('elevation',
    141142#                    0.
     143#                    filename = project.onshore_dem_name + '.pts',
    142144                    filename = project.combineddemname + '.pts',
    143145#                    filename = project.coarsedemname + '.pts',
     
    150152# Setup boundary conditions (all reflective)
    151153#-------------------------------------------------------------------------------
    152 '''
     154
    153155from pyvolution.data_manager import ferret2sww
    154156
     
    159161
    160162cache(ferret2sww,
    161       (source_dir+project.boundary_basename,
    162        project.boundary_basename),
     163      (source_dir + project.boundary_basename,
     164       source_dir + project.boundary_basename),
    163165      {'verbose': True,
    164        'minlat': south-1,
    165        'maxlat': north+1,
    166        'minlon': west-1,
    167        'maxlon': east+1,
     166#       'minlat': south - 1,
     167#       'maxlat': north + 1,
     168#       'minlon': west - 1,
     169#       'maxlon': east + 1,
     170       'minlat': south,
     171       'maxlat': north,
     172       'minlon': west,
     173       'maxlon': east,
    168174#       'origin': project.mesh_origin,
    169175       'origin': domain.geo_reference.get_origin(),
     
    175181       verbose = True)
    176182
    177 '''
     183
    178184print 'Available boundary tags', domain.get_boundary_tags()
    179185
     186Bf = File_boundary(source_dir + project.boundary_basename + '.sww',
     187                    domain, verbose = True)
    180188Br = Reflective_boundary(domain)
    181189Bd = Dirichlet_boundary([tide,0,0])
     
    186194                   f=lambda t: [(60<t<480)*6, 0, 0])
    187195
    188 domain.set_boundary( {'top': Bw, 'topleft': Br,
     196domain.set_boundary( {'top': Bf, 'topleft': Bf,
    189197                             'left': Br, 'bottom': Br,
    190198                             'bottomright': Br, 'topright': Br} )
     
    198206t0 = time.time()
    199207
    200 for t in domain.evolve(yieldstep = 100, finaltime = 100):
     208for t in domain.evolve(yieldstep = 50, finaltime = 50):
    201209    domain.write_time()
    202210    domain.write_boundary_statistics(tags = 'top')     
Note: See TracChangeset for help on using the changeset viewer.