请教一个温度在网格上分配的算法

半个月亮

如图所示，将单元所得温度值分配到1-9节点上，其中：
节点1379 分别为单元1-4温度值；节点2为单元1、2的平均值；节点4为单元1、3的平均值；节点6为单元2、4的平均值；节点8为单元3、4的平均值；
节点5为单元1234的平均值；
请教一下，下面的算法哪里有问题？



算法伪代码：

[Fortran] 纯文本查看 复制代码

tem = [6 12 24 36]！单元温度
NJoints = 6 ！节点数
NElems = 4 ！单元数
NJE = 4      ！每个单元节点数
Elems = [1 2 5 4; 
                2 3 6 5;
                4 5 8 7;
                5 6 9 8！单元节点构造
node_tem = zeros(NJoints,1）
do k = 1:NJoints
        t = 0;
        ne = 0;
     do iel = 1:NElems
         do jel =1： NJE 
             if  Elems(iel,jel) == k
                 ne = ne+1;
                 t = t + tem(iel);
             end
         enddo
     enddo
    node_tem = t/ne;
 enddo