C FLOWH     SOURCE    CB215821  16/04/21    21:16:53     8920
       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 = SQRT(FH)
         FH = SQRT(2.D0/3.D0)*FH

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

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

         FH = FH / KAPPA

        END IF

        END