Télécharger la plaquette Cast3M


fr en

Télécharger endo31.eso

Retour à la liste

Numérotation des lignes : 

endo31

  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