* boore1.f: This program is not complete (see boore2.f, 2004/8/11)
 
* boore2.f: Boore's Moment Rate Function (2004/8/12) 
     1) Boore's envelope function is applied to moment rate (velocity
        in far field), but not to acceleration -> boore-env.csv
     2) Random phases are used at the frequencies higher than Fran (Hz)
        -> phase.csv
     3) 10 times iterations are carried out to fit Boore's omega
        squared model -> fourier.csv
     4) Accelaration, velocity, displacement (moment rate funtion), and
        moment function are stored in wave.csv -> wave.csv
     5) Input data (dt, nt, SMo, fc, fmax, idum, fran, and Niter) are all
        writen in the program.
    -> fran shuld be close to fc

* boore2a.f: Boore's Moment Rate Function (2004/8/12) 
     1) Boore's envelope function is applied to moment rate (velocity
        in far field), but not to acceleration -> boore-env.csv
     2) Random phases are used at the frequencies higher than Fran (Hz)
        -> phase.csv
     3) 10 times iterations are carried out to fit Boore's omega
        squared model -> fourier.csv
     4) Accelaration, velocity, displacement (moment rate funtion), and
        moment function are stored in wave.csv -> wave.csv
     5) Input data (dt, nt, SMo, fc, fmax, idum, fran, and Niter) are all
        writen in the program.
     6) fc is calculated from Brune's formula
    -> fran shuld be close to fc

* boore3.f: Checking program in Homogeneous Full Space (2004/8/13) ***
     1) Using Green's funtion of homogeneous full-space, strong ground
        motions are simulated. This is modified from grfltHF1.f.
     2) Source slip velocity function is a triangle, but not Boore's
        source function for checking purposes.
     3) Green's function can be the complete function or the far field
        function (=1/r terms).
    -> Input data are boore3.in, and the output data is grfault.dxyz

* boore4.f: Boore's Moment Rate Function (2004/8/19) ***
     1) Using Green's funtion of homogeneous full-space, strong ground
        motions are simulated. This is modified from boore3.f
       (from grfltHF1.f) and boore2a.f.
     2) Green's function can be the complete function or the far field
        function (=1/r terms).
     3) Boore's envelope function is applied to moment rate (velocity
        in far field), but not to acceleration -> boore-env.csv
      (Modified: Actually, commented out in this program)
     4) Random phases are used at the frequencies higher than Fran (Hz)
        -> phase.csv
     5) Iterations are carried out to fit Boore's omega
        squared model -> fourier.csv
     6) Accelaration, velocity, displacement (moment rate funtion), and
        moment function are stored in wave.csv -> wave.csv
     7) Input data for Boore's source model, fc and fmax, and the other
        parameters, idum, fran, and Niter are read from input data
        'boore4.in'.
     8) Strong motions of sub faults are superposed by amplification
        factors of one at lower frequencies, and sqrt(NL*NW) at higher
        frequencies.
     9) To obtain wave forms, a fft pronram (e.g., grfftp.f) is needed.


* boore5.f : Boore's Moment Rate Function (2004/8/20) ***
    Same as boore4.f, but  
     10) Isotropic radiation patterns are used for S and P waves at
        high frequencies.