! 得到MSC.MARC的单元坐标
! by 陆新征，清华大学

module My_Var !公共数据库，用于存放公共数据
integer Elem_Nodal_No(7000,4) ! 单元节点编号, 最多3000个单元， 4个节点
real*8 Elem_Nodal_Coord(7000,4,3) ! 单元节点坐标, 最多3000个单元， 4个节点，3个坐标分量
end module

SUBROUTINE UELOOP(M,N,IL) ! 每个单元的循环，得到单元节点编号和坐标
use My_Var
IMPLICIT REAL *8 (A-H, O-Z)
include "../common/blnk"
real*8 ValNo(3)
ICOD=0; NQNCOMP=3; NQDATATYPE=0;
do I=1, 4;
Elem_Nodal_No(M, I)=IELEXT(LM(I));
NODEID=Elem_Nodal_No(M, I);
if(NodeID.ne.0) CALL NODVAR(ICOD,NODEID,VALNO,NQNCOMP,NQDATATYPE);
Elem_Nodal_Coord(M,I,:)=ValNo(:)
end do

RETURN
END