Esope source files

   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 :