Esope source files

   1 : C KBICK1    SOURCE    NC        11/03/11    21:15:02     6899           
   2 :       SUBROUTINE KBICK1(X,XKI1)
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8 (A-H,O-Z)
   5 : C    Include contenant quelques constantes dont XPI :
   6 : -INC CCREEL
   7 : C*********************************************************************
   8 : C Calcul des facteurs de forme en milieu absorbant en 2D-plan
   9 : C SP appele par KRCOA2
  10 : C    Calcul de la fonction de Bickley-Naylor à l'ordre 3 en X
  11 : C                 normalisée
  12 : C      (integrale de o à pi/2)
  13 : C     entree:
  14 : C      X       : abscisse
  15 : C      NT      : discrétisation de l'intervalle (0,pi/2)
  16 : C     resultat:
  17 : C      XKI1    : valeur de la fonction
  18 : C     
  19 : C*********************************************************************
  20 : C
  21 : C     méthode des trapèzes
  22 : C     découpage en NT segments identiques
  23 :          NT=100
  24 : 
  25 :          DTETA = XPI/2./NT
  26 : 
  27 :          G  = 0.D0
  28 :          DO  K=1,NT
  29 :          TETA = DTETA*(0.5 + (K-1))
  30 :          SINTETA = DSIN(TETA)
  31 :          G1 = DEXP (-X/SINTETA)
  32 :          G1 = G1 * SINTETA *SINTETA
  33 :          G = G + G1
  34 :          ENDDO
  35 :          XKI1= G*2.D0/NT
  36 : C        WRITE(6,*) 'BICKLEY: X KI1: ',X, XKI1
  37 :       RETURN
  38 :       END
  39 :  
  40 :