I had to remove a bad memory access that I have in for the 1D case

Author: danieljvickers

src/simulation/m_time_steppers.fpp

@@ -444,7 +444,39 @@ contains
             call s_update_lagrange_tdv_rk(stage=1)
         end if
 
-        call s_evolve_q_pb_mv(1, 1.0_wp, 0.0_wp, 1.0_wp)
+        !$acc parallel loop collapse(4) gang vector default(present)
+        do i = 1, sys_size
+            do l = 0, p
+                do k = 0, n
+                    do j = 0, m
+                        q_cons_ts(1)%vf(i)%sf(j, k, l) = &
+                            q_cons_ts(1)%vf(i)%sf(j, k, l) &
+                            + dt*rhs_vf(i)%sf(j, k, l)
+                    end do
+                end do
+            end do
+        end do
+
+        !Evolve pb and mv for non-polytropic qbmm
+        if (qbmm .and. (.not. polytropic)) then
+            !$acc parallel loop collapse(5) gang vector default(present)
+            do i = 1, nb
+                do l = 0, p
+                    do k = 0, n
+                        do j = 0, m
+                            do q = 1, nnode
+                                pb_ts(1)%sf(j, k, l, q, i) = &
+                                    pb_ts(1)%sf(j, k, l, q, i) &
+                                    + dt*rhs_pb(j, k, l, q, i)
+                                mv_ts(1)%sf(j, k, l, q, i) = &
+                                    mv_ts(1)%sf(j, k, l, q, i) &
+                                    + dt*rhs_mv(j, k, l, q, i)
+                            end do
+                        end do
+                    end do
+                end do
+            end do
+        end if
 
         if (bodyForces) call s_apply_bodyforces(q_cons_ts(1)%vf, q_prim_vf, rhs_vf, dt)