esputr.f90

Go to the documentation of this file.

!    Copyright (c) 2016  Forschungszentrum Juelich GmbH
!                        Markus Brenneis, Vladislav Kotov 
!
!    This file is part of ESPUTR.
!
!    ESPUTR 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.
!
!    ESPUTR 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 GNU General Public License
!    along with ESPUTR.  If not, see <http://www.gnu.org/licenses/>.
!
    

module esputr

!***********************************************************************************
!   DEFINITION OF PARAMETRS AND GLOBAL VARIABLES
!***********************************************************************************
    implicit none
    public numberoflines
    integer, public, save :: esputr_unit = 6
    integer, public, parameter :: esputr_dp=8
    real(ESPUTR_DP), public, parameter :: esputr_pi = 4*atan(1.)
    real(ESPUTR_DP), public, parameter :: esputr_pi2 = esputr_pi/2
    integer, public, parameter :: esputr_max_line_length = 256
    private
    character(*), public, parameter :: esputr_comment_start = '#'
    
    contains
    
!***********************************************************************************
!   UTILITY FUNCTIONS
!***********************************************************************************    
    integer function numberoflines(fileName, ignoreComments,err)

        character(*), intent(in) :: fileName
        logical, intent(in), optional :: ignoreComments
        integer, intent(out) :: err
        integer, parameter :: UNIT = 1337
        character(1) :: line
        character(100) :: msg
        intrinsic trim        

        err = 0
        open(unit, iostat=err, iomsg=msg, file=filename, status='old', action='read')
        
        numberoflines = 0
        if(err == 0) then
            do
                read(unit, *, iostat=err,end=100) line
                if(err /= 0) then
                    write(esputr_unit,*), "ERROR in numberOfLines while reading "//trim(filename)
                    close(unit)
                    err = 300
                    return
                end if
                if(line /= esputr_comment_start) then
                    numberoflines = numberoflines + 1
                else if(present(ignorecomments)) then
                    if(.not.ignorecomments) numberoflines = numberoflines + 1
                else if(.not.present(ignorecomments)) then
                    numberoflines = numberoflines + 1
                end if
            end do
        else
            write(esputr_unit,*), "ERROR in numberOfLines: cannot open file "//trim(filename)
            err = 300
            return
        end if
100     err = 0
        close(unit)
    end function numberoflines
   
end module esputr