20 integer(kind=kint) :: ib, in, ia
21 real(kind=kreal) :: ctime, dtime, qq, beta
23 type(hecmwst_matrix) :: hecMAT
25 do ib = 1, fstrheat%Q_NOD_tot
26 in = fstrheat%Q_NOD_node(ib)
27 ia = fstrheat%Q_NOD_ampl(ib)
29 hecmat%B(in) = hecmat%B(in) + qq*fstrheat%Q_NOD_val(ib)