Finalize.F90

Go to the documentation of this file.

!/*****************************************************************************/
! *
! *  Elmer, A Finite Element Software for Multiphysical Problems
! *
! *  Copyright 1st April 1995 - , CSC - IT Center for Science Ltd., Finland
! * 
! *  This library is free software; you can redistribute it and/or
! *  modify it under the terms of the GNU Lesser General Public
! *  License as published by the Free Software Foundation; either
! *  version 2.1 of the License, or (at your option) any later version.
! *
! *  This library is distributed in the hope that it will be useful,
! *  but WITHOUT ANY WARRANTY; without even the implied warranty of
! *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
! *  Lesser General Public License for more details.
! * 
! *  You should have received a copy of the GNU Lesser General Public
! *  License along with this library (in file ../LGPL-2.1); if not, write 
! *  to the Free Software Foundation, Inc., 51 Franklin Street, 
! *  Fifth Floor, Boston, MA  02110-1301  USA
! *
! *****************************************************************************/
!
!/******************************************************************************
! *
! *  ELMER/FEM Solver finalize program
! *
! ******************************************************************************
! *
! *  Authors: Juha Ruokolainen
! *  Email:   Juha.Ruokolainen@csc.fi
! *  Web:     http://www.csc.fi/elmer
! *  Address: CSC - IT Center for Science Ltd.
! *           Keilaranta 14
! *           02101 Espoo, Finland 
! *
! *  Original Date: 02 Jun 1997
! *
! *****************************************************************************/






!------------------------------------------------------------------------------
!------------------------------------------------------------------------------
   SUBROUTINE elmerfinalize(runs)
!------------------------------------------------------------------------------

     USE generalmodule

!------------------------------------------------------------------------------
     IMPLICIT NONE
!------------------------------------------------------------------------------

!     INTEGER :: Initialize

!------------------------------------------------------------------------------
!    Local variables
!------------------------------------------------------------------------------

!     INTEGER :: i,j,k,n,l,t,k1,k2,iter,Ndeg,istat,nproc,tlen,nthreads
!     CHARACTER(LEN=MAX_STRING_LEN) :: threads
!
!     REAL(KIND=dp) :: s,dt,dtfunc
!     REAL(KIND=dP), POINTER :: WorkA(:,:,:) => NULL()
!     REAL(KIND=dp), POINTER, SAVE :: sTime(:), sStep(:), sInterval(:), sSize(:), &
!           steadyIt(:),nonlinIt(:),sPrevSizes(:,:),sPeriodic(:)
!
!     TYPE(Element_t),POINTER :: CurrentElement
!
!     LOGICAL :: GotIt,Transient,Scanning,LastSaved
!
!     INTEGER :: TimeIntervals,interval,timestep, &
!       TotalTimesteps,SavedSteps,CoupledMaxIter,CoupledMinIter
!
!     INTEGER, POINTER, SAVE :: Timesteps(:),OutputIntervals(:), ActiveSolvers(:)
!     REAL(KIND=dp), POINTER, SAVE :: TimestepSizes(:,:)
!
!     INTEGER(KIND=AddrInt) :: ControlProcedure
!
!     LOGICAL :: InitDirichlet, ExecThis
!
!     TYPE(ElementType_t),POINTER :: elmt
!
!     TYPE(ParEnv_t), POINTER :: ParallelEnv
!
!     CHARACTER(LEN=MAX_NAME_LEN) :: ModelName, eq, ExecCommand
!     CHARACTER(LEN=MAX_STRING_LEN) :: OutputFile, PostFile, RestartFile, &
!                OutputName=' ',PostName=' ', When, OptionString
!
!     TYPE(Variable_t), POINTER :: Var
!     TYPE(Mesh_t), POINTER :: Mesh
!     TYPE(Solver_t), POINTER :: Solver
!
!     REAL(KIND=dp) :: RealTime,CPUTime,CT0,RT0,tt
!
!     LOGICAL :: FirstLoad = .TRUE., FirstTime=.TRUE., Found
!     LOGICAL :: Silent, Version, GotModelName
!
!     INTEGER :: iargc, NoArgs
!
!     INTEGER :: ExtrudeLevels
!     TYPE(Mesh_t), POINTER :: ExtrudedMesh
!
!     INTEGER :: omp_get_max_threads
!
!#ifdef HAVE_TRILINOS
!INTERFACE
!      SUBROUTINE TrilinosCleanup() BIND(C,name='TrilinosCleanup')
!      IMPLICIT NONE
!      END SUBROUTINE TrilinosCleanup
!END INTERFACE
!#endif
     INTEGER :: runs

     DO i=1,currentmodel % NumberOfSolvers
        solver => currentmodel % Solvers(i)
        IF ( solver % PROCEDURE == 0 ) cycle
        when = listgetstring( solver % Values, 'Exec Solver', gotit )
        IF ( gotit ) THEN
           IF ( when == 'after simulation' .OR. when == 'after all' ) THEN
              CALL solveractivate( currentmodel,solver,dt,transient )
              IF (ASSOCIATED(solver % Variable % Values) ) lastsaved = .false.
           END IF
        ELSE
           IF ( solver % SolverExecWhen == solver_exec_after_all ) THEN
              CALL solveractivate( currentmodel,solver,dt,transient )
              IF (ASSOCIATED(solver % Variable % Values) ) lastsaved = .false.
           END IF
        END IF
     END DO

     IF( .NOT. lastsaved ) THEN
        DO i=1,currentmodel % NumberOfSolvers
           solver => currentmodel % Solvers(i)
           IF ( solver % PROCEDURE == 0 ) cycle
           execthis = ( solver % SolverExecWhen == solver_exec_ahead_save)
           when = listgetstring( solver % Values, 'Exec Solver', gotit )
           IF ( gotit ) execthis = ( when == 'before saving') 
           IF( execthis ) CALL solveractivate( currentmodel,solver,dt,transient )
        END DO
     END IF


!------------------------------------------------------------------------------
!    Always save the last step to output
!------------------------------------------------------------------------------
       IF ( .NOT.lastsaved ) THEN
         DO i=1,currentmodel % NumberOfSolvers
            solver => currentmodel % Solvers(i)
            IF ( solver % PROCEDURE == 0 ) cycle
            execthis = ( solver % SolverExecWhen == solver_exec_ahead_save)
            when = listgetstring( solver % Values, 'Exec Solver', gotit )
            IF ( gotit ) execthis = ( when == 'before saving') 
            IF( execthis ) CALL solveractivate( currentmodel,solver,dt,transient )
         END DO

         CALL savetopost(0)
         CALL savecurrent(timestep)

         DO i=1,currentmodel % NumberOfSolvers
            solver => currentmodel % Solvers(i)
            IF ( solver % PROCEDURE == 0 ) cycle
            execthis = ( solver % SolverExecWhen == solver_exec_after_save)
            when = listgetstring( solver % Values, 'Exec Solver', gotit )
            IF ( gotit ) execthis = ( when == 'after saving') 
            IF( execthis ) CALL solveractivate( currentmodel,solver,dt,transient )
         END DO
       END IF

!       IF ( Initialize >= 2 ) EXIT
    

     
!------------------------------------------------------------------------------
!    THIS IS THE END (...,at last, the end, my friend,...)
!------------------------------------------------------------------------------
     IF ( initialize /= 1 ) CALL info( 'MyElmerSolver', '*** Elmer Solver: ALL DONE ***',level=3 )

     IF ( initialize <=0 ) CALL freemodel(currentmodel)

#ifdef HAVE_TRILINOS
  CALL trilinoscleanup()
#endif

     IF ( parenv % PEs>1 )  CALL parallelfinalize()
     CALL info('MyElmerSolver','The end',level=3)
 
     runs = 1

     RETURN
     

!------------------------------------------------------------------------------
  END SUBROUTINE elmerfinalize
!------------------------------------------------------------------------------