Pu_offset outlet area by a according to Schneider

@BearImp/calcCap_puchtler2025.m

@@ -30,13 +30,13 @@ function C_out = calcCap_puchtler2025(s,R_RE,R_R,R_ZL,R_ZR,B,R_L,epsilon_r,alpha
     Theta_1l    = @(phi)     atan((+B_RhL - Delta_x)./(Delta_z*cos(phi) - sign(R_L)*sqrt(R_ZL^2-Delta_z^2*sin(phi).^2)));
     Theta_1r    = @(phi)     atan((-B_RhR - Delta_x)./(Delta_z*cos(phi) - sign(R_L)*sqrt(R_ZR^2-Delta_z^2*sin(phi).^2)));
     if method.outlet == 'air'
-        Theta_0lIn  = @(phi) max(atan((R_RE*sin(alpha))./(-Delta_z*cos(phi) + sign(R_L)*sqrt(R_L^2 -Delta_z^2*sin(phi).^2))) - max(a*sqrt(max(1/R_RE^2-phi.^2/b^2,0)),phi.*(phi>0)),Theta_1l(phi));
-        Theta_0rIn  = @(phi) min(atan((R_RE*sin(alpha))./(-Delta_z*cos(phi) + sign(R_L)*sqrt(R_L^2 -Delta_z^2*sin(phi).^2))) + max(a*sqrt(max(1/R_RE^2-phi.^2/b^2,0)),phi.*(phi>0)),Theta_1r(phi));
+        Theta_0lIn  = @(phi) max(atan((R_RE*sin(alpha))./(-Delta_z*cos(phi) + sign(R_L)*sqrt(R_L^2 -Delta_z^2*sin(phi).^2))) - max(a*sqrt(max(1/R_RE^2-phi.^2/b^2,0)),(a/R_RE+phi).*(phi>0)),Theta_1l(phi));
+        Theta_0rIn  = @(phi) min(atan((R_RE*sin(alpha))./(-Delta_z*cos(phi) + sign(R_L)*sqrt(R_L^2 -Delta_z^2*sin(phi).^2))) + max(a*sqrt(max(1/R_RE^2-phi.^2/b^2,0)),(a/R_RE+phi).*(phi>0)),Theta_1r(phi));
         Theta_0lOut = @(phi)     atan((R_RE*sin(alpha))./(-Delta_z*cos(phi) + sign(R_L)*sqrt(R_L^2 -Delta_z^2*sin(phi).^2))) -     a*sqrt(max(1/R_RE^2-phi.^2/b^2,0));
         Theta_0rOut = @(phi)     atan((R_RE*sin(alpha))./(-Delta_z*cos(phi) + sign(R_L)*sqrt(R_L^2 -Delta_z^2*sin(phi).^2))) +     a*sqrt(max(1/R_RE^2-phi.^2/b^2,0));
-        phi_0 = a/R_RE/sqrt(1+a^2/b^2);
-        phi_1l = fzero(@(phi) Theta_0l(phi) - phi - Theta_1l(phi),pi/4);
-        phi_1r = fzero(@(phi) Theta_0r(phi) + phi - Theta_1r(phi),pi/4);
+        phi_0 = 0;%a/R_RE/sqrt(1+a^2/b^2);
+        phi_1l = fzero(@(phi) Theta_0l(phi) - (a/R_RE+phi) - Theta_1l(phi),pi/4);
+        phi_1r = fzero(@(phi) Theta_0r(phi) + (a/R_RE+phi) - Theta_1r(phi),pi/4);
         % fprintf('phi_0 = %d | phi_1l = %d | phi_1r = %d \n',phi_0,phi_1l,phi_1r)
     end
     Theta_2l    = @(phi)     atan((+B/2   - Delta_x)./(Delta_z*cos(phi) - sign(R_L)*sqrt(R_ZL^2-Delta_z^2*sin(phi).^2)));