Numérotation des lignes :

clebl4
C CLEBL4    SOURCE    PV        17/12/08    21:15:49     9660           C====================================================================C      SUBROUTINE CLEBL4(iwrk52,VWRK1,DTPS,VWRK2)CC====================================================================C     LB, 06/97C     =========C         Modele de Leblond et DevauxC         Un seul chgt de phase - cas du refroidissementC     Pour le moment zeq et tau sont calcules en localC      IMPLICIT INTEGER(I-N)       IMPLICIT REAL*8(A-H,O-Z)-INC DECHEC      REAL*8 VWRK1(*),VWRK2(7)C      T0  =VWRK1(1)      TP0 =VWRK1(2)      ZA0 =VWRK1(3)      ZF0 =VWRK1(4)      VMS1=VWRK1(7)        wrk52 = iwrk52      AC1 = xmat(1)      AC2 = xmat(2)      AC3 = xmat(3)      AC4 = xmat(4)      TAU1 = xmat(5)      TAU2 = xmat(6)      TAU3 = xmat(7)      TAU4 = xmat(8)CC     Calcul de ZeqCc       if (T0.LE.600.) thenc          zeq = 1.c       else if ((T0.GT.600.).AND.(T0.LT.800.)) thenc          zeq = (800. - T0) / 200.c       else if (T0.GE.800.) thenc          zeq = 0.c       endif       if (T0.LE.AC1) then          zeq = 1.       else if ((T0.GT.AC1).AND.(T0.LT.AC2).and.(AC3.ne.0.)) then          zeq = (AC2 - T0) / AC3       else if (T0.GE.AC2) then          zeq = 0.       endifCC      Calcul de tauCc       if (T0.LE.600.) thenc          tau = 100000.c       else if ((T0.GT.600.).AND.(T0.LT.620.)) thenc         tau = 100000. - (T0 - 600.) * 95000. / 20.c       else if (T0.GE.620.) thenc          tau = 500.c       endif       if (T0.LE.AC1) then          tau = TAU1       else if ((T0.GT.AC1).AND.(T0.LT.(AC1+AC4)).and.(AC4.ne.0.)) then          tau = TAU1 - (T0 - AC1) * TAU2 / AC4       else if (T0.GE.620.) then          tau = TAU3       endifCC     Calcul de zpointC      zpoint = (zeq - ZF0) / tauC        VWRK1(1)=T0        VWRK2(2)=TP0        VWRK2(3)=VWRK1(3) - zpoint * DTPS        VWRK2(4)=VWRK1(4) + zpoint * DTPS        VWRK2(5)=0.D0        VWRK2(6)=0.D0        VWRK2(7)=VWRK1(7)C      RETURN      END          

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