Télécharger la plaquette Cast3M


fr en

Télécharger convcp.eso

Retour à la liste

Numérotation des lignes : 

convcp

  1. C CONVCP    SOURCE    GOUNAND   19/07/23    21:15:06     10267          
  2.       SUBROUTINE CONVCP(DREE,DIMA,DMOD,DPHA,NP,INV)
  3.       IMPLICIT INTEGER(I-N)
  4.       IMPLICIT REAL*8(A-H,O-Z)
  5. C
  6. C=======================================================================
  7. C     INV=1
  8. C     CALCUL DU MODULE ET DE LA PHASE D'UN COMPLEXE DOUBLE PRECISION
  9. C     A PARTIR DE LA PARTIE REELLE ET DE LA PARTIE IMAGINAIRE
  10. C
  11. C                   ( ET VICE-VERSA      INV=-1 )
  12. C
  13. C     CREATION    : 15/12/87
  14. C     F.ROULLIER
  15. C=======================================================================
  16. C
  17.  
  18. -INC PPARAM
  19. -INC CCOPTIO
  20. -INC CCREEL
  21.       DIMENSION DREE(*),DIMA(*),DMOD(*),DPHA(*)
  22. C
  23.       IF (INV.EQ.1) THEN
  24.          RAP=180.D0/XPI
  25.          DO 10 I=1,NP
  26.             DR=DREE(I)
  27.             DI=DIMA(I)
  28.             DM=SQRT(DR*DR+DI*DI)
  29.             IF (ABS(DR).GT.XPETIT) THEN
  30.                DP=ATAN2(DI,DR)
  31.             ELSE
  32.                IF (DI.EQ.0.D0) DP=0.D0
  33.                IF (DI.GT.0.D0) DP=XPI/2.D0
  34.                IF (DI.LT.0.D0) DP=-XPI/2.D0
  35.             ENDIF
  36.             DMOD(I)=DM
  37.             DPHA(I)=DP*RAP
  38.  10      CONTINUE
  39.          RETURN
  40. C
  41.       ELSE
  42. C     
  43.          RAP=XPI/180.D0
  44.          DO 110 I=1,NP
  45.             DM=DMOD(I)
  46.             DP=DPHA(I)*RAP
  47.             DR=DM*COS(DP)
  48.             DI=DM*SIN(DP)
  49.             DREE(I)=DR
  50.             DIMA(I)=DI
  51.  110     CONTINUE
  52.       ENDIF
  53. C
  54.       RETURN
  55.       END
  56.  
  57.  
  58.  
  59.  
  60.  

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