Sources Esope | Cast3M

pcapfd
C PCAPFD    SOURCE    CHAT      05/01/13    02:11:09     5004
c
             FUNCTION PCAPFUNCDER(VINT)
      IMPLICIT INTEGER(I-N)
      real*8    vint,pcapfuncder
      real*8    pc0,exr,cca,phi,alp,pcc,paa,qaa,exm,eww,
     .          cpr,cpm,cfr
      common /MRScapp/ pc0,exr,cca,phi,alp,pcc,paa,qaa,exm,eww,
     .                 cpr,cpm,cfr
      if (exr.lt.0.D0) then
            pcapfuncder=0.D0
      else if (vint.eq.0.D0) then
        pcapfuncder=pc0*1.D-8**(1.D0/exr-1.D0)/exr
      else
        pcapfuncder=pc0*vint**(1.D0/exr-1.D0)/exr
      endif
      return
      end