Skip to content
GitLab
Projects
Groups
Snippets
/
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Contribute to GitLab
Sign in
Toggle navigation
Menu
Open sidebar
csst-sims
csst_mci_sim
Commits
26c6c9b3
Commit
26c6c9b3
authored
Apr 14, 2024
by
Yan Zhaojun
Browse files
test
parent
d419ba29
Pipeline
#4030
failed with stage
in 0 seconds
Changes
1
Pipelines
1
Hide whitespace changes
Inline
Side-by-side
csst_mci_sim/csst_mci_sim.py
View file @
26c6c9b3
...
...
@@ -2425,19 +2425,19 @@ class MCIsimulator():
"""
# read solar template
solar_template
=
pd
.
read_csv
(
information
[
'dir_path'
]
+
'MCI_inputData/refs/solar_spec.dat'
,
sep
=
'\s+'
,
solar_template = pd.read_csv(
self.
information['dir_path']+'MCI_inputData/refs/solar_spec.dat', sep='\s+',
header=None, comment='#')
template_wave = solar_template[0].values
template_flux = solar_template[1].values
# read earth shine surface brightness
earthshine_curve
=
pd
.
read_csv
(
information
[
'dir_path'
]
+
'MCI_inputData/refs/earthshine.dat'
,
earthshine_curve = pd.read_csv(
self.
information['dir_path']+'MCI_inputData/refs/earthshine.dat',
header=None, comment='#')
angle = earthshine_curve[0].values
surface_brightness = earthshine_curve[1].values
# read V-band throughtput
cat_filter_V
=
pd
.
read_csv
(
information
[
'dir_path'
]
+
'MCI_inputData/refs/filter_Bessell_V.dat'
,
sep
=
'\s+'
,
cat_filter_V = pd.read_csv(
self.
information['dir_path']+'MCI_inputData/refs/filter_Bessell_V.dat', sep='\s+',
header=None, comment='#')
filter_wave = cat_filter_V[0].values
filter_response = cat_filter_V[1].values
...
...
@@ -2504,7 +2504,7 @@ class MCIsimulator():
lamda = abs(lb_obj.lon.degree - lb_sun.lon.degree)
# interpolated zodical surface brightness at 0.5 um
zodi
=
pd
.
read_csv
(
information
[
'dir_path'
]
+
'MCI_inputData/refs/zodi_map.dat'
,
sep
=
'\s+'
,
header
=
None
,
comment
=
'#'
)
zodi = pd.read_csv(
self.
information['dir_path']+'MCI_inputData/refs/zodi_map.dat', sep='\s+', header=None, comment='#')
beta_angle = np.array([0, 5, 10, 15, 20, 25, 30, 45, 60, 75])
lamda_angle = np.array([0, 5, 10, 15, 20, 25, 30, 35, 40, 45,
60, 75, 90, 105, 120, 135, 150, 165, 180])
...
...
@@ -2513,7 +2513,7 @@ class MCIsimulator():
zodi_obj = f(beta, lamda) # 10^�? W m�? sr�? um�?
# read the zodical spectrum in the ecliptic
cat_spec
=
pd
.
read_csv
(
information
[
'dir_path'
]
+
'MCI_inputData/refs/solar_spec.dat'
,
sep
=
'\s+'
,
header
=
None
,
comment
=
'#'
)
cat_spec = pd.read_csv(
self.
information['dir_path']+'MCI_inputData/refs/solar_spec.dat', sep='\s+', header=None, comment='#')
wave = cat_spec[0].values # A
spec0 = cat_spec[1].values # 10^-8 W m^�? sr^�? μm^�?
zodi_norm = 252 # 10^-8 W m^�? sr^�? μm^�?
...
...
@@ -3359,7 +3359,7 @@ class MCIsimulator():
.. Warning:: This method does not work if the input file has exactly one line.
"""
#######################################################################
cosmetics
=
np
.
loadtxt
(
information
[
'dir_path'
]
+
'MCI_inputData/data/Cosmetics_g.txt'
)
cosmetics = np.loadtxt(
self.
information['dir_path']+'MCI_inputData/data/Cosmetics_g.txt')
x = np.round(cosmetics[:, 0]).astype(int) ## row number
y = np.round(cosmetics[:, 1]).astype(int) ## col number
...
...
Write
Preview
Supports
Markdown
0%
Try again
or
attach a new file
.
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment