! 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