Esope source files

   1 : C ETACFD    SOURCE    CHAT      05/01/12    23:45:46     5004           
   2 : c
   3 :              FUNCTION ETACONFUNCDER(VINT)
   4 :       IMPLICIT INTEGER(I-N)
   5 :       real*8    vint,etaconfuncder
   6 :       real*8    exm,qaa,pcc,eww,xk1,xk2,etb,exv,xep,paa,cco,
   7 :      .          exl,xnn,aaa,bbb
   8 :       common /MRScone/ exm,qaa,pcc,eww,xk1,xk2,etb,exv,xep,paa,cco,
   9 :      .                 exl,xnn,aaa,bbb
  10 :       if (exv.lt.0.D0) then
  11 :         etaconfuncder=0.D0
  12 :       else
  13 :         if (vint.lt.0.D0) then
  14 :           write(*,*)' ETACONDER var int',vint
  15 :           vint=0.D0
  16 :         endif
  17 :         etaconfuncder=-aaa*bbb*(vint+xk1)**(1.D0/exv)*dexp(-bbb*vint)
  18 :      .                -etb*xk2*vint/(xep+vint)**2
  19 :      .                +etb*xk2/(xep+vint)+aaa
  20 :      .                *(vint+xk1)**(1.D0/exv-1.D0)*dexp(-bbb*vint)/exv
  21 :         endif
  22 :       return
  23 :       end
  24 :  
  25 :  
  26 :