Homec4science
Diffusion FENNECS lite (3Dmpimode) Browse

fennecslite/src/main.f90
d24dd0989d0e3Dmpimode

Last Change
Raw File

main.f90
View Options

PROGRAM main
!
! Skeleton for a time dependent program
! Note: Even in this sequential version, MPI is required because of FUTILS
! (more specifcally because of the HASTABLE module)
!
USE basic
USE mpi
USE bsplines
USE mumps_bsplines
USE futils
USE omp_lib
IMPLICIT NONE
INTEGER:: required, provided
required=MPI_THREAD_SERIALIZED
CALL mpi_init_thread(required,provided,ierr)
IF(provided .ne. required) CALL MPI_abort(MPI_COMM_WORLD,-1,ierr)
WRITE(*,*)" MPI required provided", required,provided
IF(ierr .ne. 0) CALL MPI_abort(MPI_COMM_WORLD,-1,ierr)
CALL MPI_COMM_RANK(MPI_COMM_WORLD, mpirank, ierr)
IF(ierr .ne. 0) CALL MPI_abort(MPI_COMM_WORLD,-1,ierr)
CALL MPI_COMM_SIZE(MPI_COMM_WORLD, mpisize, ierr)
IF(ierr .ne. 0) CALL MPI_abort(MPI_COMM_WORLD,-1,ierr)
!--------------------------------------------------------------------------------
! 1. Prologue
CALL timera(0, 'Prologue')
CALL daytim('Start at')
! Define data specific to run
CALL basic_data ! definition of global variables and input paramaters loading
step=0
! Compute auxilliary values
CALL auxval ! time independent values
! Initial conditions
IF( .NOT. nlres ) THEN
CALL inital ! plasma initialisation
ELSE
CALL resume ! loads restart.h5 file
END IF
! Start or restart the run
CALL start
! Initial diagnostocs
CALL diagnose(0)
CALL timera(1, 'Prologue')
!--------------------------------------------------------------------------------
! 2. Time stepping
CALL timera(0, 'Main loop')
DO
step = step+1
cstep = cstep+1
time = time+dt*tnorm
CALL tesend
CALL stepon
IF(modulo(step,itrestart) .eq. 0 .and. mpirank .eq. 0) CALL chkrst(1)
CALL MPI_BARRIER(MPI_COMM_WORLD, ierr)
IF( nlend ) EXIT
END DO
CALL timera(1, 'Main loop')
!--------------------------------------------------------------------------------
! 9. Epilogue
CALL timera(0, 'Epilogue')
CALL diagnose(-1)
CALL endrun
IF(mpirank .eq. 0) THEN
CALL timera(1, 'Epilogue')
CALL timera(9, '')
CALL timera(-1, '')
CALL daytim('Done at ')
END IF
CALL mpi_finalize(ierr)
END PROGRAM main
c4science · Help