Sources Esope | Cast3M

ftemp
C FTEMP     SOURCE    STRU      08/12/22    21:15:23     6246
c     Sous-Programme du module MISTRAL0
C     --------------------------------------------------------------------------
      SUBROUTINE FTEMP (TT, ID, PAR, V,VP)
C     --------------------------------------------------------------------------
C     Calcul de V, et de sa derivee VP si ID = 1,
C     en fonction de TT et des parametres transmis dans PAR.
C     --------------------------------------------------------------------------
      IMPLICIT INTEGER (I-N)
      IMPLICIT REAL*8 (A-H, O-Z)
      DIMENSION PAR(1:*)
c     (NPTTM = 5)
      L = NINT(PAR(1))
      IP = 1
        IF (L.EQ.1) THEN
         V = PAR(IP+1)+(PAR(IP+2)+PAR(IP+3)*TT)*TT
       IF (ID.EQ.1) THEN
          VP = PAR(IP+2)+2.*PAR(IP+3)*TT
         END IF
        ELSE IF (L.EQ.2) THEN
         V = PAR(IP+1)*EXP(PAR(IP+2)/TT+(PAR(IP+3)+PAR(IP+4)*TT)*TT)
       IF (ID.EQ.1) THEN
          VP = (-PAR(IP+2)/TT**2+PAR(IP+3)+2.*PAR(IP+4)*TT)*V
         END IF
        ELSE IF (L.EQ.3) THEN
       FT = 1./(1.+EXP(PAR(IP+3)*(TT/PAR(IP+4)-1.)))
       V = PAR(IP+1)*FT+PAR(IP+2)*(1.-FT)
       IF (ID.EQ.1) THEN
        FTP = -FT*(1.-FT)*PAR(IP+3)/PAR(IP+4)
          VP = (PAR(IP+1)-PAR(IP+2))*FTP
         END IF
        END IF
        RETURN
        END