Télécharger ottodd.eso

Retour à la liste

Numérotation des lignes :

  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.  
  63.  

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