MODULE evol_co2_ice_s_mod

implicit none

!=======================================================================
contains
!=======================================================================

SUBROUTINE evol_co2_ice_s(qsurf,tendencies_co2_ice_phys,iim_input,jjm_input,ngrid,cell_area,nslope)

use time_evol_mod, only: dt_pem

implicit none

!=======================================================================
!
!  Routine that compute the evolution of the CO2 ice
!
!=======================================================================
!   arguments:
!   ----------
!   INPUT
integer,                intent(in) :: iim_input, jjm_input, ngrid, nslope ! # of grid points along longitude/latitude/ total
REAL, dimension(ngrid), intent(in) :: cell_area

!   OUTPUT
real, dimension(ngrid,nslope), intent(inout) :: qsurf                   ! physical point field: Previous and actual density of water ice
real, dimension(ngrid,nslope), intent(inout) :: tendencies_co2_ice_phys ! physical point field: Evolution of perennial ice over one year

!   local:
!   ----
integer :: i, j, ig0, islope ! loop variable
!=======================================================================
! Evolution of the CO2 ice for each physical point
do i=1,ngrid
    do islope=1,nslope
        qsurf(i,islope) = qsurf(i,islope) + tendencies_co2_ice_phys(i,islope)*dt_pem
        if (qsurf(i,islope) < 0) then
            qsurf(i,islope) = 0.
            tendencies_co2_ice_phys(i,islope) = 0.
        endif
    enddo
enddo

END SUBROUTINE evol_co2_ice_s

END MODULE evol_co2_ice_s_mod