Télécharger la plaquette Cast3M


fr en

Télécharger cztmax.eso

Retour à la liste

Numérotation des lignes : 

cztmax

  1. C CZTMAX    SOURCE    JK148537  23/08/21    21:15:11     11723          
  2.       SUBROUTINE CZTMAX(iwrk52,iwrk53,iele,igau)
  3.       IMPLICIT INTEGER(I-N)
  4.       IMPLICIT REAL*8(A-H,O-Z)
  5. C
  6.  
  7. -INC PPARAM
  8. -INC CCREEL
  9. -INC CCOPTIO
  10. -INC CCHAMP
  11. -INC DECHE
  12.  
  13.  
  14.        wrk52 = iwrk52
  15.        wrk53 = iwrk53
  16.        t0 = ture0(1)
  17.        TF = turef(1)
  18.  
  19.        pha10 = rhas0(1)
  20.          z2 = pha10
  21.        pha20 = rhas0(2)
  22.          z1 = pha20
  23.        viph0 = valma0(1)
  24.  
  25.        AC1 = valma0(2)
  26.        AC2 = valma0(3)
  27.        AC3 = valma0(4)
  28.        AC4 = valma0(5)
  29. c       if (iele.eq.1.and.igau.eq.1) 
  30. c     &write(6,*) 'enzt',pha10,pha20,viph0, ac3,ac4,t0,TF
  31.  
  32.        if ((z1.le.xpetit.and.TF.le.AC1).or.
  33.      & (z2.le.xpetit.and.TF.ge.AC3)) then
  34.         viph = viph0
  35.         pha1f = pha10
  36.         pha2f = pha20
  37.         goto 1000
  38.        endif
  39.  
  40.        if (viph0.le.AC1) then
  41.          if (TF.gt.AC1) then
  42.            if (ac2 .ne. ac1) then
  43.              z1 = (TF - ac1) / (ac2 - ac1)
  44.            else if (ac2.eq.ac1) then
  45.             z1 = TF - ac1
  46.            endif
  47.            viph = AC1
  48.        if (z1.ge.1.) then
  49.          pha2f = 1.d0
  50.          pha1f = 0.d0
  51.          viph = AC3
  52.        else if (z1.lt.1.) then
  53.         if (z1.ge.0.and.z1.lt.1) then
  54.            pha2f = z1
  55.            pha1f = 1.d0 - z1
  56.         else
  57.            pha2f = 0.d0
  58.            pha1f = 1.d0
  59.         endif
  60.        endif
  61.            goto 1000
  62.          else
  63.         viph = viph0
  64.         pha1f = pha10
  65.         pha2f = pha20
  66.         goto 1000
  67.          endif 
  68.        endif
  69.  
  70.        if (viph0.ge.AC3) then
  71.          if(TF.lt.AC3) then
  72.            if (ac4.ne.ac3) then
  73.              z2 = (tf - ac3)/(ac4 - ac3)
  74.            else if (ac3.eq.ac4) then
  75.              if (ac3.ge.0.) then
  76.               z2 = ac3 - TF
  77.              else
  78.               z2 = TF - ac3
  79.              endif
  80.            endif
  81.            viph = AC3
  82. c       if (iele.eq.1.and.igau.eq.1) 
  83. c     &write(6,*) 'cz2',z2
  84.           if (z2.ge.1.) then
  85.            pha2f = 0.d0
  86.            pha1f = 1.d0
  87.            VIPH = AC1
  88.           else if(z2.lt.1.and.z2.ge.0) then
  89.            pha2f = 1.d0 - z2
  90.            pha1f = z2
  91.           else
  92.            pha2f = 1.d0
  93.            pha1f = 0.d0
  94.           endif
  95.           goto 1000
  96.          else
  97.         viph = viph0
  98.         pha1f = pha10
  99.         pha2f = pha20
  100.         goto 1000         
  101.          endif
  102.        endif
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  1000  continue
  109.        rhasf(1) = pha1f
  110.        rhasf(2) = pha2f
  111.        xmatf(1) = viph
  112.        xmatf(2) = AC1
  113.        xmatf(3) = AC2
  114.        xmatf(4) = AC3
  115.        xmatf(5) = AC4
  116.        xmatf(6) = valma0(6)
  117.        xmatf(7) = valma0(7)
  118.  
  119. c       if (iele.eq.1.and.igau.eq.1) write(6,*) 'sorzt',pha1f,pha2f,viph
  120.  
  121.         RETURN
  122.         END
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.  
  130.  
  131.  
  132.  
  133.  
  134.  

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