Télécharger intfil.eso

Retour à la liste

Numérotation des lignes :

intfil
  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.  
  27. -INC PPARAM
  28. -INC CCOPTIO
  29. C
  30. IF (DK1.EQ.DK2) RETURN
  31. IF ((DK1.LT.0.D0).OR.(DK2.LT.0.D0)) RETURN
  32. C-----------------------------------------------------------------------
  33. C INITIALISATION
  34. C-----------------------------------------------------------------------
  35. SINT1=0
  36. SINT2=0
  37. SINT3=0
  38. SINT4=0
  39. C
  40. GP1=0
  41. GP2=0
  42. C
  43. DELTET=TETA2-TETA1
  44. C
  45. DO 10 IETA=-1,1,2
  46. TETKS1=TETA1+IETA*DKSI
  47. TETKS2=TETA2+IETA*DKSI
  48. C-----------------------------------------------------------------------
  49. X1=TETKS2*DK1
  50. X2=TETKS2*DK2
  51. X3=TETKS1*DK2
  52. X4=TETKS1*DK1
  53. C-----------------------------------------------------------------------
  54. CALL SINI(X1,SINT1)
  55. CALL SINI(X2,SINT2)
  56. CALL SINI(X3,SINT3)
  57. CALL SINI(X4,SINT4)
  58. C-----------------------------------------------------------------------
  59. C CALCUL DE GP1 ET GP2
  60. C
  61. GP1=GP1+SINT2-SINT1
  62. GP2=GP2+SINT2-SINT3-SINT1+SINT4
  63. 10 CONTINUE
  64. C
  65. XPI2=2.D0*XPI
  66. GP1=DELTET*GP1/XPI2
  67. GP2=GP2/XPI2
  68. RETURN
  69. END
  70.  
  71.  
  72.  
  73.  
  74.  

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