Skip to content
This repository was archived by the owner on Aug 28, 2024. It is now read-only.

Commit

Permalink
fixed type conversion error that a particular version of gfortran did…
Browse files Browse the repository at this point in the history
…n't like
  • Loading branch information
daminton committed Feb 15, 2024
1 parent 9d2e7be commit 6791265
Showing 1 changed file with 7 additions and 2 deletions.
9 changes: 7 additions & 2 deletions src/swiftest/swiftest_drift.f90
Original file line number Diff line number Diff line change
Expand Up @@ -587,13 +587,18 @@ module subroutine swiftest_drift_cb_rotphase_update(self, param, dt)
!!
!! initial 0 is set at the x-axis
!! phase is stored and calculated in radians. Converted to degrees for output

implicit none
! Arguments
class(swiftest_cb), intent(inout) :: self !! Swiftest central body data structure
class(swiftest_parameters), intent(in) :: param !! Current run configuration parameters
real(DP), intent(in) :: dt !! Stepsize

self%rotphase = MOD(self%rotphase + (.mag. self%rot(:)) * dt * param%TU2S, 2 * PI) ! phase angle calculated in radians and then scaled by 2pi to be unitless
! Internals
real(DP) :: rotmag

rotmag = (.mag. self%rot(:)) * dt * param%TU2S

self%rotphase = MOD(self%rotphase + rotmag, 2 * PI) ! phase angle calculated in radians and then scaled by 2pi to be unitless

end subroutine swiftest_drift_cb_rotphase_update

Expand Down

0 comments on commit 6791265

Please sign in to comment.