diff --git a/WENO.F90 b/WENO.F90 index 0b4bbb03b52bf977e316be6aced36717a9174d58..894dc48b269becbfaec1ffdccb57b52f673bc051 100644 --- a/WENO.F90 +++ b/WENO.F90 @@ -51,20 +51,6 @@ subroutine compute_flux_derivatives(nvars, ndim, nx, ny, nz, & delta_eta = 1.0 delta_zeta = 1.0 - !$loopy begin tagged: init_flux - do k=1,nz - do j=1,ny - do i=1,nx - do l=1,ndim - do m=1,nvars - flux_derivatives_generalized(m, l, i, j, k) = 0.0 - end do - end do - end do - end do - end do - !$loopy end tagged: init_flux - !$loopy begin tagged: to_generalized call convert_to_generalized(nvars, ndim, nx, ny, nz, & fluxes, metrics, metric_jacobians, generalized_fluxes) @@ -115,8 +101,7 @@ subroutine compute_flux_derivatives(nvars, ndim, nx, ny, nz, & do i=1,nx do v=1,nvars flux_derivatives_generalized(v, 1, i, j, k) = & - flux_derivatives_generalized(v, 1, i, j, k) & - + (weno_flux_tmp(v, i, j, k) - weno_flux_tmp(v, i-1, j, k)) / delta_xi + (weno_flux_tmp(v, i, j, k) - weno_flux_tmp(v, i-1, j, k)) / delta_xi end do end do end do