程序是通过质量守恒和能量守恒来建立模型，结合人体节点模型来得到人体的温度变化
之前的能量守恒方程是这样的
ar((I-1)*3+3)=(-CV(I,0))*TA(I,0)
-tt(i,0)+hct*(3*porosity_m(i)/rm)*(ta(i,0)-tm(i,16,0))*deltat
我要修改的部分是在这个方程后面加一项辐射热
添加的程序是这几句
qrad=5.67e-8*1.0*((60.+273.)**4-ta(npoin,0)**4)*(1.0e-4)
gama=-log(0.01)/x(npoin)
和 +(qrad*exp(-gama*xc))*deltat，这一项是加在能量守恒方程后面的
之前没加这几句的原始程序是不会报错的，加了以后不但运算中会报错而且出来的数据也不对
本来应该一直升高的温度也变成了波动的
我没学过fortran，眼看着毕设要交了，程序还改不明白，请吧友帮忙看看这几句程序有什么不对的地方，多谢各位啦
之前的能量守恒方程是这样的
ar((I-1)*3+3)=(-CV(I,0))*TA(I,0)
-tt(i,0)+hct*(3*porosity_m(i)/rm)*(ta(i,0)-tm(i,16,0))*deltat
我要修改的部分是在这个方程后面加一项辐射热
添加的程序是这几句
qrad=5.67e-8*1.0*((60.+273.)**4-ta(npoin,0)**4)*(1.0e-4)
gama=-log(0.01)/x(npoin)
和 +(qrad*exp(-gama*xc))*deltat，这一项是加在能量守恒方程后面的
之前没加这几句的原始程序是不会报错的，加了以后不但运算中会报错而且出来的数据也不对
本来应该一直升高的温度也变成了波动的
我没学过fortran，眼看着毕设要交了，程序还改不明白，请吧友帮忙看看这几句程序有什么不对的地方，多谢各位啦
