Télécharger la plaquette Cast3M


fr en

Télécharger flowh.eso

Retour à la liste

Numérotation des lignes : 

flowh

  1. C FLOWH     SOURCE    CB215821  16/04/21    21:16:53     8920
  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 = SQRT(FH)
  24.          FH = SQRT(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/(SQRT(3.D0))*TRS/FC
  31.  
  32.          IF (TRS .LE. 0.D0) THEN
  33.            KAPPA = -AH+SQRT(AH*AH-BH*KSIB+CH)
  34.          ELSE
  35.            KAPPA = -AH+SQRT(AH*AH+CH)
  36. *           KAPPA = 0.00002
  37.          END IF
  38.  
  39.          FH = FH / KAPPA
  40.  
  41.         END IF
  42.  
  43.         END
  44.  
  45.  
  46.  
  47.  
  48.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales