headfilereadermodule Module Reference

Data Types
type	headfileheadertype
type	headfilereadertype

Functions/Subroutines
subroutine	initialize (this, iu, iout)
subroutine	read_header (this, success, iout)
subroutine	read_record (this, success, iout)
subroutine	finalize (this)
character(len=:) function, allocatable	get_str (this)
	Get a string representation of the head file header. More...
subroutine	rewind (this)

Function/Subroutine Documentation

finalize()

subroutine headfilereadermodule::finalize

(

class(headfilereadertype)

this

)

Definition at line 148 of file HeadFileReader.f90.

    class(HeadFileReaderType) :: this
    close (this%inunit)
    if (allocated(this%head)) deallocate (this%head)
    if (allocated(this%header)) deallocate (this%header)
    if (allocated(this%headernext)) deallocate (this%headernext)

get_str()

character(len=:) function, allocatable headfilereadermodule::get_str ( class(headfileheadertype), intent(in)  this )

private

Definition at line 157 of file HeadFileReader.f90.

    class(HeadFileHeaderType), intent(in) :: this
    character(len=:), allocatable :: str
    character(len=LENBIGLINE) :: temp
 
    write (temp, '(*(G0))') &
      'Head file header (pos: ', this%pos, &
      ', kper: ', this%kper, &
      ', kstp: ', this%kstp, &
      ', pertim: ', this%pertim, &
      ', totim: ', this%totim, &
      ', text: ', trim(this%text), &
      ', ncol: ', this%ncol, &
      ', nrow: ', this%nrow, &
      ', ilay: ', this%ilay, &
      ')'
    str = trim(temp)

initialize()

subroutine headfilereadermodule::initialize ( class(headfilereadertype)  this, integer(i4b), intent(in)  iu, integer(i4b), intent(in)  iout  )

private

Definition at line 34 of file HeadFileReader.f90.

    ! -- dummy
    class(HeadFileReaderType) :: this
    integer(I4B), intent(in) :: iu
    integer(I4B), intent(in) :: iout
    ! -- local
    integer(I4B) :: kstp_last, kper_last
    logical :: success
    !
    this%inunit = iu
    this%nlay = 0
    call this%rewind()
    !
    ! -- Read the first head data record to set kstp_last, kstp_last
    call this%read_record(success)
    kstp_last = this%header%kstp
    kper_last = this%header%kper
    call this%rewind()
    !
    ! -- Determine number of records within a time step
    if (iout > 0) &
      write (iout, '(a)') &
      'Reading binary file to determine number of records per time step.'
    do
      call this%read_record(success, iout)
      if (.not. success) exit
      if (kstp_last /= this%header%kstp .or. kper_last /= this%header%kper) exit
      this%nlay = this%nlay + 1
    end do
    call this%rewind()
    if (iout > 0) &
      write (iout, '(a, i0, a)') 'Detected ', this%nlay, &
      ' unique records in binary file.'

read_header()

subroutine headfilereadermodule::read_header ( class(headfilereadertype), intent(inout)  this, logical, intent(out)  success, integer(i4b), intent(in), optional  iout  )

private

Definition at line 71 of file HeadFileReader.f90.

    ! -- dummy
    class(HeadFileReaderType), intent(inout) :: this
    logical, intent(out) :: success
    integer(I4B), intent(in), optional :: iout
    ! -- local
    integer(I4B) :: iostat
    integer(I8B) :: pos
    !
    success = .true.
    select type (h => this%header)
    type is (headfileheadertype)
      h%kstp = 0
      h%kper = 0
      h%text = ''
      h%ncol = 0
      h%nrow = 0
      h%ilay = 0
      inquire (unit=this%inunit, pos=h%pos)
      read (this%inunit, iostat=iostat) h%kstp, h%kper, &
        h%pertim, h%totim, h%text, h%ncol, h%nrow, h%ilay
      if (iostat /= 0) then
        success = .false.
        if (iostat < 0) this%endoffile = .true.
        return
      end if
      inquire (unit=this%inunit, pos=pos)
      this%header%size = pos - this%header%pos
    end select

read_record()

subroutine headfilereadermodule::read_record ( class(headfilereadertype), intent(inout)  this, logical, intent(out)  success, integer(i4b), intent(in), optional  iout  )

private

Definition at line 104 of file HeadFileReader.f90.

    ! -- modules
    use inputoutputmodule, only: fseek_stream
    ! -- dummy
    class(HeadFileReaderType), intent(inout) :: this
    logical, intent(out) :: success
    integer(I4B), intent(in), optional :: iout
    ! -- local
    integer(I4B) :: iout_opt
    integer(I4B) :: ncol, nrow
    !
    if (present(iout)) then
      iout_opt = iout
    else
      iout_opt = 0
    end if
    !
    call this%read_header(success, iout_opt)
    if (.not. success) return
    !
    select type (h => this%header)
    type is (headfileheadertype)
      ncol = h%ncol
      nrow = h%nrow
    end select
    !
    ! -- allocate head to proper size
    if (.not. allocated(this%head)) then
      allocate (this%head(ncol * nrow))
    else
      if (size(this%head) /= ncol * nrow) then
        deallocate (this%head)
        allocate (this%head(ncol * nrow))
      end if
    end if
    !
    ! -- read the head array
    read (this%inunit) this%head
    !
    call this%peek_record()

Here is the call graph for this function:

rewind()

subroutine headfilereadermodule::rewind ( class(headfilereadertype), intent(inout)  this )

private

Definition at line 176 of file HeadFileReader.f90.

    class(HeadFileReaderType), intent(inout) :: this
 
    rewind(this%inunit)
    this%endoffile = .false.
    if (allocated(this%header)) deallocate (this%header)
    if (allocated(this%headernext)) deallocate (this%headernext)
    allocate (headfileheadertype :: this%header)
    allocate (headfileheadertype :: this%headernext)
    this%header%pos = 1
    this%headernext%pos = 1