Skip to content

update leap second table to 28 Jun 2025 based on IERS Bulletin C 68 #96

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
yunjunz merged 1 commit into from
Oct 4, 2024
Merged

update leap second table to 28 Jun 2025 based on IERS Bulletin C 68 #96

Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
40 changes: 20 additions & 20 deletions src/pysolid/solid.for
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -3,16 +3,16 @@
*** Author: Dennis Milbert, 2018-06-01. The code is available at:
*** http://geodesyworld.github.io/SOFTS/solid.htm and can be downloaded as:
*** wget http://geodesyworld.github.io/SOFTS/solid.for.txt -O solid.for
*** The code is based on dehanttideinel.f provided by V. Dehant, S. Mathews,
*** J. Gipson and C. Bruyninx. The latest version of dehanttideinel.f and its
*** The code is based on dehanttideinel.f provided by V. Dehant, S. Mathews,
*** J. Gipson and C. Bruyninx. The latest version of dehanttideinel.f and its
*** dependent subroutines can be download from IERS conventions website as:
*** wget -r -l1 --no-parent -R "index.html*" -nH --cut-dirs=3 https://iers-conventions.obspm.fr/content/chapter7/software/dehanttideinel
*** Sep 2020: modify solid() to solid_point/grid() as subroutines, Z. Yunjun and S. Sangha.
*** Apr 2023: return numpy arrays instead of writing txt file, S. Staniewicz.

subroutine solid_grid(iyr,imo,idy,ihh,imm,iss,
* glad0,steplat,nlat,glod0,steplon,nlon,tide_e,tide_n,tide_u)

*** calculate solid earth tides (SET) for one spatial grid given the date/time
*** Arguments: iyr/imo/idy/ihh/imm/iss - int, date/time for YYYY/MM/DD/HH/MM/SS
*** glad0/glad1/steplat - float, north(Y_FIRST)/south/step(negative) in deg
Expand Down Expand Up @@ -488,9 +488,9 @@
* The code below is commented to prevent restoring deformation
* due to permanent tide. All the code above removes
* total tidal deformation with conventional Love numbers.
* The code above realizes a conventional tide free crust (i.e. ITRF).
* The code above realizes a conventional tide free crust (i.e. ITRF).
* This does NOT conform to Resolution 16 of the 18th General Assembly
* of the IAG (1983). This resolution has not been implemented by
* of the IAG (1983). This resolution has not been implemented by
* the space geodesy community in general (c.f. IERS Conventions 2003).
*-----------------------------------------------------------------------

Expand Down Expand Up @@ -1235,16 +1235,16 @@

implicit double precision(a-h,o-z)
double precision u,v,w,ra,az,va

s2=u*u+v*v
r2=s2+w*w

s =dsqrt(s2)
ra=dsqrt(r2)

az=datan2(v,u)
va=datan2(w,s)

return
end
*-----------------------------------------------------------------------
Expand Down Expand Up @@ -1290,38 +1290,38 @@
end
*-----------------------------------------------------------------------
subroutine rot1(theta,x,y,z,u,v,w)

*** rotate coordinate axes about 1 axis by angle of theta radians
*** x,y,z transformed into u,v,w

implicit double precision(a-h,o-z)
double precision theta,x,y,z,u,v,w

s=dsin(theta)
c=dcos(theta)

u=x
v=c*y+s*z
w=c*z-s*y

return
end
*-----------------------------------------------------------------------
subroutine rot3(theta,x,y,z,u,v,w)

*** rotate coordinate axes about 3 axis by angle of theta radians
*** x,y,z transformed into u,v,w

implicit double precision(a-h,o-z)
double precision theta,x,y,z,u,v,w

s=dsin(theta)
c=dcos(theta)

u=c*x+s*y
v=c*y-s*x
w=z

return
end
************************************************************************
Expand Down Expand Up @@ -1559,9 +1559,9 @@

implicit double precision(a-h,o-z)
double precision tsec
!*** upper limit, leap second table, 2024jun28
!*** upper limit, leap second table, 2025jun28
!*** lower limit, leap second table, 1972jan01
parameter(MJDUPPER=60489)
parameter(MJDUPPER=60854)
parameter(MJDLOWER=41317)

!*** leap second table limit flag
Expand Down Expand Up @@ -1642,7 +1642,7 @@
***** other leap second references at:
***** http://hpiers.obspm.fr/eoppc/bul/bulc/Leap_Second_History.dat
***** http://hpiers.obspm.fr/eoppc/bul/bulc/bulletinc.dat
***** File expires on 28 June 2024
***** File expires on 28 June 2025

*** test against newest leaps first

Expand Down