Sources Esope | Cast3M

tegral
C TEGRAL    SOURCE    CB215821  25/03/25    21:15:06     12217          
C
C*******************************************************************
C
      SUBROUTINE TEGRAL(SANGLE)
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
C
-INC CCREEL
C
C*******************************************************************
C  Calcul des intervalles d'integration par rapport a TETA
C*******************************************************************
C
C
      SEGMENT SANGLE
         REAL*8 TETM(NT4)
         REAL*8 TETI(NT4)
         REAL*8 DTE(NT4)
      ENDSEGMENT
C
C  calcul des positions et des elements d'arc de cercle pour
C  l'integration par rapport a TETA
C
      NT4 = TETM(/1)
      NT  = NT4/4
      XNT = REAL(NT)
      IF(NT.EQ.1000) THEN
         AK=1.009D0
         DT=0.001156464D0*XPI/XNT
      ENDIF
      IF(NT.EQ.2000) THEN
         AK=1.004D0
         DT=0.0027278D0*XPI/XNT
      ENDIF
      IF(NT.EQ.5000) THEN
         AK=1.0015D0
         DT=0.00417383D0*XPI/XNT
      ENDIF
      IF(NT.EQ.10000) THEN
         AK=1.00065D0
         DT=0.0098077883D0*XPI/XNT
      ENDIF
C
      DO 1 I=1,NT
         DTE(I) =DT*AK**(I-1)
         TETI(I)=DT*(AK**(I-1)-1.D0)/(AK-1.D0)
         TETM(I)=TETI(I)+0.5D0*DTE(I)
C
         DTE (2*NT+1-I)=DTE(I)
         TETM(2*NT+1-I)=2.D0*XPI-TETM(I)
         TETI(2*NT+1-I)=TETM(2*NT+1-I)-0.5D0*DTE(I)
 1    CONTINUE
C
      NT2=2*NT
C
      DO 2 I=1,NT2
         DTE(NT2+I)=DTE(I)
         TETM(NT2+I)=2.D0*XPI+TETM(I)
         TETI(NT2+I)=2.D0*XPI+TETI(I)
 2    CONTINUE
      RETURN
      END