New initialization routines

Two new initialization routines:

initialize_spinice_q0: the system is in the simplest q=0 state.

intialize_spinice_X: the system is in the state with momentum q=(1/2
1/2 0).

The second only works for L is even.

Both routines are tested for L=1 and 2 systems.