Commit 09d954aa authored by Alexis Carlotti's avatar Alexis Carlotti
Browse files

AO correction

parent f2b0dfba
......@@ -27,10 +27,10 @@ gamma = 1#%0.1;%1;
rot_flag = 1
#Flag for missing segments (1 = missing segments ; 0 = no missing segments)
miss_flag = 0
miss_flag = 1
#Flag for E-ELT M1 phase aberrations
eeltphase_flag = 0
eeltphase_flag = 1
#Flag for cross-talk and diffusion ; FWHM=1pix across x ; 2pix across lambda ; 0.5% diffusion across lambda
crosstalk_flag = 0
......@@ -91,7 +91,7 @@ T_exp = 1/60
N_exp = np.int(T_obs/T_exp)# %60 exposures per hour = 1 per minute
##############################################################################
BAND = 'TEST_K' # %Choose between the bands defined hereafter (use a valid string)
BAND = 'H' # %Choose between the bands defined hereafter (use a valid string)
### Bands ###
......@@ -349,9 +349,9 @@ logging.basicConfig(filename=path_directory + 'simulation_log.txt', filemode='w'
format='%(message)s', level=logging.DEBUG)
log = logging.getLogger('harmoni-hc')
if not log.handlers:
handler = logging.StreamHandler(sys.stdout)
log.addHandler(handler)
#if not log.handlers:
handler = logging.StreamHandler(sys.stdout)
log.addHandler(handler)
log.info('DATA MAKER')
log.info('Number of exposures: {0}'.format(N_exp))
......@@ -562,9 +562,9 @@ dsp_array= dsp
DiamFTO = 2#4
dim = 200#size_dsp/DiamFTO
fmax = size_dsp/2
fmax = (1024/1132)*size_dsp/2 # size_dsp/2
f1D,Xf,Yf,Rf,Tf = VECT(size_dsp,2*fmax)
fudge_df = (1024/1132)**2
fudge_df = 1#(1024/1132)**2
df = np.abs(f1D[1]-f1D[0])*fudge_df
# _____ _____ __ __ _ _ _ _______ _____ _ _ _____ ________ _______ ____ _____ _ _ _____ ______ _____
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment