Changeset 3112


Ignore:
Timestamp:
Jun 8, 2006, 11:56:20 AM (18 years ago)
Author:
nicholas
Message:

Submission of 2nd draft friction report

Location:
development/momentum_sink
Files:
2 added
3 edited

Legend:

Unmodified
Added
Removed
  • development/momentum_sink/scripts/loop_buildings.py

    r2938 r3112  
    2323
    2424
    25 # creates buildign scenarios from widths A to B in steps of n >> list(Numeric.arange(A,B,n))
     25# creates buildign scenarios from widths A to B in steps of n >>
     26# >> list(Numeric.arange(A,B,n))
    2627DR = list(Numeric.arange(3,7,1))
    2728#DR = [4,5] #[3,5,7,9,11,13,15,17,19,21,23]
     
    6768    Br = Reflective_boundary(domain)
    6869    Bt = Transmissive_boundary(domain)
    69     Bdw = Dirichlet_boundary([5,0,0])   # wave height   # <<<<<<<<< CHANGE WAVE DEPTH HERE >>>>>>>>>>>>
     70    # <<<<<<<<< CHANGE WAVE DEPTH HERE >>>>>>>>>>>>     
     71    Bdw = Dirichlet_boundary([5,0,0])   # inundating wave height       
    7072    Bdb = Dirichlet_boundary([0,0,0])   # rear boundary, keep at zero.
    7173    Bw = Time_boundary(domain=domain,
  • development/momentum_sink/scripts/loop_friction.py

    r2927 r3112  
     1""" File used for looping 'friction block.py' script to create many
     2    friction scenarios with variations in Manning's n """
    13
    24#Convention for strings representing files
     
    46#           #name does not have the extension
    57
    6 
    7 # set friction value for zone
    88
    99import time
     
    1717from region import Set_region
    1818from pyvolution.least_squares import fit_to_mesh_file, DEFAULT_ALPHA
    19 
    20 
    2119from friction_block import create_mesh
    22 #from building_generator import create_mesh
    23 
    2420from pmesh.mesh import importMeshFromFile
    2521
     22# generation of friction value list
    2623Friction_large = list(Numeric.arange(1,120,1))
    2724Friction = list(Numeric.arange(0.025,1,0.025))
     
    4138    print "-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-"
    4239    t0 = time.time()
    43 
    4440    #meshname, triagle_count =create_mesh(100,mesh_file=meshname,triangles_in_name=True)
    45 
    4641    meshname, triagle_count = cache(create_mesh,(100),
    4742                                    {'mesh_file':meshname,
     
    5247    #meshname = 'build.tsh'
    5348    #outputname = outputname[:-4] + '_' + str(triagle_count) + outputname[-4:]
    54 
    5549    print 'Initialising the mesh took %.2f seconds' %(time.time()-t0) 
    56 
    57 
    5850    #meshname = importMeshFromFile('build.tsh')
    59 
    6051    #Setup domain
    6152    domain = cache(pmesh_to_domain_instance, (meshname, Domain),
Note: See TracChangeset for help on using the changeset viewer.