def test_calc():
M = DotMap()
M.PC = 40
M.PEEP0 = 5
M.RR = 18
M.IT = 1
M.Pbar = 760
P1 = DotMap()
P1.H = 140
P1.W = 55
P1.Cr = 15
P1.PF = 100
P1.HCO3 = 24
P1.FGFO = 8
P1.FGFA = 2
P1.PEEP = 5
P2 = DotMap()
P2.H = 180
P2.W = 70
P2.Cr = 50
P2.PF = 300
P2.HCO3 = 24
P2.FGFO = 1
P2.FGFA = 8
P2.PEEP = 5
solver(M, P1, P2)
