Sub-terahertz (THz) communication is envisioned as one of the key components for 6G because of the abundantly available spectrum resource. Accurate and efficient channel models are prerequisites for developing sub-THz communication systems. Due to the sparsity and more ray optics propagation characteristics of the sub-THz channel, deterministic Ray-Tracing (RT) has attracted much attention for sub-THz channel modeling, which shows the potential of reducing the simulation complexity yet maintaining the accuracy. This paper presents an implementation of RT for sub-THz channel modeling and demonstrates its performance based on sub-THz channel measurements. A virtual massive multiple-inputmultiple-output (MIMO) channel operating at 100 GHz and a double-directional 300 GHz channel are considered in the RT implementation, where the RT achieves a high similarity compared to the channel measurements in terms of channel impulse response and power angular spectrum. Besides, the near-field and spatial non-stationary properties of the sub-THz massive MIMO channel and the dominant multipaths of the 300 GHz channel are accurately reconstructed in the RT simulation. This work can provide insights into deterministic sub-THz channel modeling research from the implementation, evaluation, and challenges perspectives.