Esope source files

   1 : C FONCG     SOURCE    CHAT      05/01/13    00:06:05     5004           
   2 : C
   3 : C*********************************************************************
   4 : C
   5 :       FUNCTION FONCG(Q,X,Y,Z,X0)
   6 :       IMPLICIT INTEGER(I-N)
   7 :       IMPLICIT REAL*8(A-H,O-Z) 
   8 :       EPSIL2=1.D-8
   9 :       D=X-X0
  10 :       IF(D.EQ.0.) THEN
  11 :          X=X+EPSIL2
  12 :       ENDIF
  13 :       DELTA=SQRT((X-X0)**2+Y**2+(Q*X+Z)**2)
  14 : C****
  15 :       DELTA1=MAX(DELTA-(Q*X+Z),EPSIL2)
  16 :       DELTA3=MAX(DELTA+(Q*X+Z),EPSIL2)
  17 :       DELTA2=MAX(Y+DELTA,EPSIL2)
  18 :       FONCG=0.5*Y*LOG(DELTA1/DELTA3)
  19 :      &     -(Q*X+Z)*LOG(DELTA2)
  20 :      &     +(X-X0)*ATAN(Y*(Q*X+Z)/(X-X0)/DELTA)
  21 : C****
  22 :       RETURN
  23 :       END
  24 :  
  25 :