Changeset 2645


Ignore:
Timestamp:
Mar 31, 2006, 4:20:53 PM (18 years ago)
Author:
nick
Message:

update to onslow

Location:
production/onslow_2006
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • production/onslow_2006/project.py

    r2624 r2645  
    2727offshore_name = 'onslow_offshore_points'
    2828
    29 boundary_basename = 'SU-AU'
     29boundary_basename = 'SU-AU_clip'
    3030
    3131#swollen/ all data output
     
    5050
    5151# boundary source data
    52 MOST_dir = 'f:'+sep+'3'+sep+'ehn'+sep+'users'+sep+'davidb'+sep+'tsunami'+sep+'WA_project'+sep+'SU-AU_90'+sep+'most_2'+sep+'detailed'+sep
     52#MOST_dir = 'f:'+sep+'3'+sep+'ehn'+sep+'users'+sep+'davidb'+sep+'tsunami'+sep+'WA_project'+sep+'SU-AU_90'+sep+'most_2'+sep+'detailed'+sep
    5353
    5454codedir = getcwd()+sep
     
    100100d4 = [310000, 7580000]
    101101d5 = [330000, 7610000]
     102d6 = [330000, 7620000]
    102103
    103 polyAll = [d0, d1, d2, d3, d4, d5]
     104polyAll = [d0, d1, d2, d3, d4, d5, d6]
    104105
    105106#Interior region - Onslow town
  • production/onslow_2006/run_onslow.py

    r2629 r2645  
    6161copy (project.codedir + 'run_onslow.py', project.outputdir + 'run_onslow.py')
    6262print'hello'
    63 print' most file', project.MOST_dir + project.boundary_basename+'_ha.nc'
    64 if access(project.MOST_dir + project.boundary_basename+'_ha.nc',F_OK) == 1 :
    65     print' most file', project.MOST_dir + project.boundary_basename
     63#print' most file', project.MOST_dir + project.boundary_basename+'_ha.nc'
     64#if access(project.MOST_dir + project.boundary_basename+'_ha.nc',F_OK) == 1 :
     65#    print' most file', project.MOST_dir + project.boundary_basename
    6666
    6767'''
     
    8080        use_cache=True,
    8181        verbose=True)
    82 
    83         '''
    84 '''
     82'''
     83
    8584print'create G1'
    8685G1 = Geospatial_data(file_name = project.offshore_dem_name + '.xya')
     
    9594G.new_export_points_file(project.combined_dem_name + '.pts')
    9695
    97 '''
     96
    9897#-------------------------------------------------------------------------------                                 
    9998# Create the triangular mesh based on overall clipping polygon with a tagged
     
    117116          {'boundary_tags': {'top': [0], 'topleft': [1],
    118117                             'left': [2], 'bottom': [3],
    119                              'bottomright': [4], 'topright': [5]},
     118                             'bottomright': [4], 'right': [5],
     119                             'topright':[6]},
    120120           'maximum_triangle_area': 1000000,
    121121           'filename': meshname,           
     
    162162
    163163tide = 0.
    164 '''
     164
    165165domain.set_quantity('stage', tide)
    166166domain.set_quantity('friction', 0.0)
     
    174174                    use_cache = False,
    175175                    verbose = True,
    176                     alpha = 0.1
     176                    alpha = 0.05
    177177                    )
    178 '''
     178
    179179print 'hi1'
    180180
     
    191191
    192192cache(ferret2sww,
    193 #      (source_dir + project.boundary_basename,
     193      (source_dir + project.boundary_basename,
     194       source_dir + project.boundary_basename),
     195#      (project.MOST_dir + project.boundary_basename,
    194196#       source_dir + project.boundary_basename),
    195       (project.MOST_dir + project.boundary_basename,
    196        source_dir + project.boundary_basename),
    197197      {'verbose': True,
    198198# note didn't work with the below
     
    229229domain.set_boundary( {'top': Bf, 'topleft': Bf,
    230230                             'left': Br, 'bottom': Br,
    231                              'bottomright': Br, 'topright': Bf} )
     231                             'bottomright': Br, 'right': Br, 'topright': Br} )
    232232
    233233
     
    238238t0 = time.time()
    239239
    240 for t in domain.evolve(yieldstep = 1000, finaltime = 10000):
     240for t in domain.evolve(yieldstep = 1000, finaltime = 3000):
    241241    domain.write_time()
    242242    domain.write_boundary_statistics(tags = 'top')     
    243243
    244 for t in domain.evolve(yieldstep = 50, finaltime = 10100,
     244for t in domain.evolve(yieldstep = 100, finaltime = 20000,
    245245                       skip_initial_step = True):
    246246    domain.write_time()
    247247    domain.write_boundary_statistics(tags = 'top')     
    248    
     248 
    249249print 'That took %.2f seconds' %(time.time()-t0)
    250250
Note: See TracChangeset for help on using the changeset viewer.