riq_ref = int(1000.0 * sqrt(float(tori_ref)) / float(yrange))