GwfVscInputData.f90

Go to the documentation of this file.

module gwfvscinputdatamodule
  use kindmodule, only: i4b, dp
  use constantsmodule, only: lenmodelname, lenauxname, dzero
 
  implicit none
  private
 
  !> Data structure to transfer input configuration to the
  !< VSC package, as opposed to reading from file
  type, public :: gwfvscinputdatatype
 
    ! options
    real(dp) :: viscref !< see VSC for description
    ! dim
    integer(I4B) :: nviscspecies !< see VSC for description
    ! pkg data
    integer(I4B), dimension(:), pointer :: ivisc => null() !< indicates if species uses linear or nonlinear relationship
    real(dp), dimension(:), pointer, contiguous :: dviscdc => null() !< see VSC for description
    real(dp), dimension(:), pointer, contiguous :: cviscref => null() !< see VSC for description
    character(len=LENMODELNAME), dimension(:), allocatable :: cmodelname !< see VSC for description
    character(len=LENAUXNAME), dimension(:), allocatable :: cauxspeciesname !< see VSC for description
 
  contains
    procedure, pass(this) :: construct
    procedure, pass(this) :: destruct
  end type gwfvscinputdatatype
 
contains
 
!> @brief Allocate the input data
!<
  subroutine construct(this, nviscspecies)
    class(gwfvscinputdatatype) :: this !< the input data block
    integer(I4B) :: nviscspecies !< the number of species
 
    allocate (this%ivisc(nviscspecies))
    allocate (this%dviscdc(nviscspecies))
    allocate (this%cviscref(nviscspecies))
    allocate (this%cmodelname(nviscspecies))
    allocate (this%cauxspeciesname(nviscspecies))
 
  end subroutine construct
 
  !> @brief clean up
  !<
  subroutine destruct(this)
    class(gwfvscinputdatatype) :: this !< the input data block
 
    deallocate (this%ivisc)
    deallocate (this%dviscdc)
    deallocate (this%cviscref)
    deallocate (this%cmodelname)
    deallocate (this%cauxspeciesname)
 
  end subroutine destruct
 
end module gwfvscinputdatamodule