;change units of area from arcseconds^2 to cm^2 function units, area_arcsec ;changes units to cm^2 area = area_arcsec * 72498958.3^2 return, area end