Esope source files

   1 : C FACOMP    SOURCE    CHAT      05/01/12    23:55:46     5004           
   2 :       SUBROUTINE FACOMP(PREEL,PIMAG,PMODU)                                      
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
   5 : C                                                                               
   6 : C=======================================================================        
   7 : C  =  ROUTINE FORTRAN CALCULANT LA PHASE D'UN COMPLEXE EN DEGRES      =         
   8 : C  =  ON ENTRE LA PARTIE REELLE, LA PARTIE IMAGINAIRE, ET LE MODULE   =         
   9 : C  =  LA PHASE RESSORT A LA PLACE DU MODULE                           =         
  10 : C  =                                                                  =         
  11 : C  =                                                                  =         
  12 : C  =  CREATION    : 13/04/87                                          =         
  13 : C  =  PROGRAMMEUR : BEAUFILS                                          =         
  14 : C=======================================================================        
  15 : C                                                                               
  16 : -INC CCOPTIO                                                                    
  17 : -INC CCREEL                                                                     
  18 : C                                                                               
  19 :       DATA EPS/1.E-20/                                                          
  20 : C                                                                               
  21 :       IF(ABS(PMODU).GT.EPS) THEN                                                
  22 :           PCOS=PREEL/PMODU                                                      
  23 :       ELSE                                                                      
  24 :           PMODU=0.D0                                                            
  25 :           RETURN                                                                
  26 :       ENDIF                                                                     
  27 : C                                                                               
  28 :       IF(ABS(PCOS).LT.EPS) THEN                                                 
  29 :           PMODU=90.D0*PIMAG/ABS(PIMAG)                                          
  30 :       ELSE                                                                      
  31 :           IF(ABS(PCOS).GT.1.D0) PCOS=(1.D0-EPS)*PCOS/ABS(PCOS)                  
  32 :           PMODU=180.D0/XPI*ACOS(PCOS)                                           
  33 :           IF(PIMAG.LT.0.D0) PMODU=-PMODU                                        
  34 :       ENDIF                                                                     
  35 : C                                                                               
  36 :       RETURN                                                                    
  37 :       END                                                                       
  38 :  
  39 :  
  40 :