Numérotation des lignes :

ftemp
C FTEMP     SOURCE    STRU      08/12/22    21:15:23     6246c     Sous-Programme du module MISTRAL0C     --------------------------------------------------------------------------      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  

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