What is your analysis objective? Are you trying to calculate JXB Lorentz forces in the volume of the superconductor in a device (e.g., Tokamak fusion reactor coils, MRI field coil)? Or are you trying to simulate the Meissner Effect which occurs when the superconductor is cooled so that it transitions from a normal conduction state to a superconducting one?
I have no recollection of anyone ever trying to simulate the London depth at the surface of a superconductor and I doubt our technology is capable of doing so.
We usually define a flux parallel boundary on the surface of a superconductor to model its effect on the surrounding field and don't concern ourselves with what's happening inside the superconducting domain. Actually, the superconductor is usually not modeled at all - it can be represented by a void in the modeled domain with flux parallel boundary condictions imposed on the void boundaries.
Including modeled representation of a superconducting body with properties like yours may result in a poorly conditioned system matrix - the magnitudes of the properties may be too different from those assigned to other regions in the model for machine precision to calculate accurately.
Kind regards,
Bill