!
! $Id $
!
SUBROUTINE criterion_ice_stop_slope(cell_area,ini_surf,qsurf,STOPPING,ngrid,initial_h2o_ice,global_ave_press_GCM,global_ave_press_new,nslope)

  USE temps_mod_evol, ONLY: alpha_criterion
  use comslope_mod, ONLY: subslope_dist

      IMPLICIT NONE

!=======================================================================
!
!  Routine that checks if the criterion to stop the PEM is reached
!
!=======================================================================

!   arguments:
!   ----------

!   INPUT
  INTEGER, intent(in) :: ngrid,nslope                  ! # of grid physical grid points 
  REAL,    intent(in) :: cell_area(ngrid)       ! physical point field : Area of the cells
  REAL,    intent(in) ::  qsurf(ngrid,nslope)          ! physical point field : Actual density of water ice
  REAL,    intent(in) :: ini_surf
  REAL,    intent(in) :: initial_h2o_ice(ngrid,nslope)
  REAL,    intent(in) :: global_ave_press_GCM
  REAL,    intent(in) :: global_ave_press_new


!   OUTPUT
  LOGICAL, intent(out) :: STOPPING              ! Logical : is the criterion reached?

!   local:
!   -----
  INTEGER :: i,islope                    ! Loop
  REAL :: present_surf  ! Initial/Actual surface of water ice

!=======================================================================

!   initialisation to false
    STOPPING=.FALSE.

!   computation of the actual surface
  present_surf=0.
  do i=1,ngrid
   do islope=1,nslope
      if (initial_h2o_ice(i,islope).GT.0.5 .and. qsurf(i,islope).GT.0.) then
         print *, "i", i
         print *, "initial_h2o_ice(i,islope)", initial_h2o_ice(i,islope)
         print *, "qsurf(i,:)", qsurf(i,:)
         print *, "cell_area(i)", cell_area(i)
         print *, "present_surf",present_surf
         present_surf=present_surf+cell_area(i)*subslope_dist(i,islope)
      endif
   enddo
  enddo

!  print *, "initial_h2o_ice", initial_h2o_ice
!  print *, "qsurf", qsurf

  print *, "present_surf", present_surf
  print *, "ini_surf", ini_surf
  print *, "ini_surf*0.8", ini_surf*(1-alpha_criterion)
  
!   check of the criterion
  if(present_surf.LT.ini_surf*(1-alpha_criterion) .OR. &
     present_surf.GT.ini_surf*(1+alpha_criterion)) then
  STOPPING=.TRUE. 
  endif

  if (ini_surf.LT. 1E-5 .and. ini_surf.GT. -1E-5) then
       STOPPING=.FALSE.
  endif

!  if(global_ave_press_GCM.LT.global_ave_press_new*(1-alpha_criterion) .OR. &
!     global_ave_press_GCM.GT.global_ave_press_new*(1+alpha_criterion)) then
!  STOPPING=.TRUE. 
!  endif

  if(global_ave_press_new.LT.global_ave_press_GCM*(0.9) .OR. &
     global_ave_press_new.GT.global_ave_press_GCM*(1.1)) then
  STOPPING=.TRUE. 
  endif

END SUBROUTINE criterion_ice_stop_slope