ResPart_set_AdvanceMarkers Subroutine

private pure subroutine ResPart_set_AdvanceMarkers(this, dt)

Advances markers to next timestep

Type Bound

ResPart_set

Arguments

Type IntentOptional Attributes Name
class(ResPart_set), intent(inout) :: this

Collection of Resolved Particles
real(kind=wp), intent(in) :: dt

Timestep

Calls

proc~~respart_set_advancemarkers~~CallsGraph proc~respart_set_advancemarkers ResPart_set%ResPart_set_AdvanceMarkers markers markers proc~respart_set_advancemarkers->markers particle particle proc~respart_set_advancemarkers->particle proc~cross_product~3 cross_product proc~respart_set_advancemarkers->proc~cross_product~3

Called by

proc~~respart_set_advancemarkers~~CalledByGraph proc~respart_set_advancemarkers ResPart_set%ResPart_set_AdvanceMarkers proc~cdifs_obj_advancesolutionrp cdifs_obj_AdvanceSolutionRP proc~cdifs_obj_advancesolutionrp->proc~respart_set_advancemarkers proc~grans_obj_advancesolution grans_obj_AdvanceSolution proc~grans_obj_advancesolution->proc~respart_set_advancemarkers interface~grans_obj_advancesolution grans_obj%grans_obj_AdvanceSolution interface~grans_obj_advancesolution->proc~grans_obj_advancesolution proc~cdifs_obj_advancesolution cdifs_obj_AdvanceSolution proc~cdifs_obj_advancesolution->proc~cdifs_obj_advancesolutionrp interface~cdifs_obj_advancesolution cdifs_obj%cdifs_obj_AdvanceSolution interface~cdifs_obj_advancesolution->proc~cdifs_obj_advancesolution

Source Code

    pure subroutine ResPart_set_AdvanceMarkers(this,dt)
      !> Advances markers to next timestep
      class(ResPart_set),  intent(inout) :: this                               !! Collection of Resolved Particles
      real(wp),            intent(in)    :: dt                                 !! Timestep
      ! Work variables
      integer        :: n,m
      integer(leapI8):: id
      real(wp)       :: rot_rate
      real(wp)       :: angle
      real(wp)       :: axis(3)
      real(wp)       :: disp_old(3)
      real(wp)       :: a1,a2,a3

      select type (particle =>this%p)
      type is (ResPart_obj)
        select type(markers=>this%ib%p)
        type is (marker_obj)
          do m=1,this%ib%count_
            ! Get ID of parent centroid
            id=markers(m)%s

            ! Get local index of the centroid
            n = this%lookup(id)

            ! prepare rotation coefficients
            rot_rate = norm2(particle(n)%w)
            axis     = particle(n)%w/(rot_rate+epsilon(1.0_wp))
            angle    = dt*rot_rate
            a1       = cos(angle)
            a2       = sin(angle)
            a3       = (1.0_wp-cos(angle))

            ! Old displacement
            disp_old = markers(m)%disp

            ! New rotated displacement
            markers(m)%disp = a1*disp_old                        &
                            + a2*cross_product(axis,disp_old)    &
                            + a3*dot_product(disp_old,axis)*axis

            ! Update position and velocity
            markers(m)%p = particle(n)%p + markers(m)%disp
            markers(m)%v = particle(n)%v + cross_product(particle(n)%w,markers(m)%disp)
            ! Update normals
            markers(m)%n = markers(m)%disp/norm2(markers(m)%disp)
          end do
        end select
      end select

      return
    end subroutine ResPart_set_AdvanceMarkers