Changeset 2519


Ignore:
Timestamp:
Mar 9, 2006, 5:32:48 PM (19 years ago)
Author:
nick
Message:

update to Onslow

Location:
production/onslow_2006
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • production/onslow_2006/project.py

    r2470 r2519  
    99
    1010from pmesh.create_mesh import convert_points_from_latlon_to_utm
     11
     12from coordinate_transforms.redfearn import degminsec2decimal_degrees
    1113               
    1214#Making assumptions about the location of scenario data
     
    1618# 250m data to be provided
    1719coarsename = 'onsl_bathydem250' # get from Neil/Ingo (DEM or topo data)
    18 """
     20
    1921# 30m data to be provided
    20 finename = 'onshore_30' # get from Neil/Ingo (DEM or topo data)
    21 """
     22onshore_name = 'onslow_onshore_30m_dted' # get from Neil/Ingo (DEM or topo data)
     23
     24
     25boundary_basename = 'SU-AU'
    2226
    2327# clipping region for fine elevation data
    24 west = 240000
    25 east = 340000
    26 south = 7570000
    27 north = 7645000
     28eastingmin = 250000
     29eastingmax = 330000
     30northingmin = 7580000
     31northingmax = 7635000
     32
     33south = degminsec2decimal_degrees(-22,00,0)
     34north = degminsec2decimal_degrees(-21,10,0)
     35west = degminsec2decimal_degrees(114,30,0)
     36east = degminsec2decimal_degrees(115,30,0)
     37
    2838
    2939#swollen/ all data output
     
    4555
    4656coarsedemname = datadir + coarsename
    47 """
    48 finedemname = datadir + finename
     57
     58onshore_dem_name = datadir + onshore_name
     59
    4960combineddemname = datadir + 'onslow_combined_elevation'
    50 """
    5161
    5262outputname = outputdir + basename  #Used by post processing
  • production/onslow_2006/run_onslow.py

    r2470 r2519  
    4141# filenames
    4242coarsedemname = project.coarsedemname
    43 '''
    44 #finedemname = project.finedemname
    45 '''
     43
     44onshore_dem_name = project.onshore_dem_name
     45
    4646meshname = project.meshname+'.msh'
     47
     48source_dir = project.boundarydir
    4749
    4850# coarse data
    4951convert_dem_from_ascii2netcdf(coarsedemname, use_cache=True, verbose=True)
    5052dem2pts(coarsedemname, use_cache=True, verbose=True)
    51 '''
     53
    5254# fine data (clipping the points file to smaller area)
    53 convert_dem_from_ascii2netcdf(finedemname, use_cache=True, verbose=True)
    54 dem2pts(finedemname,
    55         easting_min=project.eastingmin,
    56         easting_max=project.eastingmax,
    57         northing_min=project.northingmin,
    58         northing_max= project.northingmax,
     55convert_dem_from_ascii2netcdf(onshore_dem_name, use_cache=True, verbose=True)
     56dem2pts(onshore_dem_name,
     57#        easting_min=project.eastingmin,
     58#        easting_max=project.eastingmax,
     59#        northing_min=project.northingmin,
     60#        northing_max= project.northingmax,
    5961        use_cache=True,
    6062        verbose=True)
     
    6264
    6365# combining the coarse and fine data
    64 combine_rectangular_points_files(project.finedemname + '.pts',
     66combine_rectangular_points_files(project.onshore_dem_name + '.pts',
    6567                                 project.coarsedemname + '.pts',
    6668                                 project.combineddemname + '.pts')
    67 '''
     69
    6870
    6971#-------------------------------------------------------------------------------                                 
     
    109111domain.set_name(project.basename)
    110112domain.set_datadir(project.outputdir)
    111 domain.set_quantities_to_be_stored(['stage', 'xmomentum', 'ymomentum'])
     113domain.set_quantities_to_be_stored(['stage'])
    112114
    113115
     
    136138domain.set_quantity('stage', tide)
    137139domain.set_quantity('friction', 0.0)
     140
    138141domain.set_quantity('elevation',
    139 #                    filename = project.combineddemname + '.pts',
    140                     filename = project.coarsedemname + '.pts',
     142#                    0.
     143                    filename = project.combineddemname + '.pts',
     144#                    filename = project.coarsedemname + '.pts',
    141145                    use_cache = True,
    142146                    verbose = True
     
    147151# Setup boundary conditions (all reflective)
    148152#-------------------------------------------------------------------------------
    149 
     153'''
    150154from pyvolution.data_manager import ferret2sww
    151155
     
    163167       'minlon': west-1,
    164168       'maxlon': east+1,
    165        'origin': project.mesh_origin,
     169#       'origin': project.mesh_origin,
     170       'origin': domain.geo_reference.get_origin(),
    166171       'mean_stage': tide,
    167172       'zscale': 1,                 #Enhance tsunami
     
    171176       verbose = True)
    172177
    173 
     178'''
    174179print 'Available boundary tags', domain.get_boundary_tags()
    175180
     
    194199t0 = time.time()
    195200
    196 for t in domain.evolve(yieldstep = 50, finaltime = 100):
     201for t in domain.evolve(yieldstep = 100, finaltime = 100):
    197202    domain.write_time()
    198203    domain.write_boundary_statistics(tags = 'top')     
Note: See TracChangeset for help on using the changeset viewer.