UpdateGRK4_DGModel1D_t Subroutine

  subroutine UpdateGRK4_DGModel1D_t(this,m)
    implicit none
    class(DGModel1D_t),intent(inout) :: this
    integer,intent(in) :: m
    ! Local
    integer :: i,iEl,iVar

    do concurrent(i=1:this%solution%N+1, &
                  iel=1:this%mesh%nElem,ivar=1:this%solution%nVar)

      this%workSol%interior(i,iEl,iVar) = rk4_a(m)* &
                                          this%workSol%interior(i,iEl,iVar)+ &
                                          this%dSdt%interior(i,iEl,iVar)

      this%solution%interior(i,iEl,iVar) = &
        this%solution%interior(i,iEl,iVar)+ &
        rk4_g(m)*this%dt*this%workSol%interior(i,iEl,iVar)

    enddo

  endsubroutine UpdateGRK4_DGModel1D_t