C FLOWH     SOURCE    PV        07/10/04    11:38:12     5887           
       SUBROUTINE FLOWH(SIGE0,KH,AH,BH,CH,FC,FM,FH)

* Subroutine de calcul de h

*       IMPLICIT NONE

       REAL*8           SIGE0(6), KH, FH, FM(6)
       REAL*8           TRS, KAPPA, AH, BH, CH, KSIB, FC
       INTEGER I


       FH = 0.D0
       IF (KH .GE. 1.D0) THEN
         FH = 0.D0
       ELSE
         DO I = 1,3
           FH = FH + FM(I)*FM(I)
         ENDDO
         DO I = 4,6
           FH = FH + FM(I)*FM(I)
         ENDDO
         FH = DSQRT(FH)
         FH = DSQRT(2.D0/3.D0)*FH

         TRS = 0.D0
         DO I = 1,3
           TRS = TRS + SIGE0(I)
         END DO
         KSIB = 1.D0/(DSQRT(3.D0))*TRS/FC

         IF (TRS .LE. 0.D0) THEN
           KAPPA = -AH+DSQRT(AH*AH-BH*KSIB+CH)
         ELSE
           KAPPA = -AH+DSQRT(AH*AH+CH)
*           KAPPA = 0.00002
         END IF
         
         FH = FH / KAPPA
        
        END IF          
        
        END