Add a function calculating curl of vector variable.
Add current density calculation to mhd test program and reconnection test.