Changeset 4047 for anuga_work/development/friction_dam_2006

Timestamp:

Dec 4, 2006, 10:11:34 AM (18 years ago)

Author:

duncan

Message:

Checking in the scripts I've been using to determine the optimised mannings value from UQ and UA experiments

Location:

anuga_work/development/friction_dam_2006

Files:

• calc_norm.py (modified) (3 diffs)
• run_dam.py (modified) (3 diffs)

anuga_work/development/friction_dam_2006/calc_norm.py

@@ -1 +1 @@
 """
-
+Idea's
+have it use cache!
 """
 
@@ -97 +98 @@
                                       interpolation_points = gauge_locations,
                                       verbose = True,
-                                      use_cache = False)
+                                      use_cache = True)
         # file_instance(0.02, point_id=0) #This works
         norm_sum = 0
@@ -130 +131 @@
 if __name__ == "__main__":
     calc_norm('fromD_0.2_Slope0.05_x0.2-0.6.csv',
-              is_trial_run = True,
-              outputdir_name='friction_set')
+              is_trial_run = False,
+              outputdir_name='friction_set_A')
     

anuga_work/development/friction_dam_2006/run_dam.py

@@ -3 +3 @@
 
 Ole Nielsen and Duncan Gray, GA - 2006 
+
+
+Issues
+If running this is hand-set-up parallel, the python files are overwritten.
+
 """
 
@@ -55 +60 @@
                     dirname(project.__file__) \
                     + sep + project.__name__+'.py')
+    copy (pro_instance.codedir + 'run_dam.py',
+          pro_instance.outputdir + 'run_dam' + str(friction)+ '.py')
     copy (pro_instance.codedir + 'create_mesh.py',
           pro_instance.outputdir + 'create_mesh.py')
@@ -145 +152 @@
 #-------------------------------------------------------------
 if __name__ == "__main__":
-    for friction in [ 0.015]:
+    for friction in [ 0.0275]:
         main(friction, is_trial_run = False, outputdir_name='friction_set_A')
     
+    #    to do
+    # .075, -0-1
+    # 0.0075 - 0-16
+    #.01, - 0-2
+    #
+    #.0125, - 0-6
+    #0.015, - 0-17
+    #0.02, - 0-14
+    #$0.03]: - 0-7
+    # [0.016, - 0-11
+    #0.017, -0-8
+    #0.018, -0-9
+    #0.019]: - 0-12
+    # 0.005 - 0-20
+    # 0.02 - 4
+    # 0.016 - 8
+    # 0.017 - 9
+    # 0.018 - 1
+    # 0.019 -17
+    # calc - 7
+    # 0.0075 -12
 
+    # Let's fill some gaps
+    # 0.0225 - 7
+    # 0.025 - 12
+    # 0.0275 - 6
     #for friction in [0.0, .0025, .005, .0075, .01, .0125, 0.015, 0.02, 0.03]:
     #for friction in [0.016, 0.017, 0.018, 0.019]: