哪位大侠帮我看看这段命令流

首先看看有没有错误，再帮忙解释一下：yugg（i），ycgg（i）究竟指什么？tzx、idx_up,idy_up是怎么算的？谢谢各位啦
!对于公式T(Z，X)=T(Z，0)+DATAT(X/Xa)**3，其中t0为设计变量，数组T为各层的拱冠梁厚度，由公式t=t0+Z*13.9/90得到
!datal为公式中的DATAT的左边数组，datal为公式中的DATAT的左边数组，xtal为公式中的红XA的左边数组
!xtar为公式中的DATAT的右边数组
*set,t0,5
*dim,t,array,layer_num
*dim,yugg,array,layer_num
*dim,ycgg,array,layer_num
*dim,tzx,array,11,layer_num !各层各关键点处的厚度值
*do,i,1,layer_num,1
*do,j,1,11,1
t(i)=t0+idz(i)*13.9/90
yugg(i)=idz(i)*(0.315*(idz(i)/90)-0.378) !拱冠梁上游面Y方向坐标,也表示优化过程中拱冠梁上游面方程不变
ycgg(i)=yugg(i)+0.5*t(i) !yugg为各层拱轴线拱冠处的Y坐标，t(i)为各层拱冠梁厚度，ycgg也为拱轴线y=x**2/2*rc+ycgg
*if,j,lt,6,then
tzx(j,i)=t(i)+datatl(i)*((abs(xc(j,i)))/xatl(i))**3 !拱厚变化曲线方程T(z,x)=T(z,0)+DataT(x/xA)**3
idx_up(j,i)=xc(j,i)+0.5*tzx(j,i)*sin(atan(xc(j,i)/rclz(i)))
idy_up(j,i)=(xc(j,i))**2/(2*rclz(i))+ycgg(i)-0.5*tzx(j,i)*(cos(atan(xc(j,i)/rclz(i))))
idx_down(j,i)=xc(j,i)-0.5*tzx(j,i)*sin(atan(xc(j,i)/rclz(i)))
idy_down(j,i)=(xc(j,i))**2/(2*rclz(i))+ycgg(i)+0.5*tzx(j,i)*(cos(atan(xc(j,i)/rclz(i))))
*elseif,j,eq,6,then
tzx(j,i)=t(i)
idx_up(j,i)=xc(j,i)
idy_up(j,i)=(xc(j,i))**2/(2*rclz(i))+ycgg(i)-0.5*tzx(j,i)
idx_down(j,i)=xc(j,i)
idy_down(j,i)=(xc(j,i))**2/(2*rclz(i))+ycgg(i)+0.5*tzx(j,i)
*else
tzx(j,i)=t(i)+datatr(i)*((abs(xc(j,i))/xatr(i)))**3
idx_up(j,i)=xc(j,i)+0.5*tzx(j,i)*sin(atan(xc(j,i)/rcrz(i)))
idy_up(j,i)=(xc(j,i))**2/(2*rcrz(i))+ycgg(i)-0.5*tzx(j,i)*(cos(atan(xc(j,i)/rcrz(i))))
idx_down(j,i)=xc(j,i)-0.5*tzx(j,i)*sin(atan(xc(j,i)/rcrz(i)))
idy_down(j,i)=(xc(j,i))**2/(2*rcrz(i))+ycgg(i)+0.5*tzx(j,i)*(cos(atan(xc(j,i)/rcrz(i))))
*endif
*enddo
*enddo