This function may be performed alone by calling the script with entrystage and finalstage=`addskyarea'. At present, all that is done is that the total sky area of non-zero exposure of the mosaic of images is calculated and stored as keyword SKY_AREA in the `reciprocal sensitivity map'.
