Closed will-henney closed 11 years ago
Nahiely is now working on this. The strategy is to first create the entire PPV cube (e.g., with model_profile_fast.py), then to extract a particular aperture from that.
Now we have it. The python script is emission/extract_aperture.py.
It reads the ppv cubes and extract and aperture cube. After that, it calculates the aperture profile (saved in an special folder inside de model, each line in a each file.FITS) and the total aperture flux normalized to Hbeta (just printed).
First step would be assume no seeing, so just reject all model points whose projected position is outside the aperture.
Second step would use erf functions to simulate the seeing