!!  Copyright (C) Stichting Deltares, 2005-2018.
!!
!!  This file is part of iMOD.
!!
!!  This program is free software: you can redistribute it and/or modify
!!  it under the terms of the GNU General Public License as published by
!!  the Free Software Foundation, either version 3 of the License, or
!!  (at your option) any later version.
!!
!!  This program is distributed in the hope that it will be useful,
!!  but WITHOUT ANY WARRANTY; without even the implied warranty of
!!  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
!!  GNU General Public License for more details.
!!
!!  You should have received a copy of the GNU General Public License
!!  along with this program.  If not, see <http://www.gnu.org/licenses/>.
!!
!!  Contact: imod.support@deltares.nl   
!!  Stichting Deltares
!!  P.O. Box 177
!!  2600 MH Delft, The Netherlands.

module mod_gwf2ipsutl
use gwf2ipspar
contains 
 !##=================================    
 subroutine gwf2ipsar( )
 !##=================================
 return
 end subroutine gwf2ipsar

 !##=================================
 subroutine gwf2ipsfmwel( q, l )
 !##=================================
 ! Formulate the new value of the well
 use gwf2ipspar
 implicit none
 integer, intent(in)  :: l
 real, intent(out)    :: q

 q = q * fwell
 end subroutine gwf2ipsfmwel

 !##=================================
 subroutine gwf2ipsbd( buff )
 !##=================================
 ! Collect the cell-by-cell budget from mf2005
 use global, only: ncol,nrow,nlay
 implicit none
 double precision, intent(in) :: buff(nlay,nrow,ncol)
 integer :: ir, ic, il
 
 ipsbuff = buff
 end subroutine gwf2ipsbd

 !##=================================
 subroutine gwf2ipsda( )
 !##=================================
 ! Deallocate
 use gwf2ipspar
 implicit none
 deallocate( ipsbuff )
 return
 end subroutine gwf2ipsda
end module mod_gwf2ipsutl