H = ω_1ad_1a_1 + ω_2ad_2a_2 + EC_33(n_3-ng_3)^2 + EL_1(-φ_1-φ_2)(-0.5φ_e1) + EL_2(-φ_1+φ_2)(-0.5φ_e1) + EJ_1cos(-φ_1+φ_3) + EJ_2cos(-φ_1-φ_3) \n------------------------------------------------------------\nmode 1:       harmonic       φ_1 = zp_1(a_1+ad_1)       ω_1/2π = 3.22489       zp_1 = 2.49e+00\nmode 2:       harmonic       φ_2 = zp_2(a_2+ad_2)       ω_2/2π = 0.39497       zp_2 = 8.72e-01\nmode 3:       charge         ng_3 = 0\n------------------------------------------------------------\nparameters:       EC_33 = 0.296       EL_1 = 0.13       EL_2 = 0.13       EJ_1 = 5.0       EJ_2 = 5.0       \nloops:            φ_e1/2π = 0.5       