ind_compaction function

Calculate indicator for subsoil compaction