Sources Esope | Cast3M

ottodd
C OTTODD    SOURCE    CHAT      05/01/13    02:07:21     5004
      SUBROUTINE OTTODD(D,YOUN,XNU,G,GS,XINVL,W)
C=========================================================================
C
C  ENTREES :
C    YOUN,XNU,G,GS,MFR
C    XINVL,W
C
C  SORTIES :
C    D
C
C==========================================================================
C
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
 
-INC PPARAM
-INC CCOPTIO
C
      PARAMETER (XZER=0.D0,UN=1.D0,DEUX=2.D0)
C
      DIMENSION W(3),XINVL(3),FK(3),D(18),GS(3)
*
      FK(1)=XINVL(1)*W(1)/GS(1)+XINVL(2)*W(2)/GS(2)
      FK(2)=XINVL(1)*W(1)/GS(1)+XINVL(3)*W(3)/GS(3)
      FK(3)=XINVL(2)*W(2)/GS(2)+XINVL(3)*W(3)/GS(3)
*
      DET=(UN+XNU)*(UN+XNU)*(UN-DEUX*XNU)/(YOUN*YOUN*YOUN)
*
      D(1)=(UN-XNU*XNU)/(YOUN*YOUN)/DET
      D(2)=(UN-XNU*XNU)/(YOUN*YOUN)/DET
      D(3)=(UN-XNU*XNU)/(YOUN*YOUN)/DET
      D(4)=(XNU*(UN+XNU)/(YOUN*YOUN))/DET
      D(5)=(XNU*(UN+XNU)/(YOUN*YOUN))/DET
      D(6)=(XNU*(UN+XNU)/(YOUN*YOUN))/DET
      D(7)=G/(UN+G*FK(1))
      D(8)=G/(UN+G*FK(2))
      D(9)=G/(UN+G*FK(3))
      D(10)=XZER
      D(11)=XZER
      D(12)=XZER
      D(13)=XZER
      D(14)=XZER
      D(15)=XZER
      D(16)=XZER
      D(17)=XZER
      D(18)=XZER
*
*      IF(IIMPI.EQ.42) THEN
*         WRITE(IOIMP,77100) DET,D(1),D(6),D(5)
*77100    FORMAT(2X,'OTTODD  DET=',1PE12.5,2X,'D(1)=',1PE12.5,2X,
*     &          'D(6)=',1PE12.5,2X,'D(5)=',1PE12.5/)
*      ENDIF
*
      RETURN
      END