fix frontier bug 2

Author: ChrisZYJ

src/common/m_derived_types.fpp

@@ -233,7 +233,7 @@ module m_derived_types
         !! the partial densities, density, velocity, pressure, volume fractions,
         !! specific heat ratio function and the liquid stiffness function.
 
-        type(vec3_dt) :: B !<
+        real(wp) :: Bx, By, Bz !<
         !! Magnetic field components; B%x is not used for 1D
 
         real(wp), dimension(6) :: tau_e !<

src/pre_process/m_assign_variables.fpp

@@ -480,20 +480,20 @@ contains
         if (mhd) then
             if (n == 0) then ! 1D: By, Bz
                 q_prim_vf(B_idx%beg)%sf(j, k, l) = &
-                    eta*patch_icpp(patch_id)%B%y &
+                    eta*patch_icpp(patch_id)%By &
                     + (1._wp - eta)*orig_prim_vf(B_idx%beg)
                 q_prim_vf(B_idx%beg + 1)%sf(j, k, l) = &
-                    eta*patch_icpp(patch_id)%B%z &
+                    eta*patch_icpp(patch_id)%Bz &
                     + (1._wp - eta)*orig_prim_vf(B_idx%beg + 1)
             else ! 2D/3D: Bx, By, Bz
                 q_prim_vf(B_idx%beg)%sf(j, k, l) = &
-                    eta*patch_icpp(patch_id)%B%x &
+                    eta*patch_icpp(patch_id)%Bx &
                     + (1._wp - eta)*orig_prim_vf(B_idx%beg)
                 q_prim_vf(B_idx%beg + 1)%sf(j, k, l) = &
-                    eta*patch_icpp(patch_id)%B%y &
+                    eta*patch_icpp(patch_id)%By &
                     + (1._wp - eta)*orig_prim_vf(B_idx%beg + 1)
                 q_prim_vf(B_idx%beg + 2)%sf(j, k, l) = &
-                    eta*patch_icpp(patch_id)%B%z &
+                    eta*patch_icpp(patch_id)%Bz &
                     + (1._wp - eta)*orig_prim_vf(B_idx%beg + 2)
             end if
         end if

src/pre_process/m_check_patches.fpp

@@ -535,9 +535,9 @@ contains
         @:PROHIBIT(model_eqns == 2 .and. any(patch_icpp(patch_id)%alpha_rho(1:num_fluids) < 0._wp), &
             "Patch "//trim(iStr)//": alpha_rho(1:num_fluids) must be greater than or equal to zero when model_eqns = 2")
 
-        is_set_B(1) = .not. f_is_default(patch_icpp(patch_id)%B%x)
-        is_set_B(2) = .not. f_is_default(patch_icpp(patch_id)%B%y)
-        is_set_B(3) = .not. f_is_default(patch_icpp(patch_id)%B%z)
+        is_set_B(1) = .not. f_is_default(patch_icpp(patch_id)%Bx)
+        is_set_B(2) = .not. f_is_default(patch_icpp(patch_id)%By)
+        is_set_B(3) = .not. f_is_default(patch_icpp(patch_id)%Bz)
 
         @:PROHIBIT(.not. mhd .and. any(is_set_B), &
             "Bx, By, and Bz must not be set if MHD is not enabled")

src/pre_process/m_global_parameters.fpp

@@ -401,9 +401,9 @@ contains
             patch_icpp(i)%qv = 0._wp
             patch_icpp(i)%qvp = 0._wp
             patch_icpp(i)%tau_e = 0._wp
-            patch_icpp(i)%B%x = dflt_real
-            patch_icpp(i)%B%y = dflt_real
-            patch_icpp(i)%B%z = dflt_real
+            patch_icpp(i)%Bx = dflt_real
+            patch_icpp(i)%By = dflt_real
+            patch_icpp(i)%Bz = dflt_real
             patch_icpp(i)%a(2) = dflt_real
             patch_icpp(i)%a(3) = dflt_real
             patch_icpp(i)%a(4) = dflt_real

src/pre_process/m_mpi_proxy.fpp

@@ -86,7 +86,7 @@ contains
                 & 'length_x', 'length_y', 'length_z', 'radius', 'epsilon',     &
                 & 'beta', 'smooth_coeff', 'rho', 'p0', 'm0', 'r0', 'v0',       &
                 & 'pres', 'gamma', 'pi_inf', 'hcid', 'cv', 'qv', 'qvp',        &
-                & 'model_threshold', 'cf_val', 'B%x', 'B%y', 'B%z']
+                & 'model_threshold', 'cf_val', 'Bx', 'By', 'Bz']
                 call MPI_BCAST(patch_icpp(i)%${VAR}$, 1, mpi_p, 0, MPI_COMM_WORLD, ierr)
             #:endfor