towards production

data/tools/support/friction_stats.py

@@ -2,7 +2,7 @@ import numpy as np
 import matplotlib.pyplot as plt
 
 def friction_stats(h5file, body_ID, FINAL_TIME, patch=[], interval=[], TANGENTIAL_COORDS=1):
-    body = 'body' + str(body_ID) # 'frictionalBoundary'     'body' + str(body_ID)
+    body = 'body' + str(body_ID)   # 'frictionalBoundary'     'body' + str(body_ID)
 
     coords = np.array(h5file[body + '/coordinates'])
     if len(patch) == 0:
@@ -33,7 +33,7 @@ def friction_stats(h5file, body_ID, FINAL_TIME, patch=[], interval=[], TANGENTIA
     #ax_slip.plot(time, max_v, color='gray', linestyle='--')
     ax_slip.set_ylabel('slip rate V [m/s]')
     ax_slip.set_yscale('log')
-    ax_slip.set_ylim([1e-6,1e-2])
+    #ax_slip.set_ylim([1e-6,1e-2])
     #-------------------------
 
     print(np.min(min_v))
@@ -46,8 +46,8 @@ def friction_stats(h5file, body_ID, FINAL_TIME, patch=[], interval=[], TANGENTIA
     states_tx = states_t[:,patch]
     # statistics
     avg_states = np.average(states_tx, axis=1)
-    min_states = np.min(states_tx, axis=1)
-    max_states = np.max(states_tx, axis=1)
+    #min_states = np.min(states_tx, axis=1)
+    #max_states = np.max(states_tx, axis=1)
     # plot
     ax_state = fig.add_subplot(2, 1, 2)
     #ax_state.plot(time, min_states, color='gray', linestyle='--')