FEPG实现动态边界的要点

元计算技术

real*8 function bound(r,t,j,nf)

implicit real*8 (a-h,o-z)

dimension r(2)

pi=3.14159265358979

bound=dsin(pi*t/10.0)*100.0d0

c
write(*,*) 'bound =',bound

return

end

r ：坐标
t：当前时间自由度的序号
nf
： gcn文件中bft后的场序号

geoyang

梁老师您好，请教一下动态温度变化的边界条件咋下啊？ 比如温度 T=15t+10，那在bft中应该如何修改，修改完以后咋编译啊？

元计算技术

可以参考附件的文件