Replace intrinsics in allocate statements

src/common/m_derived_types.fpp

@@ -443,4 +443,10 @@ module m_derived_types
 
     end type bubbles_lagrange_parameters
 
+    !> Max and min number of cells in a direction of each combination of x-,y-, and z-
+    type cell_num_bounds
+        integer :: mn_max, np_max, mp_max, mnp_max
+        integer :: mn_min, np_min, mp_min, mnp_min
+    end type cell_num_bounds
+
 end module m_derived_types

src/common/m_helper_basic.fpp

@@ -16,7 +16,8 @@ module m_helper_basic
               f_is_default, &
               f_all_default, &
               f_is_integer, &
-              s_configure_coordinate_bounds
+              s_configure_coordinate_bounds, &
+              s_update_cell_bounds
 
 contains
 
@@ -146,4 +147,24 @@ contains
 
     end subroutine s_configure_coordinate_bounds
 
+    !> Updates the min and max number of cells in each set of axes
+    !! @param bounds Min ans max values to update
+    !! @param m Number of cells in x-axis
+    !! @param n Number of cells in y-axis
+    !! @param p Number of cells in z-axis
+    pure elemental subroutine s_update_cell_bounds(bounds, m, n, p)
+        type(cell_num_bounds), intent(out) :: bounds
+        integer, intent(in) :: m, n, p
+
+        bounds%mn_max = max(m, n)
+        bounds%np_max = max(n, p)
+        bounds%mp_max = max(m, p)
+        bounds%mnp_max = max(m, n, p)
+        bounds%mn_min = min(m, n)
+        bounds%np_min = min(n, p)
+        bounds%mp_min = min(m, p)
+        bounds%mnp_min = min(m, n, p)
+
+    end subroutine s_update_cell_bounds
+
 end module m_helper_basic

src/common/m_mpi_common.fpp

@@ -67,10 +67,10 @@ contains
                                          & (m + 2*buff_size + 1)* &
                                          & (n + 2*buff_size + 1)* &
                                          & (p + 2*buff_size + 1)/ &
-                                         & (min(m, n, p) + 2*buff_size + 1))
+                                         & (cells_bounds%mnp_min + 2*buff_size + 1))
             else
                 halo_size = -1 + buff_size*(v_size)* &
-                                         & (max(m, n) + 2*buff_size + 1)
+                                         & (cells_bounds%mn_max + 2*buff_size + 1)
             end if
         else
             halo_size = -1 + buff_size*(v_size)
@@ -1446,6 +1446,8 @@ contains
             end if
         end do
 
+        call s_update_cell_bounds(cells_bounds, m, n, p)
+
         ! Boundary condition at the beginning
         if (proc_coords(1) > 0 .or. (bc_x%beg == BC_PERIODIC .and. num_procs_x > 1)) then
             proc_coords(1) = proc_coords(1) - 1

src/common/m_variables_conversion.fpp

@@ -637,7 +637,7 @@
 #ifdef MFC_SIMULATION
 
         if (viscous) then
-            @:ALLOCATE(Res(1:2, 1:maxval(Re_size)))
+            @:ALLOCATE(Res(1:2, 1:Re_size_max))
             do i = 1, 2
                 do j = 1, Re_size(i)
                     Res(i, j) = fluid_pp(Re_idx(i, j))%Re(i)

src/post_process/m_global_parameters.fpp

@@ -38,6 +38,9 @@ module m_global_parameters
     integer :: p
     !> @}
 
+    !> @name Max and min number of cells in a direction of each combination of x-,y-, and z-
+    type(cell_num_bounds) :: cells_bounds
+
     integer(8) :: nGlobal ! Total number of cells in global domain
 
     !> @name Cylindrical coordinates (either axisymmetric or full 3D)
@@ -336,6 +339,8 @@ contains
 
         ! Computational domain parameters
         m = dflt_int; n = 0; p = 0
+        call s_update_cell_bounds(cells_bounds, m, n, p)
+
         m_root = dflt_int
         cyl_coord = .false.
 

src/post_process/m_start_up.f90

@@ -111,6 +111,9 @@ impure subroutine s_read_input_file
             end if
 
             close (1)
+
+            call s_update_cell_bounds(cells_bounds, m, n, p)
+
             ! Store m,n,p into global m,n,p
             m_glb = m
             n_glb = n

src/pre_process/m_global_parameters.fpp

@@ -40,6 +40,9 @@ module m_global_parameters
     integer :: n
     integer :: p
 
+    !> @name Max and min number of cells in a direction of each combination of x-,y-, and z-
+    type(cell_num_bounds) :: cells_bounds
+
     integer(8) :: nGlobal !< Global number of cells in the domain
 
     integer :: m_glb, n_glb, p_glb !< Global number of cells in each direction
@@ -306,6 +309,8 @@ contains
         ! Computational domain parameters
         m = dflt_int; n = 0; p = 0
 
+        call s_update_cell_bounds(cells_bounds, m, n, p)
+
         cyl_coord = .false.
 
         x_domain%beg = dflt_real

src/pre_process/m_start_up.fpp

@@ -169,6 +169,9 @@ contains
                                  'likely due to a datatype mismatch. Exiting.')
             end if
             close (1)
+
+            call s_update_cell_bounds(cells_bounds, m, n, p)
+
             ! Store m,n,p into global m,n,p
             m_glb = m
             n_glb = n

src/simulation/m_global_parameters.fpp

@@ -43,6 +43,9 @@
     integer :: m, n, p
     !> @}
 
+    !> @name Max and min number of cells in a direction of each combination of x-,y-, and z-
+    type(cell_num_bounds) :: cells_bounds
+
     !> @name Global number of cells in each direction
     !> @{
     integer :: m_glb, n_glb, p_glb
@@ -273,10 +276,11 @@
     !! numbers, will be non-negligible.
     !> @{
     integer, dimension(2) :: Re_size
+    integer :: Re_size_max
     integer, allocatable, dimension(:, :) :: Re_idx
     !> @}
 
-    $:GPU_DECLARE(create='[Re_size,Re_idx]')
+    $:GPU_DECLARE(create='[Re_size,Re_size_max,Re_idx]')
 
     ! The WENO average (WA) flag regulates whether the calculation of any cell-
     ! average spatial derivatives is carried out in each cell by utilizing the
@@ -525,6 +529,7 @@
 
         ! Computational domain parameters
         m = dflt_int; n = 0; p = 0
+        call s_update_cell_bounds(cells_bounds, m, n, p)
 
         cyl_coord = .false.
 
@@ -809,6 +814,7 @@
         ! of fluids for which the physical and geometric curvatures of the
         ! interfaces will be computed
         Re_size = 0
+        Re_size_max = 0
 
         ! Gamma/Pi_inf Model
         if (model_eqns == 1) then
@@ -1036,13 +1042,15 @@
             if (Re_size(1) > 0._wp) shear_stress = .true.
             if (Re_size(2) > 0._wp) bulk_stress = .true.
 
-            $:GPU_UPDATE(device='[Re_size,viscous,shear_stress,bulk_stress]')
+            Re_size_max = maxval(Re_size)
+
+            $:GPU_UPDATE(device='[Re_size,Re_size_max,viscous,shear_stress,bulk_stress]')
 
             ! Bookkeeping the indexes of any viscous fluids and any pairs of
             ! fluids whose interface will support effects of surface tension
             if (viscous) then
 
-                @:ALLOCATE(Re_idx(1:2, 1:maxval(Re_size)))
+                @:ALLOCATE(Re_idx(1:2, 1:Re_size_max))
 
                 k = 0
                 do i = 1, num_fluids
@@ -1169,8 +1177,6 @@
         if (ib) allocate (MPI_IO_IB_DATA%var%sf(0:m, 0:n, 0:p))
         Np = 0
 
-        $:GPU_UPDATE(device='[Re_size]')
-
         if (elasticity) then
             fd_number = max(1, fd_order/2)
         end if

src/simulation/m_mpi_proxy.fpp

@@ -62,10 +62,10 @@ contains
                                             & (m + 2*gp_layers + 1)* &
                                             & (n + 2*gp_layers + 1)* &
                                             & (p + 2*gp_layers + 1)/ &
-                                            & (min(m, n, p) + 2*gp_layers + 1)
+                                            & (cells_bounds%mnp_min + 2*gp_layers + 1)
                 else
                     i_halo_size = -1 + gp_layers* &
-                                            & (max(m, n) + 2*gp_layers + 1)
+                                            & (cells_bounds%mn_max + 2*gp_layers + 1)
                 end if
             else
                 i_halo_size = -1 + gp_layers

src/simulation/m_pressure_relaxation.fpp

@@ -42,7 +42,7 @@
         $:GPU_UPDATE(device='[gamma_min, pres_inf]')
 
         if (viscous) then
-            @:ALLOCATE(Res(1:2, 1:maxval(Re_size)))
+            @:ALLOCATE(Res(1:2, 1:Re_size_max))
             do i = 1, 2
                 do j = 1, Re_size(i)
                     Res(i, j) = fluid_pp(Re_idx(i, j))%Re(i)

src/simulation/m_riemann_solvers.fpp

@@ -3169,7 +3169,7 @@
         $:GPU_UPDATE(device='[Gs]')
 
         if (viscous) then
-            @:ALLOCATE(Res(1:2, 1:maxval(Re_size)))
+            @:ALLOCATE(Res(1:2, 1:Re_size_max))
         end if
 
         if (viscous) then

src/simulation/m_start_up.fpp

@@ -210,6 +210,8 @@ contains
             n_glb = n
             p_glb = p
 
+            call s_update_cell_bounds(cells_bounds, m, n, p)
+
             if (cfl_adap_dt .or. cfl_const_dt) cfl_dt = .true.
 
             if (any((/bc_x%beg, bc_x%end, bc_y%beg, bc_y%end, bc_z%beg, bc_z%end/) == -17) .or. &

src/simulation/m_viscous.fpp

@@ -35,7 +35,7 @@
 
         integer :: i, j !< generic loop iterators
 
-        @:ALLOCATE(Res_viscous(1:2, 1:maxval(Re_size)))
+        @:ALLOCATE(Res_viscous(1:2, 1:Re_size_max))
 
         do i = 1, 2
             do j = 1, Re_size(i)