'#' and format is free.     Meaning of control parameters is as follows:	ltdo, oltm, lndo,olnm	  : origin of short distance conversion 				    [deg. min.], E and S are negative	x,y			  : x,y coordinates of lower left corner				    (= origin of numerical grid)			            rel. to short distance origin 				   	baselev			  : top of model in km above sea level	h			  : grid spacing of numerical grid	xrad,yrad,zrad		  : initial search radius [km] for grid search	neqs			  : no. of earthquakes	nitmax, fatray		  : >= 1  no. inversion iterations				     0 = do only eq. relocation				    -1 = compute synthetic traveltimes only        				 if fatray = 1 fat rays are computed					 if fatray = 0 only tt-times						ifxl			  : fix earthquakes for ifxl iterations				 	vthet			  : damping of velocity parameters		delta_t			  : fat ray width [s]. should be set as				    T/2, where T is period of the dominant                     		    frequency				   	hitct			  : inversion blocks with less than hitct fat				    rays will be fixed in the inversion			    	ires			  : = 0 no resolution				    = 1 only diagonal elements of the resolution				      matrix (RDE) will be computed				    = 2 full resolution matrix is computed	rout,fout,rmscubes	  : output parameters for rays/fat rays				    rout = 1 print out ray points  				    fout = 1 print out fat rays 				    rmscubes =1 RMS volumes of earthquake				    relocation      	  STATION INPUT (file STNS)      Station input format is identical to SIMULPS13Q with one exception:      Origin for short distance conversion is given via the control file.     File must be named STNS. At first comes the number of stations (free      format) followed by a station card for each station (one station per line).     Station card consists of (format:2x,a4,i2,a1,f5.2,i4,a1,f5.2,f5.2):      stat. name, latitude, longitude, elevation [m], Vp station correction         Example station file         22       AME23S55.00  69W45.49  850 0.00       MEX23S17.08  70W 3.77 1200 0.00       PRT23S27.41  70W13.52  851 0.00       ERC23S 8.38  69W52.45 1368 0.00       CAC22S59.95  70W18.30  100 0.00       ....         EVENT INPUT (file EQKS)      Phase data input format is identical to SIMULPS13Q with no exception.     The file must be named EQKS. Each event contains:           event card (format:a4,a2,1x,a2,i2,1x,f5.2,i3,a1,f5.2,1x,i3,a1,f5.2,2f7.2)       iyrmo,iday : year, month, and day       ihr, mino, seco : hour, minute, and second       ltde, cns, eltm : latitude       lnde, cew, elnm : longitude       dep : depth [km]       rmag : magnitude      observed traveltimes to station, six sets per line (format:2a4,f6.2)       sta : station name       rmki : phase typ and observation weight       tt : travel time [s]      blank line            Example event input:         950818 0721 49.46 23S23.12  70W21.87  37.50  2.20      6      0.00      AME P 0 14.40 MEX P 0  7.86 PRT P 0  6.44 ERC P 0 10.76 ...      MJC P 0  6.72 CAC P 0  8.86 PEL P 0 10.75 QC1 P 0  7.20 ...      CRP P 0 18.88 AMA P 0 17.07 VEN P 1 13.02 CER P 0 21.62 ...      DES P 0  5.98 MJS P 0  6.55 MJP P 0  6.82 MJN P 0  8.06 ...     950818 0947 19.22 23S52.94  70W14.48  53.82  2.50     11      0.00      AME P 0 10.94 MEX P 1 13.02 PRT P 0 10.71 COL P 0  8.56 ...      PEL P 0....		--------------------------------------  OUTPUT files    printed output (file 'output')      Fatomo creates a file named 'output' for the main printout. Among the      control parameters and initial velocities, it contains the inversion       results. Model pertubations (% change of slowness), HIT-, KHIT-, DWS-      matrix, and corrected (final) velocities are arranged in planes starting      with the top plane and the UPPER LEFT corner of the model. So, print       out of the results is arranged exactly in a way as it is expected, with      north at the upper model boundary.       Model pertubations and final velocities are given in the SEISMIC grid,       HIT-, KHIT-, DWS-matrix, and RDEs (diagonal elements of the resolution      matrix) in the INVERSION grid.       Using TOMO2GMT (available via www.sg.geophys.ethz.ch/aes/husen) output of       FATOMO can be reformatted to be used with GMT-software.   final velocity model (file 'f23')      The corrected (final) velocity model after the inversion is printed in       the file 'f23'. The format is the same as the model-input file 'MOD'.      Hence, the new velocities can be used as initial ones for subsequent      inversions.   synthetic travel times (file 'f24')      Relocated hypocenter location are printed to file 'f23' in the same format      as the file 'EQKS'.      If synthetic travel times are computed (nitmax=-1) they are output      in the file 'f24' in the same format as in the file 'EQKS'.   mapping of fixed inv. blocks (file 'nodes.out')      Mapping from the full grid (all inv. blocks) to inv. blocks participating      in the inversion. Zero means inv. block was a priori fixed. This file is      needed by 'tomo2gmt' when computing spread values.        full resolution matrix (file 'resol.out')      For each inv. block participating in the inversion the complete      row of the resolution matrix (ires=2) is printed to file 'resol.out'.     ray points (file 'ray_XXXX.xyz')      Fatomo also computes rays following the steepest gradient in the travel      time field. If desired (rout=1) ray points of these rays can be output       in files named ray_XXXX.xyz, where XXXX denotes the absolute number of       observation. Ray points start at the receiver and are given in x,y,z      coordinates relative to the origin of the numerical grid (upper left      corner of the model).          fat rays (file 'fat_XXXX.xyz')      Fat rays can be output in files named 'fat_XXXX.xyz'. As for the ray      points, for each observation a new file will be opened. In FATOMO fat      rays are represented by cells of the size of the numerical grid spacing      and the i,j,k coordinates of these blocks are printed in the fat ray       files. Thereby i,j,k coordinates do not represent absolute x,y,z       coordinates, but the number of the block in x,y,z direction.      Also the absolute number of the numerical and the inversion block are      printed out.       RMS volumes after earthquake relocation (file 'rmsXXXX.xyz')        RMS values of each grid point of the fine grid used in the grid search      can be output in the file 'rmsXXXX.xyz'. For each grid point the x,y,z       position in km and the corresponding RMS value is given. A new file is       opened for each event. Plotting RMS volumes can be useful to assess the      location accuracy.