POST.f90

Go to the documentation of this file.

! (from http://patorjk.com/software/taag/ ANSI shadow)
 
!   ██████╗  ██████╗ ███████╗████████╗
!   ██╔══██╗██╔═══██╗██╔════╝╚══██╔══╝
!   ██████╔╝██║   ██║███████╗   ██║   
!   ██╔═══╝ ██║   ██║╚════██║   ██║   
!   ██║     ╚██████╔╝███████║   ██║   
!   ╚═╝      ╚═════╝ ╚══════╝   ╚═╝   
 
!------------------------------------------------------------------------------!
!>  Peeling Ballooning in 3D: postprocessing
!------------------------------------------------------------------------------!
!>  \author
!!  Toon Weyens,
!!  Contact: weyenst@gmail.com
!------------------------------------------------------------------------------!
!>  \version    2.41
!!  \date       2012-2018
!!  Copyright (C) 2019 Toon Weyens
!!
!!  This program is free software: you can redistribute it and/or modify
!!  it under the terms of the GNU General Public License as published by
!!  the Free Software Foundation, either version 3 of the License, or
!!  (at your option) any later version.
!!
!!  This program 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 General Public License for more details.
!!
!!  You should have received a copy of the License along with this program.
!!  If not, contact weyenst@gmail.com
!------------------------------------------------------------------------------!
!>  \see
!!  References:
!!  \cite weyens2014theory
!!  \cite Weyens2017PB3D
!------------------------------------------------------------------------------!
#define CHCKERR if(ierr.ne.0) then; call sudden_stop(ierr); end if
program post
    use str_utilities, only: i2str
    use num_vars, only: prog_name, prog_style, rank
    use messages
    use hdf5_vars, only: init_hdf5
    use mpi_ops, only: start_mpi, stop_mpi, broadcast_input_opts, &
        &sudden_stop
    use files_ops, only: init_files, parse_args, open_input, open_output, &
        &close_output
    use input_ops, only: read_input_opts
    use driver_post, only: run_driver_post, init_post, stop_post
    use pb3d_ops, only: reconstruct_pb3d_in
    use input_utilities, only: dealloc_in
    use eq_utilities, only: do_eq, eq_info
    use grid_vars, only: grid_type
#if ldebug
    use num_vars, only: ltest
    use test, only: generic_tests
#endif
    
    implicit none
 
    ! local variables
    integer :: ierr                                                             ! error
    
    !------------------------------!
    !   Initialize some routines   !
    !------------------------------!
    ierr = start_mpi()                                                          ! start MPI
    chckerr
    prog_name = 'POST'                                                          ! program name
    prog_style = 2                                                              ! post-processing part
    call print_hello()                                                          ! print message with time, etc
    call init_output()                                                          ! initialize output
    call init_files()                                                           ! initialize file operations
    call init_time()                                                            ! initialize time
    call init_hdf5()                                                            ! initialize HDF5
 
    !--------------------------!
    !   Read the PB3D output   !
    !--------------------------!
    call start_time
    call writo('Initialization')
    call lvl_ud(1)
    if (rank.eq.0) then
        ierr = parse_args()                                                     ! parse argument (options are used in open_input)
        chckerr
        ierr = open_input()                                                     ! open the input files
        chckerr
        ierr = reconstruct_pb3d_in('in')                                        ! reconstruct miscellaneous PB3D output variables
        chckerr
        ierr = read_input_opts()                                                ! read input options file
        chckerr
        ierr = open_output()                                                    ! open output file per alpha group
        chckerr
        call dealloc_in()                                                       ! clean up input from equilibrium codes
    end if
    ierr = broadcast_input_opts()                                               ! broadcast input options to other processors
    chckerr
    call writo('')
    call passed_time
    call writo('')
    call lvl_ud(-1)
    
#if ldebug
    !-------------------!
    !   Do some tests   !
    !-------------------!
    if (ltest) then
        call start_time
        call writo('Generic Tests')
        call lvl_ud(1)
        ierr = generic_tests()
        chckerr
        call writo('')
        call passed_time
        call writo('')
        call lvl_ud(-1)
    end if
#endif
    
    !-----------------!
    !   Main driver   !
    !-----------------!
    call start_time
    call writo('Initializing POST')
    call lvl_ud(1)
    ierr = init_post()
    chckerr
    call stop_time
    call writo('')
    call lvl_ud(-1)
    par: do while(do_eq())
        call start_time
        call writo('Main driver'//trim(eq_info()))
        call lvl_ud(1)
        ierr = run_driver_post()
        chckerr
        call writo('')
        call passed_time
        call writo('')
        call lvl_ud(-1)
    end do par
    call stop_post()
 
    !--------------!
    !   clean up   !
    !--------------!
    call writo('Clean up')
    call lvl_ud(1)
    ierr = stop_mpi()
    chckerr
    call close_output
    call lvl_ud(-1)
    
    call print_goodbye
end program post