format

Author: Anand

docs/documentation/case.md

@@ -637,6 +637,7 @@ To restart the simulation from $k$-th time step, see @ref running "Restarting Ca
 | `cf_wrt`                | Logical | Write color function field |
 | `chem_wrt_T`            | Logical | Write temperature field for chemistry output |
 | `fft_wrt`               | Logical | Enable FFT output |
+| `double_mach`           | Logical | Enable Double Mach Reflection case |
 | `sim_data`              | Logical | Write interface and energy data files (post_process) |
 | `integral_wrt`          | Logical | Write integral data |
 | `num_integrals`         | Integer | Number of integral regions |

src/common/include/2dHardcodedIC.fpp

@@ -374,7 +374,8 @@
 
             q_prim_vf(contxb)%sf(i, j, 0) = 0.5_wp*(1._wp + tanh(cf*sb))*rhoshock + 0.5_wp*(1._wp - tanh(cf*sb))*rho0_dm
             q_prim_vf(E_idx)%sf(i, j, 0) = 0.5_wp*(1._wp + tanh(cf*sb))*pshock + 0.5_wp*(1._wp - tanh(cf*sb))*p0_dm
-            q_prim_vf(momxb)%sf(i, j, 0) = 0.5_wp*(1._wp + tanh(cf*sb))*(velshock*tan(theta_dm)) + 0.5_wp*(1._wp - tanh(cf*sb))*u0_dm
+            q_prim_vf(momxb)%sf(i, j, &
+                      & 0) = 0.5_wp*(1._wp + tanh(cf*sb))*(velshock*tan(theta_dm)) + 0.5_wp*(1._wp - tanh(cf*sb))*u0_dm
             q_prim_vf(momxb + 1)%sf(i, j, 0) = 0.5_wp*(1._wp + tanh(cf*sb))*(-velshock) + 0.5_wp*(1._wp - tanh(cf*sb))*v0_dm
         end if
     case default

src/common/m_boundary_common.fpp

@@ -70,27 +70,27 @@ contains
 
         if (double_mach) then
             ! Smoothing parameter
-            cf = 20._wp 
+            cf = 20._wp
             ! Mach number
-            Mach = 10._wp 
+            Mach = 10._wp
 
-            !Pre shock
-            rho0_dm = 1.4_wp 
+            ! Pre shock
+            rho0_dm = 1.4_wp
             p0_dm = 1._wp
             u0_dm = 0._wp
             v0_dm = 0._wp
 
-            gam_dm = 1._wp + 1._wp / fluid_pp(1)%gamma
+            gam_dm = 1._wp + 1._wp/fluid_pp(1)%gamma
 
-            !Post shock
+            ! Post shock
             pshock = (2._wp*gam_dm*Mach**2 - (gam_dm - 1._wp))/(1._wp + gam_dm)
             rhoshock = ((1._wp + gam_dm)*Mach**2)*rho0_dm/((gam_dm - 1._wp)*Mach**2 + 2._wp)
-            velshock = (Mach - (Mach*rho0_dm/rhoshock))/2._wp    
+            velshock = (Mach - (Mach*rho0_dm/rhoshock))/2._wp
 
             ! Reflecting wall location
-            xr_dm = 1._wp/6._wp 
+            xr_dm = 1._wp/6._wp
             ! Shock Angle in radian
-            theta_dm = (pi/180_wp)*60_wp
+            theta_dm = (pi/180._wp)*60._wp
 
             $:GPU_UPDATE(device='[double_mach, Mach, pshock, rhoshock, velshock, rho0_dm, p0_dm, u0_dm, v0_dm, xr_dm, theta_dm, gam_dm]')
         end if
@@ -483,8 +483,9 @@ contains
                         if (igr) then
                             q_prim_vf(momxb)%sf(k, -j, l) = q_prim_vf(contxb)%sf(k, -j, l)*q_prim_vf(momxb)%sf(k, -j, l)
                             q_prim_vf(momxb + 1)%sf(k, -j, l) = q_prim_vf(contxb)%sf(k, -j, l)*q_prim_vf(momxb + 1)%sf(k, -j, l)
-                            q_prim_vf(E_idx)%sf(k, -j, l) = (gam_dm + 1._wp)*q_prim_vf(E_idx)%sf(k, -j, l) + 0.5_wp*(q_prim_vf(momxb)%sf(k, &
-                                      & -j, l)**2 + q_prim_vf(momxb + 1)%sf(k, -j, l)**2)/q_prim_vf(contxb)%sf(k, -j, l)
+                            q_prim_vf(E_idx)%sf(k, -j, l) = (gam_dm + 1._wp)*q_prim_vf(E_idx)%sf(k, -j, &
+                                      & l) + 0.5_wp*(q_prim_vf(momxb)%sf(k, -j, l)**2 + q_prim_vf(momxb + 1)%sf(k, -j, &
+                                      & l)**2)/q_prim_vf(contxb)%sf(k, -j, l)
                         end if
 #endif
                         if (.not. igr) then

src/post_process/m_global_parameters.fpp

@@ -504,9 +504,9 @@ contains
         rhoshock = dflt_real
         pshock = dflt_real
         velshock = dflt_real
-        u0_dm = dflt_real 
+        u0_dm = dflt_real
         v0_dm = dflt_real
-        p0_dm = dflt_real 
+        p0_dm = dflt_real
         rho0_dm = dflt_real
         theta_dm = dflt_real
         gam_dm = dflt_real

src/pre_process/m_global_parameters.fpp

@@ -543,16 +543,16 @@ contains
         bub_pp%cp_g = dflt_real; cp_g = dflt_real
         bub_pp%R_v = dflt_real; R_v = dflt_real
         bub_pp%R_g = dflt_real; R_g = dflt_real
-        
+
         double_mach = .false.
         xshock = dflt_real
         cf = dflt_real
         rhoshock = dflt_real
         pshock = dflt_real
         velshock = dflt_real
-        u0_dm = dflt_real 
+        u0_dm = dflt_real
         v0_dm = dflt_real
-        p0_dm = dflt_real 
+        p0_dm = dflt_real
         rho0_dm = dflt_real
         theta_dm = dflt_real
         gam_dm = dflt_real

src/simulation/m_global_parameters.fpp

@@ -498,7 +498,8 @@ module m_global_parameters
     !> @name Double Mach parameters
     logical  :: double_mach
     real(wp) :: xshock, cf, Mach, pshock, rhoshock, velshock, rho0_dm, p0_dm, u0_dm, v0_dm, xr_dm, theta_dm, gam_dm
-    $:GPU_DECLARE(create='[double_mach, xshock, cf, Mach, pshock, rhoshock, velshock, rho0_dm, p0_dm, u0_dm, v0_dm, xr_dm, theta_dm, gam_dm]')
+    $:GPU_DECLARE(create='[double_mach, xshock, cf, Mach, pshock, rhoshock, velshock, rho0_dm, p0_dm, u0_dm, v0_dm, xr_dm, &
+                  & theta_dm, gam_dm]')
 
 contains
 
@@ -805,9 +806,9 @@ contains
         rhoshock = dflt_real
         pshock = dflt_real
         velshock = dflt_real
-        u0_dm = dflt_real 
+        u0_dm = dflt_real
         v0_dm = dflt_real
-        p0_dm = dflt_real 
+        p0_dm = dflt_real
         rho0_dm = dflt_real
         theta_dm = dflt_real
         xr_dm = dflt_real