Télécharger la plaquette Cast3M


fr en

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