diff --git a/src/mesh/FEM_quadrature.f90 b/src/mesh/FEM_quadrature.f90
index ed0bbef9a..3af535c2f 100644
--- a/src/mesh/FEM_quadrature.f90
+++ b/src/mesh/FEM_quadrature.f90
@@ -20,6 +20,10 @@ module FEM_quadrature
                            -1.0_pReal,  1.0_pReal, -1.0_pReal, &
                            -1.0_pReal, -1.0_pReal,  1.0_pReal], shape=[3,4])
 
+  type :: group_float                                                                               !< variable length datatype used for storage of state
+    real(pReal), dimension(:), pointer :: p
+  end type group_float
+
   integer,             dimension(2:3,maxOrder), public, protected :: &
     FEM_nQuadrature                                                                                 !< number of quadrature points for a given spatial dimension(2-3) and interpolation order(1-maxOrder)
   type(group_float),   dimension(2:3,maxOrder), public, protected :: &
@@ -39,6 +43,9 @@ subroutine FEM_quadrature_init
 
   print'(/,a)', ' <<<+-  FEM_quadrature init  -+>>>'; flush(6)
 
+  print*, 'L. Zhang et al., Journal of Computational Mathematics 27(1):89-96, 2009'
+  print*, 'https://www.jstor.org/stable/43693493'
+
 !--------------------------------------------------------------------------------------------------
 ! 2D linear
   FEM_nQuadrature(2,1) = 1
diff --git a/src/prec.f90 b/src/prec.f90
index 7613cfe46..733286a8a 100644
--- a/src/prec.f90
+++ b/src/prec.f90
@@ -27,11 +27,6 @@ module prec
 
   real(pReal), parameter :: tol_math_check = 1.0e-8_pReal                                           !< tolerance for internal math self-checks (rotation)
 
-
-  type :: group_float                                                                               !< variable length datatype used for storage of state
-    real(pReal), dimension(:), pointer :: p
-  end type group_float
-
   type :: tState
     integer :: &
       sizeState        = 0, &                                                                       !< size of state