24 ,RR,SINK,TZERO,MM,TERM1,TERM2,NOD )
33 implicit real(kind=kreal)(a-h,o-z)
35 integer(kind=kint) :: NOD(MM)
36 real(kind=kreal) :: xx(nn), yy(nn), zz(nn), temp(nn), term1(mm*mm), term2(mm)
37 real(kind=kreal) :: xg(2), wgt(2), h(2), hr(2)
39 data xg/-0.5773502691896, 0.5773502691896/
44 else if(ltype.EQ.2)
then
47 else if(ltype.EQ.3)
then
68 gx=gx+hr(i)*xx(nod(i))
69 gy=gy+hr(i)*yy(nod(i))
71 xsum = dsqrt( gx*gx+gy*gy )
72 tt=h(1)*temp(nod(1))+h(2)*temp(nod(2))
75 rrr=(t1+t2)*(t1**2+t2**2)*rr
76 wg=xsum*wgt(lx)*rrr*thick
77 term1(ic)=term1(ic)-wg*h(ip)*h(jp)
78 if( ip.EQ.jp ) term2(ip)=term2(ip)-wg*h(jp)*sink
88 ,RR,SINK,TZERO,MM,TERM1,TERM2,NOD )
97 implicit real(kind=kreal) (a - h, o - z)
99 integer(kind=kint) :: NOD(MM)
100 real(kind=kreal) :: xx(nn), yy(nn), zz(nn), temp(nn), term1(mm*mm), term2(mm)
101 real(kind=kreal) :: xg(3), wgt(3), h(3), hr(3)
105 xg(1) = -0.7745966692
108 wgt(1) = 0.5555555555
109 wgt(2) = 0.8888888888
110 wgt(3) = 0.5555555555
116 else if(ltype.EQ.2)
then
120 else if(ltype.EQ.3)
then
135 h(1)=-0.5*(1.0-ri)*ri
136 h(2)= 0.5*(1.0+ri)*ri
144 gx=gx+hr(i)*xx(nod(i))
145 gy=gy+hr(i)*yy(nod(i))
147 xsum = dsqrt( gx*gx+gy*gy )
148 tt=h(1)*temp(nod(1))+h(2)*temp(nod(2))+h(3)*temp(nod(3))
151 rrr=(t1+t2)*(t1**2+t2**2)*rr
152 wg=xsum*wgt(lx)*rrr*thick
153 term1(ic)=term1(ic)-wg*h(ip)*h(jp)
154 if( ip.EQ.jp ) term2(ip)=term2(ip)-wg*h(jp)*sink
164 ,RR,SINK,TZERO,MM,TERM1,TERM2,NOD )
174 implicit real(kind=kreal) (a - h, o - z)
176 integer(kind=kint) :: NOD(MM)
177 real(kind=kreal) :: xx(nn), yy(nn), zz(nn), temp(nn), term1(mm*mm), term2(mm)
178 real(kind=kreal) :: xg(2), wgt(2), h(2), hr(2)
180 data xg/-0.5773502691896, 0.5773502691896/
185 else if(ltype.EQ.2)
then
188 else if(ltype.EQ.3)
then
191 else if(ltype.EQ.4)
then
212 gx=gx+hr(i)*xx(nod(i))
213 gy=gy+hr(i)*yy(nod(i))
215 xsum = dsqrt( gx*gx+gy*gy )
216 tt=h(1)*temp(nod(1))+h(2)*temp(nod(2))
219 rrr=(t1+t2)*(t1**2+t2**2)*rr
220 wg=xsum*wgt(lx)*rrr*thick
221 term1(ic)=term1(ic)-wg*h(ip)*h(jp)
222 if( ip.EQ.jp ) term2(ip)=term2(ip)-wg*h(jp)*sink
232 ,RR,SINK,TZERO,MM,TERM1,TERM2,NOD )
242 implicit real(kind=kreal) (a - h, o - z)
244 integer(kind=kint) :: NOD(MM)
245 real(kind=kreal) :: xx(nn), yy(nn), zz(nn), temp(nn), term1(mm*mm), term2(mm)
246 real(kind=kreal) :: xg(3), wgt(3), h(3), hr(3)
250 xg(1) = -0.7745966692
253 wgt(1) = 0.5555555555
254 wgt(2) = 0.8888888888
255 wgt(3) = 0.5555555555
261 else if(ltype.EQ.2)
then
265 else if(ltype.EQ.3)
then
269 else if(ltype.EQ.4)
then
284 h(1)=-0.5*(1.0-ri)*ri
285 h(2)= 0.5*(1.0+ri)*ri
293 gx=gx+hr(i)*xx(nod(i))
294 gy=gy+hr(i)*yy(nod(i))
296 xsum = dsqrt( gx*gx+gy*gy )
297 tt=h(1)*temp(nod(1))+h(2)*temp(nod(2))+h(3)*temp(nod(3))
300 rrr=(t1+t2)*(t1**2+t2**2)*rr
301 wg=xsum*wgt(lx)*rrr*thick
302 term1(ic)=term1(ic)-wg*h(ip)*h(jp)
303 if( ip.EQ.jp ) term2(ip)=term2(ip)-wg*h(jp)*sink
313 ,RR,SINK,TZERO,MM,TERM1,TERM2,NOD )
323 implicit real(kind=kreal) (a - h, o - z)
325 integer(kind=kint) :: NOD(MM)
326 real(kind=kreal) :: xx(nn), yy(nn), zz(nn), temp(nn), term1(mm*mm), term2(mm)
328 if ( ltype.EQ.1 )
then
332 else if( ltype.EQ.2 )
then
336 else if( ltype.EQ.3 )
then
340 else if( ltype.EQ.4 )
then
345 ax = xx( nod(2) ) - xx( nod(1) )
346 ay = yy( nod(2) ) - yy( nod(1) )
347 az = zz( nod(2) ) - zz( nod(1) )
348 bx = xx( nod(3) ) - xx( nod(1) )
349 by = yy( nod(3) ) - yy( nod(1) )
350 bz = zz( nod(3) ) - zz( nod(1) )
351 aa = dsqrt( (ay*bz-az*by )**2+ (az*bx-ax*bz )**2+(ax*by-ay*bx )**2 )/6.0
353 tt = ( temp(nod(1)) + temp(nod(2)) + temp(nod(3)) ) / 3.0
356 rrr = ( t1 + t2 )*( t1**2 + t2**2 ) * rr
359 term2(ip) = - aa * rrr * sink
362 term1(ic) = - aa * rrr / 3.0
371 ,RR,SINK,TZERO,MM,TERM1,TERM2,NOD )
381 implicit real(kind=kreal) (a - h, o - z)
383 integer(kind=kint) :: NOD(MM)
384 real(kind=kreal) :: xx(nn), yy(nn), zz(nn), temp(nn), term1(mm*mm), term2(mm)
385 real(kind=kreal) :: xg(3), wgt(3), h(6), hl1(6), hl2(6), hl3(6)
389 xg(1) = -0.7745966692
392 wgt(1) = 0.5555555555
393 wgt(2) = 0.8888888888
394 wgt(3) = 0.5555555555
403 else if(ltype.EQ.2)
then
410 else if(ltype.EQ.3)
then
417 else if(ltype.EQ.4)
then
438 x1=0.5*(1.0-x2)*(xl1+1.0)
477 g1x=g1x+(hl1(i)-hl3(i))*xx(nod(i))
478 g1y=g1y+(hl1(i)-hl3(i))*yy(nod(i))
479 g1z=g1z+(hl1(i)-hl3(i))*zz(nod(i))
480 g2x=g2x+(hl2(i)-hl3(i))*xx(nod(i))
481 g2y=g2y+(hl2(i)-hl3(i))*yy(nod(i))
482 g2z=g2z+(hl2(i)-hl3(i))*zz(nod(i))
487 xsum=dsqrt(g3x**2+g3y**2+g3z**2)
508 tt=h(1)*temp(nod(1)) + h(2)*temp(nod(2)) &
509 +h(3)*temp(nod(3)) + h(4)*temp(nod(4)) &
510 +h(5)*temp(nod(5)) + h(6)*temp(nod(6))
513 rrr=( t1 + t2 )*( t1**2 + t2**2 )*rr
514 wg=wgt(l1)*wgt(l2)*det*(1.0-x2)*0.25*rrr
515 term1(ic)=term1(ic)-wg*h(ip)*h(jp)
516 if( ip.EQ.jp ) term2(ip)=term2(ip)-wg*h(jp)*sink
528 ,RR,SINK,TZERO,MM,TERM1,TERM2,NOD )
539 implicit real(kind=kreal) (a - h, o - z)
541 integer(kind=kint) :: NOD(MM)
542 real(kind=kreal) :: xx(nn), yy(nn), zz(nn), temp(nn), term1(mm*mm), term2(mm)
544 real(kind=kreal) :: h(4), hr(4), hs(4), ht(4)
545 real(kind=kreal) :: xg(2), wgt(2)
549 data xg/-0.5773502691896,0.5773502691896/
553 if( ltype.EQ.1 )
then
558 else if( ltype.EQ.2 )
then
563 else if( ltype.EQ.3 )
then
569 else if( ltype.EQ.4 )
then
575 else if( ltype.EQ.5 )
then
583 ax = xx( nod(2) ) - xx( nod(1) )
584 ay = yy( nod(2) ) - yy( nod(1) )
585 az = zz( nod(2) ) - zz( nod(1) )
586 bx = xx( nod(3) ) - xx( nod(1) )
587 by = yy( nod(3) ) - yy( nod(1) )
588 bz = zz( nod(3) ) - zz( nod(1) )
589 aa = dsqrt((ay*bz-az*by)**2+(az*bx-ax*bz)**2+(ax*by-ay*bx)**2)/6.0
591 tt = ( temp(nod(1)) + temp(nod(2)) + temp(nod(3)) ) / 3.0
594 rrr = ( t1 + t2 )*( t1**2 + t2**2 ) * rr
597 term2(ip) = - aa * rrr * sink
600 term1(ic) = - aa * rrr / 3.0
603 elseif( isuf.EQ.2 )
then
614 h(1)=0.25*(1.0-ri)*(1.0-si)
615 h(2)=0.25*(1.0+ri)*(1.0-si)
616 h(3)=0.25*(1.0+ri)*(1.0+si)
617 h(4)=0.25*(1.0-ri)*(1.0+si)
633 g1x=g1x+hr(i)*xx(nod(i))
634 g1y=g1y+hr(i)*yy(nod(i))
635 g1z=g1z+hr(i)*zz(nod(i))
636 g2x=g2x+hs(i)*xx(nod(i))
637 g2y=g2y+hs(i)*yy(nod(i))
638 g2z=g2z+hs(i)*zz(nod(i))
643 xsum=dsqrt(g3x**2+g3y**2+g3z**2)
644 tt=h(1)*temp(nod(1))+h(2)*temp(nod(2)) &
645 +h(3)*temp(nod(3))+h(4)*temp(nod(4))
648 rrr=(t1+t2)*(t1**2+t2**2)*rr
649 wg=xsum*wgt(ig1)*wgt(ig2)*rrr
650 term1(ic)=term1(ic)-wg*h(ip)*h(jp)
651 if( ip.EQ.jp ) term2(ip)=term2(ip)-wg*h(jp)*sink
662 ,RR,SINK,TZERO,MM,TERM1,TERM2,NOD )
674 implicit real(kind=kreal) (a - h, o - z)
676 integer(kind=kint) :: NOD(MM)
677 real(kind=kreal) :: xx(nn), yy(nn), zz(nn), temp(nn), term1(mm*mm), term2(mm)
679 real(kind=kreal) :: h(8), hr(8), hs(8), ht(8)
680 real(kind=kreal) :: xg(3), wgt(3), hl1(6), hl2(6), hl3(6)
684 xg(1) = -0.7745966692
687 wgt(1) = 0.5555555555
688 wgt(2) = 0.8888888888
689 wgt(3) = 0.5555555555
692 if( ltype.EQ.1 )
then
700 else if( ltype.EQ.2 )
then
708 else if( ltype.EQ.3 )
then
718 else if( ltype.EQ.4 )
then
728 else if( ltype.EQ.5 )
then
753 x1=0.5*(1.0-x2)*(xl1+1.0)
792 g1x=g1x+(hl1(i)-hl3(i))*xx(nod(i))
793 g1y=g1y+(hl1(i)-hl3(i))*yy(nod(i))
794 g1z=g1z+(hl1(i)-hl3(i))*zz(nod(i))
795 g2x=g2x+(hl2(i)-hl3(i))*xx(nod(i))
796 g2y=g2y+(hl2(i)-hl3(i))*yy(nod(i))
797 g2z=g2z+(hl2(i)-hl3(i))*zz(nod(i))
802 xsum=dsqrt(g3x**2+g3y**2+g3z**2)
823 tt=h(1)*temp(nod(1)) + h(2)*temp(nod(2)) &
824 +h(3)*temp(nod(3)) + h(4)*temp(nod(4)) &
825 +h(5)*temp(nod(5)) + h(6)*temp(nod(6))
828 rrr=( t1 + t2 )*( t1**2 + t2**2 )*rr
829 wg=wgt(l1)*wgt(l2)*det*(1.0-x2)*0.25*rrr
830 term1(ic)=term1(ic)-wg*h(ip)*h(jp)
831 if( ip.EQ.jp ) term2(ip)=term2(ip)-wg*h(jp)*sink
837 elseif( isuf.EQ.2 )
then
852 h(1)=0.25*rm*sm*(-1.0-ri-si)
853 h(2)=0.25*rp*sm*(-1.0+ri-si)
854 h(3)=0.25*rp*sp*(-1.0+ri+si)
855 h(4)=0.25*rm*sp*(-1.0-ri+si)
856 h(5)=0.5*(1.0-ri*ri)*(1.0-si)
857 h(6)=0.5*(1.0-si*si)*(1.0+ri)
858 h(7)=0.5*(1.0-ri*ri)*(1.0+si)
859 h(8)=0.5*(1.0-si*si)*(1.0-ri)
860 hr(1)=-.25*sm*(-1.0-ri-si)-0.25*rm*sm
861 hr(2)= .25*sm*(-1.0+ri-si)+0.25*rp*sm
862 hr(3)= .25*sp*(-1.0+ri+si)+0.25*rp*sp
863 hr(4)=-.25*sp*(-1.0-ri+si)-0.25*rm*sp
865 hr(6)= .5*(1.0-si*si)
867 hr(8)=-.5*(1.0-si*si)
868 hs(1)=-.25*rm*(-1.0-ri-si)-0.25*rm*sm
869 hs(2)=-.25*rp*(-1.0+ri-si)-0.25*rp*sm
870 hs(3)= .25*rp*(-1.0+ri+si)+0.25*rp*sp
871 hs(4)= .25*rm*(-1.0-ri+si)+0.25*rm*sp
872 hs(5)=-.5*(1.0-ri*ri)
874 hs(7)= .5*(1.0-ri*ri)
883 g1x=g1x+hr(i)*xx(nod(i))
884 g1y=g1y+hr(i)*yy(nod(i))
885 g1z=g1z+hr(i)*zz(nod(i))
886 g2x=g2x+hs(i)*xx(nod(i))
887 g2y=g2y+hs(i)*yy(nod(i))
888 g2z=g2z+hs(i)*zz(nod(i))
893 xsum=dsqrt(g3x**2+g3y**2+g3z**2)
894 tt=h(1)*temp(nod(1)) + h(2)*temp(nod(2)) &
895 +h(3)*temp(nod(3)) + h(4)*temp(nod(4)) &
896 +h(5)*temp(nod(5)) + h(6)*temp(nod(6)) &
897 +h(7)*temp(nod(7)) + h(8)*temp(nod(8))
900 rrr=( t1 + t2 )*( t1**2 + t2**2 )*rr
901 wg=xsum*wgt(ig1)*wgt(ig2)*rrr
902 term1(ic)=term1(ic)-wg*h(ip)*h(jp)
903 if( ip.EQ.jp ) term2(ip)=term2(ip)-wg*h(jp)*sink
914 ,RR,SINK,TZERO,MM,TERM1,TERM2,NOD )
926 implicit real(kind=kreal) (a - h, o - z)
928 integer(kind=kint) :: NOD(MM)
929 real(kind=kreal) :: xx(nn), yy(nn), zz(nn), temp(nn), term1(mm*mm), term2(mm)
931 real(kind=kreal) :: h(4), hr(4), hs(4), ht(4)
932 real(kind=kreal) :: xg(2), wgt(2)
936 data xg/-0.5773502691896,0.5773502691896/
939 if( ltype.EQ.1 )
then
944 else if( ltype.EQ.2 )
then
949 else if( ltype.EQ.3 )
then
954 else if( ltype.EQ.4 )
then
959 else if( ltype.EQ.5 )
then
964 else if( ltype.EQ.6 )
then
981 h(1)=0.25*(1.0-ri)*(1.0-si)
982 h(2)=0.25*(1.0+ri)*(1.0-si)
983 h(3)=0.25*(1.0+ri)*(1.0+si)
984 h(4)=0.25*(1.0-ri)*(1.0+si)
1000 g1x=g1x+hr(i)*xx(nod(i))
1001 g1y=g1y+hr(i)*yy(nod(i))
1002 g1z=g1z+hr(i)*zz(nod(i))
1003 g2x=g2x+hs(i)*xx(nod(i))
1004 g2y=g2y+hs(i)*yy(nod(i))
1005 g2z=g2z+hs(i)*zz(nod(i))
1010 xsum=dsqrt(g3x**2+g3y**2+g3z**2)
1011 tt=h(1)*temp(nod(1)) + h(2)*temp(nod(2)) &
1012 +h(3)*temp(nod(3)) + h(4)*temp(nod(4))
1015 rrr=( t1 + t2 )*( t1**2 + t2**2 )*rr
1016 wg=xsum*wgt(ig1)*wgt(ig2)*rrr
1017 term1(ic)=term1(ic)-wg*h(ip)*h(jp)
1018 if( ip.EQ.jp ) term2(ip)=term2(ip)-wg*h(jp)*sink
1028 ,RR,SINK,TZERO,MM,TERM1,TERM2,NOD )
1040 implicit real(kind=kreal) (a - h, o - z)
1042 integer(kind=kint) :: NOD(MM)
1043 real(kind=kreal) :: xx(nn), yy(nn), zz(nn), temp(nn), term1(mm*mm), term2(mm)
1045 real(kind=kreal) :: h(8), hr(8), hs(8), ht(8)
1046 real(kind=kreal) :: xg(3), wgt(3)
1050 xg(1) = -0.7745966692
1052 xg(3) = 0.7745966692
1053 wgt(1) = 0.5555555555
1054 wgt(2) = 0.8888888888
1055 wgt(3) = 0.5555555555
1057 if( ltype.EQ.1 )
then
1066 else if( ltype.EQ.2 )
then
1075 else if( ltype.EQ.3 )
then
1084 else if( ltype.EQ.4 )
then
1093 else if( ltype.EQ.5 )
then
1102 else if( ltype.EQ.6 )
then
1127 h(1)=0.25*rm*sm*(-1.0-ri-si)
1128 h(2)=0.25*rp*sm*(-1.0+ri-si)
1129 h(3)=0.25*rp*sp*(-1.0+ri+si)
1130 h(4)=0.25*rm*sp*(-1.0-ri+si)
1131 h(5)=0.5*(1.0-ri*ri)*(1.0-si)
1132 h(6)=0.5*(1.0-si*si)*(1.0+ri)
1133 h(7)=0.5*(1.0-ri*ri)*(1.0+si)
1134 h(8)=0.5*(1.0-si*si)*(1.0-ri)
1135 hr(1)=-.25*sm*(-1.0-ri-si)-0.25*rm*sm
1136 hr(2)= .25*sm*(-1.0+ri-si)+0.25*rp*sm
1137 hr(3)= .25*sp*(-1.0+ri+si)+0.25*rp*sp
1138 hr(4)=-.25*sp*(-1.0-ri+si)-0.25*rm*sp
1140 hr(6)= .5*(1.0-si*si)
1142 hr(8)=-.5*(1.0-si*si)
1143 hs(1)=-.25*rm*(-1.0-ri-si)-0.25*rm*sm
1144 hs(2)=-.25*rp*(-1.0+ri-si)-0.25*rp*sm
1145 hs(3)= .25*rp*(-1.0+ri+si)+0.25*rp*sp
1146 hs(4)= .25*rm*(-1.0-ri+si)+0.25*rm*sp
1147 hs(5)=-.5*(1.0-ri*ri)
1149 hs(7)= .5*(1.0-ri*ri)
1158 g1x=g1x+hr(i)*xx(nod(i))
1159 g1y=g1y+hr(i)*yy(nod(i))
1160 g1z=g1z+hr(i)*zz(nod(i))
1161 g2x=g2x+hs(i)*xx(nod(i))
1162 g2y=g2y+hs(i)*yy(nod(i))
1163 g2z=g2z+hs(i)*zz(nod(i))
1168 xsum=dsqrt(g3x**2+g3y**2+g3z**2)
1169 tt=h(1)*temp(nod(1)) + h(2)*temp(nod(2)) &
1170 +h(3)*temp(nod(3)) + h(4)*temp(nod(4)) &
1171 +h(5)*temp(nod(5)) + h(6)*temp(nod(6)) &
1172 +h(7)*temp(nod(7)) + h(8)*temp(nod(8))
1175 rrr=( t1 + t2 )*( t1**2 + t2**2 )*rr
1176 wg=xsum*wgt(ig1)*wgt(ig2)*rrr
1177 term1(ic)=term1(ic)-wg*h(ip)*h(jp)
1178 if( ip.EQ.jp ) term2(ip)=term2(ip)-wg*h(jp)*sink
1188 ,RR,SINK,TZERO,TERM1,TERM2 )
1195 implicit real(kind=kreal) (a - h, o - z)
1197 real(kind=kreal) :: xx(nn), yy(nn), zz(nn), temp(nn), term1(nn*nn), term2(nn)
1205 aa=dsqrt((ay*bz-az*by)**2+(az*bx-ax*bz)**2+(ax*by-ay*bx)**2)/6.0
1207 tt = ( temp(1) + temp(2) + temp(3) ) / 3.0
1210 rrr = ( t1 + t2 )*( t1**2 + t2**2 ) * rr
1214 term2(ip) = - aa * rrr * sink
1217 term1(ic) = - aa * rrr / 3.0
1226 ,RR,SINK,TZERO,TERM1,TERM2 )
1233 implicit real(kind=kreal) (a - h, o - z)
1235 real(kind=kreal) :: xx(nn), yy(nn), zz(nn), temp(nn), term1(nn*nn), term2(nn)
1237 real(kind=kreal) :: h(4), hr(4), hs(4), ht(4)
1238 real(kind=kreal) :: xg(2), wgt(2)
1242 data xg/-0.5773502691896,0.5773502691896/
1255 h(1)=0.25*(1.0-ri)*(1.0-si)
1256 h(2)=0.25*(1.0+ri)*(1.0-si)
1257 h(3)=0.25*(1.0+ri)*(1.0+si)
1258 h(4)=0.25*(1.0-ri)*(1.0+si)
1284 xsum=dsqrt(g3x**2+g3y**2+g3z**2)
1286 tt = h(1)*temp(1) + h(2)*temp(2)+ h(3)*temp(3) + h(4)*temp(4)
1289 rrr = ( t1 + t2 ) * ( t1**2 + t2**2 ) * rr
1291 term1(ic) = term1(ic) - xsum*wgt(ig1)*wgt(ig2)*h(ip)*h(jp)*rrr
1293 term2(ip) = term2(ip) - xsum*wgt(ig1)*wgt(ig2)*h(jp)*sink*rrr
This module provides subroutines to generate heat radiate boundary.
subroutine heat_radiate_342(nn, xx, yy, zz, temp, ltype, rr, sink, tzero, mm, term1, term2, nod)
subroutine heat_radiate_361(nn, xx, yy, zz, temp, ltype, rr, sink, tzero, mm, term1, term2, nod)
subroutine heat_radiate_242(nn, xx, yy, zz, thick, temp, ltype, rr, sink, tzero, mm, term1, term2, nod)
subroutine heat_radiate_352(nn, xx, yy, zz, temp, ltype, rr, sink, tzero, mm, term1, term2, nod)
subroutine heat_radiate_231(nn, xx, yy, zz, thick, temp, ltype, rr, sink, tzero, mm, term1, term2, nod)
subroutine heat_radiate_351(nn, xx, yy, zz, temp, ltype, rr, sink, tzero, mm, term1, term2, nod)
subroutine heat_radiate_731(nn, xx, yy, zz, temp, ltype, rr, sink, tzero, term1, term2)
subroutine heat_radiate_341(nn, xx, yy, zz, temp, ltype, rr, sink, tzero, mm, term1, term2, nod)
subroutine heat_radiate_741(nn, xx, yy, zz, temp, ltype, rr, sink, tzero, term1, term2)
subroutine heat_radiate_241(nn, xx, yy, zz, thick, temp, ltype, rr, sink, tzero, mm, term1, term2, nod)
subroutine heat_radiate_362(nn, xx, yy, zz, temp, ltype, rr, sink, tzero, mm, term1, term2, nod)
subroutine heat_radiate_232(nn, xx, yy, zz, thick, temp, ltype, rr, sink, tzero, mm, term1, term2, nod)