! Type your code here, or load an example.
real function square(rhogw_vm)
    implicit none
    integer, parameter :: ADM_kmin = 2
    integer, parameter :: ADM_kmax = 41
    integer, parameter :: ADM_gall = 1156
    integer, parameter :: ADM_kall = 42
    integer, parameter :: ADM_lall = 2
    real(8), intent(inout) :: rhogw_vm (ADM_gall, ADM_kall, ADM_lall)
    integer :: i, l

    DO i = 1, 10000
      do l = 1, ADM_lall
        rhogw_vm(:,ADM_kmin  ,l) = 0.D0
        rhogw_vm(:,ADM_kmax+1,l) = 0.D0
      enddo
    enddo

    return
end function square