H = ω_1ad_1a_1 + EC_22(n_2-ng_2)^2  - EJ_1cos(-φ_1+φ_2+0.5φ_e1) - EJ_2cos(-φ_1-φ_2+0.5φ_e1)\n------------------------------------------------------------\nmode 1:       harmonic       φ_1 = zp_1(a_1+ad_1)       ω_1/2π = 3.22489       zp_1 = 2.49e+00\nmode 2:       charge         ng_2 = 0\n------------------------------------------------------------\nparameters:       EC_22 = 0.296       EJ_1 = 5.0       EJ_2 = 5.0       \nloops:            φ_e1/2π = 0.5       