Rather than "joining" the doped regions, you should approach the problems as usual when constructing a device, and set it up without doping. Then, once the device model is eastablished, incl. the DeviceCalculator, you edit the Python script and add the "charge" keywords to the electrode calculators.