def CCD_select(image,row=1,col=2):
    """Returns subset of image defined by row and column, with 0,0 representing
    the upper left CCD. Assumes 3840,3840 with 4x4 array of CCD chips. """
    return image[row*960:(row+1)*960,col*960:(col+1)*960]