Télécharger flowh.eso

Retour à la liste

Numérotation des lignes :

  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