hydrogeoscience / pygtide

A Python module and wrapper for ETERNA PREDICT to compute gravitational tides on Earth
Mozilla Public License 2.0
50 stars 19 forks source link

Wave group parameters wrong #12

Closed Jollyfant closed 4 years ago

Jollyfant commented 4 years ago

Hi, I had a problem with different results between ETERNA and this. When running the software it seems that the wave groups all have amplitude 1 and 0 phase. This is maybe related to the broken etpolut1.bin and etpolut1.dat.

  Wave groups and observed tidal parameters

    no.  from    to ampl.fac. phase lead      ampl.  WD body
                                 [deg]      [  nm/s**2 ]

       1     1     3    1.0000    0.0000     DC     35.5811    1.2165
       2     4    20    1.0000    0.0000     long    0.0001    1.0728
       3    21   213    1.0000    0.0000     SA      3.1589    1.2165
       4   247   499    1.0000    0.0000     SSA     3.4928    1.2165
       5   623   828    1.0000    0.0000     MSM     0.7610    1.2165
       6   868  1313    1.0000    0.0000     MM      4.5043    1.0728
       7  1398  1645    1.0000    0.0000     MSF     0.6600    1.2165
       8  1664  2098    1.0000    0.0000     MF      7.5327    1.2165
       9  2184  2435    1.0000    0.0000     MSTM    0.2739    1.2165
      10  2460  2830    1.0000    0.0000     MTM     1.4422    1.2165
      11  2893  3228    1.0000    0.0000     MSQM    0.2303    1.2165
      12  3249  4149    1.0000    0.0000     MQM     0.1908    1.2165
      13  4489  5557    1.0000    0.0000     SGQ1    2.2200    1.1550
      14  5672  5980    1.0000    0.0000     2Q1     7.6138    1.1550
      15  6006  6365    1.0000    0.0000     SGM1    9.1884    1.1549
      16  6490  6954    1.0000    0.0000     Q1     57.5379    1.1548
      17  7013  7354    1.0000    0.0000     RO1    10.9286    1.1548
      18  7469  7939    1.0000    0.0000     O1    300.5193    1.1546
      19  7994  8380    1.0000    0.0000     TAU1    3.9251    1.1546
      20  8513  8949    1.0000    0.0000     NO1    23.6343    1.1540
      21  9021  9305    1.0000    0.0000     CHI1    4.5200    1.1539
      22  9507  9916    1.0000    0.0000     P1    139.9121    1.1491
      23  9998 10952    1.0000    0.0000     S1K1  422.6650    1.1336
      24 11135 11365    1.0000    0.0000     TET1    4.5199    1.1569
      25 11422 11939    1.0000    0.0000     J1     23.6339    1.1567
      26 12124 12401    1.0000    0.0000     SO1     3.9251    1.1560
      27 12437 12976    1.0000    0.0000     OO1    12.9311    1.1559
      28 13151 15621    1.0000    0.0000     NU1     2.4758    1.1557
      29 16304 17474    1.0000    0.0000     EPS2    3.4744    1.1586
      30 17572 18263    1.0000    0.0000     2N2    14.3800    1.1586
      31 18385 19214    1.0000    0.0000     N2     90.0485    1.1586
      32 19347 20143    1.0000    0.0000     M2    470.3240    1.1586
      33 20246 20346    1.0000    0.0000     LAM2    3.4684    1.1586
      34 20395 20882    1.0000    0.0000     L2     13.2952    1.1586
      35 21039 21459    1.0000    0.0000     S2    218.8214    1.1586
      36 21543 22278    1.0000    0.0000     K2     59.4836    1.1586
      37 22528 23121    1.0000    0.0000     ETA2    3.3254    1.1586
      38 23265 25611    1.0000    0.0000     2K2     0.8708    1.1586
      39 26500 27914    1.0000    0.0000     M3      7.3140    1.0709
      40 28390 28754    1.0000    0.0000     M4      0.1056    1.0355
      41 28755 28806    1.0000    0.0000     M5M6    0.0015    1.0000