Hi Francesco,

Please take note that I'm just a novice in using Hypoellipse. So my description
on the way the program is used may be questionable.   Anyhow, here goes.

It should be noted the file names can be different from the one's used here.
(I'll call them the same as the manual file names.)

1)
The first thing  is to copy all the "DOS" and "Documents" files form
Dr. John Lahr's great site.

http://lahr.org/iaspei/northam/ak/s_ak/programs/hypoel/hypoel.html

Its best to place each individual section in a separate folder. Its important for this
example to copy all the "Quickrun" files to a new "quickrun folder". Make a backup
of this folder.

2) Print out the Quick-Start manual.

3) (Read chapter on "crustal.prm" )
In the "Quickrun folder" load the file called "crustal.prm" into a text editor. You
need to modify this "VELOCITY" file for your location. I'll give you the data I have
used but you will need to talk to Dr. Lahr or you local Geophysical Observatory and
get the velocity of the "P"  and "S" waves at various depths around you location.

In My area its  6.13 Km/sec for a P wave at    0 Km  Depth.
                7.14 Km/sec for a P wave at   19 km  Depth
                8.27 Km/sec for a P wave at 36.0 km  Depth

In My area its  3.62 Km/sec for a S wave at    0 Km  Depth.
                3.96 Km/sec for a S wave at   19 km  Depth
                4.75 Km/sec for a S wave at 36.0 km  Depth

Now calculate the (Vp/Vs) data for each depth.

                0 Km Depth is (6.13/3.62) = 1.69
               19 Km Depth is (7.14/3.96) = 1.80
               36 Km Depth is (8.27/4.75) = 1.74

Your model may have different depths but the principle is the same.
Using the information from the "quickrun" manual modify the opened
file "crustal.prm" to reflect the velocity model data. Note the column
spacing. It Should look someting like this for my area.


! Model 1: Western Australia
!
VELOCITY           6.13       0.0      1.69
VELOCITY           7.14      19.0      1.80
VELOCITY           8.27      36.0      1.74

Close and save the file.

4) Input Data. Its important to follow the manual and watch your column spacing.
For each seismic station you should calculate, using "Winquake" an arrival time
of the "P" wave or the "S" wave or both waves. Accuracy, depends on how good your
data is and how many stations you have. Not less than three stations.

I used my station plus two AutoDrm stations.

BICK station  P arrival is 1:46:38.5  On 99/07/13
              S arrival    1:46:56.6

FORT station  P arrival is 1:48:22.4  On 99/07/13
              S arrival    1:49:59.5  - Very  Poor data

MEEK station  P arrival is 1:47:22.9  On 99/07/13
              S arrival    1:48:14.8

I also calculated the amplitude of the wave data, but its not used in this
example. Magnitude is a bigger problem and wont be considered in this example.

"Winquake" gives a very good estimate on the magnitude when the right corrections
are entered.

Amplitude   BICK  5660
            FORT   376
            MEEK  1319


In the "Quickrun folder" load the file called "akutan.pha" into a text editor. You
need to modify this file for your Input Data. Again I'll give you the data from my
station plus the two AutoDrm stations. Armed with the arrival time data and your
"quickrun" manual;  the modified "akutan.pha" file should look something like this.

BICKIP   9907130146 3850          5660        4121
FORTEP   9907130148 2240         11950         376
MEEKEP   9907130147 2290          7480        1319
C*       YrMoDyHrMn P-Sec         S-Sec        Amp
C*(F) 960423065129i.
C*(C) Any record starting with C* is a comment.

Make sure the columns are correct. The printout here may not reflect the correct
column spacing.

Close and save the file. You have now set the velocity Model and input data for
a particular quake.

5) Station Locations.

In the "Quickrun folder" load the file called "akutan.sta" into a text editor. You
need to modify this file for your station Input Data. Again I'll give you the data
from my station plus the two AutoDrm stations.

BICK   -32.008 Latitude by 116.135E Longitude 300 Meteres above Mean sea level.
FORT   -30.778          by 128.059E           165
MEEK   -26.638          by 118.615E           530

So BICK is -32 degrees 00.48 minutes south by 114 degees 08.10 minutes east
   FORT    -30         46.68                  128        03.54
   MEEK    -26         38.28                  118        36.90


Armed with the station location's and your "quickrun" manual. The "akutan.sta" file
should look something like this.

BICK32S00.48 116E08.10  300
BICK*    18
FORT30S46.68 128E03.54  165
FORT*    18
MEEK26S38.28 118E36.90  530
MEEK*    18

Save and close the file.

6) Seismograph amplitude response.
I entered this data since I had it available, its need for magnitude calculation.
I'm going to do this later.

In the "Quickrun folder" load the file called "caldata.prm" into a text editor. You
need to modify this file for your seismograph amplitude response. Again I'll give
you the data from my station plus the two AutoDrm stations. The file should look like
this. Check the manual for column layout.

BICK D 980713 999999  1 1.00 400000.
FORT D 980713 999999  1 1.00 591716.
MEEK D 980713 999999  1 1.00 680272.

Save and close the file.

7) Options. Hypoellipse allows many, many options is the way it process data.

I changed these parameters so it approached the reduction in the way I like
it to. You will need to printout the full manual to see the options.  See
page 2-23 "Big manual". Maybe for this exercise just print out the relevant
part of the manual. (2-23 -> 2.29)

In the "Quickrun folder" load the file called "headopts.vol" into a text editor. You
need to modify this file for your various options. In this case I'll give my selection
and why. I'm a novice in all of this so shot me down, with compassion.

Options and Why

! Options - control file
header content       Bickley Western Australia    ! Run Name
reset test        1   1.69    !Ratio P/S wave velocity   **
reset test        2   6.13    !P wave velocity           **
reset test        3   -31.0   !First Trial Latitude      **
reset test        4   -117.0  !First Trial Longitude     **
reset test        5   -1.0    !First trial depth
reset test        6   -1.0    !RMS See page 2-24
reset test        8    0.0    !Top of Model, respect to sea level **
reset test        21   20.    !Maximum nuber of iteration allowed (important)
reset test        22   1.0    !Limit chang of focal depth (km)   **
reset test        31   -0.87  !C1 constant (magnitude)           **
reset test        32    2.0   !C2 constant (magnitude)           **
reset test        33    0.0035!C3 constant (magnitude)           **
! end of headopts.vol

NOTE: the Lines ending with "**" aren't really needed in this epicentre location example.

Anyhow the new "headopts.vol" should look like this.

! Options - control file
header content       Bickley Western Australia
reset test        1   1.69
reset test        2   6.13
reset test        3   -31.0
reset test        4   -117.0
reset test        5   -1.0
reset test        6   -1.0
reset test        8    0.0
reset test        21   20.
reset test        22   1.0
reset test        31   -0.87
reset test        32    2.0
reset test        33    0.0035
! end of headopts.vol

Save and close the file.



8) Input control file. (Hypo.ctl) This file controls which files are used and controls.

In the "Quickrun folder" load the file called "hypo.ctl" into a text editor. You
need to modify this file for your file name and options. (Example)
It should look something like this.

! headopts.vol contain the setup parameters
! for running HYPOELLIPSE.
jump headopts.vol
!
! crustal.prm specifies the velocity model.
jump crustal.prm
!
! caldata.prm contains the calibration parameters.
uofacal option    caldata.prm
!
! Constants noprint = 1 will cause documentation of the
! parameter values, crustal model, station locations, etc
! to be added to the output (.out) file.  Note that if the
! same parameter value it specified more than once, the last
! setting will be the one used (in this case the value will
! be set to 1).
constants noprint    0
constants noprint    1
!
! Printer option 1 adds a blow by blow description of
! every iteration step and is useful for debugging purposes.
! Reverse the order of the following records to turn this
! option on.
printer option       1
printer option       0
tabulation option    4
!
begin station list +1 19960101
jump akutan.sta
arrival times next
jump akutan.pha


Save and close the file.

You will need to consult the "Big Manual" for a full explenation of this file.

                        --------------------------------------------

The files (akutan.sta, crustal.prm, caldata.prm, and headopt.vol ) once set shouldn't
need to be changed for other events.

9)
Copy the file "Hypoel.ex" into your quickrun folder and rename it "Hypoel.exe".
At this point I recommend to close all other windows except the "quickrun" folder.


10) Lets give it a go!!

Open an "Msdos" Window. I recommend a window and not the full screen. See "MsDos" prompt,
properties. In the "MsDos" window change the directory until you are in the "quickrun"
folder. At the prompt, type.

hypoel < hypo.in                        (and hit the return key)

If all goes well you should see some new files in the "quickrun" folder. Look at the
"Hypoel.out" file with a text editor and see the program process and  reduction.
Hopefully you should see the quakes location at

99/07/13  01:46:13.6 UT  -31.21 latitude (-31 degrees 12.6 minutes)
                         117.48 East     (117 degrees 28.8 minutes)
                         3 km Depth.

Or something close to this "A" accuracy position. More data means more accuracy!

If you want to re-compute this run with some modifications then you will need to delete
the previously generated output files. Go to the "MsDos" window and repeat the line

hypoel < hypo.in                        (and hit the return key)


Francesco, I hope this is of help and I'll send you these files in a  "ZIP".

Again I wish to thank Dr. John Lahr for this great resource.


Have Fun. I did.

Arie




_____________________________________________________________________

Public Seismic Network Mailing List (PSN-L)