params_brag_init_select_m.f90 Source File


Contents


Source Code

module params_brag_init_select_m
    !! Parameters related with selection of initial state
    use error_handling_grillix_m, only: handle_error
    use status_codes_grillix_m, only : GRILLIX_ERR_NAMELIST
    use screen_io_m, only :  get_stdout
    use technical_constants_m, only : PATHLEN_MAX
    implicit none
         
    character(len=64), protected :: brag_iselect = 'NONE'
    !! Type of initial state (BLOB, CONTINUE)
    character(len=PATHLEN_MAX), protected :: brag_init_path = 'params_braginskii.in'
    !! Path to where parameters for initial state are possibly read from
    
    public :: read_params_brag_init_select
    public :: write_params_brag_init_select
    
    namelist / brag_init_select / &
        brag_iselect, &
        brag_init_path
    private brag_init_select
    
contains

    subroutine read_params_brag_init_select(filename)
        !! Reads parameters related with brag_init_select
        character(len=*), intent(in) :: filename
        !! Filename, to read from 
        
        integer :: io_unit, io_error
        character(len=256) :: io_errmsg 
         
        open(newunit=io_unit, file=filename, status='old', action='read', &
            iostat=io_error, iomsg=io_errmsg)
        if (io_error /= 0) then
            call handle_error(io_errmsg, GRILLIX_ERR_NAMELIST, &
                              __LINE__, __FILE__)
        endif
  
        read(io_unit, nml=brag_init_select, iostat=io_error, iomsg=io_errmsg)
        if (io_error /= 0) then
            call handle_error(io_errmsg, GRILLIX_ERR_NAMELIST, &
                              __LINE__, __FILE__)
        endif

        close(io_unit)

    end subroutine

    subroutine write_params_brag_init_select(filename)
        !! Reads parameters related with brag_init_select
        character(len=*), intent(in), optional :: filename
        !! If present, filename where params are written to, 
        !! if not present, writes to screen 
        
        integer :: io_unit, io_error
        character(len=256) :: io_errmsg 
         
        if (present(filename)) then
            open(newunit=io_unit, file=filename, status='unknown', &
                access='append', action='write', &
                iostat=io_error, iomsg=io_errmsg )
            if (io_error /= 0) then
                call handle_error(io_errmsg, GRILLIX_ERR_NAMELIST, &
                                  __LINE__, __FILE__)
            endif
        else
            io_unit = get_stdout()
        endif

        write(io_unit, nml=brag_init_select, iostat=io_error, iomsg=io_errmsg)
        if (io_error /= 0) then
            call handle_error(io_errmsg, GRILLIX_ERR_NAMELIST, &
                              __LINE__, __FILE__)
        endif

        if (present(filename)) then
            close(io_unit)
        endif

    end subroutine    
    

end module