Télécharger w2freq.eso

Retour à la liste

Numérotation des lignes :

w2freq
  1. C W2FREQ SOURCE CHAT 05/01/13 04:11:20 5004
  2. SUBROUTINE W2FREQ (W2CALC,W2SHIF,W2,FREQ,LIMAGE)
  3. ************************************************************************
  4. *
  5. * W 2 F R E Q
  6. * -----------
  7. * (PRONONCER "OMEGA2 FREQ")
  8. *
  9. * FONCTION:
  10. * ---------
  11. *
  12. * FOURNIR LA FREQUENCE PROPRE ET LA PULSATION AU CARRE, CONNAISSANT
  13. * LA PULSATION DE DECALAGE (AU CARRE) ET LA PULSATION PROPRE (AU
  14. * CARRE) DECALEE.
  15. * (CONTEXTE: CALCUL DE MODES PROPRES)
  16. *
  17. * MODE D'APPEL:
  18. * -------------
  19. *
  20. * CALL W2FREQ (W2CALC,W2SHIF,W2,FREQ)
  21. *
  22. * PARAMETRES: (E)=ENTREE (S)=SORTIE
  23. * -----------
  24. *
  25. * W2CALC REEL DP (E) PULSATION CALCULEE AU CARRE (DECALEE).
  26. * W2SHIF REEL DP (E) PULSATION DE DECALAGE AU CARRE.
  27. * W2 REEL DP (S) = W2CALC + W2SHIF.
  28. * FREQ REEL DP (S) FREQUENCE ASSOCIEE A LA PULSATION PROPRE
  29. * AU CARRE "W2".
  30. *
  31. * AUTEUR, DATE DE CREATION:
  32. * -------------------------
  33. *
  34. * PASCAL MANIGOT 28 NOVEMBRE 1984
  35. *
  36. * LANGAGE:
  37. * --------
  38. *
  39. * FORTRAN77
  40. *
  41. ************************************************************************
  42. *
  43. IMPLICIT INTEGER(I-N)
  44. -INC CCREEL
  45.  
  46. -INC PPARAM
  47. -INC CCOPTIO
  48. *
  49. REAL*8 W2CALC,W2SHIF,W2,FREQ,DEUXPI
  50. *
  51. LOGICAL LIMAGE
  52. *
  53. PARAMETER (DEUXPI = 2.D0*XPI)
  54. *
  55. * PULSATION AU CARRE:
  56. W2 = W2CALC + W2SHIF
  57. IF ((W2 .LT. 0.D0).AND.(.NOT.LIMAGE)) THEN
  58. REAERR(1) = REAL(SQRT(ABS(W2SHIF)) / DEUXPI)
  59. NUMERR = 142
  60. CALL ERREUR (NUMERR)
  61. W2 = -1.D0 * W2
  62. END IF
  63. *
  64. * FREQUENCE:
  65. FREQ = SQRT(ABS(W2)) / DEUXPI
  66. IF(LIMAGE) THEN
  67. FREQ=SIGN(FREQ,W2)
  68. ENDIF
  69. *
  70. END
  71.  
  72.  
  73.  
  74.  
  75.  
  76.  

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