Télécharger intfil.eso

Retour à la liste

Numérotation des lignes :

  1. C INTFIL SOURCE PV 05/04/25 21:15:30 5077
  2. SUBROUTINE INTFIL(GP1,GP2,DKSI,TETA1,TETA2,DK1,DK2)
  3. IMPLICIT INTEGER(I-N)
  4. IMPLICIT REAL*8(A-H,O-Z)
  5. C
  6. C=======================================================================
  7. C SOUS-PROGRAMME FORTRAN QUI CALCULE LES TERMES DES MATRICES ANBN
  8. C
  9. C TETA2 DK2
  10. C 1 / / SIN (K * (TETA-DKSI))
  11. C GP1(T) = ___ / / _____________________ DK DTETA
  12. C PI / / K
  13. C TETA1 DK1
  14. C TETA2 DK2
  15. C 1 / /
  16. C GP2(T) = ___ / / COS(K*(TETA-DKSI)) DK DTETA
  17. C PI / /
  18. C TETA1 DK1
  19. C
  20. C
  21. C-----------------------------------------------------------------------
  22. C CREATION: J-M GAUTIER AVRIL 87
  23. C=======================================================================
  24. C
  25. -INC CCREEL
  26. -INC CCOPTIO
  27. C
  28. IF (DK1.EQ.DK2) RETURN
  29. IF ((DK1.LT.0.D0).OR.(DK2.LT.0.D0)) RETURN
  30. C-----------------------------------------------------------------------
  31. C INITIALISATION
  32. C-----------------------------------------------------------------------
  33. SINT1=0
  34. SINT2=0
  35. SINT3=0
  36. SINT4=0
  37. C
  38. GP1=0
  39. GP2=0
  40. C
  41. DELTET=TETA2-TETA1
  42. C
  43. DO 10 IETA=-1,1,2
  44. TETKS1=TETA1+IETA*DKSI
  45. TETKS2=TETA2+IETA*DKSI
  46. C-----------------------------------------------------------------------
  47. X1=TETKS2*DK1
  48. X2=TETKS2*DK2
  49. X3=TETKS1*DK2
  50. X4=TETKS1*DK1
  51. C-----------------------------------------------------------------------
  52. CALL SINI(X1,SINT1)
  53. CALL SINI(X2,SINT2)
  54. CALL SINI(X3,SINT3)
  55. CALL SINI(X4,SINT4)
  56. C-----------------------------------------------------------------------
  57. C CALCUL DE GP1 ET GP2
  58. C
  59. GP1=GP1+SINT2-SINT1
  60. GP2=GP2+SINT2-SINT3-SINT1+SINT4
  61. 10 CONTINUE
  62. C
  63. XPI2=2.D0*XPI
  64. GP1=DELTET*GP1/XPI2
  65. GP2=GP2/XPI2
  66. RETURN
  67. END
  68.  
  69.  
  70.  
  71.  
  72.  

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