anchor.diff:          do j = 1, pasy
               b(1,j) = b(pasx,j)
               b(pasx+1,j) = b(2,j)
  +
               q(1,j) = q(pasx,j)
               q(pasx+1,j) = q(2,j)
           end do    
  
positive.code: !   . Calcul des vitesses
  !   . Sauvegarde tous les X pas de temps
  ! ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  
  module GEOPHY
  
  
negative.code:     allocate( v2(pasx+1, pasy) )
  
      ! Calcul du profil initial
      call init_xy(x, y, dx, dy, Ly)
  
      temp1 = g*h0/f/L/U
      do i = 1, pasx+1
  