Formatting and added reference to Mercury-T

src/modules/swiftest_classes.f90

@@ -121,7 +121,7 @@ module swiftest_classes
       real(DP), dimension(NDIM) :: rot     = 0.0_DP !! Body rotation vector in inertial coordinate frame (units rad / TU)
       real(DP)                  :: k2      = 0.0_DP !! Tidal Love number
       real(DP)                  :: Q       = 0.0_DP !! Tidal quality factor
-      real(DP)                  :: tlag     = 0.0_DP !! Tidal phase lag angle
+      real(DP)                  :: tlag    = 0.0_DP !! Tidal phase lag angle
       real(DP), dimension(NDIM) :: L0      = 0.0_DP !! Initial angular momentum of the central body
       real(DP), dimension(NDIM) :: dL      = 0.0_DP !! Change in angular momentum of the central body
    contains

src/tides/tides_getacch_pl.f90

@@ -5,8 +5,15 @@ module subroutine tides_getacch_pl(self, system)
       !! author: Jennifer L.L. Pouplin, Carlisle A. wishard, and David A. Minton
       !!
       !! Calculated tidal torques from central body to any planet and from any planet to central body
-      !! planet - planet interactions are considered negligeable
-      !! Adapted from Mercury-T code from Bolmont et al. 2015
+      !! planet - planet interactions are considered negligable.
+      !! This is a constant time lag model.
+      !! 
+      !! Adapted from Mercury-T code from Bolmont et al. (2015)
+      !!
+      !! Reference:
+      !!    Bolmont, E., Raymond, S.N., Leconte, J., Hersant, F., Correia, A.C.M., 2015. 
+      !!       Mercury-T : A new code to study tidally evolving multi-planet systems. 
+      !!       Applications to Kepler-62. A&A 583, A116. https://doi.org/10.1051/0004-6361/201525909
       implicit none
       ! Arguments
       class(swiftest_pl),           intent(inout) :: self   !! Swiftest massive body object
@@ -30,7 +37,7 @@ module subroutine tides_getacch_pl(self, system)
             r_unit(:) = pl%xh(:,i) / rmag
             v_unit(:) = pl%vh(:,i) / vmag
             h_unit(:) = r_unit(:) .cross. v_unit(:)
-              
+
 
             !Ftr = 
             !Pto =