Bugfix: decrease indices of up- and down-type mass matrics on r.h.s.

of constraints.  In this way, the user can write constraints of the
form

LowScaleInput = {
    {Yu[1,1], topDRbar[1,1] / (-Sqrt[2] v)},
    {Yu[2,2], topDRbar[2,2] / (-Sqrt[2] v)},
    {Yu[3,3], topDRbar[3,3] / (-Sqrt[2] v)},
    {Yd[1,1], bottomDRbar[1,1] / (Sqrt[2] v)},
    {Yd[2,2], bottomDRbar[2,2] / (Sqrt[2] v)},
    {Yd[3,3], bottomDRbar[3,3] / (Sqrt[2] v)},
    {Ye, Automatic}
};

meta/FlexibleSUSY.m

@@ -1679,7 +1679,12 @@ corresponding tadpole is real or imaginary (only in models with CP
            allParameters = ((#[[1]])& /@ Join[Join @@ susyBetaFunctions, Join @@ susyBreakingBetaFunctions]) /.
                                a_[Susyno`LieGroups`i1] :> a /.
                                a_[Susyno`LieGroups`i1,SARAH`i2] :> a;
-           allIndexReplacementRules = Parameters`CreateIndexReplacementRules[allParameters];
+           allIndexReplacementRules = Join[
+               Parameters`CreateIndexReplacementRules[allParameters],
+               {Global`topDRbar[i_,j_]      :> Global`topDRbar[i-1,j-1],
+                Global`bottomDRbar[i_,j_]   :> Global`bottomDRbar[i-1,j-1],
+                Global`electronDRbar[i_,j_] :> Global`electronDRbar[i-1,j-1]}
+           ];
            Parameters`SetModelParameters[allParameters];
 
            (* collect all phases from SARAH *)