diff --git a/tracy/tools/soltracy.cc b/tracy/tools/soltracy.cc
index 6b275d20931113aa20ac41060f7b01f64574bfe1..2b534d73c7ba3fafef98eb095659b6ebcdc10254 100644
--- a/tracy/tools/soltracy.cc
+++ b/tracy/tools/soltracy.cc
@@ -209,7 +209,7 @@ int main(int argc, char *argv[]) {
     printglob(); /* print parameter list */
   }
 
-  // add coupling by random rotating of the quadrupole magnets
+  // add coupling by random rotating of the full quadrupole magnets
   if (ErrorCouplingFlag == true) {
     SetErr(err_seed, err_rms);
     Ring_GetTwiss(chroma = true, 0.0); /* Compute and get Twiss parameters */
@@ -218,6 +218,16 @@ int main(int argc, char *argv[]) {
     GetEmittance(ElemIndex("cav"), true);
     printglob(); /* print parameter list */
   }
+  
+  // add coupling by random rotating of the half quadrupole magnets, delicated for soleil
+  if (ErrorCoupling2Flag == true) {
+    SetErr2(err_seed, err_rms);
+    Ring_GetTwiss(chroma = true, 0.0); /* Compute and get Twiss parameters */
+    printlatt(); /* dump linear lattice functions into "linlat.dat" */
+    Coupling_Edwards_Teng();
+    GetEmittance(ElemIndex("cav"), true);
+    printglob(); /* print parameter list */
+  }
 
   // WARNING Fit tunes and chromaticities before applying errors !!!!
   //set multipoles in all magnets