Télécharger endo31.eso

Retour à la liste

Numérotation des lignes :

  1. C ENDO31 SOURCE CB215821 16/04/21 21:16:36 8920
  2. SUBROUTINE ENDO31(D3,D3R,WWY,AKK,ANN,DT)
  3.  
  4. C Resolution Newton equation D3
  5. IMPLICIT INTEGER(I-N)
  6. REAL*8 D3,D3R,WWY,akk,ann,DT
  7. REAL*8 D31,PRE ,CC1,CC2,RI1,DRI1D,DD3,TEST
  8. INTEGER IK
  9.  
  10. D31= D3R
  11. PRE = 1.0D-6
  12. IK=1
  13. 31 CONTINUE
  14.  
  15. CC1 =( ABS((WWY-D31)**ann) +(WWY-D31)**ann )/2.0D0
  16. CC2 =( ABS((WWY-D31)**(ann-1)) +(WWY-D31)**(ann-1) )/2.0D0
  17. RI1= D31-D3R-AKK*DT*(CC1)
  18. DRI1D= 1.0D0 +AKK*DT*(CC2)
  19. DD3 = -RI1/ DRI1D
  20. D31=D31 + DD3
  21.  
  22. C Convergence Newton
  23.  
  24. Test = ABS(DD3)
  25. IF(Test.le.PRE) THEN
  26. D3 = D31
  27. C PRINT*,' Convergence Newton d3 Iteration',IK
  28. C PRINT*,'ENDO D3',D3
  29. RETURN
  30.  
  31. ELSE
  32.  
  33. IK =IK +1
  34.  
  35. IF(IK.EQ. 100) THEN
  36.  
  37. C Print*,'Pas de convergence IKMAX',IK
  38.  
  39. D3 = D3R
  40. RETURN
  41. ELSE
  42. C Iterations Newton
  43. GOTO 31
  44.  
  45. ENDIF
  46.  
  47. ENDIF
  48. END
  49.  
  50.  
  51.  
  52.  
  53.  

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