program sort_wf
 implicit none
 read_wf = .true.
 touch read_wf 
 call save_wavefunction

end