Theta_r = tLoc(i,r)
resTheta_r = (Theta_r*thk_r-bblTheta_r*bblEta_r)/resThk_r
Theta_r = tLoc(r,j)
_RL Theta_r, bblTheta_d, bblTheta_r, resTheta_r