Télécharger la plaquette Cast3M


fr en

Télécharger ottvai.eso

Retour à la liste

Numérotation des lignes : 

ottvai

  1. C OTTVAI    SOURCE    FD218221  21/06/10    21:15:52     11030          
  2.       SUBROUTINE OTTVAI(VAR1,XVAL,VAR2,RCZ,KV0,IERUT)
  3. *
  4.       IMPLICIT INTEGER(I-N)
  5.       IMPLICIT REAL*8 (A-H,O-Z)
  6.       DIMENSION VAR1(*),XVAL(*),VAR2(*)
  7. *
  8.         IERUT=0
  9.         RCZ=0.
  10.         GO TO  (1,2,3,4,5,6,7,8,9,10),KV0
  11. *
  12.  1      CONTINUE
  13.           RTRAC = XVAL(3)
  14.           GFTR = XVAL(4)
  15.           WO1= VAR1(2)
  16.                   FIL = VAR2(13)
  17.                   IF(FIL.EQ.0.D0) THEN
  18.                     WO0 = 0.D0
  19.                   ELSE  
  20.                     WO0  = RTRAC*RTRAC/(2.D0*GFTR*FIL)
  21.           ENDIF           
  22.           WO2 = RTRAC - WO0 * WO1
  23.           RCZ = - WO2
  24.           GO TO 20
  25. *
  26.  2      CONTINUE
  27.           RTRAC  = XVAL(3)
  28.           GFTR   = XVAL(4)
  29.           BETA   = XVAL(5)
  30.           WO3  = VAR1(1)
  31.           WO1 = VAR1(2)
  32.           IF(WO1.GT.0.) THEN
  33.                         FIL  = VAR2(13)
  34.                         WO0  = RTRAC*RTRAC/(2.D0*GFTR*FIL)
  35.                     WO4 = RTRAC - WO0 * WO1
  36.             WO4=  MAX(WO4,0.D0)
  37.             IF(WO4.GT.0.) THEN
  38.               WO2 = WO4*(BETA - WO3/WO1)/(1.-BETA)
  39.               RCZ = WO2
  40.             ELSE
  41.               RCZ = (BETA - WO3/WO1) * RTRAC
  42.             ENDIF
  43.           ELSE
  44.             RCZ = - RTRAC 
  45.           ENDIF
  46.         GO TO 20
  47. *
  48.  3      CONTINUE
  49.  
  50.           RTRAC  = XVAL(3)
  51.           WO3  = VAR1(1)
  52.           WO1 = VAR1(2)
  53.           IF(WO1.GT.0.) THEN
  54.               RCZ = (WO3/WO1 - 1.D0) * RTRAC
  55.           ELSE
  56.             RCZ = - RTRAC 
  57.           ENDIF
  58.         GO TO 20
  59. *
  60.  4      CONTINUE
  61.           RTRAC = XVAL(3)
  62.           GFTR  = XVAL(4)
  63.           WO1= VAR1(4)
  64.                   FIL = VAR2(14)
  65.                   IF(FIL.EQ.0.D0) THEN
  66.                     WO0 = 0.D0
  67.                   ELSE  
  68.                     WO0  = RTRAC*RTRAC/(2.D0*GFTR*FIL)
  69.           ENDIF 
  70.           WO2 = RTRAC - WO0 * WO1
  71.           RCZ = - WO2
  72.           GO TO 20
  73. *
  74.  5      CONTINUE
  75.           RTRAC  = XVAL(3)
  76.           GFTR   = XVAL(4)
  77.           BETA   = XVAL(5)
  78.           WO3  = VAR1(3)
  79.           WO1 = VAR1(4)
  80.           IF(WO1.GT.0.) THEN
  81.                         FIL  = VAR2(14)
  82.                         WO0  = RTRAC*RTRAC/(2.D0*GFTR*FIL)
  83.                     WO4 = RTRAC - WO0 * WO1
  84.             WO4=  MAX(WO4,0.D0)
  85.             IF(WO4.GT.0.) THEN
  86.               WO2 = WO4*(BETA - WO3/WO1)/(1.-BETA)
  87.               RCZ = WO2
  88.             ELSE
  89.               RCZ = (BETA - WO3/WO1) * RTRAC
  90.             ENDIF
  91.           ELSE
  92.             RCZ = - RTRAC 
  93.           ENDIF
  94.         GO TO 20
  95. *
  96.  6      CONTINUE
  97.           RTRAC  = XVAL(3)
  98.           WO3  = VAR1(3)
  99.           WO1 = VAR1(4)
  100.           IF(WO1.GT.0.) THEN
  101.               RCZ = (WO3/WO1 - 1.D0) * RTRAC
  102.           ELSE
  103.             RCZ = - RTRAC 
  104.           ENDIF
  105.         GO TO 20
  106. *
  107.  7      CONTINUE
  108.           RTRAC = XVAL(3)
  109.           GFTR  = XVAL(4)
  110.           WO1= VAR1(6)
  111.                   FIL = VAR2(15)
  112.                   IF(FIL.EQ.0.D0) THEN
  113.                     WO0 = 0.D0
  114.                   ELSE  
  115.                     WO0  = RTRAC*RTRAC/(2.D0*GFTR*FIL)
  116.           ENDIF
  117.           WO2 = RTRAC - WO0 * WO1
  118.           RCZ = - WO2
  119.           GO TO 20
  120. *
  121.  8      CONTINUE
  122.           RTRAC  = XVAL(3)
  123.           GFTR   = XVAL(4)
  124.           BETA   = XVAL(5)
  125.           WO3  = VAR1(5)
  126.           WO1 = VAR1(6)
  127.           IF(WO1.GT.0.) THEN
  128.                         FIL  = VAR2(15)
  129.                         WO0  = RTRAC*RTRAC/(2.D0*GFTR*FIL)
  130.                     WO4 = RTRAC - WO0 * WO1
  131.             WO4=  MAX(WO4,0.D0)
  132.             IF(WO4.GT.0.) THEN
  133.               WO2 = WO4*(BETA - WO3/WO1)/(1.-BETA)
  134.               RCZ = WO2
  135.             ELSE
  136.               RCZ = (BETA - WO3/WO1) * RTRAC
  137.             ENDIF
  138.           ELSE
  139.             RCZ = - RTRAC 
  140.           ENDIF
  141.         GO TO 20
  142. *
  143.  9      CONTINUE
  144.           RTRAC  = XVAL(3)
  145.           WO3  = VAR1(5)
  146.           WO1 = VAR1(6)
  147.           IF(WO1.GT.0.) THEN
  148.               RCZ = (WO3/WO1 - 1) * RTRAC
  149.           ELSE
  150.             RCZ = - RTRAC 
  151.           ENDIF
  152.         GO TO 20
  153. *
  154.  10     CONTINUE
  155.           RTRAC = XVAL(3)
  156.           WO2 = RTRAC
  157.           RCZ = - WO2
  158. 20    CONTINUE
  159.       RETURN
  160.       END
  161.  
  162.  
  163.  

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