Fix for benchmarking on frontier

Author: Anand Radhakrishnan

src/simulation/m_bubbles_EE.fpp

@@ -37,24 +37,23 @@ contains
 
         @:ALLOCATE(rs(1:nb))
         @:ALLOCATE(vs(1:nb))
-        if (.not. polytropic) then
-            @:ALLOCATE(ps(1:nb))
-            @:ALLOCATE(ms(1:nb))
-        end if
+        @:ALLOCATE(ps(1:nb))
+        @:ALLOCATE(ms(1:nb))
 
         do l = 1, nb
             rs(l) = bub_idx%rs(l)
             vs(l) = bub_idx%vs(l)
             if (.not. polytropic) then
                 ps(l) = bub_idx%ps(l)
                 ms(l) = bub_idx%ms(l)
+            else
+                ps(l) = rs(l)
+                ms(l) = rs(l)
             end if
         end do
 
         $:GPU_UPDATE(device='[rs, vs]')
-        if (.not. polytropic) then
-            $:GPU_UPDATE(device='[ps, ms]')
-        end if
+        $:GPU_UPDATE(device='[ps, ms]')
 
         @:ALLOCATE(divu%sf(idwbuff(1)%beg:idwbuff(1)%end, idwbuff(2)%beg:idwbuff(2)%end, idwbuff(3)%beg:idwbuff(3)%end))
         @:ACC_SETUP_SFs(divu)

src/simulation/m_riemann_solvers.fpp

@@ -2790,13 +2790,8 @@ contains
                                         $:GPU_LOOP(parallelism='[seq]')
                                         do i = 1, nb
                                             if (.not. qbmm) then
-                                                if (polytropic) then
-                                                    pbw_L(i) = f_cpbw_KM(R0(i), R0_L(i), V0_L(i), 0._wp)
-                                                    pbw_R(i) = f_cpbw_KM(R0(i), R0_R(i), V0_R(i), 0._wp)
-                                                else
-                                                    pbw_L(i) = f_cpbw_KM(R0(i), R0_L(i), V0_L(i), P0_L(i))
-                                                    pbw_R(i) = f_cpbw_KM(R0(i), R0_R(i), V0_R(i), P0_R(i))
-                                                end if
+                                                pbw_L(i) = f_cpbw_KM(R0(i), R0_L(i), V0_L(i), P0_L(i))
+                                                pbw_R(i) = f_cpbw_KM(R0(i), R0_R(i), V0_R(i), P0_R(i))
                                             end if
                                         end do