========= nodensity ========= .. f:autosrcfile:: nodensity.f90