search

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

Calculating solid tides before 1962 #23

Closed glPeng-seis closed 2 years ago

glPeng-seis commented 2 years ago

I found that we can not update the data in the file etpolut1.dat before 1962. we thus can not calculate the tides before 1962. How to solve this problem. Thanks for your help.

hydrogeoscience commented 2 years ago

This may be an overall limitation of ETERNA. I will look into this when I have some time. Before doing so, can you tell me where you get pole position data from before the year 1962?

glPeng-seis commented 2 years ago

I also do not know where to get pole position data before the year 1962. I have an idea that pygtide still calculate tides although without pole position data. It's just a little less accurate.

hydrogeoscience commented 2 years ago

Does it not work for times before 1962 right now?

GPeng00 commented 2 years ago

Yes, it do not work. The traceback as following,

Traceback (most recent call last): File "run_calc_tide.py", line 62, in ranking = calc_volume_strain_indays(date_onset, nterms, samp, vol_nm, la, lo, elev) File "/mnt/j/Vols/calc_tide.py", line 65, in calc_volume_strain_indays pt.predict(lat, lon, height, start, duration, splrate, tidalcompo=8) File "/mnt/j/Vols/pygtide/pygtide.py", line 368, in predict "Please update file '%s'." % (self.etpolut1_start, self.etpolut1_end, file)) ValueError: Dates exceed permissible range for pole/LOD tide correction (interval 1962-01-01 00:00:00 to 2023-06-10 00:00:00). Please update file 'etpolut1.dat'.

hydrogeoscience commented 2 years ago

Thanks. I have updated core.py and changed the error to a user warning. This should now work for years before 1962.