Esope source files

   1 : C OTTODD    SOURCE    CHAT      05/01/13    02:07:21     5004           
   2 :       SUBROUTINE OTTODD(D,YOUN,XNU,G,GS,XINVL,W)
   3 : C=========================================================================
   4 : C
   5 : C  ENTREES :
   6 : C    YOUN,XNU,G,GS,MFR
   7 : C    XINVL,W
   8 : C
   9 : C  SORTIES :
  10 : C    D
  11 : C
  12 : C==========================================================================
  13 : C
  14 :       IMPLICIT INTEGER(I-N)
  15 :       IMPLICIT REAL*8(A-H,O-Z)
  16 : -INC CCOPTIO
  17 : C
  18 :       PARAMETER (XZER=0.D0,UN=1.D0,DEUX=2.D0)
  19 : C
  20 :       DIMENSION W(3),XINVL(3),FK(3),D(18),GS(3)
  21 : *
  22 :       FK(1)=XINVL(1)*W(1)/GS(1)+XINVL(2)*W(2)/GS(2)
  23 :       FK(2)=XINVL(1)*W(1)/GS(1)+XINVL(3)*W(3)/GS(3)
  24 :       FK(3)=XINVL(2)*W(2)/GS(2)+XINVL(3)*W(3)/GS(3)
  25 : *
  26 :       DET=(UN+XNU)*(UN+XNU)*(UN-DEUX*XNU)/(YOUN*YOUN*YOUN)
  27 : *
  28 :       D(1)=(UN-XNU*XNU)/(YOUN*YOUN)/DET
  29 :       D(2)=(UN-XNU*XNU)/(YOUN*YOUN)/DET
  30 :       D(3)=(UN-XNU*XNU)/(YOUN*YOUN)/DET
  31 :       D(4)=(XNU*(UN+XNU)/(YOUN*YOUN))/DET
  32 :       D(5)=(XNU*(UN+XNU)/(YOUN*YOUN))/DET
  33 :       D(6)=(XNU*(UN+XNU)/(YOUN*YOUN))/DET
  34 :       D(7)=G/(UN+G*FK(1))
  35 :       D(8)=G/(UN+G*FK(2))
  36 :       D(9)=G/(UN+G*FK(3))
  37 :       D(10)=XZER
  38 :       D(11)=XZER
  39 :       D(12)=XZER
  40 :       D(13)=XZER
  41 :       D(14)=XZER
  42 :       D(15)=XZER
  43 :       D(16)=XZER
  44 :       D(17)=XZER
  45 :       D(18)=XZER
  46 : *
  47 : *      IF(IIMPI.EQ.42) THEN
  48 : *         WRITE(IOIMP,77100) DET,D(1),D(6),D(5)
  49 : *77100    FORMAT(2X,'OTTODD  DET=',1PE12.5,2X,'D(1)=',1PE12.5,2X,
  50 : *     &          'D(6)=',1PE12.5,2X,'D(5)=',1PE12.5/)
  51 : *      ENDIF
  52 : *
  53 :       RETURN      
  54 :       END
  55 :  
  56 :  
  57 :  
  58 :  
  59 :  
  60 :  
  61 :  
  62 :