From 5b29275468fc662e30797fe0560b5e812d13570c Mon Sep 17 00:00:00 2001 From: "zhaojunyan@shao.ac.cn" Date: Sun, 14 Apr 2024 17:49:09 +0800 Subject: [PATCH] test --- csst_mci_sim/csst_mci_sim.py | 17 ++++++++++++----- 1 file changed, 12 insertions(+), 5 deletions(-) diff --git a/csst_mci_sim/csst_mci_sim.py b/csst_mci_sim/csst_mci_sim.py index 085595f..23b9403 100644 --- a/csst_mci_sim/csst_mci_sim.py +++ b/csst_mci_sim/csst_mci_sim.py @@ -1923,7 +1923,7 @@ class MCIsimulator(): z_sat=float(self.orbit_pars[self.orbit_exp_num,3]) - self.log.info('load star catlog 111111111111111') + wave0, zodi0 = self.zodiacal(ra, dec, self.TianCe_day) # erg/s/cm^2/A/arcsec^2 @@ -1931,7 +1931,10 @@ class MCIsimulator(): # EarthShine from straylight sl = StrayLight(self.information['dir_path'],jtime=time_jd, sat=np.array([x_sat, y_sat, z_sat]), radec=np.array([(ra*u.degree).value, (dec*u.degree).value])) - + + self.log.info('load star catlog 3333333333333333333333333') + + earth_e = sl.caculateEarthShineFilter(filter='r') star_e = sl.caculateStarLightFilter(filter='r') angle_earth = earth_angle(time_jd, x_sat, y_sat, z_sat, ra, dec) @@ -1941,6 +1944,8 @@ class MCIsimulator(): earthshine_wave0, earthshine_flux0 = ill2flux(earth_e+star_e,self.information['dir_pah']) + + self.log.info('load star catlog 444444444') # sample as ifs wavelength wave_mci = np.linspace(2500, 11000, 8501) #np.arange(2500, 11000, 1) f1 = interp1d(wave0, zodi0) @@ -1961,7 +1966,7 @@ class MCIsimulator(): ######################################################################################## - + self.log.info('load star catlog 5555555555555') ######################################################################################### self.cal_sky_noise() @@ -1984,7 +1989,9 @@ class MCIsimulator(): ####################################################################### - if self.TianceEffect: + if self.TianceEffect: + + self.log.info('load star catlog 666666666666666') ra_list = self.star['ra_gaia'].tolist() dec_list = self.star['dec_gaia'].tolist() pmra_list = self.star['pmra_gaia'].tolist() @@ -2022,7 +2029,7 @@ class MCIsimulator(): st_magz=[] - + self.log.info('load star catlog 777777777') #################### generate star image ########## if self.debug: -- GitLab