Hi, I think the parameter tw_Col in lines 80 and 81 should change with the parameter tpz based on equation 18 in the paper 2020.