Jean-Michel, what is most accurate way to compute heat flux across lines of latitude with rStar turned on? Should I use VVELTH or VTHMASS? And if I use the latter, what vertical thickness is most appropriate at the surface: "delR(1)" or "delR(1)+ETAN"? Dimitris