subroutine fluxmethod_DGModel2D(this)
implicit none
class(DGModel2D),intent(inout) :: this
! Local
integer :: iel
integer :: i
integer :: j
real(prec) :: s(1:this%nvar),dsdx(1:this%nvar,1:2)
do iel = 1,this%mesh%nelem
do j = 1,this%solution%interp%N+1
do i = 1,this%solution%interp%N+1
s = this%solution%interior(i,j,iel,1:this%nvar)
dsdx = this%solutionGradient%interior(i,j,iel,1:this%nvar,1:2)
this%flux%interior(i,j,iel,1:this%nvar,1:2) = this%flux2d(s,dsdx)
enddo
enddo
enddo
call gpuCheck(hipMemcpy(this%flux%interior_gpu, &
c_loc(this%flux%interior), &
sizeof(this%flux%interior), &
hipMemcpyHostToDevice))
endsubroutine fluxmethod_DGModel2D