Subroutine READPROP2
Loads element properties into arrays.
Calling convention
call readprop2(elprop, elgeom, nodes, nodelist,'Y','Y','Y','Y')
Arguments | Data type | Input/Output | Description |
---|---|---|---|
elprop |
Real (10,*) array | Return | Description of internal element properties:
The second dimension of the array is equal to or greater than the value of the |
elgeom |
Real (3,4,*) array | Return |
Description of element geometry. X,Y, Z are coordinates of the geometry properties for internal element j
The third dimension of the array is equal to or greater than the value of the |
nodelist |
Real (3,*) array | Return | real(1,i) thru real(3,i) are the X,Y,Z coordinates of node i.Since nodes are referenced by external label the second dimension of this array must be at equal to, or greater than, the largest node label in the model. This value is stored in |
nodes |
Integer (8,*) array | Return | nodes(i,j) will contain the i'th node of internal element j.The second dimension of the array must be equal to or greater than the value of the |
PropFlag |
Character*1 | Call | Is set to 'Y' or 'N'. The array If set to 'N', the argument |
GeomFlag |
Character*1 | Call | Is set to 'Y' or 'N'. The array If set to 'N', the argument |
NodelistFlag |
Character*1 | Call | Is set to 'Y' or 'N'. The array If set to 'N', the argument |
NodesFlag |
Character*1 | Call | Is set to 'Y' or 'N'. The array nodes is only filled is PropFlag is set to 'Y'. If set to 'N', the argument nodes need not be an array. |
Example:
real params
common/params/ params(80000)
..
integer maxnodes
parameter (maxnodes=10000)
integer maxelem
parameter (maxelem=10000)
..
integer nodes(8,maxelem)
real elprop(10,maxelem), elgeom(3,3,maxelem),
+ nodelist(3,maxnodes)
..
save aeroheat, numaero
..
C
C Initialization code
C
if(maxelem .lt. maxno)then
print ,' **FATAL* Element data arrays are not large enough'
print *,' Minimum array size for this model is ', maxno
end if
if(maxnodes .lt. param(292))then
print ,' **FATAL* Node data array is not large enough'
print *,' Minimum array size for this model is ',
+ param(292)
end if
call readprop2(elprop, elgeom, nodes, nodelist,'Y','Y','Y','Y')
..
FORTRAN SAVE
statement to preserve the variables between calls.