diff --git a/Python/COARE3.6/coare36vnWarm_et.py b/Python/COARE3.6/coare36vnWarm_et.py
index ee31a18..b9aae07 100644
--- a/Python/COARE3.6/coare36vnWarm_et.py
+++ b/Python/COARE3.6/coare36vnWarm_et.py
@@ -431,7 +431,7 @@ def ComputeEsat(T = None,P = None):
     fnameA = os.path.join(path,'test_36_output_py_082022_withwavesinput_withwarmlayer.txt')
     A=coare36vnWarm_et(Jd, U, Zu, Tair, Zt, RH, Zq, P, Tsg, SW_dn, LW_dn, Lat, Lon, Zi, Rainrate, Ts_depth, Ss, None, None, zrf_u, zrf_t, zrf_q)
     fnameA = os.path.join(path,'test_36_output_py_082022_withnowavesinput_withwarmlayer.txt')
-    A_hdr = 'usr\ttau\thsb\thlb\thbb\thlwebb\ttsr\tqsr\tzo\tzot\tzoq\tCd\t'
+    A_hdr = 'usr\ttau\thsb\thlb\thbb\thsbb\thlwebb\ttsr\tqsr\tzo\tzot\tzoq\tCd\t'
     A_hdr += 'Ch\tCe\tL\tzeta\tdT_skinx\tdq_skinx\tdz_skin\tUrf\tTrf\tQrf\t'
     A_hdr += 'RHrf\tUrfN\tTrfN\tQrfN\tlw_net\tsw_net\tLe\trhoa\tUN\tU10\tU10N\t'
     A_hdr += 'Cdn_10\tChn_10\tCen_10\thrain\tQs\tEvap\tT10\tT10N\tQ10\tQ10N\tRH10\t'
diff --git a/Python/COARE3.6/coare36vn_zrf_et.py b/Python/COARE3.6/coare36vn_zrf_et.py
index d49f76b..babe7ea 100644
--- a/Python/COARE3.6/coare36vn_zrf_et.py
+++ b/Python/COARE3.6/coare36vn_zrf_et.py
@@ -985,7 +985,7 @@ def albedo_vector(sw_dn = None,jd = None,lon = None,lat = None,eorw = None):
     fnameA = os.path.join(path,'test_36_output_py_082022_withwavesinput.txt')
     # A=coare36vn_zrf_et(u, zu , t, zt, rh, zq, P, ts, sw_dn, lw_dn, lat, lon,jd, zi,rain, Ss, None , None, zrf_u, zrf_t, zrf_q)
     # fnameA = os.path.join(path,'test_36_output_py_082022_withnowavesinput.txt')
-    A_hdr = 'usr\ttau\thsb\thlb\thbb\thlwebb\ttsr\tqsr\tzo\tzot\tzoq\tCd\t'
+    A_hdr = 'usr\ttau\thsb\thlb\thbb\thsbb\thlwebb\ttsr\tqsr\tzo\tzot\tzoq\tCd\t'
     A_hdr += 'Ch\tCe\tL\tzeta\tdT_skinx\tdq_skinx\tdz_skin\tUrf\tTrf\tQrf\t'
     A_hdr += 'RHrf\tUrfN\tTrfN\tQrfN\tlw_net\tsw_net\tLe\trhoa\tUN\tU10\tU10N\t'
     A_hdr += 'Cdn_10\tChn_10\tCen_10\thrain\tQs\tEvap\tT10\tT10N\tQ10\tQ10N\tRH10\t'