Continuing to mark subroutines as pure or impure

Author: prathi-wind

src/common/m_variables_conversion.fpp

@@ -115,7 +115,7 @@ contains
         !! @param pres Pressure to calculate
         !! @param stress Shear Stress
         !! @param mom Momentum
-    pure subroutine s_compute_pressure(energy, alf, dyn_p, pi_inf, gamma, rho, qv, rhoYks, pres, T, stress, mom, G, pres_mag)
+    subroutine s_compute_pressure(energy, alf, dyn_p, pi_inf, gamma, rho, qv, rhoYks, pres, T, stress, mom, G, pres_mag)
 
 #ifdef _CRAYFTN
         !DIR$ INLINEALWAYS s_compute_pressure

src/post_process/m_checker.fpp

@@ -23,7 +23,7 @@ contains
 
     !> Checks compatibility of parameters in the input file.
         !! Used by the post_process stage
-    subroutine s_check_inputs
+    impure subroutine s_check_inputs
 
         call s_check_inputs_output_format
         call s_check_inputs_partial_domain

src/post_process/m_derived_variables.fpp

@@ -142,7 +142,7 @@ contains
         !!      values of the liquid stiffness, which are stored in the
         !!      derived flow quantity storage variable, q_sf.
         !!  @param q_sf Liquid stiffness
-    subroutine s_derive_liquid_stiffness(q_sf)
+    pure subroutine s_derive_liquid_stiffness(q_sf)
 
         real(wp), &
             dimension(-offset_x%beg:m + offset_x%end, &
@@ -171,7 +171,7 @@ contains
         !!      derived flow quantity storage variable, q_sf.
         !! @param q_prim_vf Primitive variables
         !! @param q_sf Speed of sound
-    subroutine s_derive_sound_speed(q_prim_vf, q_sf)
+    pure subroutine s_derive_sound_speed(q_prim_vf, q_sf)
 
         type(scalar_field), &
             dimension(sys_size), &
@@ -228,7 +228,7 @@ contains
         !!  @param i Component indicator
         !!  @param q_prim_vf Primitive variables
         !!  @param q_sf Flux limiter
-    subroutine s_derive_flux_limiter(i, q_prim_vf, q_sf)
+    pure subroutine s_derive_flux_limiter(i, q_prim_vf, q_sf)
 
         integer, intent(in) :: i
 
@@ -375,7 +375,7 @@ contains
         !!  @param i Vorticity component indicator
         !!  @param q_prim_vf Primitive variables
         !!  @param q_sf Vorticity component
-    subroutine s_derive_vorticity_component(i, q_prim_vf, q_sf)
+    pure subroutine s_derive_vorticity_component(i, q_prim_vf, q_sf)
 
         integer, intent(in) :: i
 
@@ -564,7 +564,7 @@ contains
         !!      variable, q_sf.
         !!  @param q_cons_vf Conservative variables
         !!  @param q_sf Numerical Schlieren function
-    subroutine s_derive_numerical_schlieren_function(q_cons_vf, q_sf)
+    impure subroutine s_derive_numerical_schlieren_function(q_cons_vf, q_sf)
 
         type(scalar_field), &
             dimension(sys_size), &
@@ -691,7 +691,7 @@ contains
     end subroutine s_derive_numerical_schlieren_function
 
     !>  Deallocation procedures for the module
-    subroutine s_finalize_derived_variables_module
+    impure subroutine s_finalize_derived_variables_module
 
         ! Deallocating the variable containing the gradient magnitude of the
         ! density field provided that the numerical Schlieren function was

src/post_process/m_global_parameters.fpp

@@ -324,7 +324,7 @@ contains
     !> Assigns default values to user inputs prior to reading
         !!      them in. This allows for an easier consistency check of
         !!      these parameters once they are read from the input file.
-    subroutine s_assign_default_values_to_user_inputs
+    impure subroutine s_assign_default_values_to_user_inputs
 
         integer :: i !< Generic loop iterator
 
@@ -463,7 +463,7 @@ contains
 
     !>  Computation of parameters, allocation procedures, and/or
         !!      any other tasks needed to properly setup the module
-    subroutine s_initialize_global_parameters_module
+    impure subroutine s_initialize_global_parameters_module
 
         integer :: i, j, fac
 
@@ -885,7 +885,7 @@ contains
     end subroutine s_initialize_global_parameters_module
 
     !> Subroutine to initialize parallel infrastructure
-    subroutine s_initialize_parallel_io
+    impure subroutine s_initialize_parallel_io
 
         num_dims = 1 + min(1, n) + min(1, p)
 
@@ -919,7 +919,7 @@ contains
     end subroutine s_initialize_parallel_io
 
     !> Deallocation procedures for the module
-    subroutine s_finalize_global_parameters_module
+    impure subroutine s_finalize_global_parameters_module
 
         integer :: i
 

src/post_process/m_mpi_proxy.fpp

@@ -48,7 +48,7 @@ contains
 
     !>  Computation of parameters, allocation procedures, and/or
         !!      any other tasks needed to properly setup the module
-    subroutine s_initialize_mpi_proxy_module
+    impure subroutine s_initialize_mpi_proxy_module
 
 #ifdef MFC_MPI
 
@@ -144,7 +144,7 @@ contains
         !!      these are not available to the remaining processors. This
         !!      subroutine is then in charge of broadcasting the required
         !!      information.
-    subroutine s_mpi_bcast_user_inputs
+    impure subroutine s_mpi_bcast_user_inputs
 
 #ifdef MFC_MPI
         integer :: i !< Generic loop iterator
@@ -206,7 +206,7 @@ contains
         !!      as well as recomputing some of the global parameters so
         !!      that they reflect the configuration of sub-domain that
         !!      is overseen by the local processor.
-    subroutine s_mpi_decompose_computational_domain
+    impure subroutine s_mpi_decompose_computational_domain
 
 #ifdef MFC_MPI
 
@@ -646,7 +646,7 @@ contains
         !!      cell-boundary locations is communicated.
         !!  @param pbc_loc Processor boundary condition (PBC) location
         !!  @param sweep_coord Coordinate direction normal to the processor boundary
-    subroutine s_mpi_sendrecv_grid_vars_buffer_regions(pbc_loc, sweep_coord)
+    impure subroutine s_mpi_sendrecv_grid_vars_buffer_regions(pbc_loc, sweep_coord)
 
         character(LEN=3), intent(in) :: pbc_loc
         character, intent(in) :: sweep_coord
@@ -846,7 +846,7 @@ contains
         !!  @param pbc_loc Processor boundary condition (PBC) location
         !!  @param sweep_coord Coordinate direction normal to the processor boundary
         !!  @param q_particle Projection of the lagrangian particles in the Eulerian framework
-    subroutine s_mpi_sendrecv_cons_vars_buffer_regions(q_cons_vf, pbc_loc, &
+    impure subroutine s_mpi_sendrecv_cons_vars_buffer_regions(q_cons_vf, pbc_loc, &
                                                        sweep_coord, q_particle)
 
         type(scalar_field), &
@@ -1496,7 +1496,7 @@ contains
         !!  @param spatial_extents Spatial extents for each processor's sub-domain. First dimension
         !!  corresponds to the minimum and maximum values, respectively, while
         !!  the second dimension corresponds to the processor rank.
-    subroutine s_mpi_gather_spatial_extents(spatial_extents)
+    impure subroutine s_mpi_gather_spatial_extents(spatial_extents)
 
         real(wp), dimension(1:, 0:), intent(INOUT) :: spatial_extents
 
@@ -1615,7 +1615,7 @@ contains
         !!      puts back together the grid of the entire computational
         !!      domain on the rank 0 processor. This is only done for 1D
         !!      simulations.
-    subroutine s_mpi_defragment_1d_grid_variable
+    impure subroutine s_mpi_defragment_1d_grid_variable
 
 #ifdef MFC_MPI
 
@@ -1651,7 +1651,7 @@ contains
         !!   First dimension of array corresponds to the former's minimum and
         !!  maximum values, respectively, while second dimension corresponds
         !!  to each processor's rank.
-    subroutine s_mpi_gather_data_extents(q_sf, data_extents)
+    impure subroutine s_mpi_gather_data_extents(q_sf, data_extents)
 
         real(wp), dimension(:, :, :), intent(in) :: q_sf
 
@@ -1681,7 +1681,7 @@ contains
         !!      This is only done for 1D simulations.
         !!  @param q_sf Flow variable defined on a single computational sub-domain
         !!  @param q_root_sf Flow variable defined on the entire computational domain
-    subroutine s_mpi_defragment_1d_flow_variable(q_sf, q_root_sf)
+    impure subroutine s_mpi_defragment_1d_flow_variable(q_sf, q_root_sf)
 
         real(wp), &
             dimension(0:m), &
@@ -1705,7 +1705,7 @@ contains
     end subroutine s_mpi_defragment_1d_flow_variable
 
     !> Deallocation procedures for the module
-    subroutine s_finalize_mpi_proxy_module
+    impure subroutine s_finalize_mpi_proxy_module
 
 #ifdef MFC_MPI
 

src/post_process/m_start_up.f90

@@ -48,7 +48,7 @@ module m_start_up
     !>  Reads the configuration file post_process.inp, in order
         !!      to populate parameters in module m_global_parameters.f90
         !!      with the user provided inputs
-    subroutine s_read_input_file
+    impure subroutine s_read_input_file
 
         character(LEN=name_len) :: file_loc !<
             !! Generic string used to store the address of a particular file
@@ -125,7 +125,7 @@ end subroutine s_read_input_file
         !!      individual choices are compatible with the code's options
         !!      and that the combination of these choices results into a
         !!      valid configuration for the post-process
-    subroutine s_check_input_file
+    impure subroutine s_check_input_file
 
         character(LEN=len_trim(case_dir)) :: file_loc !<
             !! Generic string used to store the address of a particular file
@@ -151,7 +151,7 @@ subroutine s_check_input_file
 
     end subroutine s_check_input_file
 
-    subroutine s_perform_time_step(t_step)
+    impure subroutine s_perform_time_step(t_step)
 
         integer, intent(inout) :: t_step
         if (proc_rank == 0) then
@@ -189,7 +189,7 @@ subroutine s_perform_time_step(t_step)
 
     end subroutine s_perform_time_step
 
-    subroutine s_save_data(t_step, varname, pres, c, H)
+    impure subroutine s_save_data(t_step, varname, pres, c, H)
 
         integer, intent(inout) :: t_step
         character(LEN=name_len), intent(inout) :: varname
@@ -671,7 +671,7 @@ subroutine s_save_data(t_step, varname, pres, c, H)
 
     end subroutine s_save_data
 
-    subroutine s_initialize_modules
+    impure subroutine s_initialize_modules
         ! Computation of parameters, allocation procedures, and/or any other tasks
         ! needed to properly setup the modules
         call s_initialize_global_parameters_module()
@@ -695,7 +695,7 @@ subroutine s_initialize_modules
         end if
     end subroutine s_initialize_modules
 
-    subroutine s_initialize_mpi_domain
+    impure subroutine s_initialize_mpi_domain
         ! Initialization of the MPI environment
         call s_mpi_initialize()
 
@@ -720,7 +720,7 @@ subroutine s_initialize_mpi_domain
 
     end subroutine s_initialize_mpi_domain
 
-    subroutine s_finalize_modules
+    impure subroutine s_finalize_modules
         ! Disassociate pointers for serial and parallel I/O
         s_read_data_files => null()
 

src/pre_process/m_assign_variables.fpp

@@ -66,7 +66,7 @@ module m_assign_variables
 
 contains
 
-    subroutine s_initialize_assign_variables_module
+    impure subroutine s_initialize_assign_variables_module
 
         allocate (alf_sum%sf(0:m, 0:n, 0:p))
 
@@ -101,7 +101,7 @@ contains
         !! @param eta pseudo volume fraction
         !! @param q_prim_vf Primitive variables
         !! @param patch_id_fp Array to track patch ids
-    subroutine s_assign_patch_mixture_primitive_variables(patch_id, j, k, l, &
+    pure subroutine s_assign_patch_mixture_primitive_variables(patch_id, j, k, l, &
                                                           eta, q_prim_vf, patch_id_fp)
         !$acc routine seq
 
@@ -281,7 +281,7 @@ contains
         !! @param eta pseudo volume fraction
         !! @param q_prim_vf Primitive variables
         !! @param patch_id_fp Array to track patch ids
-    subroutine s_assign_patch_species_primitive_variables(patch_id, j, k, l, &
+    impure subroutine s_assign_patch_species_primitive_variables(patch_id, j, k, l, &
                                                           eta, q_prim_vf, patch_id_fp)
         !$acc routine seq
 
@@ -704,7 +704,7 @@ contains
 
     end subroutine s_assign_patch_species_primitive_variables
 
-    subroutine s_finalize_assign_variables_module
+    impure subroutine s_finalize_assign_variables_module
 
         ! Nullifying procedure pointer to the subroutine assigning either
         ! the patch mixture or species primitive variables to a cell in the

src/pre_process/m_boundary_conditions.fpp

@@ -29,7 +29,7 @@ module m_boundary_conditions
  s_write_parallel_boundary_condition_files
 
 contains
-    subroutine s_line_segment_bc(patch_id, q_prim_vf, bc_type)
+    impure subroutine s_line_segment_bc(patch_id, q_prim_vf, bc_type)
 
         type(scalar_field), dimension(sys_size), intent(inout) :: q_prim_vf
         type(integer_field), dimension(1:num_dims, -1:1), intent(inout) :: bc_type
@@ -79,7 +79,7 @@ contains
 
     end subroutine s_line_segment_bc
 
-    subroutine s_circle_bc(patch_id, q_prim_vf, bc_type)
+    impure subroutine s_circle_bc(patch_id, q_prim_vf, bc_type)
 
         type(scalar_field), dimension(sys_size), intent(inout) :: q_prim_vf
         type(integer_field), dimension(1:num_dims, -1:1), intent(inout) :: bc_type
@@ -143,7 +143,7 @@ contains
 
     end subroutine s_circle_bc
 
-    subroutine s_rectangle_bc(patch_id, q_prim_vf, bc_type)
+    impure subroutine s_rectangle_bc(patch_id, q_prim_vf, bc_type)
 
         type(scalar_field), dimension(sys_size), intent(inout) :: q_prim_vf
         type(integer_field), dimension(1:num_dims, -1:1), intent(inout) :: bc_type
@@ -233,7 +233,7 @@ contains
 
     end subroutine s_rectangle_bc
 
-    subroutine s_apply_boundary_patches(q_prim_vf, bc_type)
+    impure subroutine s_apply_boundary_patches(q_prim_vf, bc_type)
 
         type(scalar_field), dimension(sys_size) :: q_prim_vf
         type(integer_field), dimension(1:num_dims, -1:1) :: bc_type
@@ -267,7 +267,7 @@ contains
 
     end subroutine s_apply_boundary_patches
 
-    subroutine s_write_serial_boundary_condition_files(q_prim_vf, bc_type, step_dirpath)
+    impure subroutine s_write_serial_boundary_condition_files(q_prim_vf, bc_type, step_dirpath)
 
         type(scalar_field), dimension(sys_size) :: q_prim_vf
         type(integer_field), dimension(1:num_dims, -1:1) :: bc_type
@@ -307,7 +307,7 @@ contains
 
     end subroutine s_write_serial_boundary_condition_files
 
-    subroutine s_write_parallel_boundary_condition_files(q_prim_vf, bc_type)
+    impure subroutine s_write_parallel_boundary_condition_files(q_prim_vf, bc_type)
 
         type(scalar_field), dimension(sys_size) :: q_prim_vf
         type(integer_field), dimension(1:num_dims, -1:1) :: bc_type
@@ -369,9 +369,9 @@ contains
 
     end subroutine s_write_parallel_boundary_condition_files
 
-    subroutine s_pack_boundary_condition_buffers(q_prim_vf)
+    impure subroutine s_pack_boundary_condition_buffers(q_prim_vf)
 
-        type(scalar_field), dimension(sys_size) :: q_prim_vf
+        type(scalar_field), dimension(sys_size), intent(in) :: q_prim_vf
         integer :: i, j, k
 
         do k = 0, p