Changeset 7555

Timestamp:

Nov 11, 2009, 9:53:07 AM (15 years ago)

Author:

jgriffin

Message:

 

Location:

misc/tools/event_selection/wave_energy

Files:

• wave_energy.doc (modified) (previous)
• wave_energy.py (modified) (11 diffs)
• wave_energy_and_momentum.py (added)

misc/tools/event_selection/wave_energy/wave_energy.py

@@ -19 +19 @@
 import glob
 import numpy
+import run_up
 
 # Change path here!!!
@@ -42 +43 @@
 max_energy_list = []
 max_energy_inst_list = []
+run_up_list = []
 counter = 0
 
@@ -69 +71 @@
     max_energy = 0
     max_energy_inst = 0
+    max_stage = 0
+    temp_max_stage = 0
     
     ###############################################################
@@ -85 +89 @@
         x_vel = x_momentum[i]/(100+stage[i])
         y_vel = y_momentum[i]/(100+stage[i])
+        velocity = numpy.sqrt(x_vel*x_vel,y_vel*y_vel)
         # Kinetic energy KE = 1/2 roh h u^2
         KE = (1./2.)*(numpy.power(x_vel,2)+numpy.power(y_vel,2))*(1000)*(100+stage[i])
@@ -105 +110 @@
             if temp_max_energy > max_energy:
                 max_energy = temp_max_energy
-                
+            if temp_max_stage > max_stage:
+                max_stage = temp_max_stage
+                max_time = time[i] - start_time   
             energy_sum = energy
             start_time = time[i]
+            max_stage_temp = 0
             
         elif stage[i] < 0 and sign != -1:
@@ -114 +122 @@
             if temp_max_energy > max_energy:
                 max_energy = temp_max_energy
-                
+            if temp_max_stage > max_stage:
+                max_stage = temp_max_stage
+                max_time = time[i] - start_time                 
             energy_sum = energy
             start_time = time[i]
@@ -120 +130 @@
         elif stage[i] == 0 and sign != 0:
             sign = 0
+            if temp_max_stage > max_stage:
+                max_stage = temp_max_stage
+                max_time = time[i] - start_time 
             energy_sum = energy
             start_time = time[i]
-            
+                      
         else:
             energy_sum += energy
+            temp_max_stage = max(temp_max_stage,stage[i])
             
     # Add maximum values to list        
@@ -130 +144 @@
     max_energy_list.append( max_energy)
 
+    # Account for both halves of the wave
+    max_time = max_time*2
+    print 'max stage',max_stage
+    print 'max time',max_time
+    #####################################################################
+    # call run-up module
+    #####################################################################
+    Run_up = run_up.Madsen(100,0.01,max_stage,max_time)
+    print 'runup', Run_up
+    run_up_list.append(Run_up)
+
+    
     counter += 1
+    
+
 
     
@@ -137 +165 @@
 total_max_inst_energy = 0
 
+
+
+
+
 # Print results for each file and calculate maximum value across all events
 for filename in filepathlist:
@@ -142 +174 @@
     print 'max crest-integrated energy:', max_energy_list[counter],'J.s/m^2'
     print 'max instantatneous energy:', max_energy_inst_list[counter],'J/m^2'
-
+    print 'Run_up:', run_up_list[counter], 'm'
+    
     if max_energy_list[counter] > total_max_energy:
         total_max_energy = max_energy_list[counter]
@@ -156 +189 @@
 print '\n File with maximum instantaneous energy is:', filepathlist[max_energy_inst_index]
 print 'Energy is', max_energy_inst_list[max_energy_inst_index], 'J/m^2'
+
+