pdisp

Fortran source

0001 C
0002 C   Calculation of the plasma dispersion function
0003 C
0004 C   Requires the use of IMSL library
0005 C
0006 C   by A.K. Ram and J. Decker 10/26/05
0007 C
0008       subroutine pdisp(rnar,cargr,cargi,zr,zi)
0009 C
0010       implicit double precision (a-h, o-z)
0011 C
0012       REAL*8 rnar,cargr(0:1000),cargi(0:1000),zr(0:1000),zi(0:1000)
0013 C
0014       double complex dcmplx
0015       double complex csqpi,zerfe
0016 C
0017       double complex carg, z
0018 C
0019       external zerfe
0020 C
0021       sqpi=1.772453850905516d0
0022         pi=3.141592653589793d0
0023 C
0024       csqpi=dcmplx(0.d0,sqpi)
0025       nar = int(rnar)
0026 C     
0027       do 1 i=0,nar
0028         carg = dcmplx(cargr(i),cargi(i))
0029         z=csqpi*zerfe(carg)      
0030         zr(i) = dreal(z)
0031         zi(i) = dimag(z)
0032     1 continue
0033 C
0034       return
0035       end