最近看到一段程序，但是不能理解，帮忙！

cbwang

原文如下：

用差值计算的方法求得节点的热流密度，具体差值计算方法是：

*do,i,1,nodesum-1

nend=nnear(nstart)

*get,ys,node,nstart,loc,y
*get,ye,node,nend,loc,y
*get,xs,node,nstart,loc,x
*get,xe,node,nend,loc,x

s=2*3.1415926*(xs+xe)/2*0.01   面积

j=ys*1000-mod(ys*1000,1)+1

*if,ys*1000,eq,j-1,then

Tis=Ti(j)

*else

Tis=((ys-(j-1)/1000)*Ti(j+1)+(j/1000-ys)*Ti(j))/(1/1000)   节点S的电流

*endif

j=ye*1000-mod(ye*1000,1)+1

*if,ye*1000,eq,j-1,then

Tie=Ti(j)

*else

Tie=((ye-(j-1)/1000)*Ti(j+1)+(j/1000-ye)*Ti(j))/(1/1000)   节点E的电流

*endif

THF=(Tie-Tis)*3300*0.03/s        电流*电压*系数/面积

FLST,2,2,1,ORDE,2
FITEM,2,nstart
FITEM,2,nend
SF,P51X,HFLUX,THF
nstart=nend

*enddo