links_circulator_update Subroutine

public subroutine links_circulator_update(me, dt)

Arguments

Type IntentOptional Attributes Name
type(circulator_t), intent(inout) :: me
real(kind=dp), intent(in) :: dt

Called by

proc~~links_circulator_update~~CalledByGraph proc~links_circulator_update links_circulator_update proc~main_loop main_loop proc~main_loop->proc~links_circulator_update program~reims_p reims_p program~reims_p->proc~main_loop

Source Code

subroutine links_circulator_update(me,dt)
    real(dp), intent(in) :: dt
    type(circulator_t), intent(inout) :: me
  
    real(dp) :: VarNCcenter
    
    ! Link IN-out
    VarNCcenter=0.0_dp
    if(R_Correction) VarNCcenter=me%br(1)%p%Prim(Pri_ro)*(me%br(1)%p%Prim(Pri_c)**2)
    me%MxInOut(:,:)=me%br(1)%p%Sgn4j*(dt/me%br(1)%p%dxLoc)*me%in_DerFlxdBr_out(:,:)
    me%MxInOut(:,Con_R)=me%MxInOut(:,Con_R)+me%br(1)%p%Sgn4j*(dt/me%br(1)%p%dxLoc)*VarNCcenter*me%in_DerVitdBr_out(:)
  
    ! Link OUT-in
    VarNCcenter=0.0_dp
    if(R_Correction) VarNCcenter=me%br(2)%p%Prim(Pri_ro)*(me%br(2)%p%Prim(Pri_c)**2)
    me%MxOutIn(:,:)=me%br(2)%p%Sgn4j*(dt/me%br(2)%p%dxLoc)*me%out_DerFlxdBr_in(:,:)
    me%MxOutIn(:,Con_R)=me%MxOutIn(:,Con_R)+me%br(2)%p%Sgn4j*(dt/me%br(2)%p%dxLoc)*VarNCcenter*me%out_DerVitdBr_in(:)    
end subroutine links_circulator_update