source: branches/GRISLIv3/SOURCES/ablation_bord.f90 @ 467

!> \file ablation_bord.f90
!! Clacul de l'ablation sur le bord
!<
      

!> SUBROUTINE: Ablation_bord
!! Calcule l'ablation pour les points du bord 
!! \author ...
!! \date ....
!! @note cas des bords de calotte : un point du bord peut perdre en plus la glace 
!! amenée par l'écoulement.
!! @note Calcule l'ablation pour les points du bord avec H=0 mais apport par l'amont
!! diagnostic pour la fermeture du bilan d'eau
!! par convention ablation negative dans la zone d'ablation 
!! @note Used modules:
!! @note  - use module3D_phy
!>

subroutine ablation_bord


  use module3D_phy, only: ice,ablbord,bmelt,bm,dt,ablbord_dtt
  use geography, only: nx,ny
  implicit none
  integer :: ip1,im1,jp1,jm1   ! pour chercher les voisins
  integer :: somm
  integer :: i,j

  ! Dans les zones ou le calcul icethick a donne un H negatif : ablbord=H

  ! pour être au bord, il suffit que la somme des ice voisins soit > 0
  do j=1,ny
     do i=1,nx
        if (ice(i,j).eq.0) then   ! seulment pour les noeuds sans glace
           im1=max(i-1,1)
           jm1=max(j-1,1)
           ip1=min(i+1,nx)
           jp1=min(j+1,ny)

           somm=ice(ip1,j)+ice(im1,j)+ice(i,jp1)+ice(i,jm1)
           if (somm.gt.0.and.(ablbord(i,j).LT.0.)) then                               ! voisins englaces
              !~ !         if (somm.gt.0.and.((ablbord(i,j).LT.0.).OR.(flot(ip1,j).OR.flot(im1,j).OR.flot(i,jp1).OR.flot(i,jm1)))) then   ! voisins englaces
              !~ !            ablbord(i,j)= ablbord(i,j) - (Acc(i,j) * dt)     !  Ablvrai=bmpdd-Acc-Hcalc/dt
              !~             if (ablbord(i,j).LT.0..and.front(i,j).LE.4.and.front(i,j).GE.0) then
              ablbord(i,j)= max(0.,ablbord(i,j) + ((Bmelt(i,j) - Bm(i,j)) * dt))  ! version Tof : ablbord = Hcons_mass - Bm - Bmelt

           else
              !~ ! version cat            ablbord(i,j)=hdot(i,j)-acc(i,j)   ! on ne fond que ce qui est tombe + existant
              ablbord(i,j)=0.
           endif
        else
           ablbord(i,j)=0.
        end if
     end do
  end do

  ablbord_dtt(:,:) = ablbord_dtt(:,:) + ablbord(:,:)

  !cdc deplace dans bilan_flux_output_mod ablbord(:,:)=0.

end subroutine ablation_bord