Top: intensity maps produced by MCFOST for our model M50 of RY Lup in the K band for the whole disk (left) and zoom in at the milliarcsecond scale of the central region (right). Middle: observed squared visibilities (red circles) with GRAVITY (left) and PIONIER (right) compared with the visibilities computed with our model M50 (black circles). Bottom: observed closure phases (red circles) with GRAVITY (left) and PIONIER (right) compared with the closure phases computed with our model M50 (black circles).