removed some unused code and added headers to modules

Author: wilfonba

src/simulation/m_monopole.f90

@@ -1,3 +1,9 @@
+!>
+!! @file m_viscous.f90
+!! @brief Contains module m_viscous
+
+!> @brief The module contains the subroutines used to
+
 module m_monopole
 
     ! Dependencies =============================================================

src/simulation/m_rhs.f90

@@ -5,7 +5,7 @@
 !> @brief The module contains the subroutines used to calculate the right-
 !!              hand-side (RHS) in the quasi-conservative, shock- and interface-
 !!              capturing finite-volume framework for the multicomponent Navier-
-!!       f       Stokes equations supplemented by appropriate advection equations
+!!              Stokes equations supplemented by appropriate advection equations
 !!              used to capture the material interfaces. The system of equations
 !!              is closed by the stiffened gas equation of state, as well as any
 !!              required mixture relationships. Capillarity effects are included

src/simulation/m_viscous.f90

@@ -1,3 +1,10 @@
+!>
+!! @file m_viscous.f90
+!! @brief Contains module m_viscous
+
+!> @brief The module contains the subroutines used to
+
+
 module m_viscous
 
     ! Dependencies =============================================================
@@ -10,7 +17,8 @@ module m_viscous
 
     private; public  s_get_viscous, &
     s_compute_viscous_stress_tensor, &
-    s_initialize_viscous_module
+    s_initialize_viscous_module, &
+    s_finalize_viscous_module
 
     type(int_bounds_info) :: iv
 
@@ -36,19 +44,14 @@ subroutine s_initialize_viscous_module()
         end do
 !$acc update device(gammas, pi_infs)
 
-        if (any(Re_size > 0)) then
-            allocate (Res(1:2, 1:maxval(Re_size)))
-        end if
+        allocate (Res(1:2, 1:maxval(Re_size)))
 
-        if (any(Re_size > 0)) then
-            do i = 1, 2
-                do j = 1, Re_size(i)
-                    Res(i, j) = fluid_pp(Re_idx(i, j))%Re(i)
-                end do
+        do i = 1, 2
+            do j = 1, Re_size(i)
+                Res(i, j) = fluid_pp(Re_idx(i, j))%Re(i)
             end do
+        end do
 !$acc update device(Res, Re_idx, Re_size)
-        end if
-
 
         momxb = mom_idx%beg
         momxe = mom_idx%end
@@ -85,27 +88,6 @@ subroutine s_compute_viscous_stress_tensor(q_prim_vf, grad_x_vf, grad_y_vf, grad
 
         type(int_bounds_info) :: ix, iy, iz
 
-        !ix%beg = -buff_size; iy%beg = 0; iz%beg = 0
-        !if (n > 0) iy%beg = -buff_size; if (p > 0) iz%beg = -buff_size
-        !ix%end = m - ix%beg; iy%end = n - iy%beg; iz%end = p - iz%beg
-
-        !$acc update device(ix, iy, iz)
-
-        do i = 1, num_dims
-            allocate (tau_Re_vf(cont_idx%end + i)%sf(ix%beg:ix%end, &
-                                                     iy%beg:iy%end, &
-                                                     iz%beg:iz%end))
-!$acc enter data create(tau_Re_vf(cont_idx%end + i)%sf(ix%beg:ix%end, &
-!$acc                                                  iy%beg:iy%end, &
-!$acc                                                  iz%beg:iz%end))
-        end do
-        allocate (tau_Re_vf(E_idx)%sf(ix%beg:ix%end, &
-                                      iy%beg:iy%end, &
-                                      iz%beg:iz%end))
-!$acc enter data create (tau_Re_vf(E_idx)%sf(ix%beg:ix%end, &
-!$acc                                        iy%beg:iy%end, &
-!$acc                                        iz%beg:iz%end))
-
     !$acc parallel loop collapse(3) gang vector default(present)
         do l = iz%beg, iz%end
             do k = iy%beg, iy%end
@@ -1390,4 +1372,9 @@ 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)
+        allocate (Res(1:2, 1:maxval(Re_size)))
+    end subroutine s_finalize_viscous_module
+
 end module m_viscous

src/simulation/p_main.fpp

@@ -356,6 +356,10 @@ program p_main
     call s_finalize_mpi_proxy_module()
     call s_finalize_global_parameters_module()
 
+    if (any(Re_size > 0)) then
+        call s_finalize_viscous_module()
+    end if
+
     ! Terminating MPI execution environment
     call s_mpi_finalize()