                                            2                                  ↪
H = ω₁⋅ad₁⋅a₁ + ω₂⋅ad₂⋅a₂ + EC₃₃⋅(-ng₃ + n₃)  - EJ₁⋅cos(φₑ₁/2 - φ₁ + φ₃) - EJ₂ ↪

↪                       
↪ ⋅cos(-φₑ₁/2 + φ₁ + φ₃)
------------------------------------------------------------
mode 1:       harmonic       φ₁ = zp₁⋅(a₁ + ad₁)       ω₁/2⋅π = 3.22489       zp₁ = 2.49
mode 2:       harmonic       φ₂ = zp₂⋅(a₂ + ad₂)       ω₂/2⋅π = 0.39497       zp₂ = 0.872
mode 3:       charge       ng₃ = 0
------------------------------------------------------------
parameters:       EC₃₃ = 0.296       EJ₁ = 5.0       EJ₂ = 5.0
loops:       φₑ₁/2⋅π = 0.5