Télécharger kaltet.eso

Retour à la liste

Numérotation des lignes :

  1. C KALTET SOURCE CHAT 05/01/13 00:51:01 5004
  2. SUBROUTINE KALTET(TPLUS,PRANDT)
  3. C
  4. C PROGRAMMA PER LA TABULAZIONE DI TPLUS
  5. C PER UN ASSEGNATO VALORE DI PRANDTL
  6. C
  7. IMPLICIT INTEGER(I-N)
  8. IMPLICIT REAL*8 (A-H,O-Z)
  9. DIMENSION TPLUS(101)
  10. APLUS=26.0
  11. BPLUS=34.96+28.79*(LOG10(PRANDT))+33.95*(LOG10(PRANDT))**2
  12. *+6.33*(LOG10(PRANDT))**3-1.186*(LOG10(PRANDT))**4
  13. TPLUS(1)=0.0D0
  14. YPLUS=0.0D0
  15. ZL=ZLM(YPLUS,APLUS)
  16. UP=UPRIMO(ZL)
  17. ST=SIGMAT(PRANDT,YPLUS,APLUS,BPLUS)
  18. ZETAV=PRANDT/(1.0+ZL*UP*PRANDT/ST)
  19. DO 10 I = 2,101
  20. YPLUS=YPLUS+1.0D0
  21. ZL=ZLM(YPLUS,APLUS)
  22. UP=UPRIMO(ZL)
  23. ST=SIGMAT(PRANDT,YPLUS,APLUS,BPLUS)
  24. ZETAN=PRANDT/(1.0+ZL*UP*PRANDT/ST)
  25. TPLUS(I)=TPLUS(I-1)+(ZETAV+ZETAN)/2.0
  26. ZETAV=ZETAN
  27. 10 CONTINUE
  28. RETURN
  29. END
  30.  
  31.  
  32.  
  33.  

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