grfltHF1.f : Seismic Wave using Homogeneous Green's function and randam Vs. 
             Gassian integral formula is used for the fault integration.
grfltHF2.f : Seismic Wave using Homogeneous Green's function and randam Vs. 
             The reqtangular formula is used for the fault integration.
	     Time Delay =dis/VR, and VR=VRavg+VRstd*dev
grfltHF3.f : Same as grfltHF2.f, but except 
	     Time Delay =delay+dis/VRavg+TMstd*dev
grfltHF4.f : Same as grfltHF3.f, but except corrections for the horizontal 
             and veritical amplitudes to adjust with those of excact half-space
             Green's functions
grfltHF5.f : Seismic Wave using Homogeneous Green's function and randam Vs.
             For randam Vs, the local Vs is introduced on sub sub-faults. 
grfltHF6.f : Distributions of Boore's point source.

For the 85 Mexico, 85 Chile, and Tokai quakes grfltHF5.f and grfltHF6.f for low
and high frequencies, respectively.          

grfltHF1a.f : Seismic Wave using Homogeneous Green's function and randam Vs. 
             Gassian integral formula is used for the fault integration.
           * Complete Green's Function or 1/R terms can be chosen (2004/8/16).