Toolbox.f90

MODULE Toolbox

  IMPLICIT NONE

CONTAINS

  SUBROUTINE CheckFileOpen(ios,fileName,openUnit)
  
    IMPLICIT NONE

    INTEGER, INTENT(IN) :: ios, openUnit
    CHARACTER(*), INTENT(IN) :: fileName

    IF (ios < 0) then
      WRITE(*,*) "END-OF-FILE OR END-OF-RECORD OCCURED OPENING", trim(adjustl(fileName)), "ON UNIT", openUnit
    ELSE IF (ios > 0) then
      WRITE(*,*) "ERROR CONDITION OCCURED OPENING", trim(adjustl(fileName)) ,"ON UNIT", openUnit
    ENDIF

  END SUBROUTINE CheckFileOpen

!*

  SUBROUTINE CheckFileClose(ios,fileName,unit)
  
    IMPLICIT NONE

    INTEGER, INTENT(IN) :: ios, unit
    CHARACTER(*), INTENT(IN) :: fileName

    IF (ios < 0) then
      WRITE(*,*) "END-OF-FILE OR END-OF-RECORD OCCURED CLOSING", trim(adjustl(fileName)), "ON UNIT", unit
    ELSE IF (ios > 0) then
      WRITE(*,*) "ERROR CONDITION OCCURED CLOSING", trim(adjustl(fileName)) ,"ON UNIT", unit
    ENDIF

  ENDSUBROUTINE CheckFileClose

!*

  PURE REAL(8) FUNCTION CleanTrigArgument(x) result(x_clean); IMPLICIT NONE

    real(8), intent(in) :: x

    IF (abs(x - 1.0D0) < 1.0D-7) THEN
      x_clean = 1.0D0
    ELSE IF (abs(x+1.0D0) < 1.0D-7) THEN
      x_clean = -1.0D0
    ELSE
      x_clean = x
    ENDIF

  END FUNCTION CleanTrigArgument

END MODULE Toolbox