diff --git a/module_definitions.f90 b/module_definitions.f90
index c815102ac012b8d8ce6ff59635414887730452cb..ed28c2c6d089db631a0ec6d6ce26831b149894b0 100644
--- a/module_definitions.f90
+++ b/module_definitions.f90
@@ -285,33 +285,39 @@ module definitions
 ! following is a general interface to sort an array of numbers
 
       interface qsort
+
         subroutine iqsort (array,n,perm)
         integer,intent(in) :: n
-        integer,dimension(n),intent(out)::array
-        integer,dimension(n),intent(out)::perm
+        integer,dimension(n),intent(inout)::array
+        integer,dimension(n),intent(inout)::perm
         end subroutine iqsort
+
         subroutine rqsort (array,n,perm)
         integer,intent(in) :: n
-        real,dimension(n),intent(out)::array
-        integer,dimension(n),intent(out)::perm
+        real,dimension(n),intent(inout)::array
+        integer,dimension(n),intent(inout)::perm
         end subroutine rqsort
+
         subroutine dpqsort (array,n,perm)
         integer,intent(in) :: n
-        real*8,dimension(n),intent(out)::array
-        integer,dimension(n),intent(out)::perm
+        real*8,dimension(n),intent(inout)::array
+        integer,dimension(n),intent(inout)::perm
         end subroutine dpqsort
+
         subroutine iqsort_s (array,n)
         integer,intent(in) :: n
-        integer,dimension(n),intent(out)::array
+        integer,dimension(n),intent(inout)::array
         end subroutine iqsort_s
+
         subroutine rqsort_s (array,n)
         integer,intent(in) :: n
-        real,dimension(n),intent(out)::array
+        real,dimension(n),intent(inout)::array
         end subroutine rqsort_s
-        subroutine dpqsort_s (array,n)
-        integer,intent(in) :: n
-        real*8,dimension(n),intent(out)::array
-        end subroutine dpqsort_s
+!        subroutine dpqsort_s (array,n)
+!        integer,intent(in) :: n
+!        real*8,dimension(n),intent(inout)::array
+!        end subroutine dpqsort_s
+
       end interface