Correction AO radius

hrm-hc/MAIN.py

@@ -8,7 +8,7 @@ import os
 import sys
 from time import localtime, strftime
 import time
-#import matplotlib.pyplot as plt
+import matplotlib.pyplot as plt
 from astropy.io import fits
 from utils_basic import minterp2d, DISPER, VECT, VECT1, BEAMSHIFT, ELEVATION, PARANG, ft, MAKE_ELT, PSF, LD2MAS, MAS2LD, M4Magic, DSP, EELTNLS, APPLY_CROSSTALK, DISPER_ADC
 
@@ -16,8 +16,8 @@ start_time = time.time()
 
 #HARMONI_DIR = '/Users/avigan/data/HARMONI-HC/'
 #COM_DIR     = '/Users/avigan/data/HARMONI-HC/COMMON_FILES/'
-HARMONI_DIR = '/Users/carlotal/harmoni/'
-COM_DIR     = '/Users/carlotal/harmoni/COMMON_FILES/'
+HARMONI_DIR = '/Users/carlotal/harmoni-hc/hrm-hc'
+COM_DIR     = '/Users/carlotal/harmoni-hc/hrm-hc/COMMON_FILES/'
 
 # Artifical dispersion corrector (dispersion is gamma times what it should be)
 gamma = 1#%0.1;%1;
@@ -50,7 +50,7 @@ z_max = 50
 seeing = 'MED' 
 
 #Sensing wavelength, ZELDA [m]
-lambda_wfs = 1.2e-6
+lambda_wfs = 1.175e-6
 
 #CONFIG='Ideal';
 #CONFIG='Relay_PDR_ZELDA_REQUIREMENT';
@@ -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 = 'TEST_H'  # %Choose between the bands defined hereafter (use a valid string)
 
 ### Bands ###
 
@@ -115,7 +115,7 @@ K1_high  = [2.024, 2.113, 2.202, 17000]
 K2_high  = [2.197, 2.294, 2.390, 17000]
 
 if BAND=='TEST_H':
-    N_LD = 3
+    N_LD = 3#48
     LBD_min = H_band[0]*1e-6
     LBD_max = H_band[2]*1e-6
     if APODIZER == 'SP_1':
@@ -559,9 +559,9 @@ dsp = dsp*DSP_gamma**2
 size_dsp = len(dsp[0])
 dsp_array= dsp
 
-DiamFTO = 4
-dim = size_dsp/DiamFTO
-fmax     = size_dsp/2
+DiamFTO = 2#4
+dim = 200#size_dsp/DiamFTO
+fmax     = size_dsp/2 
 f1D,Xf,Yf,Rf,Tf = VECT(size_dsp,2*fmax)
 fudge_df = (1024/1132)**2
 df = np.abs(f1D[1]-f1D[0])*fudge_df