states_L_and_R_channel Subroutine

public subroutine states_L_and_R_channel(me, i, Cell_L, Cell_R)

Arguments

Type IntentOptional Attributes Name
type(channel_t), intent(in) :: me
integer, intent(in) :: i
type(cell_state_t), intent(out) :: Cell_L
type(cell_state_t), intent(out) :: Cell_R

Called by

proc~~states_l_and_r_channel~~CalledByGraph proc~states_l_and_r_channel states_L_and_R_channel proc~he_riemann_solver_channel He_Riemann_solver_channel proc~he_riemann_solver_channel->proc~states_l_and_r_channel proc~main_loop main_loop proc~main_loop->proc~he_riemann_solver_channel program~reims_p reims_p program~reims_p->proc~main_loop

Source Code

subroutine states_L_and_R_channel(me,i,Cell_L,Cell_R)  

    type(channel_t), intent(in) :: me
    type(cell_state_t), intent(out) :: Cell_L,Cell_R
    integer, intent(in) :: i

    ! Internal boundary cells
    Cell_L%Prim(:)=me%StVar%He_pr(:,i)
    Cell_R%Prim(:)=me%StVar%He_pr(:,i+1)

    Cell_L%Cons(:)=me%StVar%He_cs(:,i)
    Cell_R%Cons(:)=me%StVar%He_cs(:,i+1)

    Cell_L%Jcb(:,:)=me%big%Jacob(:,:,i)
    Cell_R%Jcb(:,:)=me%big%Jacob(:,:,i+1)

    Cell_L%SuFrct=me%srcPP%Su(i)
    Cell_R%SuFrct=me%srcPP%Su(i+1)
    Cell_L%DerSuFrct(:)=me%srcPP%DerSu(:,i)
    Cell_R%DerSuFrct(:)=me%srcPP%DerSu(:,i+1)

    Cell_L%SrFrct=me%srcPP%Sr(i)
    Cell_R%SrFrct=me%srcPP%Sr(i+1)
    Cell_L%DerSrFrct(:)=me%srcPP%DerSr(:,i)
    Cell_R%DerSrFrct(:)=me%srcPP%DerSr(:,i+1)

    Cell_L%dx=me%HeProp%dxLoc(i)
    Cell_R%dx=me%HeProp%dxLoc(i+1)
end subroutine states_L_and_R_channel