format

Author: wilfonba

examples/3D_IGR_33jet/case.py

@@ -3,18 +3,18 @@
 import json
 
 # Domain parameters
-D = 2.5       # Jet diameter
+D = 2.5  # Jet diameter
 Nd = 711  # Cells per jet diameter
 
-x0 = 0    # x_beg coordinate
-x1 = 19*D   # x_end coordinate
-y0 = -22*D/2   # y_beg coordinate
-y1 = 22*D/2    # y_end coordinate
-z0 = -22*D/2   # y_beg coordinate
-z1 = 22*D/2   # y_end coordinate
-Nx = int(Nd*(x1 - x0)/D) - 1
-Ny = int(Nd*(y1 - y0)/D) - 1
-Nz = int(Nd*(z1 - z0)/D) - 1
+x0 = 0  # x_beg coordinate
+x1 = 19 * D  # x_end coordinate
+y0 = -22 * D / 2  # y_beg coordinate
+y1 = 22 * D / 2  # y_end coordinate
+z0 = -22 * D / 2  # y_beg coordinate
+z1 = 22 * D / 2  # y_end coordinate
+Nx = int(Nd * (x1 - x0) / D) - 1
+Ny = int(Nd * (y1 - y0) / D) - 1
+Nz = int(Nd * (z1 - z0) / D) - 1
 
 time_end = 5
 igrIters = 5
@@ -42,8 +42,8 @@
             "p": int(Nz),
             "dt": dt,
             "t_step_start": 0,
-            "t_step_stop": 10, #Nt,
-            "t_step_save": 10, #int(Nt/50),
+            "t_step_stop": 10,  # Nt,
+            "t_step_save": 10,  # int(Nt/50),
             # Simulation Algorithm Parameters
             "num_patches": 1,
             "num_bc_patches": 0,
@@ -74,12 +74,12 @@
             "down_sample": "T",
             # Patch
             "patch_icpp(1)%geometry": 9,
-            "patch_icpp(1)%x_centroid": (x1 + x0)/2,
-            "patch_icpp(1)%y_centroid": (y1 + y0)/2,
-            "patch_icpp(1)%z_centroid": (z1 + z0)/2,
-            "patch_icpp(1)%length_x": 2*(x1 - x0),
-            "patch_icpp(1)%length_y": 2*(y1 - y0),
-            "patch_icpp(1)%length_z": 2*(z1 - z0),
+            "patch_icpp(1)%x_centroid": (x1 + x0) / 2,
+            "patch_icpp(1)%y_centroid": (y1 + y0) / 2,
+            "patch_icpp(1)%z_centroid": (z1 + z0) / 2,
+            "patch_icpp(1)%length_x": 2 * (x1 - x0),
+            "patch_icpp(1)%length_y": 2 * (y1 - y0),
+            "patch_icpp(1)%length_z": 2 * (z1 - z0),
             "patch_icpp(1)%hcid": 303,
             "patch_icpp(1)%vel(1)": 0.0e00,
             "patch_icpp(1)%vel(2)": 0.0e00,
@@ -112,6 +112,7 @@
             "fluid_pp(1)%pi_inf": 0.0,
             "viscous": "T",
             "fluid_pp(1)%Re(1)": 5e5,
-        }, indent=4
+        },
+        indent=4,
     )
 )

src/common/include/3dHardcodedIC.fpp

@@ -132,19 +132,19 @@
         rcut = f_cut_on(r - r_th, eps_smooth)
         xcut = f_cut_on(x_cc(i), eps_smooth)
 
-        q_prim_vf(momxb)%sf(i, j, k) = ux_th * rcut * xcut + ux_am
+        q_prim_vf(momxb)%sf(i, j, k) = ux_th*rcut*xcut + ux_am
         q_prim_vf(momxb + 1)%sf(i, j, k) = 0._wp
         q_prim_vf(momxe)%sf(i, j, k) = 0._wp
 
         if (num_fluids == 1) then
-            q_prim_vf(contxb)%sf(i, j, k) = (rho_th - rho_am) * rcut * xcut + rho_am
+            q_prim_vf(contxb)%sf(i, j, k) = (rho_th - rho_am)*rcut*xcut + rho_am
         else
-            q_prim_vf(advxb)%sf(i, j, k) = (1._wp - 2._wp*eps) * rcut * xcut + eps
+            q_prim_vf(advxb)%sf(i, j, k) = (1._wp - 2._wp*eps)*rcut*xcut + eps
             q_prim_vf(contxb)%sf(i, j, k) = rho_th*q_prim_vf(advxb)%sf(i, j, k)
             q_prim_vf(contxe)%sf(i, j, k) = rho_am*(1._wp - q_prim_vf(advxb)%sf(i, j, k))
         end if
 
-        q_prim_vf(E_idx)%sf(i, j, k) = p_th * rcut * xcut + p_am
+        q_prim_vf(E_idx)%sf(i, j, k) = p_th*rcut*xcut + p_am
 
     case (303) ! 3D Multijet
 
@@ -160,19 +160,19 @@
         rcut = rcut_arr(j, k)
         xcut = f_cut_on(x_cc(i), eps_smooth)
 
-        q_prim_vf(momxb)%sf(i, j, k) = ux_th * rcut * xcut + ux_am
+        q_prim_vf(momxb)%sf(i, j, k) = ux_th*rcut*xcut + ux_am
         q_prim_vf(momxb + 1)%sf(i, j, k) = 0._wp
         q_prim_vf(momxe)%sf(i, j, k) = 0._wp
 
         if (num_fluids == 1) then
-            q_prim_vf(contxb)%sf(i, j, k) = (rho_th - rho_am) * rcut * xcut + rho_am
+            q_prim_vf(contxb)%sf(i, j, k) = (rho_th - rho_am)*rcut*xcut + rho_am
         else
-            q_prim_vf(advxb)%sf(i, j, k) = (1._wp - 2._wp*eps) * rcut * xcut + eps
+            q_prim_vf(advxb)%sf(i, j, k) = (1._wp - 2._wp*eps)*rcut*xcut + eps
             q_prim_vf(contxb)%sf(i, j, k) = rho_th*q_prim_vf(advxb)%sf(i, j, k)
             q_prim_vf(contxe)%sf(i, j, k) = rho_am*(1._wp - q_prim_vf(advxb)%sf(i, j, k))
         end if
 
-        q_prim_vf(E_idx)%sf(i, j, k) = p_th * rcut * xcut + p_am
+        q_prim_vf(E_idx)%sf(i, j, k) = p_th*rcut*xcut + p_am
 
     case (370)
         ! This hardcoded case extrudes a 2D profile to initialize a 3D simulation domain

src/simulation/m_ibm.fpp

@@ -912,9 +912,9 @@ contains
                             do q = 1, nb
                                 do l = 1, nnode
                                     presb_IP((q - 1)*nnode + l) = presb_IP((q - 1)*nnode + l) + &
-                                                                    coeff*real(pb_in(i, j, k, l, q), kind=wp)
+                                                                  coeff*real(pb_in(i, j, k, l, q), kind=wp)
                                     massv_IP((q - 1)*nnode + l) = massv_IP((q - 1)*nnode + l) + &
-                                                                    coeff*real(mv_in(i, j, k, l, q), kind=wp)
+                                                                  coeff*real(mv_in(i, j, k, l, q), kind=wp)
                                 end do
                             end do
                         end if