program pouet
implicit none
read_wf = .True.
touch read_wf
call routine
end
subroutine routine
provide integrated_delta_rho_all_points