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:
  • elprop(1,j) is the area of one side of shell element j, or the volume of solid element j.
  • elprop(2,j) is the emissivity of element j.
  • elprop(3,j) is the solar surface properties of element j combined into a single TSA value.
  • elprop(4,j) is the conductance parameter, K, of element j.

(K * cross sectional area) for a beam

(K * thickness) for a solid shell

  • elprop(5,j) is the capacitance parameter of element j.

(C*ρ * cross sectional area) for a beam

(KC*ρ* thickness) for a solid shell

  • elprop(6,j) is the hydraulic diameter of a fluid element.
  • elprop(7,j) is the number of nodes in element j topology.

The second dimension of the array is equal to or greater than the value of the USER1 argument, maxno.

elgeom Real (3,4,*) array Return

Description of element geometry. X,Y, Z are coordinates of the geometry properties for internal element j

  • elgeom(1,1,j) thru elgeom(3,1,j) define the coordinates of the centroid of element j.
  • elgeom(1,2,j) thru elgeom(3,2,j) define the coordinates of the element calculation point of element j.
  • elgeom(1,3,j) thru elgeom(3,3,j) define a vector corresponding to the surface normal of the shell element j.

The third dimension of the array is equal to or greater than the value of the USER1 argument, maxno.

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 PARAM(292). For more information, see Arguments and Common Blocks.

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 USER1 argument, MAXNO.

PropFlag Character*1 Call Is set to 'Y' or 'N'.

The array elprop is only filled is PropFlag is set to 'Y'.

If set to 'N', the argument elprop need not be an array.

GeomFlag Character*1 Call Is set to 'Y' or 'N'.

The array elgeom is only filled is PropFlag is set to 'Y'.

If set to 'N', the argument elgeom need not be an array.

NodelistFlag Character*1 Call Is set to 'Y' or 'N'.

The array nodelist is only filled is PropFlag is set to 'Y'.

If set to 'N', the argument nodelist need not be an array.

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')
        ..
Note:
The routine should only be called once. It is therefore necessary to use the FORTRAN SAVE statement to preserve the variables between calls.