anchor.diff:      dt    = dx*dx/2.0D0
       dy    = 2.0D0*Ly/(pasy-1)
       kmoy  = 0.0D0
  -    save  = 1
  -
  -    write(*,*) 'Ly/L =', Ly, '(', (3.0D0/4.0D0 + g*h0/(2.0D0*f*U*L)), ')'
  -    write(*,*) 'Bu   =', Bu,'Ro   =', Ro
  -    write(*,*) 'dx   =', dx, ' ; dy   =', dy 
  -    write(*,*) 'pasx =', pasx,' ; pasy =', pasy
  -    write(*,*) 'T    =', T, 'dH   =', f*L*U/g
  -    write(*,*) 'dt (max) =', dt, 'sec.' 
  -    write(*,*) 'Le canal fait ', 2.0D0*Lx*L, 'mtres de long'
  -    write(*,*) '           et ', 2.0D0*Ly*L, 'mtres de large'
  -    write(*,*) 'tol (1D-6)?'
  +    nsave  = 1
  +
  +    print *, 'Ly/L =', Ly, '(', (3.0D0/4.0D0 + g*h0/(2.0D0*f*U*L)), ')'
  +    print *, 'Bu   =', Bu, ' ; Ro   =', Ro
  +    print *, 'dx   =', dx, ' ; dy   =', dy 
  +    print *, 'pasx =', pasx,' ; pasy =', pasy
  +    print *, 'T    =', T, ' ; dH   =', f*L*U/g
  +    print *, 'dt (max) =', dt, 'sec.' 
  +    print *, 'Le canal fait ', 2.0D0*Lx*L, 'm de long'
  +    print *, '           et ', 2.0D0*Ly*L, 'm de large'
  +    print *, 'tol (1D-6)?'
       read(*,*) tol
  -    write(*,*) 'w (1.825)?'
  +    print *, 'w (1.825)?'
       read(*,*) w
  -    write(*,*) 'dt (2D-2)?'
  +    print *, 'dt (2D-2)?'
       read(*,*) dt
  -    write(*,*) 'Nmax (100)?'
  +    print *, 'Nmax (100)?'
       read(*,*) Nmax
  -    write(*,*) 'fsave (Nmax)?'
  +    print *, 'fsave (Nmax)?'
       read(*,*) fsave
  -    write(*,*) 'save U (1 = oui)?'
  +    print *, 'save U (1 = oui)?'
       read(*,*) saveu
   
  -    if(fsave.GT.Nmax) then
  +    if(fsave > Nmax) then
           fsave = Nmax
       endif
   
       write(*,*)
  -    write(*,*) 'Un pas de temps =', T*dt, 'sec.'
  +    print *, 'Un pas de temps =', T*dt, 'sec.'
   
       ! Calcul du profil initial
         
  
positive.code: 
      !  Commentaires finaux
  
      write(*,*) 'Calcul effectu sur ', Nmax*T*dt, 'sec.'
      write(*,*) 'Nb itrations (moy) ', kmoy/Nmax
      write(*,*) 'Rsultats dans RES.M'
        
  end program NIHOUL
  
negative.code: 
          end do  
         
          print *, '    no it. =', compt
          kmoy = kmoy + compt
  
          ! Copie de eta (futur calcul des vitesses)
  
          if(n == nsave*fsave-1) then
              do i = 1, pasx+1
                  do j = 1, pasy
                      leta(i,j) = eta(i,j)
                  end do
  