CalculiX  2.13
A Free Software Three-Dimensional Structural Finite Element Program
networkinum.f File Reference

Go to the source code of this file.

Functions/Subroutines

subroutine networkinum (ipkon, inum, kon, lakon, ne, itg, ntg, network)
 

Function/Subroutine Documentation

◆ networkinum()

subroutine networkinum ( integer, dimension(*)  ipkon,
integer, dimension(*)  inum,
integer, dimension(*)  kon,
character*8, dimension(*)  lakon,
integer  ne,
integer, dimension(*)  itg,
integer  ntg,
integer  network 
)
21 !
22 ! assigns a negative sign to the inum values corresponding to
23 ! network nodes
24 !
25  implicit none
26 !
27  character*8 lakon(*),lakonl
28 !
29  integer ipkon(*),inum(*),kon(*),ne,indexe,i,node1,node2,node3,
30  & itg(*),ntg,network
31 !
32 ! changing the sign of inum for all network nodes
33 !
34  do i=1,ntg
35  if(inum(itg(i)).gt.0) inum(itg(i))=-inum(itg(i))
36  enddo
37 !
38 ! changing the sign of inum for all nodes belonging to
39 ! "D "-elements (network elements without type in purely
40 ! thermal networks)
41 !
42  do i=1,ne
43 !
44  if(ipkon(i).lt.0) cycle
45  lakonl=lakon(i)
46  if((lakonl(1:1).ne.'D ').and.
47  & ((lakonl(1:1).ne.'D').or.(network.ne.1))) cycle
48 !
49  indexe=ipkon(i)
50  if(kon(indexe+1).ne.0) then
51  node1=kon(indexe+1)
52  if(inum(node1).gt.0) inum(node1)=-inum(node1)
53  endif
54  node2=kon(indexe+2)
55  if(inum(node2).gt.0) inum(node2)=-inum(node2)
56  if(kon(indexe+3).ne.0) then
57  node3=kon(indexe+3)
58  if(inum(node3).gt.0) inum(node3)=-inum(node3)
59  endif
60  enddo
61 !
62  return
Hosted by OpenAircraft.com, (Michigan UAV, LLC)