readme

这是一层析静校正的源码

     '#' 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.