写一段使用jplephem计算地球在太阳系质心位置坐标的代码

JPL Ephemeris库是一个用于天文学计算的Python模块，它提供了一种方便的方式来获取行星和其他天体的位置信息。以下是一个简单的示例，展示了如何使用`jplephem`计算地球在太阳系质心的位置坐标：

from jplephem.spk import SPK
from jplephem import Ephemeris

# 加载SPK数据文件，包含天体的位置数据
spk = SPK.open('de421.bsp')  # 使用De421 ephemerides 数据

# 创建Ephemeris对象，指定地球作为目标天体
earth = 'Earth'
eph = Ephemeris(spk, earth)

# 获取当前日期，这里假设是UTC时间
import datetime
date_now = datetime.datetime.utcnow()

# 计算地球在太阳系质心的位置
location = eph.at(date_now)
position = location['x'], location['y'], location['z']  # 返回的是AU（天文单位）

print(f"地球在太阳系质心的坐标 (AU): ({position[0]:.6f}, {position[1]:.6f}, {position[2]:.6f})")

注意：这个例子需要首先安装`jplephem`库，你可以通过`pip install jplephem`来安装。