Esope source files

   1 : C MISIMP    SOURCE    AF221230  13/08/01    21:15:05     7808           
   2 : C MISIMP    SOURCE              
   3 :       SUBROUTINE MISIMP(MTABI,JIMPD,NF,TITR,IL,IC)
   4 : C
   5 : C=======================================================================  
   6 : C      ECRITURE DES IMPEDANCES MISS3D DANS UNE EVOLUTION COMPLEXE
   7 : C
   8 : C      Appelle par la routine MISL                        
   9 : C=======================================================================
  10 : C
  11 :       IMPLICIT INTEGER(I-N)
  12 :       IMPLICIT REAL*8(A-H,O-Z)
  13 :       CHARACTER*4 TITR
  14 :       CHARACTER*72 lemot
  15 :       LOGICAL OK1
  16 :       PARAMETER(DEUX=2.0D0)
  17 : 
  18 : -INC SMEVOLL
  19 : -INC SMLREEL
  20 :       SEGMENT MATIMPD
  21 :         COMPLEX*16 IMPD(6,6,NFR)
  22 :       ENDSEGMENT
  23 :          MATIMPD=JIMPD
  24 :          SEGACT MATIMPD
  25 : 
  25 :          CALL ACCTAB(MTABI,'MOT',0,0.0D0,TITR,.TRUE.,0,
  26 :      &        'EVOLUTIO',IP,RR,lemot,OK1,ICO)
  27 :          MEVOLL=ICO
  28 :          SEGACT MEVOLL*MOD
  29 :          IEVTEX='Impedance '//TITR//'       Rouge : partie reelle -- '//
  30 :      &          'Vert : partie imaginaire'
  31 :          KEVOL1=IEVOLL(1)
  32 :          KEVOL2=IEVOLL(2)
  33 :          SEGACT KEVOL1*MOD,KEVOL2*MOD
  34 :          MLREE1=KEVOL1.IPROGY
  35 :          MLREE2=KEVOL2.IPROGY
  36 :          SEGACT MLREE1*MOD,MLREE2*MOD
  37 :          KEVOL1.NUMEVX=2
  38 :          KEVOL2.NUMEVX=4
  39 :          KEVOL1.KEVTEX='Partie reelle'
  40 :          KEVOL2.KEVTEX='Partie imaginaire'
  41 :          
  42 :          IF(IL.EQ.IC)THEN
  43 :            DO JF=1, NF
  44 :              MLREE1.PROG(JF)=DREAL(IMPD(IL,IC,JF))
  45 :              MLREE2.PROG(JF)=-DIMAG(IMPD(IL,IC,JF))
  46 :            ENDDO
  47 :          ELSE
  48 :            DO JF=1, NF
  49 :              MLREE1.PROG(JF)=(DREAL(IMPD(IL,IC,JF))+
  50 :      &                        DREAL(IMPD(IC,IL,JF)))/DEUX
  51 :              MLREE2.PROG(JF)=-(DIMAG(IMPD(IL,IC,JF))+
  52 :      &                        DIMAG(IMPD(IC,IL,JF)))/DEUX
  53 :            ENDDO
  54 :          ENDIF
  55 :          SEGDES MLREE1,MLREE2
  56 :          SEGDES KEVOL1,KEVOL2
  57 :          SEGDES MEVOLL
  58 : 
  59 :        SEGDES MATIMPD
  60 : 
  61 :       END
  62 :