iunitmodule Module Reference

– Module to manage unit numbers. Allows for multiple unit numbers – assigned to a single package type, as shown below. – row(i) cunit(i) iunit(i)nval iunit(i)iunit iunit(i)ipos – 1 BCF6 1 (1000) (1) – 2 WEL 3 (1001,1003,1005) (2,5,7) – 3 GHB 1 (1002) (4) – 4 EVT 2 (1004,1006) (6,10) – 5 RIV 0 () () – ...

Data Types
type	iunitrowtype
type	iunittype

Functions/Subroutines
subroutine	init (this, niunit, cunit)
	Allocate the cunit and iunit entries of this object, and copy cunit into the object. More...
subroutine	addfile (this, ftyp, iunit, ipos, namefilename)
	Add an ftyp and unit number. More...
subroutine	getunitnumber (this, ftyp, iunit, iremove)
	Get the last unit number for type ftyp or return 0 for iunit. More...

Function/Subroutine Documentation

addfile()

subroutine iunitmodule::addfile ( class(iunittype), intent(inout)  this, character(len=*), intent(in)  ftyp, integer(i4b), intent(in)  iunit, integer(i4b), intent(in)  ipos, character(len=*), intent(in)  namefilename  )

private

Find the row for the ftyp and store another iunit value.

Definition at line 63 of file Iunit.f90.

    ! -- dummy
    class(IunitType), intent(inout) :: this
    character(len=*), intent(in) :: ftyp
    integer(I4B), intent(in) :: iunit
    integer(I4B), intent(in) :: ipos
    character(len=*), intent(in) :: namefilename
    ! -- local
    character(len=LINELENGTH) :: errmsg
    integer(I4B), allocatable, dimension(:) :: itemp
    integer(I4B) :: i, irow
    !
    ! -- Find the row containing ftyp
    irow = 0
    do i = 1, this%niunit
      if (this%cunit(i) == ftyp) then
        irow = i
        exit
      end if
    end do
    if (irow == 0) then
      write (errmsg, '(a,a)') 'Package type not supported: ', ftyp
      call store_error(errmsg)
      call store_error_filename(namefilename, terminate=.true.)
    end if
    !
    ! -- Store the iunit number for this ftyp
    if (this%iunit(irow)%nval == 0) then
      allocate (this%iunit(irow)%iunit(1))
      allocate (this%iunit(irow)%ipos(1))
      this%iunit(irow)%nval = 1
    else
      !
      ! -- Increase size of iunit
      allocate (itemp(this%iunit(irow)%nval))
      itemp(:) = this%iunit(irow)%iunit(:)
      deallocate (this%iunit(irow)%iunit)
      this%iunit(irow)%nval = this%iunit(irow)%nval + 1
      allocate (this%iunit(irow)%iunit(this%iunit(irow)%nval))
      this%iunit(irow)%iunit(1:this%iunit(irow)%nval - 1) = itemp(:)
      !
      ! -- Increase size of ipos
      itemp(:) = this%iunit(irow)%ipos(:)
      deallocate (this%iunit(irow)%ipos)
      allocate (this%iunit(irow)%ipos(this%iunit(irow)%nval))
      this%iunit(irow)%ipos(1:this%iunit(irow)%nval - 1) = itemp(:)
      !
      ! -- Cleanup temp
      deallocate (itemp)
    end if
    this%iunit(irow)%iunit(this%iunit(irow)%nval) = iunit
    this%iunit(irow)%ipos(this%iunit(irow)%nval) = ipos

Here is the call graph for this function:

getunitnumber()

subroutine iunitmodule::getunitnumber ( class(iunittype), intent(inout)  this, character(len=*), intent(in)  ftyp, integer(i4b), intent(inout)  iunit, integer(i4b), intent(in)  iremove  )

private

If iremove is 1, then remove this unit number. Similar to a list.pop().

Definition at line 121 of file Iunit.f90.

    ! -- dummy
    class(IunitType), intent(inout) :: this
    character(len=*), intent(in) :: ftyp
    integer(I4B), intent(inout) :: iunit
    integer(I4B), intent(in) :: iremove
    ! -- local
    integer(I4B) :: i, irow, nval
    !
    ! -- Find the row
    irow = 0
    do i = 1, this%niunit
      if (this%cunit(i) == ftyp) then
        irow = i
        exit
      end if
    end do
    !
    ! -- Find the unit number
    iunit = 0
    if (irow > 0) then
      nval = this%iunit(irow)%nval
      if (nval > 0) then
        iunit = this%iunit(irow)%iunit(nval)
        if (iremove > 0) then
          this%iunit(irow)%iunit(nval) = 0
          this%iunit(irow)%nval = nval - 1
        end if
      else
        iunit = 0
      end if
    end if

init()

subroutine iunitmodule::init ( class(iunittype), intent(inout)  this, integer(i4b), intent(in)  niunit, character(len=*), dimension(niunit), intent(in)  cunit  )

private

Definition at line 43 of file Iunit.f90.

    ! -- dummy
    class(IunitType), intent(inout) :: this
    integer(I4B), intent(in) :: niunit
    character(len=*), dimension(niunit), intent(in) :: cunit
    ! -- local
    integer(I4B) :: i
    !
    allocate (this%cunit(niunit))
    allocate (this%iunit(niunit))
    this%niunit = niunit
    do i = 1, niunit
      this%cunit(i) = cunit(i)
    end do