Obtaining the current from CHI
During pw_jobs nmr_shift the magnetic susceptibility is written to a file CHI. During pw_jobs nmr_plot CHI is reread and the induced current calculated. For plotting perposes the magnetic field is 1Tesla in the direction of ket eg ket 0 0 1 . Plots along a line or plane can be defined see section 5.10. As the default the current will be plotted as a 3 component vector. To plot the projection of the current on the direction bra set output_flags_curproj and, for example bra 1 0 0.
Obtaining the current directly
By setting output_flags_current the file JGS is written during pw_jobs nmr_shift. If output_flags_current is present during pw_jobs nmr_plot plots can be made as in the previous section. (This is only implemented for the crystal method, only trivial modifications are needed to make it work in eqn3 and eqn8)