Télécharger la plaquette Cast3M


fr en

Télécharger clebl4.eso

Retour à la liste

Numérotation des lignes : 

clebl4

  1. C CLEBL4    SOURCE    PV        17/12/08    21:15:49     9660           
  2. C====================================================================
  3. C
  4.       SUBROUTINE CLEBL4(iwrk52,VWRK1,DTPS,VWRK2)
  5. C
  6. C====================================================================
  7. C     LB, 06/97
  8. C     =========
  9. C         Modele de Leblond et Devaux
  10. C         Un seul chgt de phase - cas du refroidissement
  11. C     Pour le moment zeq et tau sont calcules en local
  12. C
  13.       IMPLICIT INTEGER(I-N)
  14.        IMPLICIT REAL*8(A-H,O-Z)
  15. -INC DECHE
  16. C
  17.       REAL*8 VWRK1(*),VWRK2(7)
  18. C
  19.       T0  =VWRK1(1)
  20.       TP0 =VWRK1(2)
  21.       ZA0 =VWRK1(3)
  22.       ZF0 =VWRK1(4)
  23.       VMS1=VWRK1(7)
  24.  
  25.        wrk52 = iwrk52
  26.       AC1 = xmat(1)
  27.       AC2 = xmat(2)
  28.       AC3 = xmat(3)
  29.       AC4 = xmat(4)
  30.       TAU1 = xmat(5)
  31.       TAU2 = xmat(6)
  32.       TAU3 = xmat(7)
  33.       TAU4 = xmat(8)
  34. C
  35. C     Calcul de Zeq
  36. C
  37. c       if (T0.LE.600.) then
  38. c          zeq = 1.
  39. c       else if ((T0.GT.600.).AND.(T0.LT.800.)) then
  40. c          zeq = (800. - T0) / 200.
  41. c       else if (T0.GE.800.) then
  42. c          zeq = 0.
  43. c       endif
  44.        if (T0.LE.AC1) then
  45.           zeq = 1.
  46.        else if ((T0.GT.AC1).AND.(T0.LT.AC2).and.(AC3.ne.0.)) then
  47.           zeq = (AC2 - T0) / AC3
  48.        else if (T0.GE.AC2) then
  49.           zeq = 0.
  50.        endif
  51. C
  52. C      Calcul de tau
  53. C
  54. c       if (T0.LE.600.) then
  55. c          tau = 100000.
  56. c       else if ((T0.GT.600.).AND.(T0.LT.620.)) then
  57. c         tau = 100000. - (T0 - 600.) * 95000. / 20.
  58. c       else if (T0.GE.620.) then
  59. c          tau = 500.
  60. c       endif
  61.        if (T0.LE.AC1) then
  62.           tau = TAU1
  63.        else if ((T0.GT.AC1).AND.(T0.LT.(AC1+AC4)).and.(AC4.ne.0.)) then
  64.           tau = TAU1 - (T0 - AC1) * TAU2 / AC4
  65.        else if (T0.GE.620.) then
  66.           tau = TAU3
  67.        endif
  68. C
  69. C     Calcul de zpoint
  70. C
  71.       zpoint = (zeq - ZF0) / tau
  72. C
  73.         VWRK1(1)=T0
  74.         VWRK2(2)=TP0
  75.         VWRK2(3)=VWRK1(3) - zpoint * DTPS
  76.         VWRK2(4)=VWRK1(4) + zpoint * DTPS
  77.         VWRK2(5)=0.D0
  78.         VWRK2(6)=0.D0
  79.         VWRK2(7)=VWRK1(7)
  80. C
  81.       RETURN
  82.       END
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  

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