added diagnostic 'bwprel'=P/P_o
added diagnostic 'effbwp' = P_o - P