interfacemapmodule Module Reference

Data Types
type	interfacemaptype

Functions/Subroutines
subroutine	init (this, nr_models, nr_exchanges)
subroutine	add (this, map_to_add)
	@ Adds a map, either by extending the existing map for a certain model or exchange, or by assigning the map to an empty slot. More...
type(indexmaptype) function, pointer	get_node_map (this, model_id)
type(indexmaptype) function, pointer	get_connection_map (this, model_id)
subroutine	print_interface (this, outunit)
	Dumps interface data to the screen. More...
subroutine	destroy (this)

Function/Subroutine Documentation

add()

subroutine interfacemapmodule::add ( class(interfacemaptype)  this, class(interfacemaptype)  map_to_add  )

private

The map to which is added, should be properly

Definition at line 63 of file InterfaceMap.f90.

    class(InterfaceMapType) :: this
    class(InterfaceMapType) :: map_to_add
    ! local
    integer(I4B) :: im, ie
    integer(I4B) :: m_id, m_index
    integer(I4B) :: e_id, e_index
 
    ! add models
    do im = 1, map_to_add%nr_models
      m_id = map_to_add%model_ids(im)
      m_index = ifind(this%model_ids, m_id)
      if (m_index > 0) then
        ! extend existing index map
        call this%node_maps(m_index)%add(map_to_add%node_maps(im))
        call this%conn_maps(m_index)%add(map_to_add%conn_maps(im))
      else
        ! place in first empty spot
        m_index = ifind(this%model_ids, -1)
        this%model_ids(m_index) = m_id
        this%model_names(m_index) = map_to_add%model_names(im)
        call this%node_maps(m_index)%copy(map_to_add%node_maps(im))
        call this%conn_maps(m_index)%copy(map_to_add%conn_maps(im))
      end if
    end do
 
    ! add exchanges
    do ie = 1, map_to_add%nr_exchanges
      e_id = map_to_add%exchange_ids(ie)
      e_index = ifind(this%exchange_ids, e_id)
      if (e_index > 0) then
        ! extend existing index map
        call this%exchange_maps(e_index)%add(map_to_add%exchange_maps(ie))
      else
        ! place in first empty spot
        e_index = ifind(this%exchange_ids, -1)
        this%exchange_ids(e_index) = e_id
        this%exchange_names(e_index) = map_to_add%exchange_names(ie)
        call this%exchange_maps(e_index)%copy(map_to_add%exchange_maps(ie))
      end if
    end do
 

destroy()

subroutine interfacemapmodule::destroy ( class(interfacemaptype)  this )

private

Definition at line 180 of file InterfaceMap.f90.

    class(InterfaceMapType) :: this
    ! local
    integer(I4B) :: i
 
    do i = 1, this%nr_models
      if (this%model_ids(i) == -1) cycle
      deallocate (this%node_maps(i)%src_idx)
      deallocate (this%node_maps(i)%tgt_idx)
      deallocate (this%conn_maps(i)%src_idx)
      deallocate (this%conn_maps(i)%tgt_idx)
    end do
    deallocate (this%node_maps)
    deallocate (this%conn_maps)
 
    do i = 1, this%nr_exchanges
      if (this%exchange_ids(i) == -1) cycle
      deallocate (this%exchange_maps(i)%src_idx)
      deallocate (this%exchange_maps(i)%tgt_idx)
      deallocate (this%exchange_maps(i)%sign)
    end do
    deallocate (this%exchange_maps)
 
    deallocate (this%model_ids)
    deallocate (this%model_names)
    deallocate (this%exchange_ids)
    deallocate (this%exchange_names)
 

Here is the caller graph for this function:

get_connection_map()

type(indexmaptype) function, pointer interfacemapmodule::get_connection_map	(	class(interfacemaptype)	this,
		integer(i4b)	model_id
	)

Definition at line 125 of file InterfaceMap.f90.

    class(InterfaceMapType) :: this
    integer(I4B) :: model_id
    type(IndexMapType), pointer :: connection_map
    ! local
    integer(I4B) :: m_idx
 
    connection_map => null()
    m_idx = ifind(this%model_ids, model_id)
    if (m_idx > 0) then
      connection_map => this%conn_maps(m_idx)
    end if
 

get_node_map()

type(indexmaptype) function, pointer interfacemapmodule::get_node_map ( class(interfacemaptype)  this, integer(i4b)  model_id  )

private

Definition at line 107 of file InterfaceMap.f90.

    use simmodule, only: ustop
    class(InterfaceMapType) :: this
    integer(I4B) :: model_id
    type(IndexMapType), pointer :: node_map
    ! local
    integer(I4B) :: m_idx
 
    node_map => null()
    m_idx = ifind(this%model_ids, model_id)
    if (m_idx > 0) then
      node_map => this%node_maps(m_idx)
    else
      call ustop("Internal error. Can't find node map in interface")
    end if
 

Here is the call graph for this function:

init()

subroutine interfacemapmodule::init	(	class(interfacemaptype)	this,
		integer(i4b)	nr_models,
		integer(i4b)	nr_exchanges
	)

Definition at line 34 of file InterfaceMap.f90.

    class(InterfaceMapType) :: this
    integer(I4B) :: nr_models
    integer(I4B) :: nr_exchanges
 
    this%nr_models = nr_models
    this%nr_exchanges = nr_exchanges
 
    allocate (this%model_ids(nr_models))
    allocate (this%model_names(nr_models))
    allocate (this%exchange_ids(nr_exchanges))
    allocate (this%exchange_names(nr_exchanges))
 
    allocate (this%node_maps(nr_models))
    allocate (this%conn_maps(nr_models))
    allocate (this%exchange_maps(nr_exchanges))
 
    ! model id == -1 when not set
    this%model_ids = -1
    this%exchange_ids = -1
 

print_interface()

subroutine interfacemapmodule::print_interface ( class(interfacemaptype)  this, integer(i4b)  outunit  )

private

Definition at line 142 of file InterfaceMap.f90.

    class(InterfaceMapType) :: this
    integer(I4B) :: outunit
    ! local
    integer(I4B) :: i, n
 
    write (outunit, '(a,i0)') "nr. models: ", this%nr_models
    write (outunit, '(a,i0)') "nr. exchanges: ", this%nr_exchanges
    do i = 1, this%nr_models
      if (this%model_ids(i) == -1) cycle
      write (outunit, '(3a,i0,a)') "model: ", trim(this%model_names(i)), &
        "[", this%model_ids(i), "]"
      write (outunit, *) "node map:"
      do n = 1, size(this%node_maps(i)%src_idx)
        write (outunit, '(i7,a,i7)') this%node_maps(i)%src_idx(n), &
          " ", this%node_maps(i)%tgt_idx(n)
      end do
      write (outunit, *) "connection map:"
      do n = 1, size(this%conn_maps(i)%src_idx)
        write (outunit, '(i7,a,i7)') this%conn_maps(i)%src_idx(n), &
          " ", this%conn_maps(i)%tgt_idx(n)
      end do
    end do
 
    do i = 1, this%nr_exchanges
      if (this%exchange_ids(i) == -1) cycle
      write (outunit, '(3a,i0,a)') "exchange: ", trim(this%exchange_names(i)), &
        "[", this%exchange_ids(i), "]"
      write (outunit, *) "exchange map:"
      do n = 1, size(this%exchange_maps(i)%src_idx)
        write (outunit, '(i7,a,i7,a,i7)') this%exchange_maps(i)%src_idx(n), &
          " ", this%exchange_maps(i)%tgt_idx(n), &
          " ", this%exchange_maps(i)%sign(n)
      end do
    end do