Download flowh.eso

Back to the list

   1 : C FLOWH     SOURCE    PV        07/10/04    11:38:12     5887           
   2 :        SUBROUTINE FLOWH(SIGE0,KH,AH,BH,CH,FC,FM,FH)
   3 : 
   4 : * Subroutine de calcul de h
   5 : 
   6 : *       IMPLICIT NONE
   7 : 
   8 :        REAL*8           SIGE0(6), KH, FH, FM(6)
   9 :        REAL*8           TRS, KAPPA, AH, BH, CH, KSIB, FC
  10 :        INTEGER I
  11 : 
  12 : 
  13 :        FH = 0.D0
  14 :        IF (KH .GE. 1.D0) THEN
  15 :          FH = 0.D0
  16 :        ELSE
  17 :          DO I = 1,3
  18 :            FH = FH + FM(I)*FM(I)
  19 :          ENDDO
  20 :          DO I = 4,6
  21 :            FH = FH + FM(I)*FM(I)
  22 :          ENDDO
  23 :          FH = DSQRT(FH)
  24 :          FH = DSQRT(2.D0/3.D0)*FH
  25 : 
  26 :          TRS = 0.D0
  27 :          DO I = 1,3
  28 :            TRS = TRS + SIGE0(I)
  29 :          END DO
  30 :          KSIB = 1.D0/(DSQRT(3.D0))*TRS/FC
  31 : 
  32 :          IF (TRS .LE. 0.D0) THEN
  33 :            KAPPA = -AH+DSQRT(AH*AH-BH*KSIB+CH)
  34 :          ELSE
  35 :            KAPPA = -AH+DSQRT(AH*AH+CH)
  36 : *           KAPPA = 0.00002
  37 :          END IF
  38 :          
  39 :          FH = FH / KAPPA
  40 :         
  41 :         END IF          
  42 :         
  43 :         END 
  44 :  
  45 :  

© Cast3M 2003 - All rights reserved.
Disclaimer