Commit 08d92562 authored by Femi Kadri's avatar Femi Kadri

Modifications made

parent 5f96039f
......@@ -46,21 +46,21 @@ SUBROUTINE INIT_GEOMETRY_backwards_step(&
INTEGER :: I,J,K
IF(debug) WRITE (9,*) My_MPI_Process_ID,'INIT_GEOMETRY_backwards_step: entered subroutine'
IF(debug) WRITE (0,*) My_MPI_Process_ID,'INIT_GEOMETRY_backwards_step: entered subroutine'
global_STEP_x_max_real = INT( 0.1 * REAL(global_x_max_real) ) !USER INPUT (CONTROLS STEP WIDTH RELATIVE TO DOMAIN WIDTH)
IF(debug) WRITE (0,*) My_MPI_Process_ID,'INIT_GEOMETRY_backwards_step: entered subroutine'
global_STEP_x_max_real = INT( 0.1 * REAL(global_x_max_real) ) !USER INPUT (CONTROLS STEP WIDTH RELATIVE TO DOMAIN WIDTH)
global_STEP_y_max_real = INT( 0.5 * REAL(global_y_max_real) ) !USER INPUT (CONTROLS STEP HEIGHT RELATIVE TO DOMAIN HEIGHT)
IF(debug) WRITE (0,*) My_MPI_Process_ID,'INIT_GEOMETRY_backwards_step: global_STEP_x_max_real',global_STEP_x_max_real
IF(debug) WRITE (0,*) My_MPI_Process_ID,'INIT_GEOMETRY_backwards_step: global_STEP_y_max_real',global_STEP_y_max_real
IF(debug) WRITE (0,*) My_MPI_Process_ID,'INIT_GEOMETRY_backwards_step: global_STEP_x_max_real',global_STEP_x_max_real
IF(debug) WRITE (0,*) My_MPI_Process_ID,'INIT_GEOMETRY_backwards_step: global_STEP_y_max_real',global_STEP_y_max_real
IF(debug) WRITE (0,*) My_MPI_Process_ID,'INIT_GEOMETRY_backwards_step: local_z_min_ghost',local_z_min_ghost
IF(debug) WRITE (0,*) My_MPI_Process_ID,'INIT_GEOMETRY_backwards_step: local_z_max_ghost',local_z_max_ghost
IF(debug) WRITE (0,*) My_MPI_Process_ID,'INIT_GEOMETRY_backwards_step: local_y_min_ghost',local_y_min_ghost
IF(debug) WRITE (0,*) My_MPI_Process_ID,'INIT_GEOMETRY_backwards_step: local_y_max_ghost',local_y_max_ghost
IF(debug) WRITE (0,*) My_MPI_Process_ID,'INIT_GEOMETRY_backwards_step: local_x_min_ghost',local_x_min_ghost
IF(debug) WRITE (0,*) My_MPI_Process_ID,'INIT_GEOMETRY_backwards_step: local_x_max_ghost',local_x_max_ghost
IF(debug) WRITE (0,*) My_MPI_Process_ID,'INIT_GEOMETRY_backwards_step: local_z_min_ghost',local_z_min_ghost
IF(debug) WRITE (0,*) My_MPI_Process_ID,'INIT_GEOMETRY_backwards_step: local_z_max_ghost',local_z_max_ghost
IF(debug) WRITE (0,*) My_MPI_Process_ID,'INIT_GEOMETRY_backwards_step: local_y_min_ghost',local_y_min_ghost
IF(debug) WRITE (0,*) My_MPI_Process_ID,'INIT_GEOMETRY_backwards_step: local_y_max_ghost',local_y_max_ghost
IF(debug) WRITE (0,*) My_MPI_Process_ID,'INIT_GEOMETRY_backwards_step: local_x_min_ghost',local_x_min_ghost
IF(debug) WRITE (0,*) My_MPI_Process_ID,'INIT_GEOMETRY_backwards_step: local_x_max_ghost',local_x_max_ghost
IF(debug) WRITE (0,*) My_MPI_Process_ID,'INIT_GEOMETRY_backwards_step: global_z_min_real',global_z_min_real
IF(debug) WRITE (0,*) My_MPI_Process_ID,'INIT_GEOMETRY_backwards_step: global_z_max_real',global_z_max_real
......@@ -80,13 +80,13 @@ SUBROUTINE INIT_GEOMETRY_backwards_step(&
local_walls(I,J,K) = .TRUE.
ENDIF !!(J.LE.global_STEP_y_max_real).AND.(I.LE.global_STEP_x_max_real)
IF( (J.EQ.global_y_min_real).OR.(J.EQ.global_y_max_real) ) THEN
local_walls(I,J,K) = .TRUE.
ENDIF !!J.EQ.global_y_min_real).AND.(J.EQ.global_y_max_real)
ENDIF !!(J.LE.global_STEP_y_max_real).AND.(I.LE.global_STEP_x_max_real)
IF( (J.EQ.global_y_min_real).OR.(J.EQ.global_y_max_real) ) THEN
local_walls(I,J,K) = .TRUE.
ENDIF !!J.EQ.global_y_min_real).AND.(J.EQ.global_y_max_real)
END DO !k = 1, nz
......
......@@ -17,27 +17,27 @@ SUBROUTINE INIT_GEOMETRY_fiber_layers_MAIN(&
& local_y_min_ghost,&
& local_y_max_ghost,&
& local_z_min_ghost,&
& local_z_max_ghost,&
& local_x_min_real,&
& local_x_max_real,&
& local_y_min_real,&
& local_y_max_real,&
& local_z_min_real,&
& local_z_max_real,&
& divide_x,&
& divide_y,&
& divide_z,&
& local_x_nodes_ghost,&
& local_y_nodes_ghost,&
& local_z_nodes_ghost,&
& local_x_nodes_real,&
& local_y_nodes_real,&
& local_z_nodes_real,&
& IN,&
& OUT,&
& WEST,&
& EAST,&
& UP,&
& local_z_max_ghost,&
& local_x_min_real,&
& local_x_max_real,&
& local_y_min_real,&
& local_y_max_real,&
& local_z_min_real,&
& local_z_max_real,&
& divide_x,&
& divide_y,&
& divide_z,&
& local_x_nodes_ghost,&
& local_y_nodes_ghost,&
& local_z_nodes_ghost,&
& local_x_nodes_real,&
& local_y_nodes_real,&
& local_z_nodes_real,&
& IN,&
& OUT,&
& WEST,&
& EAST,&
& UP,&
& DOWN,&
& local_walls)
......@@ -57,23 +57,23 @@ USE debug_module
INTEGER,INTENT(IN) :: local_x_min_ghost,local_x_max_ghost
INTEGER,INTENT(IN) :: local_y_min_ghost,local_y_max_ghost
INTEGER,INTENT(IN) :: local_z_min_ghost,local_z_max_ghost
INTEGER,INTENT(IN) :: local_x_min_real,local_x_max_real
INTEGER,INTENT(IN) :: local_y_min_real,local_y_max_real
INTEGER,INTENT(IN) :: local_z_min_real,local_z_max_real
INTEGER,INTENT(IN) :: IN,OUT,EAST,WEST,UP,DOWN
INTEGER,INTENT(IN) :: local_x_nodes_ghost
INTEGER,INTENT(IN) :: local_y_nodes_ghost
INTEGER,INTENT(IN) :: local_z_nodes_ghost
INTEGER,INTENT(IN) :: local_x_nodes_real
INTEGER,INTENT(IN) :: local_y_nodes_real
INTEGER,INTENT(IN) :: local_z_min_ghost,local_z_max_ghost
INTEGER,INTENT(IN) :: local_x_min_real,local_x_max_real
INTEGER,INTENT(IN) :: local_y_min_real,local_y_max_real
INTEGER,INTENT(IN) :: local_z_min_real,local_z_max_real
INTEGER,INTENT(IN) :: IN,OUT,EAST,WEST,UP,DOWN
INTEGER,INTENT(IN) :: local_x_nodes_ghost
INTEGER,INTENT(IN) :: local_y_nodes_ghost
INTEGER,INTENT(IN) :: local_z_nodes_ghost
INTEGER,INTENT(IN) :: local_x_nodes_real
INTEGER,INTENT(IN) :: local_y_nodes_real
INTEGER,INTENT(IN) :: local_z_nodes_real
INTEGER,INTENT(IN) :: ROOT,cartesian_communicator
INTEGER,INTENT(IN) :: ROOT,cartesian_communicator
LOGICAL*4,INTENT(IN) :: divide_x,divide_y,divide_z
LOGICAL*4,PARAMETER :: debug = GLOBAL_DEBUG
......@@ -529,8 +529,8 @@ USE debug_module
END DO !j=1,...
END DO !i=1,...
END DO !i=1,...
!WRITE(1234,*)layer,1.0d0-REAL(count_grid)/REAL(layer_length_in_lattice_units*layer_width_in_lattice_units*layer_thickness_in_lattice_units)
......@@ -539,21 +539,21 @@ USE debug_module
END DO ! total_number_of_layers
IF(debug) WRITE (0,*) My_MPI_Process_ID,&
& 'INIT_GEOMETRY_fiber_layers_MAIN: about to CALL INIT_GEOMETRY_fiber_layers_MAIN: local_walls'
IF(debug) WRITE (9,*) My_MPI_Process_ID,&
& 'INIT_GEOMETRY_fiber_layers_MAIN: about to CALL INIT_GEOMETRY_fiber_layers_MAIN: local_walls'
!UPDATE GHOST BOUNDARIES OF local_walls
CALL update_ghost_boundaries_logical(local_x_min_ghost,local_x_max_ghost,local_y_min_ghost,local_y_max_ghost,&
& local_z_min_ghost,local_z_max_ghost,local_x_min_real,local_x_max_real,local_y_min_real,local_y_max_real,local_z_min_real,&
& local_z_max_real,My_MPI_Process_ID,IN,OUT,WEST,EAST,UP,DOWN,cartesian_communicator,local_x_nodes_ghost,&
& local_y_nodes_ghost,local_z_nodes_ghost,local_x_nodes_real,local_y_nodes_real,local_z_nodes_real,&
& divide_x,divide_y,divide_z,local_walls)
IF(debug) WRITE (0,*) My_MPI_Process_ID,&
& 'INIT_GEOMETRY_fiber_layers_MAIN: finished CALL INIT_GEOMETRY_fiber_layers_MAIN: local_walls'
IF(debug) WRITE (9,*) My_MPI_Process_ID,&
IF(debug) WRITE (0,*) My_MPI_Process_ID,&
& 'INIT_GEOMETRY_fiber_layers_MAIN: about to CALL INIT_GEOMETRY_fiber_layers_MAIN: local_walls'
IF(debug) WRITE (9,*) My_MPI_Process_ID,&
& 'INIT_GEOMETRY_fiber_layers_MAIN: about to CALL INIT_GEOMETRY_fiber_layers_MAIN: local_walls'
!UPDATE GHOST BOUNDARIES OF local_walls
CALL update_ghost_boundaries_logical(local_x_min_ghost,local_x_max_ghost,local_y_min_ghost,local_y_max_ghost,&
& local_z_min_ghost,local_z_max_ghost,local_x_min_real,local_x_max_real,local_y_min_real,local_y_max_real,local_z_min_real,&
& local_z_max_real,My_MPI_Process_ID,IN,OUT,WEST,EAST,UP,DOWN,cartesian_communicator,local_x_nodes_ghost,&
& local_y_nodes_ghost,local_z_nodes_ghost,local_x_nodes_real,local_y_nodes_real,local_z_nodes_real,&
& divide_x,divide_y,divide_z,local_walls)
IF(debug) WRITE (0,*) My_MPI_Process_ID,&
& 'INIT_GEOMETRY_fiber_layers_MAIN: finished CALL INIT_GEOMETRY_fiber_layers_MAIN: local_walls'
IF(debug) WRITE (9,*) My_MPI_Process_ID,&
& 'INIT_GEOMETRY_fiber_layers_MAIN: finished CALL INIT_GEOMETRY_fiber_layers_MAIN: local_walls'
......
SUBROUTINE INIT_GEOMETRY_narrowing_slit(&
& global_x_min_real,&
& global_x_max_real,&
& global_x_min_ghost,&
& global_x_max_ghost,&
& global_y_min_real,&
& global_y_max_real,&
& global_y_min_ghost,&
& global_y_max_ghost,&
& global_z_min_real,&
& global_z_max_real,&
& My_MPI_Process_ID,&
& local_x_min_ghost,&
& local_x_max_ghost,&
& local_y_min_ghost,&
& local_y_max_ghost,&
& local_z_min_ghost,&
& local_z_max_ghost,&
& local_walls)
!THIS SUBROUTINE CREATES A NARROWING SLIT GEOMETRY
!(TOP WALL AND BOTTOM WALL COMING TOGETHER TOWARDS THE END OF SLIT)
!WARNING, TAKE CARE WHEN DOING THIS WITH PERIODIC BOUNDARY CONDITIONS
USE debug_module
IMPLICIT NONE
!-----------INTENT(IN)-----------------------------------------------------------------------------------------------
INTEGER,INTENT(IN) :: global_x_min_real,global_x_max_real,global_x_min_ghost,global_x_max_ghost
INTEGER,INTENT(IN) :: global_y_min_real,global_y_max_real,global_y_min_ghost,global_y_max_ghost
INTEGER,INTENT(IN) :: global_z_min_real,global_z_max_real,My_MPI_Process_ID
INTEGER,INTENT(IN) :: local_x_min_ghost,local_x_max_ghost
INTEGER,INTENT(IN) :: local_y_min_ghost,local_y_max_ghost
INTEGER,INTENT(IN) :: local_z_min_ghost,local_z_max_ghost
LOGICAL*4,PARAMETER :: debug = GLOBAL_DEBUG
!-----------INTENT(INOUT)-----------------------------------------------------------------------------------------------
LOGICAL*4,INTENT(INOUT),DIMENSION(local_x_min_ghost:local_x_max_ghost, &
& local_y_min_ghost:local_y_max_ghost,local_z_min_ghost:local_z_max_ghost) &
& :: local_walls
!-----------NO INTENT-----------------------------------------------------------------------------------------------
INTEGER :: I,J,K
REAL :: L, L_plus_one
REAL :: h_0,h_L,height,h_L_per_L
IF(debug) WRITE (9,*) My_MPI_Process_ID,'INIT_GEOMETRY_narrowing_slit: entered subroutine'
IF(debug) WRITE (0,*) My_MPI_Process_ID,'INIT_GEOMETRY_narrowing_slit: entered subroutine'
IF(debug) WRITE (0,*) My_MPI_Process_ID,'entered SUBROUTINE INIT_GEOMETRY_narrowing_slit'
! IF(debug) WRITE (0,*) My_MPI_Process_ID,'INIT_GEOMETRY_narrowing_slit: global_z_min_ghost',global_z_min_ghost
! IF(debug) WRITE (0,*) My_MPI_Process_ID,'INIT_GEOMETRY_narrowing_slit: global_z_max_ghost',global_z_max_ghost
IF(debug) WRITE (0,*) My_MPI_Process_ID,'INIT_GEOMETRY_narrowing_slit: global_y_min_ghost',global_y_min_ghost
IF(debug) WRITE (0,*) My_MPI_Process_ID,'INIT_GEOMETRY_narrowing_slit: global_y_max_ghost',global_y_max_ghost
IF(debug) WRITE (0,*) My_MPI_Process_ID,'INIT_GEOMETRY_narrowing_slit: global_x_min_ghost',global_x_min_ghost
IF(debug) WRITE (0,*) My_MPI_Process_ID,'INIT_GEOMETRY_narrowing_slit: global_x_max_ghost',global_x_max_ghost
IF(debug) WRITE (0,*) My_MPI_Process_ID,'INIT_GEOMETRY_narrowing_slit: global_z_min_real',global_z_min_real
IF(debug) WRITE (0,*) My_MPI_Process_ID,'INIT_GEOMETRY_narrowing_slit: global_z_max_real',global_z_max_real
IF(debug) WRITE (0,*) My_MPI_Process_ID,'INIT_GEOMETRY_narrowing_slit: global_y_min_real',global_y_min_real
IF(debug) WRITE (0,*) My_MPI_Process_ID,'INIT_GEOMETRY_narrowing_slit: global_y_max_real',global_y_max_real
IF(debug) WRITE (0,*) My_MPI_Process_ID,'INIT_GEOMETRY_narrowing_slit: global_x_min_real',global_x_min_real
IF(debug) WRITE (0,*) My_MPI_Process_ID,'INIT_GEOMETRY_narrowing_slit: global_x_max_real',global_x_max_real
L = REAL(global_x_max_real - global_x_min_real) / 2.0d0
L_plus_one = L+1
h_0 = REAL(global_y_max_real - global_y_min_real)
h_L = h_0 / 2.0d0
h_L_per_L = h_L/L
DO K=local_z_min_ghost,local_z_max_ghost
DO J=local_y_min_ghost,local_y_max_ghost
DO I=local_x_min_ghost,local_x_max_ghost
IF ( (REAL(I).LT.(L_plus_one)).OR.&
& (REAL(I).EQ.(L_plus_one)) ) THEN
height = (h_L_per_L) * REAL(1-I) + h_0
IF(J.GT.height.OR.&
& J.EQ.global_y_min_real.OR.&
& K.EQ.global_z_min_real.OR.&
& K.EQ.global_z_max_real) THEN
local_walls(I,J,K) = .TRUE.
ENDIF !!J.EQ.global_y_min_real.OR.J.EQ.global_y_max_real
ELSE
height = (h_L_per_L) * REAL(I-1)
IF(J.GT.height.OR.&
& J.EQ.global_y_min_real.OR.&
& K.EQ.global_z_min_real.OR.&
& K.EQ.global_z_max_real) THEN
local_walls(I,J,K) = .TRUE.
ENDIF !!J.EQ.global_y_min_real.OR.J.EQ.global_y_max_real
END IF
END DO !k = 1, nz
END DO !j = 1, ny
END DO !i = 1, nx
IF(debug) WRITE (9,*) My_MPI_Process_ID,'INIT_GEOMETRY_narrowing_slit: about to return from subroutine'
IF(debug) WRITE (0,*) My_MPI_Process_ID,'INIT_GEOMETRY_narrowing_slit: about to return from subroutine'
RETURN
SUBROUTINE INIT_GEOMETRY_narrowing_slit(&
& global_x_min_real,&
& global_x_max_real,&
& global_x_min_ghost,&
& global_x_max_ghost,&
& global_y_min_real,&
& global_y_max_real,&
& global_y_min_ghost,&
& global_y_max_ghost,&
& global_z_min_real,&
& global_z_max_real,&
& My_MPI_Process_ID,&
& local_x_min_ghost,&
& local_x_max_ghost,&
& local_y_min_ghost,&
& local_y_max_ghost,&
& local_z_min_ghost,&
& local_z_max_ghost,&
& local_walls)
!THIS SUBROUTINE CREATES A NARROWING SLIT GEOMETRY
!(TOP WALL AND BOTTOM WALL COMING TOGETHER TOWARDS THE END OF SLIT)
!WARNING, TAKE CARE WHEN DOING THIS WITH PERIODIC BOUNDARY CONDITIONS
USE debug_module
IMPLICIT NONE
!-----------INTENT(IN)-----------------------------------------------------------------------------------------------
INTEGER,INTENT(IN) :: global_x_min_real,global_x_max_real,global_x_min_ghost,global_x_max_ghost
INTEGER,INTENT(IN) :: global_y_min_real,global_y_max_real,global_y_min_ghost,global_y_max_ghost
INTEGER,INTENT(IN) :: global_z_min_real,global_z_max_real,My_MPI_Process_ID
INTEGER,INTENT(IN) :: local_x_min_ghost,local_x_max_ghost
INTEGER,INTENT(IN) :: local_y_min_ghost,local_y_max_ghost
INTEGER,INTENT(IN) :: local_z_min_ghost,local_z_max_ghost
LOGICAL*4,PARAMETER :: debug = GLOBAL_DEBUG
!-----------INTENT(INOUT)-----------------------------------------------------------------------------------------------
LOGICAL*4,INTENT(INOUT),DIMENSION(local_x_min_ghost:local_x_max_ghost, &
& local_y_min_ghost:local_y_max_ghost,local_z_min_ghost:local_z_max_ghost) &
& :: local_walls
!-----------NO INTENT-----------------------------------------------------------------------------------------------
INTEGER :: I,J,K
REAL :: L, L_plus_one
REAL :: h_0,h_L,height,h_L_per_L
IF(debug) WRITE (9,*) My_MPI_Process_ID,'INIT_GEOMETRY_narrowing_slit: entered subroutine'
IF(debug) WRITE (0,*) My_MPI_Process_ID,'INIT_GEOMETRY_narrowing_slit: entered subroutine'
IF(debug) WRITE (0,*) My_MPI_Process_ID,'entered SUBROUTINE INIT_GEOMETRY_narrowing_slit'
! IF(debug) WRITE (0,*) My_MPI_Process_ID,'INIT_GEOMETRY_narrowing_slit: global_z_min_ghost',global_z_min_ghost
! IF(debug) WRITE (0,*) My_MPI_Process_ID,'INIT_GEOMETRY_narrowing_slit: global_z_max_ghost',global_z_max_ghost
IF(debug) WRITE (0,*) My_MPI_Process_ID,'INIT_GEOMETRY_narrowing_slit: global_y_min_ghost',global_y_min_ghost
IF(debug) WRITE (0,*) My_MPI_Process_ID,'INIT_GEOMETRY_narrowing_slit: global_y_max_ghost',global_y_max_ghost
IF(debug) WRITE (0,*) My_MPI_Process_ID,'INIT_GEOMETRY_narrowing_slit: global_x_min_ghost',global_x_min_ghost
IF(debug) WRITE (0,*) My_MPI_Process_ID,'INIT_GEOMETRY_narrowing_slit: global_x_max_ghost',global_x_max_ghost
IF(debug) WRITE (0,*) My_MPI_Process_ID,'INIT_GEOMETRY_narrowing_slit: global_z_min_real',global_z_min_real
IF(debug) WRITE (0,*) My_MPI_Process_ID,'INIT_GEOMETRY_narrowing_slit: global_z_max_real',global_z_max_real
IF(debug) WRITE (0,*) My_MPI_Process_ID,'INIT_GEOMETRY_narrowing_slit: global_y_min_real',global_y_min_real
IF(debug) WRITE (0,*) My_MPI_Process_ID,'INIT_GEOMETRY_narrowing_slit: global_y_max_real',global_y_max_real
IF(debug) WRITE (0,*) My_MPI_Process_ID,'INIT_GEOMETRY_narrowing_slit: global_x_min_real',global_x_min_real
IF(debug) WRITE (0,*) My_MPI_Process_ID,'INIT_GEOMETRY_narrowing_slit: global_x_max_real',global_x_max_real
L = REAL(global_x_max_real - global_x_min_real) / 2.0d0
L_plus_one = L+1
h_0 = REAL(global_y_max_real - global_y_min_real)
h_L = h_0 / 2.0d0
h_L_per_L = h_L/L
DO K=local_z_min_ghost,local_z_max_ghost
DO J=local_y_min_ghost,local_y_max_ghost
DO I=local_x_min_ghost,local_x_max_ghost
IF ( (REAL(I).LT.(L_plus_one)).OR.&
& (REAL(I).EQ.(L_plus_one)) ) THEN
height = (h_L_per_L) * REAL(1-I) + h_0
IF(J.GT.height.OR.&
& J.EQ.global_y_min_real.OR.&
& K.EQ.global_z_min_real.OR.&
& K.EQ.global_z_max_real) THEN
local_walls(I,J,K) = .TRUE.
ENDIF !!J.EQ.global_y_min_real.OR.J.EQ.global_y_max_real
ELSE
height = (h_L_per_L) * REAL(I-1)
IF(J.GT.height.OR.&
& J.EQ.global_y_min_real.OR.&
& K.EQ.global_z_min_real.OR.&
& K.EQ.global_z_max_real) THEN
local_walls(I,J,K) = .TRUE.
ENDIF !!J.EQ.global_y_min_real.OR.J.EQ.global_y_max_real
END IF
END DO !k = 1, nz
END DO !j = 1, ny
END DO !i = 1, nx
IF(debug) WRITE (9,*) My_MPI_Process_ID,'INIT_GEOMETRY_narrowing_slit: about to return from subroutine'
IF(debug) WRITE (0,*) My_MPI_Process_ID,'INIT_GEOMETRY_narrowing_slit: about to return from subroutine'
RETURN
END SUBROUTINE INIT_GEOMETRY_narrowing_slit
\ No newline at end of file
rm -irf debug/* backup/* results_LBM/* results_LST/*
mkdir results_LBM/stresses
rm -irf LBM_Time_Steps_Prev_Completed.txt LBM_previously_converged.txt wrote_HARD_resume_files.txt *LBM*out.txt *LBM*err.txt LBM_currently_converged.dat
\ No newline at end of file
This diff is collapsed.
head -2 LBM_Time_Steps_Prev_Completed.txt 6000 clear all file *.dat
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.