diff --git a/tmd/utils/trace.py b/tmd/utils/trace.py
index 67507932..ddad925a 100644
--- a/tmd/utils/trace.py
+++ b/tmd/utils/trace.py
@@ -33,18 +33,38 @@ def trace1(loc: str, vdf: pd.DataFrame) -> None:
     wght_min = wght.min()
     wght_max = wght.max()
     wght_results = (
-        f">{loc} weights all,puf,cps (#M)= "
+        f">{loc} WEIGHT all,puf,cps (#M)= "
         f"{wtot:.3f} {wpuf:.3f} [160.8] {wcps:.3f}"
     )
     print(wght_results)
-    print(f">{loc} weights all_min,all_max (#)= {wght_min:.1f} {wght_max:.1f}")
+    print(f">{loc} WEIGHT all_min,all_max= {wght_min:.1f} {wght_max:.1f}")
+    # population tabulations
+    people = vdf.XTOT
+    pop_tot = (wght * people).sum() * 1e-6
+    pop_puf = (wght * filer * people).sum() * 1e-6
+    pop_cps = (wght * ~filer * people).sum() * 1e-6
+    pop_results = (
+        f">{loc} POPULATION all,puf,cps (#M)= "
+        f"{pop_tot:.3f} [334.181] {pop_puf:.3f} {pop_cps:.3f}"
+    )
+    print(pop_results)
+    # pencon tabulations
+    pencon = vdf.pencon_p + vdf.pencon_s
+    pencon_amt = (wght * filer * pencon).sum() * 1e-9
+    pencon_num = (wght * filer * (pencon > 0)).sum() * 1e-6
+    print(f">{loc} weighted puf PENCON ($B)= {pencon_amt:.3f}")
+    print(f">{loc} weighted puf PENCON>0 (#M)= {pencon_num:.3f}")
+    pencon_amt = (wght * ~filer * pencon).sum() * 1e-9
+    pencon_num = (wght * ~filer * (pencon > 0)).sum() * 1e-6
+    print(f">{loc} weighted cps PENCON ($B)= {pencon_amt:.3f}")
+    print(f">{loc} weighted cps PENCON>0 (#M)= {pencon_num:.3f}")
     # CTC tabulations
     if "ctc_total" in vdf:
         ctc = vdf.ctc_total * filer
         ctc_amt = (wght * ctc).sum() * 1e-9
         ctc_num = (wght * (ctc > 0)).sum() * 1e-6
         print(f">{loc} weighted puf CTC ($B)= {ctc_amt:.3f} [124.6]")
-        print(f">{loc} weighted puf CTC (#M)= {ctc_num:.3f} [36.5...47.4]")
+        print(f">{loc} weighted puf CTC>0 (#M)= {ctc_num:.3f} [36.5...47.4]")
     else:
         print(f">{loc} CTC not in DataFrame")
     # SALT tabulations
@@ -52,14 +72,14 @@ def trace1(loc: str, vdf: pd.DataFrame) -> None:
     salt_amt = (wght * salt).sum() * 1e-9
     salt_num = (wght * (salt > 0)).sum() * 1e-6
     print(f">{loc} weighted puf SALT ($B)= {salt_amt:.3f} [?]")
-    print(f">{loc} weighted puf SALT (#M)= {salt_num:.3f} [14.3...27.1]")
+    print(f">{loc} weighted puf SALT>0 (#M)= {salt_num:.3f} [14.3...27.1]")
     # PT_binc_w2_wages tabulations
     w2wages = vdf.PT_binc_w2_wages * filer
     wages_min = w2wages.min()
     wages_max = w2wages.max()
     wages_wtot = (wght * w2wages).sum() * 1e-9
-    print(f">{loc} W2_wages min,max ($)= {wages_min:.0f} {wages_max:.0f}")
-    print(f">{loc} weighted puf W2_wages ($B)= {wages_wtot:.3f}")
+    print(f">{loc} W2_WAGES min,max ($)= {wages_min:.0f} {wages_max:.0f}")
+    print(f">{loc} weighted puf W2_WAGES ($B)= {wages_wtot:.3f}")
     # QBID tabulations
     if "qbided" in vdf:
         qbid = vdf.qbided * filer
@@ -67,3 +87,18 @@ def trace1(loc: str, vdf: pd.DataFrame) -> None:
         print(f">{loc} weighted puf QBID ($B)= {qbid_wtot:.3f} [205.8]")
     else:
         print(f">{loc} QBID not in DataFrame")
+    # IITAX tabulations
+    if "iitax" in vdf:
+        wiitax = wght * vdf.iitax
+        itax_tot = (wiitax).sum() * 1e-9
+        itax_puf = (wiitax * filer).sum() * 1e-9
+        itax_cps = (wiitax * ~filer).sum() * 1e-9
+        itax_results = (
+            f">{loc} weighted all,puf,cps IITAX ($B)= "
+            f"{itax_tot:.3f} {itax_puf:.3f} {itax_cps:.3f}"
+        )
+        print(itax_results)
+        pos_itax_cps = (wght * (vdf.iitax > 0) * ~filer).sum() * 1e-6
+        print(f">{loc} weighted cps IITAX>0 (#M)= {pos_itax_cps:.3f}")
+    else:
+        print(f">{loc} IITAX not in DataFrame")