From de7ef430958183843867ea2b42430f16919af16f Mon Sep 17 00:00:00 2001
From: Martin Diehl <mail@martin-diehl.net>
Date: Wed, 21 Jul 2021 23:43:55 +0200
Subject: [PATCH] documenting/silencing gfortran

---
 src/mesh/FEM_quadrature.f90      | 6 +++---
 src/mesh/discretization_mesh.f90 | 6 ++----
 src/quit.f90                     | 5 ++++-
 3 files changed, 9 insertions(+), 8 deletions(-)

diff --git a/src/mesh/FEM_quadrature.f90 b/src/mesh/FEM_quadrature.f90
index bf1f1567f..0c470ca6b 100644
--- a/src/mesh/FEM_quadrature.f90
+++ b/src/mesh/FEM_quadrature.f90
@@ -124,7 +124,7 @@ subroutine FEM_quadrature_init()
   allocate(FEM_quadrature_weights(3,2)%p(FEM_nQuadrature(3,2)))
   FEM_quadrature_weights(3,2)%p(1:4) = 0.25_pReal
 
-  FEM_quadrature_points (3,2)%p = permutationStar31([0.1381966011250105_pReal])
+  FEM_quadrature_points (3,2)%p = permutationStar31([1.3819660112501052e-1_pReal])
 
 !--------------------------------------------------------------------------------------------------
 ! 3D cubic
@@ -141,7 +141,7 @@ subroutine FEM_quadrature_init()
     permutationStar22([4.5503704125649649e-2_pReal]) ]
 
 !--------------------------------------------------------------------------------------------------
-! 3D quartic
+! 3D quartic (lower precision/unknown source)
   FEM_nQuadrature(3,4) = 35
 
   allocate(FEM_quadrature_weights(3,4)%p(FEM_nQuadrature(3,4)))
@@ -159,7 +159,7 @@ subroutine FEM_quadrature_init()
     permutationStar4([0.25_pReal]) ]
 
 !--------------------------------------------------------------------------------------------------
-! 3D quintic
+! 3D quintic (lower precision/unknown source)
   FEM_nQuadrature(3,5) = 56
 
   allocate(FEM_quadrature_weights(3,5)%p(FEM_nQuadrature(3,5)))
diff --git a/src/mesh/discretization_mesh.f90 b/src/mesh/discretization_mesh.f90
index 49eaf0e18..422f4ca85 100644
--- a/src/mesh/discretization_mesh.f90
+++ b/src/mesh/discretization_mesh.f90
@@ -71,16 +71,14 @@ subroutine discretization_mesh_init(restart)
 
   logical, intent(in) :: restart
 
-  integer, allocatable, dimension(:) :: chunkPos
   integer :: dimPlex, &
     mesh_Nnodes, &                                                                                  !< total number of nodes in mesh
-    j, l, &
+    j, &
     debug_element, debug_ip
   PetscSF :: sf
   DM :: globalMesh
   PetscInt :: nFaceSets
   PetscInt, pointer, dimension(:) :: pFaceSets
-  character(len=pStringLen), dimension(:), allocatable :: fileContent
   IS :: faceSetIS
   PetscErrorCode :: ierr
   integer, dimension(:), allocatable :: &
@@ -88,7 +86,7 @@ subroutine discretization_mesh_init(restart)
   class(tNode), pointer :: &
     num_mesh
   integer :: integrationOrder                                                                       !< order of quadrature rule required
-  type(tvec)  :: coords_node0
+  type(tvec) :: coords_node0
 
   print'(/,a)',   ' <<<+-  discretization_mesh init  -+>>>'
 
diff --git a/src/quit.f90 b/src/quit.f90
index 786241d85..6361ff783 100644
--- a/src/quit.f90
+++ b/src/quit.f90
@@ -6,7 +6,10 @@
 !--------------------------------------------------------------------------------------------------
 subroutine quit(stop_id)
 #include <petsc/finclude/petscsys.h>
-  use PetscSys
+  use PETScSys
+#if (PETSC_VERSION_MAJOR==3 && PETSC_VERSION_MINOR>14) && !defined(PETSC_HAVE_MPI_F90MODULE_VISIBILITY)
+  use MPI_f08
+#endif
   use HDF5
 
   implicit none