call Heat (integer(IN) :: blockCount, integer(IN) :: blockList(blockCount), real(IN) :: dt, real(IN) :: time)
Apply the stat+gauss source term operator to a block of zones. The energy generation rate is used to update the internal energy in the zone. The phonomenological heating rate is described as a 3-D Gauss function. After we call stat+gauss, call the eos to update the pressure and temperature based on the phenomenological heating.
blockCount : number of blocks to operate on blockList : list of blocks to operate on dt : current timestep time : current time