Esope source files

   1 : C TEGRAL    SOURCE    CHAT      05/01/13    03:35:53     5004           
   2 : C
   3 : C*******************************************************************
   4 : C
   5 :       SUBROUTINE TEGRAL(SANGLE)
   6 :       IMPLICIT INTEGER(I-N)
   7 :       IMPLICIT REAL*8(A-H,O-Z)
   8 : C
   9 : -INC CCREEL
  10 : C
  11 : C*******************************************************************
  12 : C  Calcul des intervalles d'integration par rapport a TETA
  13 : C*******************************************************************
  14 : C
  15 : C
  16 :       SEGMENT SANGLE
  17 :          REAL*8 TETM(NT4)
  18 :          REAL*8 TETI(NT4)
  19 :          REAL*8 DTE(NT4)
  20 :       ENDSEGMENT
  21 : C
  22 : C  calcul des positions et des elements d'arc de cercle pour
  23 : C  l'integration par rapport a TETA
  24 : C
  25 :       NT4=TETM(/1)
  26 :       NT=NT4/4 
  27 :       IF(NT.EQ.1000) THEN
  28 :          AK=1.009
  29 :          DT=0.001156464*XPI/NT
  30 :       ENDIF
  31 :       IF(NT.EQ.2000) THEN
  32 :          AK=1.004
  33 :          DT=0.0027278*XPI/NT
  34 :       ENDIF
  35 :       IF(NT.EQ.5000) THEN
  36 :          AK=1.0015
  37 :          DT=0.00417383*XPI/NT
  38 :       ENDIF
  39 :       IF(NT.EQ.10000) THEN
  40 :          AK=1.00065
  41 :          DT=0.0098077883*XPI/NT
  42 :       ENDIF
  43 : C
  44 :       DO 1 I=1,NT
  45 :          DTE(I)=DT*AK**(I-1)
  46 :          TETI(I)=DT*(AK**(I-1)-1.)/(AK-1.)
  47 :          TETM(I)=TETI(I)+0.5*DTE(I)
  48 : C
  49 :          DTE(2*NT+1-I)=DTE(I)
  50 :          TETM(2*NT+1-I)=2.*XPI-TETM(I)
  51 :          TETI(2*NT+1-I)=TETM(2*NT+1-I)-0.5*DTE(I)
  52 :  1    CONTINUE
  53 : C
  54 :       NT2=2*NT
  55 : C
  56 :       DO 2 I=1,NT2
  57 :          DTE(NT2+I)=DTE(I)
  58 :          TETM(NT2+I)=2.*XPI+TETM(I)
  59 :          TETI(NT2+I)=2.*XPI+TETI(I)
  60 :  2    CONTINUE
  61 :       RETURN
  62 :       END
  63 :  
  64 :  
  65 :  
  66 :