Télécharger ottofl.eso

Retour à la liste

Numérotation des lignes :

  1. C OTTOFL SOURCE CHAT 05/01/13 02:07:35 5004
  2. SUBROUTINE OTTOFL(NCA,NN,XINVL,PENTE,SIGMA,GS,SMAX,XLTR,
  3. & DFDS,DGDS,HDFDQ,XCOMP,XLAMC,PRECIE,
  4. & PRECIZ,KERRE)
  5. C----------------------------------------------------------------------
  6. C
  7. C ENTREES : NCA,NN,XINVL,PENTE,SIGMA,GS,SMAX
  8. C
  9. C SORTIES : DFDS,DGDS,HDFDQ
  10. C
  11. C-----------------------------------------------------------------------
  12. C
  13. IMPLICIT INTEGER(I-N)
  14. IMPLICIT REAL*8(A-H,O-Z)
  15. -INC CCOPTIO
  16. *
  17. PARAMETER (XZER=0.D0)
  18. *
  19. DIMENSION NN(*),DFDS(6,*),DGDS(6,*),HDFDQ(*)
  20. DIMENSION XINVL(3),PENTE(*),SIGMA(*),GS(*),XLTR(*)
  21. DIMENSION SS(6)
  22. DIMENSION SMAX(*),XCOMP(*)
  23. DIMENSION IV1(3),IV2(3)
  24.  
  25. DIMENSION DFF(6),DGG(6)
  26. *
  27. DATA IV1/4,4,5/
  28. DATA IV2/5,6,6/
  29. *
  30. KERRE=0
  31. *
  32. *
  33. DO IC=1,NCA
  34. JC=NN(IC)
  35. DO I=1,6
  36. DFDS(I,JC)=XZER
  37. DGDS(I,JC)=XZER
  38. ENDDO
  39. HDFDQ(JC)=XZER
  40. *
  41. GO TO (1,1,1,4),JC
  42. *
  43. KERRE=99
  44. RETURN
  45. *
  46. 1 CONTINUE
  47. *
  48. IC1=IV1(JC)
  49. IC2=IV2(JC)
  50. *
  51. IF (XINVL(JC).NE.XZER) THEN
  52. IF(IIMPI.EQ.42) THEN
  53. WRITE(IOIMP,77005) XINVL(JC),PENTE(JC)
  54. 77005 FORMAT(2X,'XINVL=',1PE12.5,2X,'PENTE=',1PE12.5/)
  55. ENDIF
  56. *
  57. *
  58. H =PENTE(JC)/XINVL(JC)
  59. IF(IIMPI.EQ.42) THEN
  60. WRITE(IOIMP,77009) H
  61. 77009 FORMAT(2X,' MODULE D ECROUISSAGE H= ',1PE12.5/)
  62. ENDIF
  63. HDFDQ(JC)=-H
  64. DFDS(JC,JC)=1.D0
  65. DGDS(JC,JC)=1.D0
  66. DGDS(IC1,JC)=SIGMA(IC1)/GS(JC)/2.D0
  67. DGDS(IC2,JC)=SIGMA(IC2)/GS(JC)/2.D0
  68. ENDIF
  69. GO TO 100
  70. *
  71. 4 CONTINUE
  72.  
  73.  
  74. CALL OTTOCP(SIGMA,FCR4,XLTR,DFF,DGG,H4,
  75. & PRECIE,PRECIZ,XCOMP,XLAMC,KERRE)
  76.  
  77. DO I=1,6
  78. DFDS(I,4)=DFF(I)
  79. DGDS(I,4)=DGG(I)
  80. ENDDO
  81. *
  82. HDFDQ(4)=H4
  83.  
  84.  
  85.  
  86. *
  87. * DFDS(1,4)=0.D0
  88. * DGDS(1,4)=0.D0
  89. * HDFDQ(4)=0.D0
  90. GO TO 100
  91. *
  92. 100 CONTINUE
  93. ENDDO
  94.  
  95. RETURN
  96. END
  97.  
  98.  
  99.  
  100.  
  101.  

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