Télécharger ftemp.eso

Retour à la liste

Numérotation des lignes :

  1. C FTEMP SOURCE STRU 08/12/22 21:15:23 6246
  2. c Sous-Programme du module MISTRAL0
  3. C --------------------------------------------------------------------------
  4. SUBROUTINE FTEMP (TT, ID, PAR, V,VP)
  5. C --------------------------------------------------------------------------
  6. C Calcul de V, et de sa derivee VP si ID = 1,
  7. C en fonction de TT et des parametres transmis dans PAR.
  8. C --------------------------------------------------------------------------
  9. IMPLICIT INTEGER (I-N)
  10. IMPLICIT REAL*8 (A-H, O-Z)
  11. DIMENSION PAR(1:*)
  12. c (NPTTM = 5)
  13. L = NINT(PAR(1))
  14. IP = 1
  15. IF (L.EQ.1) THEN
  16. V = PAR(IP+1)+(PAR(IP+2)+PAR(IP+3)*TT)*TT
  17. IF (ID.EQ.1) THEN
  18. VP = PAR(IP+2)+2.*PAR(IP+3)*TT
  19. END IF
  20. ELSE IF (L.EQ.2) THEN
  21. V = PAR(IP+1)*EXP(PAR(IP+2)/TT+(PAR(IP+3)+PAR(IP+4)*TT)*TT)
  22. IF (ID.EQ.1) THEN
  23. VP = (-PAR(IP+2)/TT**2+PAR(IP+3)+2.*PAR(IP+4)*TT)*V
  24. END IF
  25. ELSE IF (L.EQ.3) THEN
  26. FT = 1./(1.+EXP(PAR(IP+3)*(TT/PAR(IP+4)-1.)))
  27. V = PAR(IP+1)*FT+PAR(IP+2)*(1.-FT)
  28. IF (ID.EQ.1) THEN
  29. FTP = -FT*(1.-FT)*PAR(IP+3)/PAR(IP+4)
  30. VP = (PAR(IP+1)-PAR(IP+2))*FTP
  31. END IF
  32. END IF
  33. RETURN
  34. END
  35.  
  36.  

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