Esope source files

   1 : C W2FRQC    SOURCE    BP208322  09/03/20    21:15:26     6331           
   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 : *     
  19 : *     CALL W2FRQC (RXVP, IXVP, W2SHIF, XRW2,  XIW2, XRFREQ, XIFREQ)
  20 : *     
  21 : * PARAMETRES:   (E)=ENTREE   (S)=SORTIE
  22 : * -----------
  23 : *     
  24 : *     XRVP    REEL DP   (E)  PARTIE REELLE DE PULSATION CALCULEE AU CARRE (DECALEE).
  25 : *     XIVP    REEL DP   (E)  PARTIE IMGINAIRE DE PULSATION CALCULEE
  26 : *     W2SHIF  REEL DP   (E)  PULSATION REELLE DE DECALAGE AU CARRE.
  27 : *     XRW2    REEL DP   (S)  = XRVP + W2SHIF.
  28 : *     XIW2    REEL DP   (S)  = XIVP
  29 : *     XRFREQ  REEL DP   (S)  FREQUENCE REELLE ASSOCIEE A LA PULSATION PROPRE
  30 : *     XIFREQ  REEL DP   (S)  FREQUANCE IMAGINAIRE ASSOCIEE A LA PULSATION PROPRE
  31 : *     
  32 : * AUTEUR, DATE DE CREATION:
  33 : * -------------------------
  34 : *     
  35 : *     C. LE BIDEAU     JUILLET 2001
  36 : *     
  37 : * LANGAGE:
  38 : * --------
  39 : *     
  40 : *     FORTRAN77 + ESOPE
  41 : *     
  42 : ************************************************************************
  43 : *     
  44 :       IMPLICIT INTEGER(I-N)
  45 :       IMPLICIT real*8 (a-h,o-z)
  46 : -INC CCREEL
  47 : -INC CCOPTIO
  48 : *     
  49 :       REAL*8 XRVP, XIVP, W2SHIF, XRW2, XIW2 ,XRFREQ, XIFREQ, DEUXPI, TETA,
  50 :      &  XNW2, XRW, XIW, ALPHA
  51 : *
  52 :       PARAMETER (DEUXPI = 2.D0*XPI)
  53 : *     
  54 : *     PULSATION AU CARRE
  55 : *
  56 :       XRW2 = XRVP + W2SHIF
  57 :       XIW2 = XIVP
  58 : *
  59 : *     CALCUL DE LA FREQUENCE
  60 : *
  61 : *     ARGUMENT DE W2
  62 : *
  63 :       TETA = DATAN2 (XIW2,XRW2)
  64 : *
  65 : *     MODULE DE W2
  66 : *
  67 :       XNW2 = SQRT( XRW2**2 + XIW2**2)
  68 : *     RACINE DE W2 (W2=XRVP + i*XIVP)
  69 : *
  70 :       XNW = SQRT (XNW2)
  71 :       XRW = XNW * COS(TETA/2)
  72 :       XIW = XNW * SIN(TETA/2)
  73 : *     FREQUENCE
  74 : *
  75 : *     LA FREQUENCE CORRESPODANTE EST A 
  76 : *     PARTIE IMAGINAIRE NEGATIVE
  77 : *     BP208322 (mars 2009) : ne semble pas justifié....
  78 :       XRFREQ = XRW/DEUXPI
  79 : *      XIFREQ = -1.D0*XIW/DEUXPI
  80 :       XIFREQ = XIW/DEUXPI
  81 : *
  82 :       END
  83 :  
  84 :  
  85 :  
  86 :