anchor.diff:      kmoy  = 0.
       save  = 1
   
  -    write(*,*) 'Ly/L =',Ly,'(',(3./4.+g*h0/(2.*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*Lx*L,'mtres de long'
  -    write(*,*) '           et ',2*Ly*L,'mtres de large'
  +    write(*,*) 'Ly/L =', Ly, '(', (3./4. + g*h0/(2.*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*Lx*L, 'mtres de long'
  +    write(*,*) '           et ', 2*Ly*L, 'mtres de large'
       write(*,*) 'tol (1D-6)?'
       read(*,*) tol
       write(*,*) 'w (1.825)?'
  
positive.code:     endif
  
      write(*,*)
      write(*,*)'Un pas de temps =',T*dt,'sec.'
  
      ! ------------------------------------------------------------------      
      !                      Calcul du profil initial
  
negative.code:     !            Boucle principale (progression dans le temps)
      ! ------------------------------------------------------------------
        
      do n = 1, Nmax
          write(*,*) 'Pas de temps n :', n
  
      !        ------------------------------------
      !        -- Calcul de Beta et du laplacien --
  