Télécharger la plaquette Cast3M


fr en

Télécharger w2frqc.eso

Retour à la liste

Numérotation des lignes : 

w2frqc

  1. C W2FRQC    SOURCE    CB215821  16/04/21    21:18:41     8920
  2.       SUBROUTINE W2FRQC (XRVP, XIVP, W2SHIF, XRW2,XIW2, XRFREQ, XIFREQ)
  3. ************************************************************************
  4. *
  5. *                             W 2 F R Q C
  6. *                             -----------
  7. *                       (PRONONCER "OMEGA2 FREQ C")
  8. *
  9. * FONCTION:
  10. * ---------
  11. *
  12. *     FOURNIR LA FREQUENCE PROPRE COMPEXEET LA PULSATION AU CARRE, CONNAISSANT
  13. *     LA PULSATION DE DECALAGE (AU CARRE ET REELLE) ET LA PULSATION PROPRE (AU
  14. *     CARRE COMPLEXE) DECALEE. ON PASSE PAR LE CALCUL DU MODULE ET DE L'ARGUMENT
  15. *     (CONTEXTE: CALCUL DE MODES PROPRES COMPLEXES)
  16. *
  17. * MODE D'APPEL:
  18. * -------------
  19. *
  20. *     CALL W2FRQC (RXVP, IXVP, W2SHIF, XRW2,  XIW2, XRFREQ, XIFREQ)
  21. *
  22. * PARAMETRES:   (E)=ENTREE   (S)=SORTIE
  23. * -----------
  24. *
  25. *     XRVP    REEL DP   (E)  PARTIE REELLE DE PULSATION CALCULEE AU CARRE (DECALEE).
  26. *     XIVP    REEL DP   (E)  PARTIE IMGINAIRE DE PULSATION CALCULEE
  27. *     W2SHIF  REEL DP   (E)  PULSATION REELLE DE DECALAGE AU CARRE.
  28. *     XRW2    REEL DP   (S)  = XRVP + W2SHIF.
  29. *     XIW2    REEL DP   (S)  = XIVP
  30. *     XRFREQ  REEL DP   (S)  FREQUENCE REELLE ASSOCIEE A LA PULSATION PROPRE
  31. *     XIFREQ  REEL DP   (S)  FREQUANCE IMAGINAIRE ASSOCIEE A LA PULSATION PROPRE
  32. *
  33. * AUTEUR, DATE DE CREATION:
  34. * -------------------------
  35. *
  36. *     C. LE BIDEAU     JUILLET 2001
  37. *
  38. * LANGAGE:
  39. * --------
  40. *
  41. *     FORTRAN77 + ESOPE
  42. *
  43. ************************************************************************
  44. *
  45.       IMPLICIT INTEGER(I-N)
  46.       IMPLICIT real*8 (a-h,o-z)
  47. -INC CCREEL
  48.  
  49. -INC PPARAM
  50. -INC CCOPTIO
  51. *
  52.       REAL*8 XRVP, XIVP, W2SHIF, XRW2, XIW2 ,XRFREQ, XIFREQ, DEUXPI, TETA,
  53.      &  XNW2, XRW, XIW, ALPHA
  54. *
  55.       PARAMETER (DEUXPI = 2.D0*XPI)
  56. *
  57. *     PULSATION AU CARRE
  58. *
  59.       XRW2 = XRVP + W2SHIF
  60.       XIW2 = XIVP
  61. *
  62. *     CALCUL DE LA FREQUENCE
  63. *
  64. *     ARGUMENT DE W2
  65. *
  66.       TETA = ATAN2(XIW2,XRW2)
  67. *
  68. *     MODULE DE W2
  69. *
  70.       XNW2 = SQRT( XRW2**2 + XIW2**2)
  71. *     RACINE DE W2 (W2=XRVP + i*XIVP)
  72. *
  73.       XNW = SQRT (XNW2)
  74.       XRW = XNW * COS(TETA/2)
  75.       XIW = XNW * SIN(TETA/2)
  76. *     FREQUENCE
  77. *
  78. *     LA FREQUENCE CORRESPODANTE EST A
  79. *     PARTIE IMAGINAIRE NEGATIVE
  80. *     BP208322 (mars 2009) : ne semble pas justifié....
  81.       XRFREQ = XRW/DEUXPI
  82. *      XIFREQ = -1.D0*XIW/DEUXPI
  83.       XIFREQ = XIW/DEUXPI
  84. *
  85.       END
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  

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