Homec4science
Diffusion FENNECS: Finite Element Non-Neutral Electron Cloud Simulator (master) Browse

fennecs/src/main.f90
dcbebf0c577cmaster

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
!
IF( .NOT. nlres ) THEN
CALL newrun !not implemented yet
ELSE
CALL restart !not implemented yet
END IF
!
! 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 !not implemented yet
!
! 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