Esope source files

   1 : C ENDO31    SOURCE    CHAT      05/01/12    23:38:53     5004           
   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 =( DABS((WWY-D31)**ann) +(WWY-D31)**ann )/2.0D0 
  16 :         CC2 =( DABS((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 = DABS(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 :