Use the crude relation
LogLbol = -0.4*M_i + 36 = LogL3000 + 0.7
to link Mi with L3000 as the default, which could still be over-ridden by a direct input.
The new hierarchy would instead be (1) Mi and L3000 both directly input (and hence controlled separately) by the user, (2) Mi is tied to L3000, where L3000 is either (i) set by the user, or (ii) inferred from the input redshift by assuming some typical redshift-luminosity relation. This would also involve changing the default L3000 behaviour such that L3000 varies with redshift (it's currently held fixed at 10^46 unless otherwise specified).
Update associated documentation in python comments and README.
Use the crude relation LogLbol = -0.4*M_i + 36 = LogL3000 + 0.7 to link Mi with L3000 as the default, which could still be over-ridden by a direct input. The new hierarchy would instead be (1) Mi and L3000 both directly input (and hence controlled separately) by the user, (2) Mi is tied to L3000, where L3000 is either (i) set by the user, or (ii) inferred from the input redshift by assuming some typical redshift-luminosity relation. This would also involve changing the default L3000 behaviour such that L3000 varies with redshift (it's currently held fixed at 10^46 unless otherwise specified).
Update associated documentation in python comments and README.