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

珊瑚虫

建议参照有限元后处理应力磨平操作进行计算，即对每个单元进行循环，循环的过程中累积节点的温度，并且记录每个节点所在的单元数，最后用累积温度除以每个节点所在的单元数即可：
do ie=1,ne       ！单元总数循环
do inode=1,nd  !每个单元的节点数
T(ELEMS(inode.ie))=T(ELEMS(inode.ie))+每个单元计算出来的节点温度  ！累积单元节点温度
nn(ELEMS(inode.ie))=nn(ELEMS(inode.ie))+1                                    ！累积每个节点所在的单元数
enddo
enddo
! 再平均
T(1:NNODE)=T(1:NNODE)/NN(1:NNODE)