extrapolated galaxy models

For HSC simulations, Francois ( @EiffL ) has been working on a way to easily extrapolate the COSMOS I<25.2 sample to fainter magnitudes and make a deeper parametric galaxy sample. Doing this extrapolation isn't likely good enough for precise tests of shear calibration for deeper samples, but it's likely good enough to just test for issues like "my real galaxies brighter than 25.2 are blended with faint galaxies some fraction of the time; what's the impact of that?" My question is whether there's interest in having a routine that will make synthetic extrapolated Sersic samples in GalSim? If there is interest then Francois would be happy to incorporate his code and make a routine for this in scene.py.

Start by reviewing scene.py and the proposed extrapolation code from Francois, then read the discussion to determine whether there is agreement on adding the routine. A complete contribution would define and implement a way to generate synthetic extrapolated Sersic samples for the described HSC simulation use case, with project tests covering the resulting behavior.