15 type(hecmwst_matrix) :: hecMAT
16 integer(kind=kint) :: ib, ii, id
17 real(kind=kreal) :: ctime, dtime, qq, beta
20 do ib = 1, fstrheat%T_FIX_tot
21 ii = fstrheat%T_FIX_node(ib)
22 id = fstrheat%T_FIX_ampl(ib)
27 call hecmw_mat_ass_bc(hecmat, ii, 1, fstrheat%T_FIX_VAL(ib)*qq)