added calculate_F. no effective speedup over calculate_age so still using calculate_age for flux