@@ -0,0 +1,658 @@
+from __future__ import print_function%0Aimport os%0Aimport numpy as np%0A%0Afrom legacypipe.forced_photom_decam import main as forced_photom_main%0Afrom legacypipe.forced_photom_decam import get_parser as forced_photom_parser%0Afrom legacypipe.runcosmos import CosmosDecals%0A    %0Adef main():%0A    parser = forced_photom_parser()%0A    parser.add_argument('--subset', type=int, help='COSMOS subset number %5B0 to 4%5D', default=0)%0A    opt = parser.parse_args()%0A    decals = CosmosDecals(subset=opt.subset)%0A    %0A    np.random.seed(1000000 + opt.subset)%0A    %0A    return forced_photom_main(decals=decals, opt=opt)%0A    %0Aif __name__ == '__main__':%0A    import sys%0A    sys.exit(main())%0A

add script for performing forced photom of DECam test images in COSMOS region