Modify the projection formula, add tolerance, and add tests.

[Chuan1937]

When dip=90° and rank=180°, tan(90°)×sin(180°)=Inf×0, causing an error.

Change the formula in the utilmeca_proj_radius function from atan(tand×sind) to atan2(sin×sin, cos) and add an epsilon tolerance.

Fixes #8053

Test results:

[Esteban82]

In Code Validator it says that there is an error with the script of this PR.

Following bash scripts may not have execute permission:
1: ./test/seis/psmeca_radius.sh
Error: Process completed with exit code 1.

[Chuan1937]

In Code Validator it says that there is an error with the script of this PR.

ok，I'll change the execute permission.