Esope source files

   1 : C LEMAIT    SOURCE    CHAT      05/01/13    01:15:12     5004           
   2 :       SUBROUTINE LEMAIT(SIG,X,DX,DT,XMAT,DEPCT)                                 
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL *8 (A-H,O-Z)                                                
   5 :       DIMENSION XMAT(*)                                                         
   6 :       DX= (SIG/XMAT(9 ))**(XMAT(8)/XMAT(7))*XMAT(11)**(1.0D0/XMAT(7))           
   7 :       DY= (SIG/XMAT(10))**XMAT(8)*XMAT(11)                                      
   8 :       XX=X                                                                      
   9 :       T1=0.5D0*DX*DT                                                            
  10 :       IF(XX.LT.T1) XX=T1                                                        
  11 :       DEPCT=XMAT(6)*(XMAT(7)*XX**(XMAT(7)-1.0D0)*DX+DY)                         
  12 :       RETURN                                                                    
  13 :       END                                                                       
  14 :