p_wfs0.set_type("sh")
p_wfs0.set_nxsub(20)
p_wfs0.set_npix(4)
p_wfs0.set_pixsize(1.5)
p_wfs0.set_fracsub(0.5)
p_wfs0.set_xpos(0.)
p_wfs0.set_ypos(0.)
p_wfs0.set_Lambda(0.647) #CHECK
p_wfs0.set_gsmag(13.)
p_wfs0.set_optthroughput(0.3)
p_wfs0.set_zerop(1.e11)
p_wfs0.set_noise(5.11)
p_wfs0.set_atmos_seen(1)
#dm
p_dm0 = conf.Param_dm()
p_dm1 = conf.Param_dm()
p_dms = [p_dm0, p_dm1]
p_dm0.set_type("pzt")
nact = p_wfs0.nxsub + 1
p_dm0.set_nact(nact)
p_dm0.set_alt(0.)
p_dm0.set_thresh(0.14)
p_dm0.set_coupling(0.1458)
p_dm0.set_unitpervolt(0.4) # From Sam
p_dm0.set_push4imat(1.) # From Sam
# Remove DTT mirror for initial simulation (but the TTM is there for Mode 1 Keck)
p_dm1.set_type("tt")
p_dm1.set_alt(0.)
p_dm1.set_unitpervolt(0.016) # unedited
p_wfs4.set_optthroughput(0.5)
p_wfs4.set_zerop(1.e11)
p_wfs4.set_noise(1.)
p_wfs4.set_atmos_seen(1)
# lgs parameters
# p_wfs0.set_gsalt(90*1.e3)
# p_wfs0.set_lltx(0.)
# p_wfs0.set_llty(0.)
# p_wfs0.set_laserpower(10)
# p_wfs0.set_lgsreturnperwatt(1.e3)
# p_wfs0.set_proftype("Exp")
# p_wfs0.set_beamsize(0.8)
# dm
p_dm0 = conf.Param_dm()
p_dm1 = conf.Param_dm()
p_dm2 = conf.Param_dm()
p_dm3 = conf.Param_dm()
p_dms = [p_dm0, p_dm1, p_dm2, p_dm3]
p_dm0.set_type("pzt")
p_dm0.set_nact(41)
p_dm0.set_alt(0.)
p_dm0.set_thresh(0.3)
p_dm0.set_coupling(0.2)
p_dm0.set_unitpervolt(0.01)
p_dm0.set_push4imat(100.)
p_dm1.set_type("pzt")
p_dm1.set_nact(51)
p_wfs0.set_optthroughput(0.5)
p_wfs0.set_zerop(1.e11)
p_wfs0.set_noise(-1.)
p_wfs0.set_atmos_seen(1)
# lgs parameters
# p_wfs0.set_gsalt(90*1.e3)
# p_wfs0.set_lltx(0)
# p_wfs0.set_llty(0)
# p_wfs0.set_laserpower(10)
# p_wfs0.set_lgsreturnperwatt(1.e3)
# p_wfs0.set_proftype("Exp")
# p_wfs0.set_beamsize(0.8)
# dm
p_dm0 = conf.Param_dm()
p_dms = [p_dm0]
p_dm0.set_type("kl")
nact = p_wfs0.nxsub + 1
p_dm0.set_nkl(30)
p_dm0.set_alt(0.)
p_dm0.set_thresh(0.3)
p_dm0.set_coupling(0.2)
p_dm0.set_unitpervolt(0.001)
p_dm0.set_push4imat(1.)
# centroiders
p_centroider0 = conf.Param_centroider()
p_centroiders = [p_centroider0]
p_centroider0.set_nwfs(0)
p_wfss[i].set_noise(3)
p_wfss[i].set_atmos_seen(1)
#lgs parameters
#p_wfss[0].set_gsalt(90*1.e3)
#p_wfss[0].set_lltx(0)
#p_wfss[0].set_llty(0)
#p_wfss[0].set_laserpower(10)
#p_wfss[0].set_lgsreturnperwatt(1.e3)
#p_wfss[0].set_proftype("Exp")
#p_wfss[0].set_beamsize(0.8)
#dm
#p_dm0=conf.Param_dm()
#p_dm1=conf.Param_dm()
p_dms = [conf.Param_dm(), conf.Param_dm()]
p_dms[0].set_type("pzt")
nact = p_wfss[0].nxsub + 1
p_dms[0].set_nact(nact)
p_dms[0].set_alt(0.)
p_dms[0].set_thresh(0.3)
p_dms[0].set_coupling(0.2)
p_dms[0].set_unitpervolt(1.)
p_dms[0].set_push4imat(1.)
p_dms[0].set_influType("radialSchwartz")
p_dms[1].set_type("tt")
p_dms[1].set_alt(0.)
p_dms[1].set_unitpervolt(1.)
p_dms[1].set_push4imat(1.)