Télécharger la plaquette Cast3M


fr en

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