# I'd like to photograph NEOWISE from Los Angeles, what is the optimal time?

I live in Los Angeles and have a camera adapter for my telescope -- nothing fancy just an Astromaster 114. What is the optimal time, in terms of dark sky and bright comet, to take a photo of Comet NEOWISE? I.e., when will it be brightest relative to the surround sky? When will it be largest?

And, lastly, what is the best time tonight to view it?

EDIT: I've spent the last two nights at sunset trying to locate it, but it doesn't really get dark until about an hour after sunset around 8pm. I also spent over an hour on the roof this morning from about 4am til 5:15am. I couldn't see it despite clear skies (or what passes for clear in LA).

You may try few hours before the sun rise tomorrow (12). Something to get started is here. https://theskylive.com/c2020f3-info

The bright start Capella, which is relatively easy to spot, could guide you to the comet.

Now that it has passed the perihelion of its orbit, it is receding from the Sun. So the brightness will gradually drop. Some details that might interest you are here.

http://astro.vanbuitenen.nl/comet/2020F3

• That first link is very helpful once I enter my own location, especially the section Comet C/2020 F3 (NEOWISE) rise and set times which tells me when it rises and sets and provides alt-az coordinates. – S. Imp Jul 12 '20 at 23:05
• I have marked this as the answer because the skylive link allows you to enter your geographic location and get information specific to your location on the globe. Additionally, it provides easy-to-use azimuth and altitude details which are easy to comprehend and interpret with just a compass. I am delighted to report that I saw the comet last night. I was unable to photograph it, but really enjoyed having a look. – S. Imp Jul 18 '20 at 19:51
• great, so glad that you did see it in the end. – mysterium Jul 19 '20 at 20:56

I don't know if this is helpful but using the Python package Skyfield one can calculate the positions of things including comet C/2020 F3 (NEOWISE) and the times that they rise and set for the month of July:

Here are the rise and set times of the Sun and the comet for reference for the next week

sunset:  2020-07-12T03:06:17Z
sunrise:  2020-07-12T12:51:20Z
sunset:  2020-07-13T03:05:55Z
sunrise:  2020-07-13T12:51:56Z
sunset:  2020-07-14T03:05:32Z
sunrise:  2020-07-14T12:52:33Z
sunset:  2020-07-15T03:05:07Z
sunrise:  2020-07-15T12:53:10Z
sunset:  2020-07-16T03:04:41Z
sunrise:  2020-07-16T12:53:48Z
sunset:  2020-07-17T03:04:13Z
sunrise:  2020-07-17T12:54:27Z
sunset:  2020-07-18T03:03:43Z
sunrise:  2020-07-18T12:55:06Z

cometset:  2020-07-12T04:15:12Z
cometrise:  2020-07-12T10:46:43Z
cometset:  2020-07-13T04:34:10Z
cometrise:  2020-07-13T10:48:01Z
cometset:  2020-07-14T04:54:27Z
cometrise:  2020-07-14T10:51:29Z
cometset:  2020-07-15T05:15:33Z
cometrise:  2020-07-15T10:57:33Z
cometset:  2020-07-16T05:36:48Z
cometrise:  2020-07-16T11:06:45Z
cometset:  2020-07-17T05:57:17Z
cometrise:  2020-07-17T11:19:29Z
cometset:  2020-07-18T06:16:07Z
cometrise:  2020-07-18T11:36:02Z


Here's the python script:

from skyfield.api import Topos
from skyfield.data import mpc
from skyfield import almanac
import numpy as np
import matplotlib.pyplot as plt
from skyfield.constants import GM_SUN_Pitjeva_2005_km3_s2 as GM_SUN

ts = loaddata.timescale() # include builtin=True if you want to use older files (you may miss some leap-seconds)

comets = comets.set_index('designation', drop=False)

row = comets.loc['C/2020 F3 (NEOWISE)']
print(row)

earth, sun, moon = [eph[x] for x in ('earth', 'sun', 'moon')]

comet = sun + mpc.comet_orbit(row, ts, GM_SUN)

times = ts.utc(2020, 7, range(1, 32))

RA, Dec = [thing._degrees for thing in (ra, dec)]

plt.plot(RA, Dec, 'ok')
plt.gca().set_aspect('equal')
plt.show()

t0 = ts.utc(2020, 7, 12)
t1 = ts.utc(2020, 7, 19)
LA = Topos('34.05 N', '118.25 W') # 34.05, -118.25

times, is_sunrises = almanac.find_discrete(t0, t1, almanac.sunrise_sunset(eph, LA))
for time, is_sunrise in zip(times, is_sunrises):
if is_sunrise:
print('sunrise: ', time.utc_iso())
else:
print('sunset: ', time.utc_iso())

times, is_rises = almanac.find_discrete(t0, t1, almanac.risings_and_settings(eph, comet, LA))
for time, is_sunrise in zip(times, is_rises):
if is_sunrise:
print('cometrise: ', time.utc_iso())
else:
print('cometset: ', time.utc_iso())

• I love this answer so much, but my understanding of RA and Dec is woefully insufficient to translate that chart into any me-centric coordinates. Those times are UTC time, yes? So I'd need to subtract 7 hours to convert those to LA local time, yes? – S. Imp Jul 12 '20 at 23:15
• @S.Imp yes subtract 7 hours to convert to Pacific Daylight time i.stack.imgur.com/lx4aE.png Okay this is new feature to me too and I need to double check it but I can see if I can give you altitude and azimuth instead. I'll ping you again when I've updated it. – uhoh Jul 13 '20 at 1:08
• I'm guessing that the me-centric coordinates will be a lot more elaborate because of the earth rotating, etc. I.e, the coordinates will change greatly throughout the day -- or even while the comet is visible at morning/night. – S. Imp Jul 13 '20 at 1:31

Nothing is more helpful to an aspiring astronomer than a good planetarium program.