Numérotation des lignes :

fvdi
C FVDI      SOURCE    STRU      08/12/22    21:15:24     6246c     Sous-Programme du module MISTRAL0C     --------------------------------------------------------------------------      SUBROUTINE FVDI (PSI,TT, ID, PAR, V,DVDPSI,DVDTT)C     --------------------------------------------------------------------------C     Calcul de V, et de ses derivees DVDPSI et DVDTT si ID = 1,C     en fonction de PSI, TT et des parametres transmis dans PAR.C     --------------------------------------------------------------------------      IMPLICIT INTEGER (I-N)      IMPLICIT REAL*8 (A-H, O-Z)      DIMENSION PAR(1:*)c     (NPPSIM = 6)      L = NINT(PAR(1))      IP = 1        IF (L.EQ.1) THEN         V = LOG(1.+PAR(IP+1)*PSI)       IF (ID.EQ.1) THEN          DVDPSI = PAR(IP+1)/(1.+PAR(IP+1)*PSI)        DVDTT = 0.         END IF        ELSE IF (L.EQ.2) THEN         V = 1.-EXP(-PAR(IP+1)*PSI)       IF (ID.EQ.1) THEN          DVDPSI = PAR(IP+1)*EXP(-PAR(IP+1)*PSI)        DVDTT = 0.         END IF        ELSE IF (L.EQ.3) THEN         V = EXP(PAR(IP+1)*PSI**PAR(IP+2))       IF (ID.EQ.1) THEN          DVDPSI = PAR(IP+1)*PAR(IP+2)*PSI**(PAR(IP+2)-1.)*V        DVDTT = 0.         END IF        ELSE IF (L.EQ.4) THEN         VP = PAR(IP+1)*PSI**PAR(IP+2)         VPT = PAR(IP+3)*EXP(PAR(IP+4)*PSI**PAR(IP+5))/TT         V = EXP(VP+VPT)       IF (ID.EQ.1) THEN        DVPDP = PAR(IP+1)*PAR(IP+2)*PSI**(PAR(IP+2)-1.)        DVPTDP = PAR(IP+4)*PAR(IP+5)*PSI**(PAR(IP+5)-1.)*VPT          DVDPSI = (DVPDP+DVPTDP)*V        DVDTT = -VPT/TT*V         END IF        END IF        RETURN        END  

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