HELP! This list is far from complete. Some items have no description and others are missing all together. Please send us an email if you can provide info.
Here is the NU table.
Be sure to check the Notes and Obsolete sections at the end.
UNKNOWN!
A(ns,ns) 
? 
B(ns,ns) ? 

COVA(NEMX,NRMX) 
used in data/graphics operations? 
EIGV(NOMX,NTMX) 
? 
ER(nv) 
? 
MXI(2,NEMX) 
used for graphics? 
NLCHOR(0:10,NRMX) 
used for BEM problems? 
NPB(1..6,0) 
? 
NPB(1..6,np) 
? 
NTCOVA(NEMX) 
used in data/graphics operations? 
XA(na,nj,nq) 
nodes.global_param_XA 
XF(ns,nj) 
? 
Y 
? 
ZA(na,nh,ne) 
nodes.deform_global_param_ZA 
ZC(nj,ne) 
used in chnode.f, drelem.f, and moelem.f 
ZD2(nj,nd) 
used in write.f calcnods.f 
ZDD(nd,nj) 
data.zdd? 
ZF(ns,nh) 
? 
A 

ACOEFF(nactv) 
nactv=1,NTACTV are coefficients for linear dynamic terms 

ALFA(nactv) 
nactv=1,NTACTV are time constants for linear dynamic terms 

AXL,AXU 
are deformed metric tensor components wrt undeformed Nucoords (KTYP53>1) or wrt undeformed theta coords (KTYP53=1) in zere50.f. SeeRI1,RI2,RI3. 

B 

colloc.bezier_control_pnts 
BCP(4,4,NJMX,NEMX) 
Bezier control points for rendering [collocation] 

C 

C(NYMX,1) 
for BEM problems (solve4.f) 

CC(NYMX,1) 
for BEM problems (solve4.f) 

CCOMP(NCYMX,NCYMX) 
for BEM problems (solve4.f) 

prob.const_mat_param 
CE(nm,ne) 
is value of piecewise constant material parameters nm in element ne 

prob.gauss_mat_param 
CG(nm,ng) 
is value of material parameter nm at Gauss point ng 

CGA(3,ng,nb) 
accumulated coordinate transformation matrices [collocation] 

CGB(3,3,ng,nb) 
accumulated coordinate transformation matrices [collocation] 

colloc.integr_scaling_factor 
CGDA(ng,nb,ne) 
integration scaling factor at the gauss points of basis nb [collocation] 

colloc.mat_param 
CGMAT(5,ng,nb,ne) 
material parameters [collocation] 

data.opt_coefs 
COYY(noy,ny) 
coefficient for optimization DOF NOYY(noy,ny) 

prob.linear_mat_const 
CP(nm,np) 
is value of piecewise linear material constant at node np 

CSEG(nosg) 
is string defining segment number NOSG as follows: 

D 

DATYPE(220) 
character string of data type drawn in graphics window (set in fe24/drdata.f). Is this ever used anywhere? 

DEL_T 
is the time step used (taken from load stepping loop in FE07) 

elem.face_area 
DF 
is the logarithmic frequency increment in LPT 

DF(nf) 
is area of face segment nf 

nodes.arc_length_deriv 
DL(1..3,nl) 
are arclength derivatives of line segment nl 

nodes.old_scaling_factors 
DLL(3,nl) 
holds old scaling factors DL when rescaling nodes (fe21/upnode.f) 

DRDN(NGMX) 
for BEM problems (solve4.f) 

DRDNS(NGMX) 
for BEM problems (solve4.f) 

DXIX 
are partial derivatives of Xi wrt Xj coords (if KTYP53=1) or wrt Nu coords (if KTYP53>1) 

DZDX 
are components of the deformation gradient tensor 

E 

colloc.elem_damping 
ED(nv,nv) 
is element damping matrix (1st order time derivatives) 

data.diff_XID_ZD 
EDD(nd) 
difference between element point at XID and data point ZD, weighted by WD (see fe03/fitfld.f) 

EG 
are physical components of Green’s strain 

march1.elem_mass 
EM(nv,nv) 
is element mass matrix (2nd order time derivatives) 

fit.er 
ER 
is element right hand side; computed in fe03/zder.f 

ERRMF 
is solution tolerance (epsilon) for the iterative solver 

colloc.elem_stiffness 
ES(nv,nv) 
is element stiffness matrix (0th order time derivatives) 

ESE(nv,nv,ne) 
is element stiffness matrix for element ne (replaces ES when using iterative solver) 

ETYP(ie) 
ie=1,12 is .TRUE. if element type ie (see TITLE2) used. 

EV1RE,EV1IM 
are the minimum (EV1*) and maximum (EV2*) Real and Imaginary parts of eigenvalues for the constraintreduced global system of equations (used in some iterative solvers) 

EXR 
are extension ratios wrt COORDS 

F 

prob.fiber_exten 
FEXT(8,NGMX,NEMX) 
from zgtg5a.f:


FLOW_COEFFS(no_coeffs) 
no_coeffs=1,nt_coeffs are Fourier coeffs for flow 

data.dof_data 
FIX(ny,5) 
FIX(ny,1)FIX(ny,2)FIX(ny,3)


prob.applied_pi 
FIXP(2,ne) 
is TRUE if a pressure incrementPEis applied to either the Xi_3=0 (1) or Xi_3=1 (2) face of element ne 

G 

colloc.global_damping 
GD(ny,ny) 
is global damping matrix (1st order time derivatives) 

colloc.global_stiffness 
GK(ny,ny) 
is global stiffness matrix (0th order time derivatives) 

GKC(NCYMX,NCYMX) 
for BEM problems (solve4.f) 

colloc.reduced_stiffness 
GKK(nz) 
is contraintreduced global stiffness matrix in 1D (vector) form 

march1.global_mass 
GM(ny,ny) 
is global mass matrix (2nd order time derivatives) 

GMC(NCYMX,NCYMX) 
for BEM problems (solve4.f) 

march1.gr1_rhs_vector 
GR1(ny) 
is RHS vector for system of global equations (with constraints) 

nonlin.gr2_rhs_vector 
GR2(ny) 
is RHS vector for system of global equations (with constraints) (BFGS method) 

GRC(NCYMX) 
is the vector resulting from BEM domain integrals (fe90/solve4.f) 

fit.grr_rhs_vector 
GRR(no) 
is RHS vector for system of constraintreduced global equations 

GRRC(NCYMX) 
for BEM problems (solve4.f) 

GXL,GXU 
are undeformed metric tensor components wrt Xicoords (zere50.f). 

GZ 
is the determinant of GZL 

GZL,GZU 
are the covariant & contravariant components 

H 

fit.smoothing_matrix 
HK(ny,ny) 
smoothing matrix (KTYP12 > 0) for fitting with smoothing constraints 

I 

elem.index_basis_type 
IBT(1,ni,nb) 
is index for basis type nb in Xi direction ni: 

IBT(2,ni,nb) 
=1 for Lagrange linear 

elem.index_deriv_order 
IDO(nk,0,nb) 
is index for derivative order: NU partial derivative number to derivative nk (seeNU table) of basis nb 

IDO(nk,ni,nb) 
is an index for derivative order: 1=zeroth order, 2=first order 

colloc.extern_face 
IECG(NK,NN,NH,NE) 
“external face” tables to figure out which element parameters are handled by the Galerkin method. [collocation] 

colloc.global_equa 
IGCG(NYMX) 
vector indicates which global equations are from collocation (=1) and which are from Galerkin conditions (>1). [collocation] 

elem.index_nodal_pos 
INP(nn,ni,nb) 
gives the index for element node nn in each Xi direction. 

IOTYPE 
1 = prompt, 2 = read, 3 = write, 4 = read & list 

ISAXES(iw) 
is segment number of axes 

ISBASE(nb) 
is segment number of basis function type nb 

ISCONO(nh,ne) 
is segment number of contour numbers in element ne 

ISCONT(nh,ne,nocont) 
is segment number of contour nocont of variable nh in element ne 

ISDANO(iw,ne) 
is segment number of data point numbers 

ISDAPR(iw,ne) 
is segment number of data point projections 

ISDATA(iw,nodata) 
is segment number of data points at set nodata 

ISDATR(iw,ne) 
is segment number of data point trace in element ne 

ISEG(nosg) 
is 0,1,2 if segment not yet created / created but not visible / created and visible 

ISELNO(iw,ne) 
is segment number of element numbers 

ISFACE(iw,nf) 
is segment number of face nf 

ISFANO(iw,nf) 
is segment number of face numbers 

ISFIBR(iw,ne,nofibr) 
is segment number of fibres in element ne at set nofibr 

ISGAUS(iw) 
is segment number of Gauss points 

ISGRID(iw) 
is segment number of grid 

ISHIST(0) 
is segment number of time history axes & labels 

ISHIST(np) 
is segment number of time history at node np 

ISINCR(iw) 
is segment number of increments 

ISISOC(iw,noisoc) 
is segment number of isochrones 

ISL2BE(nl) 
is segment number of Bezier tangent line 1 on nl 

ISL3BE(nl) 
is segment number of Bezier tangent line 2 on nl 

ISLINE(iw,noline) 
is segment number of lines at set noline 

ISLINO(iw) 
is segment number of line numbers 

ISMAP(nomap) 
is segment number of map 

ISMATL(iw,ne) 
is segment number of material in element ne 

ISN2BE(nl) 
is segment number of Bezier control pt 1 on nl 

ISN3BE(nl) 
is segment number of Bezier control pt 2 on nl 

ISNONO(iw,np) 
is segment number of node numbers 

ISREAC(iw) 
is segment number of reactions 

ISSECT(nosect) 
is segment number of section 

ISSTRE(ne,nostre) 
is segment number of principal stresses at set nostre 

ISSTRM(ne,nostrm) 
is segment number of streamline 

ISSURF(ne) 
is segment number of surface grid in element ne 

ISVELO(ne,novelo) 
is segment number of velocity field 

data.i_threshold 
ITHRES(NGMX,NEMX) 
used in threshold modelling activation pattern computations (fe30). Value at Gauss point ng of element ne is…


ITYP1 
is 3,4,5 or 9 for use of FE30,FE40,FE50 or FE90 

ITYP2 
is equation type [formerly KTYP1 ? FJV] 

ITYP3 
is equation type qualifier [formerly KTYP11 ? FJV] 

IWKDEF(0) 
is number of open windows 

IWKDEF(noiw) 
noiw=1,IWKDEF(0) is list of open windows 

IWKG(iw) 
is 0 for nongraphics window (eg menu) 

IWKS(iw) 
= 0 workstation (window) iw is not defined (ie, not open) 

IWKT(iw) 
= 1 for GKS workstation 

IWRIT1 
controls output printing frequency. 

IWRIT2 
is 1,2 for equilibrium solution only / intermediate solutions also 

IWRIT3 
is 1,2 for solution vectors only / residual vectors also 

IWRIT4 
is 0..4 for output from linear / nonlinear solver: 

IWRIT5 
is 0..5 for output from iterative solver: 

J 

data.elem_def 
JTYP1 
is 1,2 for elements defined by user / chosen from menu 

data.ens_elem_map 
JTYP2 
is 0,1 for ensembletoelement map for nodal derivatives is standard / nonstandard 

data.coor_sys 
JTYP3 
is 1..5 for coordinate system: rect. cartesian / cylindrical polar / spherical polar / prolate spheroidal / oblate spheroidal 

data.sym_type 
JTYP4 
is 1..3 for geometry unsymmetric / cyl. symm. / sph. symm.


data.basis_rep 
JTYP5 
is 1,2 for basis functions in Lagrange or Hermite / monomial format 

data.coor_const 
JTYP6 
is 1,2 for global coordinate system constant / specified by elements 

data.deformed_coor_sys 
JTYP7 
is 1..5 for dependent variable coordinate system rect. cartesian / cylindrical polar / spherical polar / prolate spheroidal / oblate spheroidal 

data.basis_output 
JTYP8 
is 0,1 for no output / output of basis functions 

data.fiber_field_type 
JTYP9 
is 0,1,2 for fiber direction field not defined / defined / sheet direction field defined 

data.isochoric_interp 
JTYP10 
is 1,2,3 for type of ‘radial’ interpolation (JTYP3>1 only) 

data.num_of_field_vars 
JTYP11 
is the number of additional geometric / field variables 

data.fiber_ref_axis 
JTYP12 
is 1,2 for fibres defined wrt Xi1 or Xi2 coordinates 

data.fiber_ang_units 
JTYP13 
is 1,2 for fibers entered in degrees / radians 

data.special_mesh_type 
JTYP14 
is mesh type for specialized meshes as follows: 1 & 2) Fractal tree with branch parameters: 

data.extern_stim 
JTYP15 
0,1 : external stimulus not defined/defined (collocation only) 

elem.num_of_mFHN_params 
JTYP16 
number of material parameters for (m)FHN equations (collocation only) 

data.scalar_field 
JTYP17 
0,1 : scalar field (e.g. consistent strains) not defined/defined (collocation only) 

K 

data.prob_type 
KTYP1 
is 1..15 for problem type [or ITYP2 ? — FJV]


data.FE_method 
KTYP2 
is 1..5 for Galerkin finite elements / direct boundary elements / indirect boundary elements / orthogonal collocation / finite element collocation (formerly arrayITYP4) 

data.time_domain_type 
KTYP3 
is 1..5 for static / time integration / modal analysis / Fourier analysis / buckling analysis (formerly arrayITYP5) 

data.prob_linear 
KTYP4 
is 1..2 for linear / nonlinear problem (formerly arrayITYP6) 

data.init_state 
KTYP5 
is 1..3 for initial solution zero / read in / restarted 

data.gauss_pnt_fit 
KTYP6 
is 1 for Gauss point fitting; is number of boundary integral equation domains 

data.equa_params 
KTYP7 
is 1..3 for equation parameters constant wrt time / user defined in subroutine USER / read from file IPC at each time step 

data.fitting_type 
KTYP8 
is 1..5 for geometry / fibre or field / motion / optical flow / Fourier motion fitting 

data.nonli_solver 
KTYP9 
is 1..4 for full Newton / modified Newton / BFGS inverse / elementbyelement method 

data.line_search 
KTYP10 
is 1..2 for solution with no search / linear search 

data.opt_prob 
KTYP11 
is 1..5 for additional options on problem type [or ITYP3 ? — FJV]


data.constraint_type 
KTYP12 
is 0..3 for fitting without / with constraints: 

data.ode_sovler_type 
KTYP13 
is 1 if pressure read from file (PRESS.VSAERO) 

data.inc_params 
KTYP14 
is > 0 if material parameter is incremented 

KTYP15 
is equation type parameter 

KTYP16 
is 1..2 for lowest / highest eigenvalues required 

KTYP17 
is number of eigenvalue pairs required 

KTYP18 
is number of subspace iteration vectors 

KTYP19 
is number of starting vectors 

KTYP22 
is 1..3 for time integration algorithm linear / quadratic / cubic 

KTYP23 
is fixed time step / automatic stepping / read from file 

KTYP25 
is type of driving function in Fourier analysis 1..3 for impulse / step / sine wave 

KTYP26 
is 1..2 for optimization of material params / geometric params 

KTYP27 
is 1..9 for type of minimization objective function 

KTYP28 
is number of sets of measurements in fit 

KTYP31 
is 1..2 for Cardaic activation model implemented forwards / backwards 

KTYP43 
is 0..3 for thermal strains not included / included as fixed initial strain / constrained by displacement b.c.s / fully coupled 

KTYP45 
is 1..4 for type of beam crosssection 

data.elast_prob 
KTYP51 
is 1..6 for plain stress / plain strain / 3D / membrane/ thin shell / thick shell 

data.compressibility 
KTYP52 
is 1..3 for compressible / incompressible / incomp with fluid perfusate 

data.anisotropy 
KTYP53 
is 1..3 for isotropic / aeleotropic / aeleotropic + active fibres 

data.hyperelasticity 
KTYP54 
is 1..3 for hyperelastic / Cauchyelasticity / creep 

data.kinematic_params 
KTYP55 
is 1..3 for princ. strain invariants / extension ratios / fibre strains 

data.energy_function_form 
KTYP56 
is 1..3 for polynomial / special function / exponential strain energy function 

data.press_bc_type 
KTYP57 
is no pressure bc / press.incr / press.read / volume incr / vol.computed 

data.isochoric 
KTYP58 
is 1,2 for conventional / isochoric element 

data.active_stress_terms 
KTYP59 
is number of terms in active fibre stress relation 

data.SL_parameter 
KTYP5A 
is material parameter number of stressfree SL distribution 

data.time_delay_param 
KTYP5B 
is material parameter number of timedelay variable 

data.Vcf_param 
KTYP5C 
is material parameter number of capillary volume fraction 

KTYP5D 
is unused 

KTYP71 
is 1 if pressure loads read from file PRESS.VSAERO (ID=14) after flow solution by VSAERO 

KTYP90 
is saturatedunsaturated / heartbody/ ? coupling 

L 

data.associated_data_point 
LD(nd) 
line or face l number associated with data point nd 

LDR(nd) 
is unused and should be removed 

LFR(no) 
is unused and should be removed 

prob.elem_location 
LGE(nv) 
is location of element variable nv in global system. LGE(nv) is negative if NE is the last element to reference that variable. 

LGEE(nv,ne) 
is location of element variable nv in global system for element ne (replaces LGE when using iterative solver) 

data.fitting_elem 
LN(0) 
number of elements in fitting 

LN(l) 
element number for l=1..LN(0) 

prob.boundary_region_var 
LRE(nv) 
is for the shell/fluid interface case. It is the boundary region variable LRE of the local finite element variable NVE. (fe02/melge.f). This may be unused. 

LUMP 
is .TRUE. if mass lumping is used 

M 

prob.last_influenced_elem 
ME(np) 
is the last element to be influenced by global node NP. WHERE IS THIS SET??? 

MOTION_TYPE 
is 1,2 for Fourier coeffs / Spreadsheet column 

N 

NA 
auxiliary element variable (na=1,NAT(nb)) 

dims.num_aux_var 
NAMX 
is maximum number of auxiliary parameters 

elem.poly_degree 
NAN(ni,na,nb) 
is polynomial degree in Xi(ni) direction for basis nb for auxiliary variable na 

NAT(nb) 
is number of auxiliary or spline basis functions 

NB 
basis function type (nb=1,NBT) 

NBC(nb) 
is basis function type choice 

elem.basis_num 
NBH(nh,ne) 
is basis number for dependent variable nh in element ne 

NBI(nb) 
is 1..5 for unit scale factors / elem scale factors read in / global scale factors read in / arclength / anglechange 

elem.basis_type_num 
NBJ(nj,ne) 
is basis function type number for geometric variable nj in element ne 

NBL(0,0,nb) 
is number of boundary line / face segments for blending functions 

NBL(0,ls,nb) 
is nb number of boundary line / face segment ls 

NBL(1..4,ls,nb) 
are element node numbers of line / face ls (used for blending function interpolants) 

colloc.num_basis 
NBMX 
maximum number of basis functions 

NBT 
number of basis function types 

NB_MOTION 
Fourier basis number 

data.data_in_elem 
NCMX 
maximum number of equations for a given dependent var (the old meaning of NCMX was the max number of data points in an element) 

elem.coor_type 
NCO(ne) 
is the coordinate type used in element ne 

dims.BE_array_size 
NCYMX 
maximum size of complex arrays (for BEM problems) 

ND 
data point (nd=1,NDT) 

NDAL(0:nde) 
is unused and should be removed 

data.global_data_pts 
NDDL(ne,nde) 
global data pt no. of local data point nde 

NDEMX 
number of data points in one element 

data.num_of_data_pts 
NDLT(ne) 
number of data points within element ne 

data.data_points 
NDMX 
number of data points total 

data.data_pts_index 
NDP(NDMX) 
index of data points that project onto an element 

NDT 
is total number of data points 

NE 
an element (ne=1,NET) 

nodes.shared_corner 
NECOR(NPMX,0:4) 
stores the number of elements which share the current corner or edge node NP. NECOR(NP,0) is the total number of elements sharing the corner node (BEM problems, fe90/solve4.f) 

NEELEM(0) 
total number of elements. 

NEELEM(noelem) 
noelem=1..NEELEM(0) are the element numbers. 

colloc.elements 
NEMX 
maximum number of elements 

NEIBS(NI,IN,NE) 
Element neighbors of NE on INth face in NIth direction [collocation] 

colloc.local_node 
NEP(1,IN,np) 
the INth local node number coincident with global node np [collocation] 

NEP(2,IN,np) 
the element the INth local node belongs to. [collocation] 

NET 
is the highest element number 

NF 
global face (nf=1,NFT) 

NFE(nb) 
number of faces for element basis type nb. 

elem.global_face_num 
NFF(nf,ne) 
are the global face numbers of side nf of element ne 

dims.global_face_seg 
NFMX 
maximum number of global face segments 

NFT 
total number global element faces 

NG 
Gaussian quadrature point (ng=1,NGT(nb)) 

elem.num_gauss_pts 
NGAP(ni,nb) 
is number of Gauss points in Xi direction ni for basis nb 

colloc.gaus_per_elem 
NGMX 
maximum number Gauss points per element 

NGT(nb) 
number of Gauss points per element 

NH 
dependent variable (nh=1,NHP(np)) 

elem.num_dep_var_elem 
NHE(ne) 
number of dependent variables defined in element ne 

colloc.depen_var 
NHMX 
maximum number of dependent variablesNote:NHMX must be equal to NJMX for problems where the dependent variable array carries deformed coordinates. 

NHO 
is Gauss variable to be fitted in Gauss point fitting 

nodes.num_dep_var_node 
NHP(np) 
is number of dependent variables defined at node np 

NHT(nje,nve,ie) 
is number of global variables required for element type ie when number of local variables is nve & number of dimensions nje.


NHV(nv,ie) 
is global variable number, 1 .. NHE(ne), corresponding to dependent variable nv in element type ie. 

NI 
Xicoordinate (ni=1,NIT(nb)) 

colloc.local_xi 
NIMX 
maximum number of local Xicoordinates 

NIT(nb) 
is number of local Xicoordinates for basis nb 

NJ 
Xjcoordinate (nj=1,NJT) 

elem.num_xj_coor_in_elem 
NJE(ne) 
number of Xjcoords defined in element ne (excluding JTYP9) 

NJG 
is geometric variable number in linear field fitting. 

colloc.global_ref_coor 
NJMX 
maximum number of global reference coordinates 

NJO 
is field variable number in linear field fitting (eg =NJT+1). 

nodes.num_xj_coor 
NJP(np) 
is number of Xjcoordinates (ie, geometric variables) defined at node np 

NK 
derivative number (nk=1,NKT(nb)) 

nodes.nodal_derv_dep_var 
NKH(nh,np) 
is number of nodal derivatives for dependent variable nh at node np (also called nodal_derv_dep_var in some places) 

nodes.nodal_derv_geom_var 
NKJ(nj,np) 
is number of nodal derivatives for geometric variable nj at node np 

colloc.deriv_per_var 
NKMX 
maximum number of derivatives per variable 

elem.nodal_spline_term 
NKT(nb) 
is number of nodal derivatives or spline polynomial terms for basis nb 

NL 
global line (nl=1,NLT) 

NLE(nb) 
is number of element line segments for basis nb 

elem.line_num_of_arcface 
NLF(naf,nf) 
are global line numbers of local arc naf of face nf 

elem.line_num_of_arcelem 
NLL(nae,ne) 
are global line numbers of local arc nae of element ne 

dims.global_line_seg 
NLMX 
maximum number of global line segments 

NLT 
total global line segments 

NLV(ie,njt) 
is number of local variables possible in element type ie. 

NM 
material parameter (nm=1,NMT) 

dims.materials 
NMMX 
maximum number of material parameters 

NMOPTI(noopti) 
noopti=NMOPTI(0)+1,NTOPTI is list of displacements in fit 

NN 
element node (nn=1,NNT(nb)) 

elem.num_of_elem_nodes_in_face 
NNF(0,nf,nb) 
is number of element nodes in face NF of element with basis NB 

NNF(1 ,nf,nb) 
is Xidirection normal to face 

NNF(11.14,nf,nb) 
are element derivative numbers in face 

NNF(2..10,nf,nb) 
are the element node numbers in face 

elem.elemnode_face_num 
NNL(1..4,nae,nb) 
are the element / face node numbers along local arc nae 

colloc.nodes_pre_elem 
NNMX 
maximum number of nodes per element 

colloc.coincident_nodes 
NNP(np) 
the number of element nodes coincident with global node np [collocation] 

elem.elem_nodes 
NNT(nb) 
is number of element nodes for basis nb 

NO 
a d.o.f in the set of global equations (no=1,NOT) 

data.dof 
NOMX 
maximum number of optimization DOFs. 

NORMX 
maximum number of size of fractal tree order arrays 

NOT 
total number of optimization DOFs 

data.mesh_opt_dof 
NOYMX 
maximum number of optimization DOFs attached to one mesh DOF 

data.num_of_optimization_dof 
NOYT(ny) 
number of optimization DOFs coupled to mesh variable ny 

data.optimization_dof_num 
NOYY(noy,ny) 
optimization DOF number of variable noy for mesh variable ny 

NP 
global node (np=1,NPT) 

NPCOR(0:NPMX,2) 
NPCOR(i,1)=0 if node i is not on a corner =column number of GM matrix for second dH/dn contribution if node i is at some corner or edge (3D) (the column number corresponds to the contribution from the element at the corner with the higher element number if it is a 2d corner or 3d edge, otherwise it corresponds to the middle corner element number). 

dims.domain_nodes 
NPDMX 
maximum number of domain nodes for BEM problems 

nodes.global_node_num 
NPE(nn,nb,ne) 
is global node number of local node nn of element ne for basis nb 

NPE(nn,nb,nf) 
are global node numbers of face nf, where nb is basis type number for first geometric variable. 

nodes.nodal_basis_type 
NPF(1,nf) 
is the 1st Xidirection of face segment NF 

NPF(2,nf) 
is the basis fn type for 1st Xidir (1,2,3 or 4) 

NPF(3,nf) 
is the 2nd Xidirection 

NPF(4,nf) 
is the basis function type for 2nd Xidir 

NPF(5..8,nf) 
are the basis function type (nb) numbers for nj=1..4 

NPF(9,nf) 
is number of elements adjoining face (1,2 for external / internal) 

nodes.connect 
NPL(1,nl) 
is the Xidirection of line segment NL 

NPL(2..5,nl) 
are the basis function types for nj=1..4 : 

NPL(6..9,nl) 
are the global nodes along line NL in the direction of Xi (8 & 9 are NK numbers of 1st derivatives wrt Xi for any geometric variable which is cubic Hermite) 

NPL(10,nl) 
is the number of elements adjoining line NL 

NPL(11..18,nl) 
are the element numbers of elements adjoining line NL 

colloc.global_nodes 
NPMX 
maximum number of global nodes 

NPNODE(0) 
is the total number of nodes. 

NPNODE(nonode) 
nonode=1..NPNODE(0) are the node numbers. 

data.global_node_fitting 
NPO(0) 
number of global nodes in data fitting 

NPO(np) 
global node number corresponding to n=1..NPO(0) 

NPT 
is the highest node number. 

NP_MOTION 
is node number for applying motion 

NQ 
a general DOF (nq=1,NQT) 

nodes.global_DOF_num 
NQE(ns,nb,ne) 
is global DOF number of local DOF ns in element ne for basis nb. 

dims.global_dof 
NQMX 
maximum number of global DOFs 

NQT 
total number of global DOFs / geometric variables in XA 

dims.contours_per_elem 
NRMX 
maximum number of contours per element & various segs 

NS 
ns=nk+(nn1)*NKT(nb) 

colloc.elem_dof 
NSMX 
maximum number of element DOFs per variable ( <= NNMX*NKMX ) 

NST(nb) 
total number of element parameters for basis nb [ =NKT(nb)*NNT(nb)+NAT(nb) ] 

NTACTV 
is number of dynamic terms in the material response function 

dims.eigenvalues 
NTMX 
maximum number of eigenvalues 

NTSG 
is current number of segments 

NTYP1 
is 1..16 for type of iterative solver: 

NTYP2 
is the number of Krylov subspace basis vectors (used for NTYP1 = 4, 6, 7, or 8)NOTE:max is 250 for RGMRES & RGMRESEV. 

NTYP3 
is the stopping criterion for the iterative solver 

NTYP4 
is the maximum number of iterative solver restarts or iterations 

NTYP5 
is the type of preconditioning for the iterative solver: 

colloc.deriv_terms 
NUMX 
maximum number of derivative terms up to 2nd order 

NU 
is index of Xicoordinate derivative
e.g.: u1 = partial(u)/partial(Xi_1) 

NUT(nb) 
is number of Xicoordinate derivatives up to 2nd order for basis nb 

NVE(ie) 
is number of local variables used in element type ie 

colloc.total_elem_dof 
NVMX 
maximum number of element DOFs ( = NHMX*NSMX ) 

elem.type_num 
NW(ne) 
is the type number (= 1,2..12) for element ne. 

dims.workstations 
NWMX 
maximum number of workstations for segment arrays (e.g. ISAXES, ISFIBR…) 

elem.elem_adjacency 
NXI(ni:ni,0:ne) 
Filled in fe02/nenxi.f. Element number adjacent to element ne: 

dims.image_cell 
NXMX 
maximum number of image cell array dimension (ASSERTed in deimag.f) 

NY 
a finite element mesh DOF (ny=1,NYT) 

colloc.mesh_dof 
NYMX 
maximum number of mesh DOFs. 

prob.ny_ne_mapping 
NYNE(na,nh,ne) 
is mapping from na,nh,ne to ny (set in ipini5.f) 

prob.ny_np_mapping 
NYNP(nk,nh,np) 
is mapping from nk,nh,np to ny (set in ipini5.f) 

NYT 
total finite element mesh DOFs (initialized in ipini5.f) 

NZ 
a stiffness matrix component in a 1D vector (nz=1,NZT) 

colloc.stiffness_vector 
NZMX 
maximum number of coefficients in 1D global stiffness vector GKK. 

NZT 
total number of nonzero coefficients in 1D global stiffness vector GKK. 

O 

P 

prob.pressure_incr 
PE(1..2,ne) 
the pressure increment applied to the Xi_3=0 face (1) or Xi_3=1 face (2) of element ne 

prob.pressure_load 
PF(1..2,ne) 
the pressure load applied to the Xi_3=0 face (1) or Xi_3=1 face(2) of element ne 

PG(ns,nu,ng,nb) 
are basis function values for element DOF ns at Gauss point ng, basis number nb. Set in fe02/gauss1.f when bases are defined. 

PGF(ns,nu,ng,nf,nb) 
basis function values (as in PG) on face nf. [collocation] 

PHI 
are the Euler angles wrt COORDS of the principal extensions 

PMAX(noopti) 
maximum parameter values allowed 

PMIN(noopti) 
minimum parameters value allowed 

PROMPT 
is .TRUE. if computations await prompt after IWRIT1 steps. 

PST 
are principal strains / stresses 

Q 

R 

R 
is the orthogonal rotation tensor 

RAD(NGMX) 
for BEM problems (solve4.f) 

RADS(NGMX) 
for BEM problems (solve4.f) 

prob.elem_resid_1 
RE1(ns,nh) 
element residual array 

prob.elem_resid_2 
RE2(ns,nh) 
element residual array 

elem.jacobian 
RG(ng) 
is the Jacobian for a length, area, or volume integral (computed as RG and RGX in fe02/xgmg.f) 

RGS(NGMX) 
for BEM problems (solve4.f) 

RHO(NNMX,NJMX,2*NNMX) 
are the (rho,phi) coordinates of a 2D boundary element which has been subdivided with a local polar coordinate system. NNMIN is set up in solve4.f and is the local node at which the origin of the local polar coordinate scheme is located. BEM 

RI1,RI2,RI3 
are principal invariants of AZL 

RM 
is the modal matrix whose columns are the eigenvectors assoc with PST 

prob.global_resid_vector 
RP1(ny) 
global residual vector 

RWG(ng) 
is the square root of the determinant of the undeformed metric tensor multiplied by the quadratrure weight at gauss point ng (this variable isRGXin fe50/zeex50.f) 

S 

data.dof_opti_scale 
SCALE(no) 
scale associated with optimization DOF no for geometric fitting (see fe03/ipfit.f and fitgeo.f). Currently unused and should be removed. 

elem.scale_factor 
SE(ns,nb,ne) 
are scaling factors for basis nb of element ne. Arclength scale factors are set in fe02/lincal.f and fe02/dlse.f; face scaling factors are set in fe02/facseg.f 

SNLPA 
is static nonlinearity parameter “a” 

data.squared_dist 
SQ(nd) 
square of distance from mesh to data point nd 

SS 
total sum of squared distances SQ(nd) 

T 

TC 
are physical components of Cauchy stress 

TG 
are tensor components of 2nd PiolaKirchhoff stresses 

data.threshold 
THRES(3,ng,ne) 
used in threshold modelling activation pattern computations (fe30). Value at Gauss point ng of element ne is…


TN 
are physical components of Nominal stresses 

TV_SLO 
is slope of the force / velocity relation in stretching (before yield) 

U 

U 
is the right stretch tensor 

V 

elem.b_spline 
VE(ns,nk,ne) 
are Bspline polynomial coefficients for each DOF ns of element ne 

W 

data.weighting_factor 
WD(nj,nd) 
weighting factor for data point nd (this is gamma in Hashima et.al. 1993) 

fit.line_waiting_factor 
WDL(nj,nde) 
weighting factor for line data point nde 

elem.gaus_weight 
WG(ng,nb) 
are Gauss point weights at Gauss point ng for basis number nb 

WGF(ng,nf,nb) 
gaussian integration weights (as in WG) on face nf. [collocation] 

WS(5*NYMX) 
is unused and should be removed 

data.Sobolev_smoothing_weights 
WU(nu,ne) 
Sobolev smoothing weights for DOF nu (these are alpha & beta in Hashima et.al. 1993) 

X 

XB(1..2,NJMX,NLMX) 
is first/second Bezier slope control point for line NL (fe21/crherm.f) 

elem.xj_geom_pos_XE 
XE(ns,nj) 
is Xj geometric position or derivative for element DOF ns, coordinate nj 

XD(NKMX,NPDMX) 
for BEM problems (solve4.f) 

XDC(NKMX,NPDMX) 
for BEM problems (solve4.f) 

elem.xj_geom_pos_XG 
XG 
are undeformed theta coords and derivatives wrt Xi 

XG(nj,nu) 
is Xj geometric position (nu=1) or derivative (nu>1) for coord nj within an element 

XG1(NJMX,NUMX,NGMX) 
for BEM problems (solve4.f) 

XGS(NJMX,NUMX,NGMX) 
for BEM problems (solve4.f) 

data.xi_coor_data 
XID(ni,nd) 
Xicoordinate of data point nd 

fit.xi_of_line 
XIDL(ni,nde) 
Xicoordinate of line data point nde 

XIDR(ni,nd) 
is unused and should be removed 

elem.xi_gaus 
XIG(ni,ng,nb) 
is Xicoordinate (ni) at Gauss point ng for basis nb 

XIGF(ni,ng,nf,nb) 
gauss point locations (as in XIG) on face nf. [collocation] 

XMGF(nfg,ng,njmat) 
Material parameter njmat evaluated at gauss point ng of face nfg [collocation] 

XN(NJMX,NGMX) 
for BEM problems (solve4.f) 

XNS(NJMX,NGMX) 
for BEM problems (solve4.f) 

data.opti_var 
XO(no) 
optimization variable 

XOC(NCYMX) 
for BEM problems (solve4.f) 

nodes.global_param_XP 
XP(nk,nj,np) 
is Xj geometric position (nk=1) or derivative (nk>1) for coordinate nj at global node np 

XPD(NKMX,NJMX,NPDMX) 
for BEM problems (solve4.f) 

Y 

YD(NKMX) 
calculated BEM solution at the point XPFP (fe90/domsol.f) 

YDC(NKMX) 
calculated BEM solution at the point XPFP (fe90/domsol.f) 

YIELDR 
is ratio of yield tension to isometric tension 

data.data_at_guass 
YG(ng,nj,ne) 
is data defined at Gauss point ng (used in fe03/yger.f and fe03/fitgau.f)
– or – 

nodes.mesh_DOF 
YP(ny,1..16) 
is value of mesh DOF ny


Z 

data.rect_cart_coord_data 
ZD(nj,nd) 
rectangular cartesian coords of data point nd 

data.rect_cart_coord 
ZDL(nj,nde) 
rectangular cartesian coords of line data point nde 

data.deformed_param_array 
ZE(ns,nh) 
local deformed element parameter array [FJV 9 April 1996] 

ZEC(ns,nh) 
used for BEM problems. This is currently unused and should be removed. 

ZFC(ns,nh) 
used for BEM problems. This is currently unused and should be removed. 

data.deformed_theta_coor 
ZG(nh,nu) 
are deformed theta coordinates and derivatives wrt undeformed coordinates (computed in fe02/zezg.f):


ZK 
is unused and should be removed 

nodes.deform_global_param_ZP 
ZP(nk,nh,np) 
global deformed element parameter array [FJV 9 April 1996] 
Notes
 Material Nucoordinates (reference for stresses): are orthogonal and (Nu1,Nu2) lie in the (Xi1Xi2) plane such that Nu(1) is aligned with the `fibres’ to which material aeolotropy is referred; The undeformed base vectors are defined such that the undeformed metric tensors wrt the Nu are delta(i,j).

The fibre coord Nu(1) lies in the Xi1Xi2 plane Nu(2) is orthog to the fibre coord & lies in the Xi1Xi2 plane & the remaining Nu coord is orthog to this plane. The fibre angle eta is eta(1), the angle between the fibre coord & the Xi1 coord. The Nu coords are stress coords & are orthonormal with metric a(i,j)=Kronecker delta.
 KTYP26 = 3 KTYP27 = 1 is used for stripe intersection calculation.

If KTYP57>1 a pressure b.c. is applied to the Xi(3)=0,1 faces of elements with NW(ne)=2,3, respectively, or both if NW(ne)=4.
 NHMX must be equal to NJMX for problems where the dependent variable array carries deformed coordinates.

NZMX must be > NOMX^{2 and NYMX}2 for some problems and large enough to take the fillin during sparse matrix solution in others.
 If JTYP2=1 NQE stores the global NK number of local DOFs ns for element ne.
For cylindrical symmetry (JTYP4=2):
if JTYP3=1
radius is nj=2 coord.
if JTYP3=2,3
radius is nj=1 coord.
if JTYP3=4
radius=focus*sinhX1*sinX2
 Dependent and Geometric Variables counterparts
Dep Var Fortran Name
Dep Var Python Name
Geometric Var Fortran Name
Geom Var Python Name
nh
nj
zp
xp
ze
xe
zg
xg
nbh
elem.basis_num
nbj
elem.basis_type_num
nhe
nje
nhp
njp
nht
njt
nkh
nkj
Obsolete
GB(no) 
is RHS vector for system of reduced global equations 
GS(nz) 
is reduced global solution matrix in 1D form 
INDEX_PLIN 
is the index for a particular polyline. 
INDEX_PLIN_TYPE(index_plin) 
indicates the type (piecewise linear / Bezier). 
ITYP4 
is Galerkin FEM / Direct bem / Indirect bem / Orthogonal collocation (replaced by KTYP2) 
ITYP5 
is static / time integration / modal / Fourier / Laplace / buckling (replaced by KTYP3) 
ITYP6 
is linear / nonlinear problem (replaced by KTYP4) 
LXI( 0,nl) 
replaced by NXI 
NCNP(nh,np) 
is the mapping between nc and nh, and np which is set up in IPEQUA or DEEQUA and possibly modified in DECORN. 
NKE(nk,nb,ne) 
is global derivative number of local derivative nk of element ne for basis nb 
NO_BRANCH_PATTERN(no_gen,inumber_of_branches) 
is number of elements which have inumber_of_branches at generation no_gen. 
NT_PLIN 
is the total number of polylines. 
NT_PLIN_SECTIONS(index_plin) 
is the no of sections to the polyline. 
PLIN_DATA(nj,no_point,index_plin) 
are the coordinates of each point of the polyline. 
RP2(ny) 
formerly passed to GENSOL, now unused. 
YNP(nk,nh,nc,np,0) 
is the global mapping between nk,nh,nc,np and ny or use in coupled problems (yet to be set up). 