对流换热子程序film

在这个子程序中，可以设置对流以及辐射换热。我们知道，对流和辐射的系数都是与温度相关的。尤其是辐射的系数，与温度的三次方相关，因此当温度较高时，工件的换热主要通过辐射的方式。具体的公式大家找一本传热学的书看一下就知道了。
今天介绍的子程序比较简单，类似于一个水火弯板或者随焊激冷过程，对流换热系数认为是一个常数。代码如下：
subroutine film(h,tinf,ts,n,time)
implicit logical (a-z)
real*8 h
integer n
real*8 time, tinf, ts
dimension h(2),n(7),ts(6)
tinf=20
if (ts(4).gt.20.and.ts(4).lt.30.and.ts(5).gt.(4*time)
$ .and.ts(5).lt.(4*time+10)) then
h(1)=10
else
h(1)=0
end if
return
end
解释一下，模型中，工件的初始温度设置为1000，温度为20度的冷源沿着y轴移动，移动速度为4mm/s，冷源的范围为x方向10mm，y方向10mm。对流换热系数为10000（随便设置的）。模拟结果如图中所示。

模型及子程序见附件。