The routine calc_physical_current of src/system/magnetic.F90 should be removed and replaced by the use of st%current