Télécharger la plaquette Cast3M


fr en

Télécharger hbmordo.eso

Retour à la liste

Numérotation des lignes : 

hbmordo

  1. C HBMORDO   SOURCE    BP208322  20/09/18    21:16:48     10718          
  2.  
  3.         SUBROUTINE HBMORDO(NM,NDDL,OMEG,EXPIM,INDEIG) 
  4.  
  5. *=======================================================================  
  6. *  tri des valeurs propres par partie imaginaire croisante en module
  7. *=======================================================================  
  8.  
  9.       IMPLICIT INTEGER(I-N)
  10.       IMPLICIT REAL*8(A-H,O-Z)
  11.  
  12.         REAL*8 EXPIM(NM),OMEG,EXPIW(NM)
  13.         INTEGER NM,NDDL,INDEIG(2*NM),MINT,I
  14.  
  15. c         DO I=1,2*NDDL
  16. c             INDEIG(I)=0
  17. c         ENDDO
  18.          DO I=1,NM
  19.            INDEIG(I)=I
  20.          ENDDO
  21. c          
  22. c c  Copier EXPIM dans le tableau de travail EXPIW
  23. *          DO I = 1,NM
  24. *            EXPIW(I) = ABS(EXPIM(I))
  25. *          ENDDO
  26. c          
  27. c c    Pour l'instant, on prend simplement les 2*NDDL VP avec les parties 
  28. c c    imaginaires les plus petites en module
  29. *          DO I=1,2*NDDL
  30. *            MINT = MINLOC(EXPIW,DIM=1)
  31. *            INDEIG(I) = MINT 
  32. *            EXPIW(MINT) = 1.E6
  33. *          ENDDO
  34.  
  35. c   bp: on remplace ci-dessus par le classement de la liste
  36. c       on recuperera seulement les 1eres valeurs
  37.         CALL ORDO13(EXPIM,NM,.TRUE.,INDEIG)
  38.  
  39.         END
  40.  
  41.  
  42.  

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