C DYFVIS    SOURCE    CHAT      05/01/12    23:09:47     5004           
      SUBROUTINE DYFVIS(XDH,XDEP,XVIT,XJEU,XCV,XFV,iannul)                             
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)                                                  
*--------------------------------------------------------------------*          
*                                                                    *          
*     Op{rateur DYNE : algorithme de Fu - de Vogelaere               *          
*     ________________________________________________               *          
*                                                                    *          
*     Calcul de la force de viscosit{.                               *          
*                                                                    *          
*     Param}tres:                                                    *          
*                                                                    *          
* e   XDEP    Valeur du d{placement                                  *          
* e   XVIT    Valeur de la vitesse                                   *          
* e   XJEU    Valeur du jeu                                          *          
* e   XCV     Valeur du coefficient de viscosit{                     *          
*  s  XFV     Valeur de la force de viscosit{                        *          
*                                                                    *          
*     Auteur, date de cr{ation:                                      *          
*                                                                    *          
*     Lionel VIVAN, le 20 mars 1991.                                 *          
*                                                                    *          
*--------------------------------------------------------------------*          
*     
      IF (iannul.EQ.0) THEN                                                                          
         XDH3 = XDH * XDH * XDH                                                    
         XFV = -XCV * XVIT / XDH3 
      ELSE
         XFV = 0.D0
      ENDIF                                                 
*                                                                               
      END