Sources Esope | Cast3M

tegral
C TEGRAL    SOURCE    CHAT      05/01/13    03:35:53     5004
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
      IF(NT.EQ.1000) THEN
         AK=1.009
         DT=0.001156464*XPI/NT
      ENDIF
      IF(NT.EQ.2000) THEN
         AK=1.004
         DT=0.0027278*XPI/NT
      ENDIF
      IF(NT.EQ.5000) THEN
         AK=1.0015
         DT=0.00417383*XPI/NT
      ENDIF
      IF(NT.EQ.10000) THEN
         AK=1.00065
         DT=0.0098077883*XPI/NT
      ENDIF
C
      DO 1 I=1,NT
         DTE(I)=DT*AK**(I-1)
         TETI(I)=DT*(AK**(I-1)-1.)/(AK-1.)
         TETM(I)=TETI(I)+0.5*DTE(I)
C
         DTE(2*NT+1-I)=DTE(I)
         TETM(2*NT+1-I)=2.*XPI-TETM(I)
         TETI(2*NT+1-I)=TETM(2*NT+1-I)-0.5*DTE(I)
 1    CONTINUE
C
      NT2=2*NT
C
      DO 2 I=1,NT2
         DTE(NT2+I)=DTE(I)
         TETM(NT2+I)=2.*XPI+TETM(I)
         TETI(NT2+I)=2.*XPI+TETI(I)
 2    CONTINUE
      RETURN
      END