Changeset 3112

Timestamp:

Jun 8, 2006, 11:56:20 AM (19 years ago)

Author:

nicholas

Message:

Submission of 2nd draft friction report

Location:

development/momentum_sink

Files:

• CART_scans/CART_tree_rot.jpg (added)
• CART_scans/Thumbs.db (modified) (previous)
• Documents/FinalFrictionReport_2nd_draft.doc (added)
• scripts/loop_buildings.py (modified) (2 diffs)
• scripts/loop_friction.py (modified) (5 diffs)

development/momentum_sink/scripts/loop_buildings.py

@@ -23 +23 @@
 
 
-# creates buildign scenarios from widths A to B in steps of n >> list(Numeric.arange(A,B,n))
+# creates buildign scenarios from widths A to B in steps of n >>
+# >> list(Numeric.arange(A,B,n))
 DR = list(Numeric.arange(3,7,1)) 
 #DR = [4,5] #[3,5,7,9,11,13,15,17,19,21,23]
@@ -67 +68 @@
     Br = Reflective_boundary(domain)
     Bt = Transmissive_boundary(domain)
-    Bdw = Dirichlet_boundary([5,0,0])   # wave height   # <<<<<<<<< CHANGE WAVE DEPTH HERE >>>>>>>>>>>> 
+    # <<<<<<<<< CHANGE WAVE DEPTH HERE >>>>>>>>>>>>     
+    Bdw = Dirichlet_boundary([5,0,0])   # inundating wave height        
     Bdb = Dirichlet_boundary([0,0,0])   # rear boundary, keep at zero.
     Bw = Time_boundary(domain=domain,

development/momentum_sink/scripts/loop_friction.py

@@ +1 @@
+""" File used for looping 'friction block.py' script to create many
+    friction scenarios with variations in Manning's n """
 
 #Convention for strings representing files
@@ -4 +6 @@
 #           #name does not have the extension
 
-
-# set friction value for zone
 
 import time
@@ -17 +17 @@
 from region import Set_region
 from pyvolution.least_squares import fit_to_mesh_file, DEFAULT_ALPHA 
-
-
 from friction_block import create_mesh
-#from building_generator import create_mesh
-
 from pmesh.mesh import importMeshFromFile
 
+# generation of friction value list
 Friction_large = list(Numeric.arange(1,120,1))
 Friction = list(Numeric.arange(0.025,1,0.025))
@@ -41 +38 @@
     print "-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-"
     t0 = time.time()
-
     #meshname, triagle_count =create_mesh(100,mesh_file=meshname,triangles_in_name=True)
-
     meshname, triagle_count = cache(create_mesh,(100),
                                     {'mesh_file':meshname,
@@ -52 +47 @@
     #meshname = 'build.tsh'
     #outputname = outputname[:-4] + '_' + str(triagle_count) + outputname[-4:]
-
     print 'Initialising the mesh took %.2f seconds' %(time.time()-t0)  
-
-
     #meshname = importMeshFromFile('build.tsh')
-
     #Setup domain
     domain = cache(pmesh_to_domain_instance, (meshname, Domain),