handling both superconducting and normal conducting cases

example/maxwell_sibc/sphere.dat

 // Geometry //
 DefineConstant[R  = { 100e-3, Name "Input/00Geometry/00Radius [m]" }];
-DefineConstant[K1 = { 1/R,    Name "Input/00Geometry/01Curvature 1 [m^-1]",
-                              ReadOnly 1}];
-DefineConstant[K2 = { 1/R,    Name "Input/00Geometry/02Curvature 2 [m^-1]",
-                              ReadOnly 1}];
 
 // Mesh //
 DefineConstant[md = { 10,     Name "Input/01Mesh/00Density [Radius^-1]" }];
 cl = R/md;
 
 // Material //
-DefineConstant[Sigma = { 1e15,
-                         Name "Input/02Material/00Conductivity [S m^-1]" }];
+DefineConstant[isSC   = { 0,
+                          Name "Input/02Material/00Is superconductor?",
+                          GmshOption "Reset",
+                          Choices {0 = "No", 1 = "Yes"} }];
+If(isSC == 0)
+DefineConstant[Sigma  = { 1e15,
+                          Name "Input/02Material/01Conductivity [S m^-1]" }];
+EndIf
+If(isSC == 1)
+DefineConstant[Sigma  = { 1e9,
+                          Name "Input/02Material/01Conductivity [S m^-1]" }];
+DefineConstant[Lambda = { 1e-8,
+                          Name "Input/02Material/02London depth [m]"      }];
+EndIf

example/maxwell_sibc/sphere.pro

@@ -17,7 +17,12 @@ Function{
   Eps0 = 1/(Mu0 * C0^2); // [F/m]
 
   // Conductivity (coming from sphere.dat) //
-  Sigma[] = Sigma; // [S/m]
+  If(isSC == 1)                                            // Superconductor
+    Sigma[] = Sigma - J[] / (Mu0 * Lambda^2 * Re[Puls[]]); // [S/m]
+  EndIf
+  If(isSC == 0)                                            // Normal conductor
+    Sigma[] = Sigma;                                       // [S/m]
+  EndIf
 
   // Leontovich SIBC (H-Formulation) //
   SL[] = Sqrt[(J[] * Re[Puls[]] * Mu0) / Sigma[]]; // H-Formulation