added runtime folder

Author: wilfonba

RunTimeTests/MeanTimeData

@@ -0,0 +1,29 @@
+Viscous
+0.00018558695652173907
+0.00020577826086956504
+0.0035785804347826084
+0.00452261956521739
+0.0037711695652173907
+0.0046880260869565215
+0.09709701956521741
+0.13853727391304346
+0.153668252173913
+0.21701702391304342
+Monopole
+0.0006033239130434785
+0.0006059804347826085
+0.0003294130434782607
+0.0007314065217391305
+0.005524604347826086
+0.010737452173913042
+0.0073501434782608675
+0.007414117391304346
+0.004163260869565218
+0.009179697826086953
+0.05305613913043479
+0.09925884130434783
+0.10973742826086959
+0.06566204130434782
+0.142803347826087
+1.3211132347826084
+0.7312944891304348

RunTimeTests/testRuntime.py

@@ -0,0 +1,62 @@
+import os
+
+ViscousJobs = [ ['C79E1D3C', 'CD9D3050'], \
+                ['9CB03CEF', 'D6BAC936', 'FB822062', 'B3AAC9C8'], \
+                ['1C0780C8', '2060F55A', '07C33719', '939D6718']]
+
+MonopoleJobs = [['EF54219C', '7FC6826B', '59D05DE9', '6B22A317', 'AF46C382', '46AA7AF8'], \
+                ['34DBFE14', '14B6198D', 'AE37D842', 'CC4F7C44', '66CFF8CC', '303B925A'], \
+                ['AD63A4A5', '8A341282', '622DEC78', 'AB04C64D', '63850240' ]]
+
+
+itr = 10
+vTimes = [[0 for x in range(10)] for y in range(itr)]
+mTimes = [[0 for x in range(17)] for y in range(itr)]  
+
+# Viscous Tests
+for k in range(0,itr):
+    os.system("rm -rf RunTimeTests/ViscousTests/*")
+    for i in ViscousJobs:
+        for j in i:
+            os.system("cp tests/" + j + "/case.py RunTimeTests/ViscousTests/case" + j + ".py")
+            os.system("./mfc.sh run RunTimeTests/ViscousTests/case" + j + ".py")
+    
+    f = open("RunTimeTests/ViscousTests/time_data.dat")
+    print(f)
+    for m in range(0, 10):
+        vTimes[k][m] = float(f.readline().split()[1])
+
+meanVTimes = []
+for m in range(0,10):
+    testSum = 0
+    for j in range(0, itr):
+        testSum = testSum + vTimes[j][m]
+    meanVTimes.append(testSum/itr)
+
+# Monopole Tests
+for k in range(0,itr):
+    os.system("rm -rf RunTimeTests/MonopoleTests/*")
+    for i in MonopoleJobs:
+        for j in i:
+            os.system("cp tests/" + j + "/case.py RunTimeTests/MonopoleTests/case" + j + ".py")
+            os.system("./mfc.sh run RunTimeTests/MonopoleTests/case" + j + ".py")
+    
+    f = open("RunTimeTests/MonopoleTests/time_data.dat")
+    for m in range(0, 17):
+        mTimes[k][m] = float(f.readline().split()[1])   
+
+meanMTimes = []
+for m in range(0,17):
+    testSum = 0
+    for j in range(0, itr):
+        testSum = testSum + mTimes[j][m]    
+    meanMTimes.append(testSum/itr)
+        
+f = open("RunTimeTests/MeanTimeData", 'x')
+f.write("Viscous\n")
+for i in range(0,10):
+    f.write(str(meanVTimes[i]) + "\n")
+f.write("Monopole\n")
+for i in range(0,17):
+    f.write(str(meanMTimes[i]) + "\n")
+f.close()

jobs.txt

@@ -0,0 +1,30 @@
+Viscous:
+ C79E1D3C
+ CD9D3050
+ 9CB03CEF
+ D6BAC936
+ FB822062
+ B3AAC9C8
+ 1C0780C8
+ 2060F55A
+ 07C33719
+ 939D6718
+
+Monopole:
+ EF54219C
+ 7FC6826B
+ 59D05DE9
+ 6B22A317
+ AF46C382
+ 46AA7AF8
+ 34DBFE14
+ 14B6198D
+ AE37D842
+ CC4F7C44
+ 66CFF8CC
+ 303B925A
+ AD63A4A5
+ 8A341282
+ 622DEC78
+ AB04C64D
+ 63850240

src/simulation/m_global_parameters.fpp

@@ -30,6 +30,8 @@ module m_global_parameters
 
     implicit none
 
+    real(kind(0d0)) :: time = 0
+
     ! Logistics ================================================================
     integer :: num_procs             !< Number of processors
     integer, parameter :: num_stcls_min = 5     !< Mininum # of stencils

src/simulation/m_time_steppers.f90

@@ -231,6 +231,8 @@ subroutine s_1st_order_tvd_rk(t_step, time_avg) ! ------------------------------
 
         call cpu_time(finish)
 
+        time = time + (finish - start)
+
         if (t_step >= 4) then
             time_avg = (abs(finish - start) + (t_step - 4)*time_avg)/(t_step - 3)
         else
@@ -313,6 +315,8 @@ subroutine s_2nd_order_tvd_rk(t_step, time_avg) ! ------------------------------
 
         call cpu_time(finish)
 
+        time = time + (finish - start)
+
         if (t_step >= 4) then
             time_avg = (abs(finish - start) + (t_step - 4)*time_avg)/(t_step - 3)
         else
@@ -422,6 +426,8 @@ subroutine s_3rd_order_tvd_rk(t_step, time_avg) ! ------------------------------
 
         call cpu_time(finish)
 
+        time = time + (finish - start)
+
         if (t_step >= 4) then
             time_avg = (abs(finish - start) + (t_step - 4)*time_avg)/(t_step - 3)
         else

src/simulation/m_viscous.f90

@@ -1376,8 +1376,8 @@ subroutine s_reconstruct_cell_boundary_values_visc(v_vf, vL_x_flat, vL_y_flat, v
     end subroutine s_reconstruct_cell_boundary_values_visc ! --------------------
 
     subroutine s_finalize_viscous_module()
-        deallocate (gammas, pi_infs)
-        deallocate (Res(1:2, 1:maxval(Re_size)))
+        !deallocate (gammas, pi_infs)
+        !deallocate (Res(1:2, 1:maxval(Re_size)))
     end subroutine s_finalize_viscous_module
 
 end module m_viscous

src/simulation/p_main.fpp

@@ -256,7 +256,7 @@ program p_main
             call s_mpi_barrier()
 
             if (num_procs > 1) then
-                call mpi_bcast_time_step_values(proc_time, time_avg)
+                call mpi_bcast_time_step_values(proc_time, time)
 
                 call mpi_bcast_time_step_values(io_proc_time, io_time_avg)
             end if