#tp.servo_error = [0.005,1] ##only relevant if trying AO!
tp.piston_error = False #piston error off/on
tp.use_coron = True #coronagraph off/on
tp.occulter_type = 'Vortex' #occulter type - vortex, none, gaussian, 8th_Order, None (Lyot Stop)
tp.detector = 'ideal'
tp.beam_ratio = 25/64./2
tp.check_args()
##set simulation parameters
sp.show_wframe = False
sp.save_obs = False
sp.num_processes = 4
sp.save_locs = np.array([['add_atmos',], ['coronagraph',]])
sp.return_E = True
sp.show_cube = False
#atmos params -- all currently set to defaults
cp.model = 'single' # single|hcipy_standard|evolving
cp.show_caosparams= True # for control over all other variables
cp.cn = 1e-14 #lower cn2 should mean less turbulence - this value is equivalent to r0~30 (>>D)
cp.L0 = 10 #longer L0 (coherence length) means less turbulent, I think
cp.v = np.asarray([5,0])
cp.h = 100
iop.update("complex-fields/")
###SET FILE SAVE NAME
iop.fields = os.path.join(iop.testdir, 'MEDIS_1sec_Aug2020.h5')
if __name__ == '__main__':
gpd.run_medis()
"""
import numpy as np
from medis.params import iop, sp, ap, tp, cdip
from medis.utils import dprint
import medis.optics as opx
from medis.plot_tools import view_spectra, view_timeseries, quick2D, plot_planes
import medis.medis_main as mm
#################################################################################################
#################################################################################################
#################################################################################################
# testname = input("Please enter test name: ")
testname = 'SCExAO-test1'
iop.update(testname)
# dprint(f"iop name = {iop.testname}")
iop.makedir()
# Telescope
tp.prescription = 'Subaru_SCExAO'
tp.entrance_d = 7.9716
tp.flen_primary = tp.entrance_d * 13.612
# Simulation & Timing
sp.numframes = 1
sp.closed_loop = False
# Grid Parameters
sp.focused_sys = True
sp.beam_ratio = 0.2 # parameter dealing with the sampling of the beam in the pupil/focal plane