program print_var_energy
implicit none
read_wf = .True.
touch read_wf
call routine
end
subroutine routine
print*,'psi_energy = ',psi_energy