Changeset 2956


Ignore:
Timestamp:
May 24, 2006, 11:54:27 AM (18 years ago)
Author:
nick
Message:

update to Onslow

Location:
production/onslow_2006
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • production/onslow_2006/export_results.py

    r2955 r2956  
    2323#which_var = int(raw_input('Stage = 0, Absolute Momentum = 1, Depth = 2, Speed = 3  '  ))
    2424which_var = 2
    25 sys.stderr.write(sys.stdout.data)
     25#sys.stderr.write(sys.stdout.data)
    2626if which_var == 0:  # Stage
    2727    outname = name + '_stage'
     
    4747
    4848print 'start sww2dem'
    49 sys.stderr.write(sys.stdout.data)
     49#sys.stderr.write(sys.stdout.data)
    5050sww2dem(name, basename_out = outname,
    5151            quantity = quantityname,
     
    6060            format = 'asc')
    6161
    62 sys.stderr.write(sys.stdout.data)
     62#sys.stderr.write(sys.stdout.data)
    6363
    6464#Check
  • production/onslow_2006/run_onslow.py

    r2902 r2956  
    2525from pyvolution.combine_pts import combine_rectangular_points_files
    2626from pyvolution.pmesh2domain import pmesh_to_domain_instance
    27 from geospatial_data import add_points_files
     27#from geospatial_data import add_points_files
    2828
    2929# Application specific imports
     
    8888#    print' most file', project.MOST_dir + project.boundary_basename
    8989
    90 
     90"""
    9191# fine data (clipping the points file to smaller area)
    9292# creates DEM from asc data
     
    113113print'export G'
    114114G.export_points_file(project.combined_dem_name + '.pts')
    115 
     115"""
    116116
    117117#-------------------------------------------------------------------------------                                 
     
    122122
    123123from pmesh.mesh_interface import create_mesh_from_regions
    124 '''
    125 # original
    126 interior_res = 5000
    127 high_res = 1500
    128 interior_regions = [[project.poly_onslow, high_res],
    129                     [project.poly_thevenard, interior_res],
    130                     [project.poly_coast, interior_res]]
    131 '''
     124
    132125#new
    133126region_res = 50000
     
    135128onslow_res = 500
    136129interior_regions = [[project.poly_onslow, onslow_res],
    137                     [project.poly_coast, coast_res],
    138                     [project.poly_region, region_res]]
     130                    [project.poly_coast, coast_res]]
     131#,
     132#                    [project.poly_region, region_res]]
    139133
    140134print 'number of interior regions', len(interior_regions)
     
    147141                             'bottom': [4], 'bottomright': [5],
    148142                             'topright':[6]},
    149            'maximum_triangle_area': 200000,
     143           'maximum_triangle_area': 1000000,
    150144           'filename': meshname,           
    151145           'interior_regions': interior_regions},
Note: See TracChangeset for help on using the changeset viewer.