Esope source files

   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 : *     
  19 : *     CALL W2FREQ (W2CALC,W2SHIF,W2,FREQ)
  20 : *     
  21 : * PARAMETRES:   (E)=ENTREE   (S)=SORTIE
  22 : * -----------
  23 : *     
  24 : *     W2CALC  REEL DP   (E)  PULSATION CALCULEE AU CARRE (DECALEE).
  25 : *     W2SHIF  REEL DP   (E)  PULSATION DE DECALAGE AU CARRE.
  26 : *     W2      REEL DP   (S)  = W2CALC + W2SHIF.
  27 : *     FREQ    REEL DP   (S)  FREQUENCE ASSOCIEE A LA PULSATION PROPRE
  28 : *                            AU CARRE "W2".
  29 : *     
  30 : * AUTEUR, DATE DE CREATION:
  31 : * -------------------------
  32 : *     
  33 : *     PASCAL MANIGOT     28 NOVEMBRE 1984
  34 : *     
  35 : * LANGAGE:
  36 : * --------
  37 : *     
  38 : *     FORTRAN77
  39 : *     
  40 : ************************************************************************
  41 : *     
  42 :       IMPLICIT INTEGER(I-N)
  43 : -INC CCREEL
  44 : -INC CCOPTIO
  45 : *     
  46 :       REAL*8 W2CALC,W2SHIF,W2,FREQ,DEUXPI
  47 : *     
  48 :       LOGICAL LIMAGE
  49 : *     
  50 :       PARAMETER (DEUXPI = 2.D0*XPI)
  51 : *     
  52 : *     PULSATION AU CARRE:
  53 :       W2 = W2CALC + W2SHIF
  54 :       IF ((W2 .LT. 0.D0).AND.(.NOT.LIMAGE)) THEN
  55 :          REAERR(1) = REAL(SQRT(ABS(W2SHIF)) / DEUXPI)
  56 :          NUMERR = 142
  56 :          CALL ERREUR (NUMERR)
  57 :          W2 = -1.D0 * W2
  58 :       END IF
  59 : *     
  60 : *     FREQUENCE:
  61 :       FREQ = SQRT(ABS(W2)) / DEUXPI
  62 :       IF(LIMAGE) THEN
  63 :          FREQ=SIGN(FREQ,W2)
  64 :       ENDIF
  65 : *     
  66 :       END
  67 : 
  68 : 
  69 :  
  70 :  
  71 :