If you don't want to simulate as lossless metal, have a look at the Sonnet EM solver. In Sonnet's loss model "general" you can define the superconductor's surface resistance and surface inductance (based on temperature dependent London penetration depth).