Commit 5cf72f42 authored by Femi Kadri's avatar Femi Kadri

updated for pgi compiler

parent 3c7d677d
......@@ -18,6 +18,7 @@
USE debug_module
USE ieee_arithmetic
!USE precision_detector
IMPLICIT NONE
......@@ -132,7 +133,7 @@ USE debug_module
!REAL,DIMENSION(:,:,:),ALLOCATABLE :: local_density_rho,local_velocity_x,local_velocity_y,local_velocity_z
REAL,DIMENSION(:,:,:),ALLOCATABLE :: local_velocity_x
LOGICAL :: ISNAN
!LOGICAL :: ISNAN
INCLUDE "mpif.h"
......@@ -772,9 +773,9 @@ ENDIF !! DEBUG
!CHECK FOR DIVERGENCE (INFINITY) AND STOP IF OCCURED
IF ( ISNAN(global_Ux_max).OR.ISNAN(global_Ux_avg).OR.ISNAN(global_Ux_min) &
& .OR.ISNAN(global_Uy_max).OR.ISNAN(global_Uy_avg).OR.ISNAN(global_Uy_min) &
& .OR.ISNAN(global_Uz_max).OR.ISNAN(global_Uz_avg).OR.ISNAN(global_Uz_min) ) THEN
IF ( IEEE_IS_NAN(global_Ux_max).OR.IEEE_IS_NAN(global_Ux_avg).OR.IEEE_IS_NAN(global_Ux_min) &
& .OR.IEEE_IS_NAN(global_Uy_max).OR.IEEE_IS_NAN(global_Uy_avg).OR.IEEE_IS_NAN(global_Uy_min) &
& .OR.IEEE_IS_NAN(global_Uz_max).OR.IEEE_IS_NAN(global_Uz_avg).OR.IEEE_IS_NAN(global_Uz_min) ) THEN
PRINT*, 'LBM_equilibrate_rhovel ERROR: SIMULATION DIVERGED (TRY CHANGING TAU OR REDUCING Forcing_Factor_x, SUCH THAT Re <10)'
WRITE (0,*) 'LBM_equilibrate_rhovel ERROR: SIMULATION DIVERGED (TRY CHANGING TAU OR REDUCING Forcing_Factor_x, SUCH THAT Re <10)'
......
This diff is collapsed.
......@@ -2964,7 +2964,7 @@ WRITE (0,*) My_MPI_Process_ID,'LST_LAGRANGIAN_TRACKING: FINISHED READING FLUID S
& WRITE(76,*) periodic_marker_old_position_x * dx_in_microns,&
& periodic_marker_old_position_y * dx_in_microns,&
& periodic_marker_old_position_z * dx_in_microns,&
& & LST_this_marker_avgrxn_const_now
& LST_this_marker_avgrxn_const_now
ENDIF !!write_LST_trajectories.AND.My_MPI_Process_ID.EQ.ANTIROOT.AND..NOT.marker_ID_is_non_brownian
......@@ -4179,10 +4179,6 @@ 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 / &
......@@ -4204,7 +4200,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) &
......@@ -4214,16 +4210,10 @@ 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
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
IF(My_MPI_Process_ID.EQ.ROOT) THEN
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]'
......@@ -4555,10 +4545,10 @@ ELSE !.NOT.LST_previously_converged
IF(write_resume_files) THEN
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'
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'
CALL LST_soft_binary_REAL_4D_write(cartesian_communicator,&
& rxn_avg,local_x_min_ghost,&
& rxn_avg+1,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,&
......@@ -4566,27 +4556,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'
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'
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'
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')
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: 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: 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'
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'
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')
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'
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'
ENDIF !!write_resume_files
......
......@@ -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
......
......@@ -422,7 +422,7 @@ USE debug_module
num_COLLSNS_PRD_bin_inst_VOLMASK(LST_reaction_bin_number) = &
& num_COLLSNS_PRD_bin_inst_VOLMASK(LST_reaction_bin_number) + &
& & global_wall_collision_INST_count(LST_marker_ID)
& global_wall_collision_INST_count(LST_marker_ID)
IF(THIS_maker_is_in_PERIODIC_domain) THEN !MARKER IS ABSOLUTELY IN THE VOLUME MASK
......@@ -669,4 +669,4 @@ USE debug_module
IF(debug) WRITE (9,*) My_MPI_Process_ID,'about to END SUBROUTINE LST_calc_abs_position_averages'
RETURN
END SUBROUTINE LST_calc_abs_position_averages
\ No newline at end of file
END SUBROUTINE LST_calc_abs_position_averages
This diff is collapsed.
......@@ -12,6 +12,7 @@ SUBROUTINE LST_time_incriment(local_velocity_field,cartesian_communicator,ROOT,S
USE debug_module
USE precision_detector
USE ieee_arithmetic
IMPLICIT NONE
......@@ -196,7 +197,7 @@ USE precision_detector
& number_of_std_dev_squared * kinematic_viscosity*Schmidt_Number*global_velocity_magnitude_max) ) / &
& (2.0d0 * Schmidt_Number * global_velocity_magnitude_max**2)
IF( LST_time_increment_max.LE.0.0.OR.ISNAN(LST_time_increment_max) ) THEN
IF( LST_time_increment_max.LE.0.0.OR.IEEE_IS_NAN(LST_time_increment_max) ) THEN
!THIS WAS THE SECOND (LESS CONSERVATIVE) SOLUTION, SINCE ITS A QUADRATIC EQUATION
!(ONLY USE THIS IF THE FIRST SOLUTION GIVES NONPHYSICAL NUMBERS FOR LST_time_increment_max)
......@@ -221,7 +222,7 @@ USE precision_detector
WRITE (29,*) My_MPI_Process_ID,'LST_time_incriment ERROR: --> LST_time_increment_max.LE.0.0 <--'
STOP
ELSEIF (ISNAN(LST_time_increment_max)) THEN
ELSEIF (IEEE_IS_NAN(LST_time_increment_max)) THEN
PRINT*, My_MPI_Process_ID,'LST_time_incriment ERROR: --> Variable LST_time_increment_max contains a NaN value <--'
IF(debug) WRITE (0,*) My_MPI_Process_ID,'LST_time_incriment ERROR: --> Variable LST_time_increment_max contains a NaN value <--'
IF(debug) WRITE (9,*) My_MPI_Process_ID,'LST_time_incriment ERROR: --> Variable LST_time_increment_max contains a NaN value <--'
......@@ -299,4 +300,4 @@ USE precision_detector
IF(debug) WRITE (9,*) My_MPI_Process_ID,'about to END SUBROUTINE LST_time_incriment'
RETURN
END SUBROUTINE LST_time_incriment
\ No newline at end of file
END SUBROUTINE LST_time_incriment
This diff is collapsed.
......@@ -9,6 +9,7 @@ SUBROUTINE numerical_2nd_deriv_2nd_order_1D(&
!NUMERICAL APPROXIMATION OF THE 2ND DERIVATIVE USING 4TH ORDER CENTRAL DIFFERENCE FORMULA
USE debug_module
USE ieee_arithmetic
!USE precision_detector
IMPLICIT NONE
......@@ -25,7 +26,7 @@ USE debug_module
second_derivative_F = ( forward_F - REAL(2) * F + backward_F )
IF( ISNAN(second_derivative_F) ) THEN
IF( IEEE_IS_NAN(second_derivative_F) ) THEN
IF(debug) WRITE (0,*) My_MPI_Process_ID,'numerical_2nd_deriv_2nd_order_1D: forward_F',forward_F
IF(debug) WRITE (9,*) My_MPI_Process_ID,'numerical_2nd_deriv_2nd_order_1D: forward_F',forward_F
......
SUBROUTINE numerical_2nd_deriv_2nd_order_1D(&
& My_MPI_Process_ID,&
& backward_F,&
& F,&
& forward_F,&
& second_derivative_F)
!NUMERICAL APPROXIMATION OF THE 2ND DERIVATIVE USING 4TH ORDER CENTRAL DIFFERENCE FORMULA
USE debug_module
USE ieee_arithmetic
!USE precision_detector
IMPLICIT NONE
!REAL,PARAMETER :: Inverse_3 = 1.0 / 3.0
!-----------INTENT(IN)-----------------------------------------------------------------------------------------------
INTEGER,INTENT(IN) :: My_MPI_Process_ID
REAL,INTENT(IN) :: forward_F,backward_F,F
LOGICAL*4,PARAMETER :: debug = GLOBAL_DEBUG
!-----------INTENT(OUT)-----------------------------------------------------------------------------------------------
REAL,INTENT(OUT) :: second_derivative_F
IF(debug) WRITE (0,*) My_MPI_Process_ID,'entered SUBROUTINE numerical_2nd_deriv_2nd_order_1D'
IF(debug) WRITE (9,*) My_MPI_Process_ID,'entered SUBROUTINE numerical_2nd_deriv_2nd_order_1D'
second_derivative_F = ( forward_F - REAL(2) * F + backward_F )
IF( IEEE_IS_NAN(second_derivative_F) ) THEN
IF(debug) WRITE (0,*) My_MPI_Process_ID,'numerical_2nd_deriv_2nd_order_1D: forward_F',forward_F
IF(debug) WRITE (9,*) My_MPI_Process_ID,'numerical_2nd_deriv_2nd_order_1D: forward_F',forward_F
IF(debug) WRITE (0,*) My_MPI_Process_ID,'numerical_2nd_deriv_2nd_order_1D: backward_F',backward_F
IF(debug) WRITE (9,*) My_MPI_Process_ID,'numerical_2nd_deriv_2nd_order_1D: backward_F',backward_F
ENDIF !!ISNAN(second_derivative_F)
IF(debug) WRITE (9,*) My_MPI_Process_ID,'numerical_2nd_deriv_2nd_order_1D: about to return from SUBROUTINE'
IF(debug) WRITE (0,*) My_MPI_Process_ID,'numerical_2nd_deriv_2nd_order_1D: about to return from SUBROUTINE'
RETURN
END SUBROUTINE numerical_2nd_deriv_2nd_order_1D
......@@ -11,6 +11,7 @@ SUBROUTINE numerical_2nd_deriv_4th_order_1D(&
!NUMERICAL APPROXIMATION OF THE 2ND DERIVATIVE USING 4TH ORDER CENTRAL DIFFERENCE FORMULA
USE debug_module
USE ieee_arithmetic
!USE precision_detector
IMPLICIT NONE
......@@ -31,7 +32,7 @@ USE debug_module
& - REAL(30) * F + REAL(16) * backward_F - backward_backward_F ) * Inverse_12
IF( ISNAN(second_derivative_F) ) THEN
IF( IEEE_IS_NAN(second_derivative_F) ) THEN
IF(debug) WRITE (0,*) My_MPI_Process_ID,'numerical_2nd_deriv_4th_order_1D: forward_forward_F',forward_forward_F
IF(debug) WRITE (9,*) My_MPI_Process_ID,'numerical_2nd_deriv_4th_order_1D: forward_forward_F',forward_forward_F
......
SUBROUTINE numerical_2nd_deriv_4th_order_1D(&
& My_MPI_Process_ID,&
& backward_backward_F,&
& backward_F,&
& F,&
& forward_F,&
& forward_forward_F,&
& second_derivative_F)
!NUMERICAL APPROXIMATION OF THE 2ND DERIVATIVE USING 4TH ORDER CENTRAL DIFFERENCE FORMULA
USE debug_module
!USE precision_detector
IMPLICIT NONE
REAL,PARAMETER :: Inverse_12 = 1.0 / 12.0
!-----------INTENT(IN)-----------------------------------------------------------------------------------------------
INTEGER,INTENT(IN) :: My_MPI_Process_ID
REAL,INTENT(IN) :: forward_forward_F,forward_F,backward_backward_F,backward_F,F
LOGICAL*4,PARAMETER :: debug = GLOBAL_DEBUG
!-----------INTENT(OUT)-----------------------------------------------------------------------------------------------
REAL,INTENT(OUT) :: second_derivative_F
IF(debug) WRITE (0,*) My_MPI_Process_ID,'entered SUBROUTINE numerical_2nd_deriv_4th_order_1D'
IF(debug) WRITE (9,*) My_MPI_Process_ID,'entered SUBROUTINE numerical_2nd_deriv_4th_order_1D'
second_derivative_F = ( - forward_forward_F + REAL(16) * forward_F &
& - REAL(30) * F + REAL(16) * backward_F - backward_backward_F ) * Inverse_12
IF( ISNAN(second_derivative_F) ) THEN
IF(debug) WRITE (0,*) My_MPI_Process_ID,'numerical_2nd_deriv_4th_order_1D: forward_forward_F',forward_forward_F
IF(debug) WRITE (9,*) My_MPI_Process_ID,'numerical_2nd_deriv_4th_order_1D: forward_forward_F',forward_forward_F
IF(debug) WRITE (0,*) My_MPI_Process_ID,'numerical_2nd_deriv_4th_order_1D: forward_F',forward_F
IF(debug) WRITE (9,*) My_MPI_Process_ID,'numerical_2nd_deriv_4th_order_1D: forward_F',forward_F
IF(debug) WRITE (0,*) My_MPI_Process_ID,'numerical_2nd_deriv_4th_order_1D: backward_F',backward_F
IF(debug) WRITE (9,*) My_MPI_Process_ID,'numerical_2nd_deriv_4th_order_1D: backward_F',backward_F
IF(debug) WRITE (0,*) My_MPI_Process_ID,'numerical_2nd_deriv_4th_order_1D: backward_backward_F',backward_backward_F
IF(debug) WRITE (9,*) My_MPI_Process_ID,'numerical_2nd_deriv_4th_order_1D: backward_backward_F',backward_backward_F
IF(debug) WRITE (0,*) My_MPI_Process_ID,'numerical_2nd_deriv_4th_order_1D: forward_F',F
IF(debug) WRITE (9,*) My_MPI_Process_ID,'numerical_2nd_deriv_4th_order_1D: forward_F',F
ENDIF !!ISNAN(second_derivative_F)
IF(debug) WRITE (9,*) My_MPI_Process_ID,'numerical_2nd_deriv_4th_order_1D: about to return from SUBROUTINE'
IF(debug) WRITE (0,*) My_MPI_Process_ID,'numerical_2nd_deriv_4th_order_1D: about to return from SUBROUTINE'
RETURN
END SUBROUTINE numerical_2nd_deriv_4th_order_1D
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