*** grflt01p.f + grflt00p1.f (grflt01p.in as the input data) (Aug.22, 2000) *** This code is modified from grflt00p.f as follows 1) The exponential source is included (when, Source Type =3) 2) Integration points are added around the branch points 3) The integration range are changed to the two from the original four 4) Phinney's integration method are added by including the imaginay omega 5) For r*k > 8, the approximate Bessel function are improved including higher order expansion functions 6) A dynamic distribution of source points is used on the fault planes, in which the points per wavelength are automatically assigned. For FFT, grfftp.f must be used to obtain correct time-domain amplitude *** grflt00p.f + grflt00p1.f (grflt00p.in as the input data) (Aug.11, 2000) *** This code is modified from grflt95.f as follows 1) The exponential source is included (when, Source Type =3) 2) Integration points are added around the branch points 3) The integration range are changed to the two from the original four 4) Phinney's integration method are added by including the imaginay omega 5) For r*k > 8, the approximate Bessel function are improved including higher order expansion functions For FFT, grfftp.f must be used to obtain correct time-domain amplitude *** grflt00.f + grflt002.f (grflt00.in as the input data) (Aug.10, 2000) *** This code is modified from grflt95.f as follows 1) The exponential source is included (when, Source Type =3) 2) Integration points are added around the branch points 3) The integration range are changed to the two from the original four 4) For r*k > 8, the approximate Bessel function are improved including higher order expansion functions *** grflt95p.f + grflt95p1.f (grflt95e.in as the input data) (Aug.9, 2000) *** This code is modified from grflt95.f as follows 1) The exponential source is included (when, Source Type =3) 2) Phinney's integration method are added by including the imaginay omega 3) For FFT, grfftp.f must be used to obtain correct time-domain amplitud *** grflt95e.f + grflt952.f (grflt95p.in as the input data) (Aug.4, 2000) *** This code is modified from grflt95.f as follows 1) The exponential source is included (when, Source Type =3)