Numérotation des lignes :

C OTTODD    SOURCE    CHAT      05/01/13    02:07:21     5004      SUBROUTINE OTTODD(D,YOUN,XNU,G,GS,XINVL,W)C=========================================================================CC  ENTREES :C    YOUN,XNU,G,GS,MFRC    XINVL,WCC  SORTIES :C    DCC==========================================================================C      IMPLICIT INTEGER(I-N)      IMPLICIT REAL*8(A-H,O-Z)-INC CCOPTIOC      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

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