budgetfilereadermodule Module Reference

Data Types
type	budgetfileheadertype
type	budgetfilereadertype

Functions/Subroutines
subroutine	initialize (this, iu, iout, ncrbud)
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 budget file header. More...
subroutine	rewind (this)

Function/Subroutine Documentation

finalize()

subroutine budgetfilereadermodule::finalize

(

class(budgetfilereadertype)

this

)

Definition at line 248 of file BudgetFileReader.f90.

    class(BudgetFileReaderType) :: this
    close (this%inunit)
    if (allocated(this%flowja)) deallocate (this%flowja)
    if (allocated(this%nodesrc)) deallocate (this%nodesrc)
    if (allocated(this%nodedst)) deallocate (this%nodedst)
    if (allocated(this%flow)) deallocate (this%flow)
    if (allocated(this%auxvar)) deallocate (this%auxvar)
    if (allocated(this%header)) deallocate (this%header)
    if (allocated(this%headernext)) deallocate (this%headernext)

get_str()

character(len=:) function, allocatable budgetfilereadermodule::get_str ( class(budgetfileheadertype), intent(in)  this )

private

Definition at line 261 of file BudgetFileReader.f90.

    class(BudgetFileHeaderType), intent(in) :: this
    character(len=:), allocatable :: str
    character(len=LENHUGELINE) :: temp
 
    write (temp, '(*(G0))') &
      'Budget file header (pos: ', this%pos, &
      ', kper: ', this%kper, &
      ', kstp: ', this%kstp, &
      ', delt: ', this%delt, &
      ', pertim: ', this%pertim, &
      ', totim: ', this%totim, &
      ', budtxt: ', trim(this%budtxt), &
      ', nval: ', this%nval, &
      ', idum1: ', this%idum1, &
      ', idum2: ', this%idum2, &
      ', imeth: ', this%imeth, &
      ', srcmodel: ', trim(this%srcmodelname), &
      ', srcpackage: ', trim(this%srcpackagename), &
      ', dstmodel: ', trim(this%dstmodelname), &
      ', dstpackage: ', trim(this%dstpackagename), &
      ', ndat: ', this%ndat, &
      ', naux: ', this%naux, &
      ', nlist: ', this%nlist, &
      ')'
    str = trim(temp)

initialize()

subroutine budgetfilereadermodule::initialize ( class(budgetfilereadertype)  this, integer(i4b), intent(in)  iu, integer(i4b), intent(in)  iout, integer(i4b), intent(out)  ncrbud  )

private

Definition at line 50 of file BudgetFileReader.f90.

    ! -- dummy
    class(BudgetFileReaderType) :: this
    integer(I4B), intent(in) :: iu
    integer(I4B), intent(in) :: iout
    integer(I4B), intent(out) :: ncrbud
    ! -- local
    integer(I4B) :: ibudterm
    integer(I4B) :: kstp_last, kper_last
    integer(I4B) :: maxaux
    logical :: success
    !
    this%inunit = iu
    this%nbudterms = 0
    ncrbud = 0
    maxaux = 0
    call this%rewind()
    !
    ! -- Determine number of budget terms within a time step
    if (iout > 0) &
      write (iout, '(a)') &
      'Reading budget file to determine number of terms per time step.'
    !
    ! -- Read through the first set of data for time step 1 and stress period 1
    do
      call this%read_record(success)
      if (.not. success) exit
      this%nbudterms = this%nbudterms + 1
      select type (h => this%header)
      type is (budgetfileheadertype)
        if (h%naux > maxaux) maxaux = h%naux
      end select
 
      if (this%endoffile) exit
      if (this%header%kstp /= this%headernext%kstp .or. &
          this%header%kper /= this%headernext%kper) &
        exit
    end do
    kstp_last = this%header%kstp
    kper_last = this%header%kper
    allocate (this%budtxtarray(this%nbudterms))
    allocate (this%imetharray(this%nbudterms))
    allocate (this%dstpackagenamearray(this%nbudterms))
    allocate (this%nauxarray(this%nbudterms))
    allocate (this%auxtxtarray(maxaux, this%nbudterms))
    this%auxtxtarray(:, :) = ''
    call this%rewind()
    !
    ! -- Now read through again and store budget text names
    do ibudterm = 1, this%nbudterms
      call this%read_record(success, iout)
      if (.not. success) exit
      select type (h => this%header)
      type is (budgetfileheadertype)
        this%budtxtarray(ibudterm) = h%budtxt
        this%imetharray(ibudterm) = h%imeth
        this%dstpackagenamearray(ibudterm) = h%dstpackagename
        this%nauxarray(ibudterm) = h%naux
        if (h%naux > 0) then
          this%auxtxtarray(1:h%naux, ibudterm) = h%auxtxt(:)
        end if
        if (h%srcmodelname == h%dstmodelname) then
          if (allocated(this%nodesrc)) ncrbud = max(ncrbud, maxval(this%nodesrc))
        end if
      end select
    end do
    call this%rewind()
    if (iout > 0) &
      write (iout, '(a, i0, a)') 'Detected ', this%nbudterms, &
      ' unique flow terms in budget file.'

read_header()

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

private

Definition at line 124 of file BudgetFileReader.f90.

    ! -- dummy
    class(BudgetFileReaderType), intent(inout) :: this
    logical, intent(out) :: success
    integer(I4B), intent(in), optional :: iout
    ! -- local
    integer(I4B) :: iostat, pos
    character(len=LINELENGTH) :: errmsg
    !
    success = .true.
    select type (h => this%header)
    type is (budgetfileheadertype)
      h%kstp = 0
      h%kper = 0
      h%budtxt = ''
      h%nval = 0
      h%idum1 = 0
      h%idum2 = 0
      h%imeth = 0
      h%srcmodelname = ''
      h%srcpackagename = ''
      h%dstmodelname = ''
      h%dstpackagename = ''
      h%ndat = 0
      h%naux = 0
      h%nlist = 0
      if (allocated(h%auxtxt)) deallocate (h%auxtxt)
 
      inquire (unit=this%inunit, pos=h%pos)
      read (this%inunit, iostat=iostat) h%kstp, h%kper, &
        h%budtxt, h%nval, h%idum1, h%idum2
      if (iostat /= 0) then
        success = .false.
        if (iostat < 0) this%endoffile = .true.
        return
      end if
      read (this%inunit) h%imeth, h%delt, h%pertim, h%totim
      if (h%imeth == 6) then
        read (this%inunit) h%srcmodelname
        read (this%inunit) h%srcpackagename
        read (this%inunit) h%dstmodelname
        read (this%inunit) h%dstpackagename
        read (this%inunit) h%ndat
        h%naux = h%ndat - 1
        if (allocated(h%auxtxt)) deallocate (h%auxtxt)
        allocate (h%auxtxt(h%naux))
        read (this%inunit) h%auxtxt
        read (this%inunit) h%nlist
      elseif (h%imeth /= 1) then
        write (errmsg, '(a, a)') 'ERROR READING: ', trim(h%budtxt)
        call store_error(errmsg)
        write (errmsg, '(a, i0)') 'INVALID METHOD CODE DETECTED: ', h%imeth
        call store_error(errmsg)
        call store_error_unit(this%inunit)
      end if
    end select
    inquire (unit=this%inunit, pos=pos)
    this%header%size = pos - this%header%pos

Here is the call graph for this function:

read_record()

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

private

Definition at line 186 of file BudgetFileReader.f90.

    ! -- modules
    use inputoutputmodule, only: fseek_stream
    ! -- dummy
    class(BudgetFileReaderType), intent(inout) :: this
    logical, intent(out) :: success
    integer(I4B), intent(in), optional :: iout
    ! -- local
    integer(I4B) :: i, n, iout_opt
    !
    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 (budgetfileheadertype)
      if (h%imeth == 1) then
        if (trim(adjustl(h%budtxt)) == 'FLOW-JA-FACE') then
          if (allocated(this%flowja)) deallocate (this%flowja)
          allocate (this%flowja(h%nval))
          read (this%inunit) this%flowja
          this%hasimeth1flowja = .true.
        else
          h%nval = h%nval * h%idum1 * abs(h%idum2)
          if (allocated(this%flow)) deallocate (this%flow)
          allocate (this%flow(h%nval))
          if (allocated(this%nodesrc)) deallocate (this%nodesrc)
          allocate (this%nodesrc(h%nval))
          read (this%inunit) this%flow
          do i = 1, h%nval
            this%nodesrc(i) = i
          end do
        end if
      elseif (h%imeth == 6) then
        if (allocated(this%nodesrc)) deallocate (this%nodesrc)
        allocate (this%nodesrc(h%nlist))
        if (allocated(this%nodedst)) deallocate (this%nodedst)
        allocate (this%nodedst(h%nlist))
        if (allocated(this%flow)) deallocate (this%flow)
        allocate (this%flow(h%nlist))
        if (allocated(this%auxvar)) deallocate (this%auxvar)
        allocate (this%auxvar(h%naux, h%nlist))
        read (this%inunit) (this%nodesrc(n), this%nodedst(n), this%flow(n), &
                            (this%auxvar(i, n), i=1, h%naux), n=1, h%nlist)
      end if
 
      if (iout_opt > 0) then
        write (iout_opt, '(1pg15.6, a, 1x, a)') h%totim, h%budtxt, &
          h%dstpackagename
      end if
    end select
    !
    call this%peek_record()

Here is the call graph for this function:

rewind()

subroutine budgetfilereadermodule::rewind ( class(budgetfilereadertype), intent(inout)  this )

private

Definition at line 289 of file BudgetFileReader.f90.

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