Colorize bars of a faceted graph of two categorical variables

[Felix Kaysers]

I am using Stata 17 Standard Edition on Windows.

I have created a faceted graph of two categorical variables using:

use https://www.stata-press.com/data/r17/nlsw88.dta
twoway (histogram industry), by(occupation)

I would now like to colorize the bars standing for each industry in the same color across all graphs and add a legend for industry.
I have found different questions that potentially relate to my question on Statalist (i.e. here and here) and I have tried out several of them, but I am fairly inexperienced in Stata and it does not seem to do the trick for me.

Stackoverflow link to my question.

[Nick Cox]

Comments at the Stack Overflow link.