Commit 46358cb7 authored by Femi Kadri's avatar Femi Kadri

modified some files.debug is set to true in most of them.Need to change that…

modified some files.debug is set to true in most of them.Need to change that when result files are getting too large
parent 3d272e8b
......@@ -488,7 +488,7 @@ USE precision_detector
LOGICAL*4 :: LST_marker_new_pos_in_TARGETMASK
LOGICAL*4 :: LST_marker_new_pos_passed_exit_x,LST_marker_old_pos_passed_exit_x
LOGICAL*4,PARAMETER :: LST_use_marker_TARGET_user_mask = .true.
LOGICAL*4,PARAMETER :: LST_use_marker_TARGET_user_mask = .false.
LOGICAL,DIMENSION(:),ALLOCATABLE :: global_VOLMSK_marker_list_INST !INSTANTENOUS LIST OF MARKERS CURRENTLY WITHIN THE USER DEFINED VOLUME MASK
LOGICAL,DIMENSION(:),ALLOCATABLE :: global_TARGET_marker_list_INST !INSTANTENOUS LIST OF MARKERS THAT HAVE HIT THE TARGET (NOT NECESSARILY FOR THE FIRST TIME)
LOGICAL,DIMENSION(:),ALLOCATABLE :: global_TARGET_VIRGIN_list_INST !INSTANTENOUS LIST OF MARKERS THAT HAVE HIT THE TARGET FOR THE 1ST TIME EVER
......@@ -4179,6 +4179,10 @@ WRITE (0,*) My_MPI_Process_ID,'LST_LAGRANGIAN_TRACKING: FINISHED READING FLUID S
WRITE (0,*) My_MPI_Process_ID,'LST_LAGRANGIAN_TRACKING: END MAIN LST TIME STEP LOOP'
WRITE (9,*) My_MPI_Process_ID,'LST_LAGRANGIAN_TRACKING: END MAIN LST TIME STEP LOOP'
WRITE (0,*) My_MPI_Process_ID,'LST_LAGRANGIAN_TRACKING: num_unRXN_PRD_bin_inst_ALLSYSTEM(total_number_of_LST_rxn_bins+1)',num_unRXN_PRD_bin_inst_ALLSYSTEM(total_number_of_LST_rxn_bins+1)
WRITE (9,*) My_MPI_Process_ID,'LST_LAGRANGIAN_TRACKING: num_unRXN_PRD_bin_inst_ALLSYSTEM(total_number_of_LST_rxn_bins+1)',num_unRXN_PRD_bin_inst_ALLSYSTEM(total_number_of_LST_rxn_bins+1)
IF(num_unRXN_PRD_bin_inst_ALLSYSTEM(total_number_of_LST_rxn_bins+1).GT.1) THEN
inv_tot_num_LST_markers_lessone = 1.0 / &
......@@ -4200,7 +4204,7 @@ WRITE (0,*) My_MPI_Process_ID,'LST_LAGRANGIAN_TRACKING: FINISHED READING FLUID S
ENDIF !!num_unRXN_PRD_bin_inst_ALLSYSTEM(total_number_of_LST_rxn_bins+1).EQ. ...
mass_diffusivity = Kinematic_Viscosity / Schmidt_Number
mass_diffusivity_physical_units = mass_diffusivity * ( dx**2 / dt )
mass_diffusivity_physical_units = mass_diffusivity * ( dx**2 / dt )
! effective_mass_diffusivity = &
! & global_mean_square_disp_init_pos(total_number_of_LST_rxn_bins+1) &
......@@ -4210,10 +4214,16 @@ WRITE (0,*) My_MPI_Process_ID,'LST_LAGRANGIAN_TRACKING: FINISHED READING FLUID S
& global_mean_square_disp_mean_pos(total_number_of_LST_rxn_bins+1) &
& * inv_tot_num_LST_markers_lessone / (6.0 * total_LST_runtime_lattice_units)
WRITE (0,*) My_MPI_Process_ID,'LST_LAGRANGIAN_TRACKING: mass_diffusivity',mass_diffusivity
WRITE (9,*) My_MPI_Process_ID,'LST_LAGRANGIAN_TRACKING: mass_diffusivity',mass_diffusivity
IF(My_MPI_Process_ID.EQ.ROOT) THEN
WRITE (0,*) My_MPI_Process_ID,'LST_LAGRANGIAN_TRACKING: mass_diffusivity_physical_units',mass_diffusivity_physical_units
WRITE (9,*) My_MPI_Process_ID,'LST_LAGRANGIAN_TRACKING: mass_diffusivity_physical_units',mass_diffusivity_physical_units
WRITE (0,*) My_MPI_Process_ID,'LST_LAGRANGIAN_TRACKING: effective_mass_diffusivity',effective_mass_diffusivity
WRITE (9,*) My_MPI_Process_ID,'LST_LAGRANGIAN_TRACKING: effective_mass_diffusivity',effective_mass_diffusivity
IF(My_MPI_Process_ID.EQ.ROOT) THEN
PRINT*
PRINT*,'mass_diffusivity',mass_diffusivity,'[lattice units]'
......@@ -4545,10 +4555,10 @@ ELSE !.NOT.LST_previously_converged
IF(write_resume_files) THEN
IF(debug) WRITE (0,*) My_MPI_Process_ID,'LST_LAGRANGIAN_TRACKING: about to CALL LST_soft_binary_REAL_4D_write: LST_local_fluid_surf_rxn_min_max'
IF(debug) WRITE (9,*) My_MPI_Process_ID,'LST_LAGRANGIAN_TRACKING: about to CALL LST_soft_binary_REAL_4D_write: LST_local_fluid_surf_rxn_min_max'
WRITE (0,*) My_MPI_Process_ID,'LST_LAGRANGIAN_TRACKING: about to CALL LST_soft_binary_REAL_4D_write: LST_local_fluid_surf_rxn_min_max'
WRITE (9,*) My_MPI_Process_ID,'LST_LAGRANGIAN_TRACKING: about to CALL LST_soft_binary_REAL_4D_write: LST_local_fluid_surf_rxn_min_max'
CALL LST_soft_binary_REAL_4D_write(cartesian_communicator,&
& rxn_avg+1,local_x_min_ghost,&
& rxn_avg,local_x_min_ghost,&
& local_x_max_ghost,local_y_min_ghost,local_y_max_ghost,&
& local_z_min_ghost,local_z_max_ghost,Print_Results_MPI_ID,&
& My_MPI_Process_ID_Digits,My_MPI_Process_ID,Number_Of_Processes,&
......@@ -4556,27 +4566,27 @@ ELSE !.NOT.LST_previously_converged
IF(debug) WRITE (0,*) My_MPI_Process_ID,'LST_LAGRANGIAN_TRACKING: finished to CALL LST_soft_binary_REAL_4D_write: LST_local_fluid_surf_rxn_min_max'
IF(debug) WRITE (9,*) My_MPI_Process_ID,'LST_LAGRANGIAN_TRACKING: finished to CALL LST_soft_binary_REAL_4D_write: LST_local_fluid_surf_rxn_min_max'
IF(debug) WRITE (0,*) My_MPI_Process_ID,'LST_LAGRANGIAN_TRACKING: about to CALL LST_soft_binary_INTEGER_3D_write: LST_local_fluid_surf_rxn_count'
IF(debug) WRITE (9,*) My_MPI_Process_ID,'LST_LAGRANGIAN_TRACKING: about to CALL LST_soft_binary_INTEGER_3D_write: LST_local_fluid_surf_rxn_count'
WRITE (0,*) My_MPI_Process_ID,'LST_LAGRANGIAN_TRACKING: about to CALL LST_soft_binary_INTEGER_3D_write: LST_local_fluid_surf_rxn_count'
WRITE (9,*) My_MPI_Process_ID,'LST_LAGRANGIAN_TRACKING: about to CALL LST_soft_binary_INTEGER_3D_write: LST_local_fluid_surf_rxn_count'
CALL LST_soft_binary_INTEGER_3D_write(cartesian_communicator,&
& local_x_min_ghost,local_x_max_ghost,&
& local_y_min_ghost,local_y_max_ghost,&
& local_z_min_ghost,local_z_max_ghost,Print_Results_MPI_ID,&
& My_MPI_Process_ID_Digits,My_MPI_Process_ID,Number_Of_Processes,&
& local_product_xyz_nodes_ghost,LST_local_fluid_surf_rxn_count,37,'backup/LST_local_fluid_surf_rxn_count')
IF(debug) WRITE (0,*) My_MPI_Process_ID,'LST_LAGRANGIAN_TRACKING: finished to CALL LST_soft_binary_INTEGER_3D_write: LST_local_fluid_surf_rxn_count'
IF(debug) WRITE (9,*) My_MPI_Process_ID,'LST_LAGRANGIAN_TRACKING: finished to CALL LST_soft_binary_INTEGER_3D_write: LST_local_fluid_surf_rxn_count'
WRITE (0,*) My_MPI_Process_ID,'LST_LAGRANGIAN_TRACKING: finished to CALL LST_soft_binary_INTEGER_3D_write: LST_local_fluid_surf_rxn_count'
WRITE (9,*) My_MPI_Process_ID,'LST_LAGRANGIAN_TRACKING: finished to CALL LST_soft_binary_INTEGER_3D_write: LST_local_fluid_surf_rxn_count'
IF(debug) WRITE (0,*) My_MPI_Process_ID,'LST_LAGRANGIAN_TRACKING: about to CALL LST_soft_binary_INTEGER_3D_write: LST_local_fluid_surf_rxn_t_count'
IF(debug) WRITE (9,*) My_MPI_Process_ID,'LST_LAGRANGIAN_TRACKING: about to CALL LST_soft_binary_INTEGER_3D_write: LST_local_fluid_surf_rxn_t_count'
WRITE (0,*) My_MPI_Process_ID,'LST_LAGRANGIAN_TRACKING: about to CALL LST_soft_binary_INTEGER_3D_write: LST_local_fluid_surf_rxn_t_count'
WRITE (9,*) My_MPI_Process_ID,'LST_LAGRANGIAN_TRACKING: about to CALL LST_soft_binary_INTEGER_3D_write: LST_local_fluid_surf_rxn_t_count'
CALL LST_soft_binary_INTEGER_3D_write(cartesian_communicator,&
& local_x_min_ghost,local_x_max_ghost,&
& local_y_min_ghost,local_y_max_ghost,&
& local_z_min_ghost,local_z_max_ghost,Print_Results_MPI_ID,&
& My_MPI_Process_ID_Digits,My_MPI_Process_ID,Number_Of_Processes,&
& local_product_xyz_nodes_ghost,LST_local_fluid_surf_rxn_t_count,39,'backup/LST_local_fluid_surf_rxn_t_count')
IF(debug) WRITE (0,*) My_MPI_Process_ID,'LST_LAGRANGIAN_TRACKING: finished to CALL LST_soft_binary_INTEGER_3D_write: LST_local_fluid_surf_rxn_t_count'
IF(debug) WRITE (9,*) My_MPI_Process_ID,'LST_LAGRANGIAN_TRACKING: finished to CALL LST_soft_binary_INTEGER_3D_write: LST_local_fluid_surf_rxn_t_count'
WRITE (0,*) My_MPI_Process_ID,'LST_LAGRANGIAN_TRACKING: finished to CALL LST_soft_binary_INTEGER_3D_write: LST_local_fluid_surf_rxn_t_count'
WRITE (9,*) My_MPI_Process_ID,'LST_LAGRANGIAN_TRACKING: finished to CALL LST_soft_binary_INTEGER_3D_write: LST_local_fluid_surf_rxn_t_count'
ENDIF !!write_resume_files
......
......@@ -139,6 +139,9 @@ USE precision_detector
! REAL,INTENT(OUT) :: global_avg_absolute_position_x,global_avg_absolute_position_y,global_avg_absolute_position_z
!-----------NO INTENT------------------------------------------
LOGICAL*4,DIMENSION(:,:,:), ALLOCATABLE :: NOT_LST_local_walls
LOGICAL*4,DIMENSION(:,:,:), ALLOCATABLE :: NOT_LST_global_walls_release_plane
INTEGER :: LST_marker_loop_counter,LST_marker_loop_index
INTEGER :: vertex000_x,vertex100_x,vertex010_x,vertex001_x,vertex101_x,vertex011_x,vertex110_x,vertex111_x
INTEGER :: vertex000_y,vertex100_y,vertex010_y,vertex001_y,vertex101_y,vertex011_y,vertex110_y,vertex111_y
......@@ -152,8 +155,6 @@ USE precision_detector
INTEGER :: mpi_abort_error,mpi_error_code
CHARACTER(LEN = My_MPI_Process_ID_Digits) Print_Results_MPI_ID
IF(debug) WRITE (0,*) My_MPI_Process_ID,'entered SUBROUTINE LST_initialize_marker_positions'
......@@ -425,7 +426,6 @@ USE precision_detector
IF(.NOT.spread_LST_markers_in_a_plane) THEN !SPREAD MARKERS UNIFORMLY (NOT IN A PLANE AND NOT A POINT RELEASE)
IF(surface_release_of_LST_markers) THEN !!surface_release_of_LST_markers (UNIFORM/SURFACE RELEASE)
IF(debug) WRITE (0,*) My_MPI_Process_ID,'LST_initialize_marker_positions: about to CALL LST_uniform_marker_release'
......@@ -449,7 +449,7 @@ USE precision_detector
& local_x_nodes_real,local_y_nodes_real,local_z_nodes_real,&
& local_x_nodes_ghost,local_y_nodes_ghost,local_z_nodes_ghost,&
& global_x_nodes_real,global_y_nodes_real,global_z_nodes_real,&
& global_x_nodes_ghost,global_y_nodes_ghost,global_z_nodes_ghost,&
& global_x_nodes_ghost,global_y_nodes_ghost,global_z_nodes_ghost,&
& LST_apply_volume_user_mask,&
& My_MPI_Process_ID_Digits,LST_local_walls, &
& Print_Results_MPI_ID,my_mpi_real_data_type,&
......@@ -525,6 +525,25 @@ USE precision_detector
!SO THE SUBROUTINE WILL TREAT EVERYTHING THAT IS NOT WALLS AS A "SURFACE"
!AND WILL INSERT LST MARKERS ONLY INTO THE SURFACE SITES (I.E. FLUID NODES)
IF(debug) WRITE (0,*) My_MPI_Process_ID,'LST_initialize_marker_positions: about to ALLOCATE NOT_LST_local_walls'
IF(debug) WRITE (9,*) My_MPI_Process_ID,'LST_initialize_marker_positions: about to ALLOCATE NOT_LST_local_walls'
ALLOCATE ( NOT_LST_local_walls(local_x_min_ghost:local_x_max_ghost,&
& local_y_min_ghost:local_y_max_ghost, local_z_min_ghost:local_z_max_ghost),STAT=memory_status )
IF (memory_status /= memory_success) THEN
WRITE (0,*) My_MPI_Process_ID, &
& "LST_initialize_marker_positions ERROR: can't allocate NOT_LST_local_walls of dimensions ", &
& local_x_min_ghost, ":", local_x_max_ghost, ",", &
& local_y_min_ghost, ":", local_y_max_ghost, ",", &
& local_z_min_ghost, ":", local_z_max_ghost, ")"
CALL MPI_Abort(cartesian_communicator,mpi_abort_error,mpi_error_code)
END IF !! (memory_status /= memory_success)
IF(debug) WRITE (0,*) My_MPI_Process_ID,'LST_initialize_marker_positions: ALLOCATED NOT_LST_local_walls'
IF(debug) WRITE (9,*) My_MPI_Process_ID,'LST_initialize_marker_positions: ALLOCATED NOT_LST_local_walls'
NOT_LST_local_walls = .NOT.LST_local_walls
IF(debug) WRITE (0,*) My_MPI_Process_ID,'LST_initialize_marker_positions: about to CALL LST_uniform_marker_release'
IF(debug) WRITE (9,*) My_MPI_Process_ID,'LST_initialize_marker_positions: about to CALL LST_uniform_marker_release'
......@@ -547,16 +566,17 @@ USE precision_detector
& local_x_nodes_real,local_y_nodes_real,local_z_nodes_real,&
& local_x_nodes_ghost,local_y_nodes_ghost,local_z_nodes_ghost,&
& global_x_nodes_real,global_y_nodes_real,global_z_nodes_real,&
& global_x_nodes_ghost,global_y_nodes_ghost,global_z_nodes_ghost,&
& global_x_nodes_ghost,global_y_nodes_ghost,global_z_nodes_ghost,&
& LST_apply_volume_user_mask,&
& My_MPI_Process_ID_Digits,LST_local_walls, &
& Print_Results_MPI_ID,my_mpi_real_data_type,&
& global_initial_marker_IDs,global_initial_positions_lu,&
& .NOT.LST_local_walls)
& NOT_LST_local_walls)
IF(debug) WRITE (0,*) My_MPI_Process_ID,'LST_initialize_marker_positions: FINISHED to CALL LST_uniform_marker_release'
IF(debug) WRITE (9,*) My_MPI_Process_ID,'LST_initialize_marker_positions: FINISHED to CALL LST_uniform_marker_release'
DEALLOCATE ( NOT_LST_local_walls )
ENDIF !! surface_release_of_LST_markers
......@@ -577,6 +597,26 @@ USE precision_detector
IF(.NOT.surface_release_of_LST_markers) THEN !REGULAR PLANE RELEASE
IF(debug) WRITE (0,*) My_MPI_Process_ID,'LST_initialize_marker_positions: about to ALLOCATE NOT_LST_global_walls_release_plane'
IF(debug) WRITE (9,*) My_MPI_Process_ID,'LST_initialize_marker_positions: about to ALLOCATE NOT_LST_global_walls_release_plane'
ALLOCATE ( NOT_LST_global_walls_release_plane(0:1,global_plane_dim2_min_ghost:global_plane_dim2_max_ghost,&
& global_plane_dim3_min_ghost:global_plane_dim3_max_ghost),STAT=memory_status )
IF (memory_status /= memory_success) THEN
WRITE (0,*) My_MPI_Process_ID, &
& "LST_initialize_marker_positions ERROR: can't allocate NOT_LST_global_walls_release_plane of dimensions ", &
& 0, ":", 1, ",", &
& global_plane_dim2_min_ghost, ":", global_plane_dim2_max_ghost, ",", &
& global_plane_dim3_min_ghost, ":", global_plane_dim3_max_ghost, ")"
CALL MPI_Abort(cartesian_communicator,mpi_abort_error,mpi_error_code)
END IF !! (memory_status /= memory_success)
IF(debug) WRITE (0,*) My_MPI_Process_ID,'LST_initialize_marker_positions: ALLOCATED NOT_LST_global_walls_release_plane'
IF(debug) WRITE (9,*) My_MPI_Process_ID,'LST_initialize_marker_positions: ALLOCATED NOT_LST_global_walls_release_plane'
NOT_LST_global_walls_release_plane = .NOT.LST_global_walls_release_plane
IF(debug) WRITE (0,*) My_MPI_Process_ID,'LST_initialize_marker_positions: about to CALL LST_uniform_plane_marker_release'
IF(debug) WRITE (9,*) My_MPI_Process_ID,'LST_initialize_marker_positions: about to CALL LST_uniform_plane_marker_release'
......@@ -604,11 +644,13 @@ USE precision_detector
& global_plane_dim2_min_real,global_plane_dim2_max_real,&
& global_plane_dim3_min_real,global_plane_dim3_max_real,&
& global_walls_surf_rel_dim1_min,global_walls_surf_rel_dim1_max,&
& .NOT.LST_global_walls_release_plane )
& NOT_LST_global_walls_release_plane)
IF(debug) WRITE (0,*) My_MPI_Process_ID,'LST_initialize_marker_positions: FINISHED to CALL LST_uniform_plane_marker_release'
IF(debug) WRITE (9,*) My_MPI_Process_ID,'LST_initialize_marker_positions: FINISHED to CALL LST_uniform_plane_marker_release'
DEALLOCATE ( NOT_LST_global_walls_release_plane )
ELSE !!surface_release_of_LST_markers (PLANE/SURFACE RELEASE)
IF(debug) WRITE (0,*) My_MPI_Process_ID,'LST_initialize_marker_positions: about to CALL LST_uniform_plane_marker_release'
......@@ -677,4 +719,4 @@ USE precision_detector
RETURN
END SUBROUTINE LST_initialize_marker_positions
\ No newline at end of file
END SUBROUTINE LST_initialize_marker_positions
......@@ -30,7 +30,8 @@ SUBROUTINE LST_soft_binary_REAL_4D_write(cartesian_communicator,&
& local_x_min_ghost:local_x_max_ghost,&
& local_y_min_ghost:local_y_max_ghost,&
& local_z_min_ghost:local_z_max_ghost) :: real_4D_array
LOGICAL*4,PARAMETER :: debug = GLOBAL_DEBUG
!LOGICAL*4,PARAMETER :: debug = GLOBAL_DEBUG
LOGICAL*4,PARAMETER :: debug = .true.
!-----------INTENT(OUT)------------------------------------------
!-----------NO INTENT--------------------------------------------
INTEGER :: K,mpi_error_code
......@@ -69,4 +70,4 @@ SUBROUTINE LST_soft_binary_REAL_4D_write(cartesian_communicator,&
IF(debug) WRITE (0,*) My_MPI_Process_ID,'LST_soft_binary_REAL_4D_write: after barrier'
RETURN
END SUBROUTINE LST_soft_binary_REAL_4D_write
\ No newline at end of file
END SUBROUTINE LST_soft_binary_REAL_4D_write
......@@ -240,8 +240,8 @@ USE debug_module
LOGICAL,INTENT(IN),DIMENSION(tot_num_markers_plus_nonbrownian) :: global_TARGET_marker_list_CUMLTV
LOGICAL,INTENT(IN) :: LST_apply_volume_user_mask,LST_use_marker_TARGET_user_mask
LOGICAL*4,PARAMETER :: debug = GLOBAL_DEBUG
!LOGICAL*4,PARAMETER :: debug = .TRUE.
!LOGICAL*4,PARAMETER :: debug = GLOBAL_DEBUG
LOGICAL*4,PARAMETER :: debug = .TRUE.
!-----------INTENT(OUT)------------------------------------------
!-----------INTENT(INOUT)------------------------------------------
INTEGER,INTENT(INOUT) :: tot_reacted_inall_bins_increment
......@@ -965,7 +965,8 @@ USE debug_module
num_COLNS_PRD_bin_AVG_VRGN_TRGT(LST_reaction_bin_number) = &
& REAL(num_COLNS_PRD_bin_CUMM_VRGN_TRGT(LST_reaction_bin_number)) * Inv_LST_Print_Results_Incriment
IF(debug) WRITE (0,*) My_MPI_Process_ID,'LST_Time_Steps_Prev_Completed',LST_Time_Steps_Prev_Completed,'LST_write_output: num_unRXN_ABS_bin_AVG_ALLSYSTEM(LST_reaction_bin_number)',num_unRXN_ABS_bin_AVG_ALLSYSTEM(LST_reaction_bin_number)
IF(debug) WRITE (9,*) My_MPI_Process_ID,'LST_Time_Steps_Prev_Completed',LST_Time_Steps_Prev_Completed,'LST_write_output: num_unRXN_ABS_bin_AVG_ALLSYSTEM(LST_reaction_bin_number)',num_unRXN_ABS_bin_AVG_ALLSYSTEM(LST_reaction_bin_number)
IF(num_unRXN_ABS_bin_AVG_ALLSYSTEM(LST_reaction_bin_number).GT.0) THEN
......@@ -977,6 +978,8 @@ USE debug_module
ENDIF !num_unRXN_ABS_bin_AVG_ALLSYSTEM(LST_reaction_bin_number).GT.0
IF(debug) WRITE (0,*) My_MPI_Process_ID,'LST_Time_Steps_Prev_Completed',LST_Time_Steps_Prev_Completed,'LST_write_output: num_unRXN_PRD_bin_AVG_ALLSYSTEM(LST_reaction_bin_number)',num_unRXN_PRD_bin_AVG_ALLSYSTEM(LST_reaction_bin_number)
IF(debug) WRITE (9,*) My_MPI_Process_ID,'LST_Time_Steps_Prev_Completed',LST_Time_Steps_Prev_Completed,'LST_write_output: num_unRXN_PRD_bin_AVG_ALLSYSTEM(LST_reaction_bin_number)',num_unRXN_PRD_bin_AVG_ALLSYSTEM(LST_reaction_bin_number)
IF(num_unRXN_PRD_bin_AVG_ALLSYSTEM(LST_reaction_bin_number).GT.0) THEN
......@@ -988,6 +991,8 @@ USE debug_module
ENDIF !num_unRXN_PRD_bin_AVG_ALLSYSTEM(LST_reaction_bin_number).GT.0
IF(debug) WRITE (0,*) My_MPI_Process_ID,'LST_Time_Steps_Prev_Completed',LST_Time_Steps_Prev_Completed,'LST_write_output: num_unrxn_ABS_bin_AVG_VOLMASK(LST_reaction_bin_number)',num_unrxn_ABS_bin_AVG_VOLMASK(LST_reaction_bin_number)
IF(debug) WRITE (9,*) My_MPI_Process_ID,'LST_Time_Steps_Prev_Completed',LST_Time_Steps_Prev_Completed,'LST_write_output: num_unrxn_ABS_bin_AVG_VOLMASK(LST_reaction_bin_number)',num_unrxn_ABS_bin_AVG_VOLMASK(LST_reaction_bin_number)
IF(num_unrxn_ABS_bin_AVG_VOLMASK(LST_reaction_bin_number).GT.0) THEN
......@@ -999,6 +1004,9 @@ USE debug_module
ENDIF !num_unrxn_ABS_bin_AVG_VOLMASK(LST_reaction_bin_number).GT.0
IF(debug) WRITE (0,*) My_MPI_Process_ID,'LST_Time_Steps_Prev_Completed',LST_Time_Steps_Prev_Completed,'LST_write_output: num_unrxn_PRD_bin_AVG_VOLMASK(LST_reaction_bin_number)',num_unrxn_PRD_bin_AVG_VOLMASK(LST_reaction_bin_number)
IF(debug) WRITE (9,*) My_MPI_Process_ID,'LST_Time_Steps_Prev_Completed',LST_Time_Steps_Prev_Completed,'LST_write_output: num_unrxn_PRD_bin_AVG_VOLMASK(LST_reaction_bin_number)',num_unrxn_PRD_bin_AVG_VOLMASK(LST_reaction_bin_number)
IF(num_unrxn_PRD_bin_AVG_VOLMASK(LST_reaction_bin_number).GT.0) THEN
INV_num_unrxn_PRD_bin_AVG_VMASK = 1 / (num_unrxn_PRD_bin_AVG_VOLMASK(LST_reaction_bin_number))
......@@ -1009,6 +1017,8 @@ USE debug_module
ENDIF !num_unrxn_PRD_bin_AVG_VOLMASK(LST_reaction_bin_number).GT.0
IF(debug) WRITE (0,*) My_MPI_Process_ID,'LST_Time_Steps_Prev_Completed',LST_Time_Steps_Prev_Completed,'LST_write_output: num_unrxn_ABS_bin_AVG_AT_TARGET(LST_reaction_bin_number)',num_unrxn_ABS_bin_AVG_AT_TARGET(LST_reaction_bin_number)
IF(debug) WRITE (9,*) My_MPI_Process_ID,'LST_Time_Steps_Prev_Completed',LST_Time_Steps_Prev_Completed,'LST_write_output: num_unrxn_ABS_bin_AVG_AT_TARGET(LST_reaction_bin_number)',num_unrxn_ABS_bin_AVG_AT_TARGET(LST_reaction_bin_number)
IF(num_unrxn_ABS_bin_AVG_AT_TARGET(LST_reaction_bin_number).GT.0) THEN
......@@ -1020,6 +1030,9 @@ USE debug_module
ENDIF !num_unrxn_ABS_bin_AVG_AT_TARGET(LST_reaction_bin_number).GT.0
IF(debug) WRITE (0,*) My_MPI_Process_ID,'LST_Time_Steps_Prev_Completed',LST_Time_Steps_Prev_Completed,'LST_write_output: num_unrxn_PRD_bin_AVG_AT_TARGET(LST_reaction_bin_number)',num_unrxn_PRD_bin_AVG_AT_TARGET(LST_reaction_bin_number)
IF(debug) WRITE (9,*) My_MPI_Process_ID,'LST_Time_Steps_Prev_Completed',LST_Time_Steps_Prev_Completed,'LST_write_output: num_unrxn_PRD_bin_AVG_AT_TARGET(LST_reaction_bin_number)',num_unrxn_PRD_bin_AVG_AT_TARGET(LST_reaction_bin_number)
IF(num_unrxn_PRD_bin_AVG_AT_TARGET(LST_reaction_bin_number).GT.0) THEN
INV_num_unrxn_PRD_bin_AVG_TARGET = 1 / (num_unrxn_PRD_bin_AVG_AT_TARGET(LST_reaction_bin_number))
......@@ -1030,6 +1043,9 @@ USE debug_module
ENDIF !num_unrxn_PRD_bin_AVG_AT_TARGET(LST_reaction_bin_number).GT.0
IF(debug) WRITE (0,*) My_MPI_Process_ID,'LST_Time_Steps_Prev_Completed',LST_Time_Steps_Prev_Completed,'LST_write_output: num_unrxn_ABS_bin_AVG_VRGN_TRGT(LST_reaction_bin_number)',num_unrxn_ABS_bin_AVG_VRGN_TRGT(LST_reaction_bin_number)
IF(debug) WRITE (9,*) My_MPI_Process_ID,'LST_Time_Steps_Prev_Completed',LST_Time_Steps_Prev_Completed,'LST_write_output: num_unrxn_ABS_bin_AVG_VRGN_TRGT(LST_reaction_bin_number)',num_unrxn_ABS_bin_AVG_VRGN_TRGT(LST_reaction_bin_number)
IF(num_unrxn_ABS_bin_AVG_VRGN_TRGT(LST_reaction_bin_number).GT.0) THEN
INV_num_unrxn_ABS_AVG_VRGN_TRGT = 1 / (num_unrxn_ABS_bin_AVG_VRGN_TRGT(LST_reaction_bin_number))
......@@ -1040,6 +1056,9 @@ USE debug_module
ENDIF !num_unrxn_ABS_bin_AVG_VRGN_TRGT(LST_reaction_bin_number).GT.0
IF(debug) WRITE (0,*) My_MPI_Process_ID,'LST_Time_Steps_Prev_Completed',LST_Time_Steps_Prev_Completed,'LST_write_output: num_unrxn_PRD_bin_AVG_VRGN_TRGT(LST_reaction_bin_number)',num_unrxn_PRD_bin_AVG_VRGN_TRGT(LST_reaction_bin_number)
IF(debug) WRITE (9,*) My_MPI_Process_ID,'LST_Time_Steps_Prev_Completed',LST_Time_Steps_Prev_Completed,'LST_write_output: num_unrxn_PRD_bin_AVG_VRGN_TRGT(LST_reaction_bin_number)',num_unrxn_PRD_bin_AVG_VRGN_TRGT(LST_reaction_bin_number)
IF(num_unrxn_PRD_bin_AVG_VRGN_TRGT(LST_reaction_bin_number).GT.0) THEN
INV_num_unrxn_PRD_AVG_VRGN_TRGT = 1 / (num_unrxn_PRD_bin_AVG_VRGN_TRGT(LST_reaction_bin_number))
......@@ -1429,4 +1448,4 @@ USE debug_module
IF(debug) WRITE (9,*) My_MPI_Process_ID,'about to END SUBROUTINE LST_write_output'
RETURN
END SUBROUTINE LST_write_output
\ No newline at end of file
END SUBROUTINE LST_write_output
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment