adding first code for subparcellation from bernard