Changes in clustering code required for compatibility with fortran scripts

[basaks]

Clustering code should be flexible to be able to deal with different cell sizes in different zones, i.e. we want to be able to specify different grid sizes for the regional and global zones.

Best option is probably to work of the same parameter file used during inversion.

Something like this for the global grid size: 72 36 16 5. 5.

And this for the regional grid: 100. 190. -54. 0. 45 27 22

Source: param2x2 parameter file passed on by @alexgorb .

[basaks]

This is required as cluster.f passed on by @alexgorb implemented uniform grids everywhere. We need very different grids based on the parameter file for inversrion. Here is an example of the relevant part of the parameter file looks like for regional and global grid, with the first block corresponding to global:

72 36 16 5. 5.
      0.
    110.
    280.
    410.
    660.
    840.
   1020.
   1250.
   1400.
   1600.
   1850.
   2050.
   2250.
   2450.
   2600.
   2750.
   2889.
100. 190. -54. 0. 45 27  22
      0.
     35.
     70.
    110.
    160.
    210.
    260.
    310.
    360.
    410.
    460.
    510.
    560.
    610.
    660.
    710.
    810.
    910.
   1010.
   1110.
   1250.
   1400.
   1600.

[basaks]

The above implies non-trivial changes in the Grid class and the _find_block functions.

[basaks]

The above mentioned changes may not be required. However, a different set of changes can achieve the output @alexgorb 's scripts expect. These are:

• Changing the zone latitude outputs to co-latitude.
• Converting station and source block number, and event number to 4 byte integers instead of 8 byte integers so @alexgorb 's fortran code can read it unchanged. (Warning: Using bigger block sizes will likely reduce the number of arrivals in the final output due to cluster grouping and filtering.) If that is a concern, consider changing fortran code to read 8 byte integers instead.
• changing zone outputs to space separated format and limit floating point precision (helps output ascii file size too).

[basaks]

All of the above is done, @alexgorb has run inversion code with these changes. Closed.