Esope source files

   1 : C PARTAN    SOURCE    CHAT      05/01/13    02:09:34     5004           
   2 :       SUBROUTINE PARTAN(IVLAMB,IVPR,IVGE,IVGM) 
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8 ( A-H,O-Z)                                                                                                                    
   5 : -INC CCOPTIO                                                                    
   6 : -INC SMLREEL
   7 :       POINTEUR MLREE4.MLREEL,MLREE5.MLREEL    
   8 :       MLREEL=IVLAMB
   9 :       MLREE1=IVPR
  10 : 
  11 :       MLREE5=IVGE
  12 :       SEGACT MLREE1,MLREEL,MLREE5*mod
  13 :       JG = PROG(/1)                               
  14 :       DO 2 I=1,JG
  15 :         MLREE5.PROG(I)=PROG(I) - MLREE1.PROG(I) 
  16 :    2  CONTINUE                  
  17 : *                                                                               
  18 : *  CALCUL DE LA DIRECTION DE MONTEE on verifie en fait que la 
  19 : * direction calculee precedemment est compatible avec le fait que 
  20 : * lambda doit rester positif
  21 : * 
  22 :       MLREE4 = IVGM
  23 :       SEGACT MLREE4*mod
  24 :       DO 20 I=1,JG
  25 :          MLREE4.PROG(I) = MLREE5.PROG(I)
  26 :    20 CONTINUE                                                                                                                               
  27 :         DO 3 I=1,JG                                                                                                     
  28 :             IF(MLREE5.PROG(I).LT.0.D0) THEN                                     
  29 :               IF(PROG(I).EQ.0.) THEN                                            
  30 :                 MLREE4.PROG(I)=0.D0                                             
  31 :               ENDIF                                                             
  32 :             ENDIF                                                               
  33 :     3   CONTINUE
  34 :          SEGSUP MLREEL
  35 :          IVLAMB=MLREE1                                                 
  36 :       RETURN                                                                    
  37 :       END                                                                       
  38 :  
  39 :