Commit 0c3b0333 authored by rfetick's avatar rfetick
Browse files

update example to match galpak PSF centering

parent 2216fc66
......@@ -63,8 +63,8 @@ class PsfaoGalpak(Psfao, PointSpreadFunction):
ww = wvl_um/self.wvl_ref_um # wavelength ratio
r0,C,A,alpha,ratio,theta,beta = self.param
p = [r0*ww**(6./5.),C*ww**(-2),A*ww**(-2),alpha,ratio,theta,beta]
dx = 0
dy = 0
dx = -(1-self.npix[0]%2)/2 # between 4 pixels if even nb of pixels
dy = -(1-self.npix[1]%2)/2 # between 4 pixels if even nb of pixels
return self.__call__(p, dx=dx, dy=dy)
......@@ -74,7 +74,7 @@ class PsfaoGalpak(Psfao, PointSpreadFunction):
if __name__ == "__main__":
### PARAMETERS
### PSF PARAMETERS
r0 = 0.10 # Fried parameter [m]
C = 1e-2
A = 0.8 # AO variance [rad²]
......@@ -95,9 +95,9 @@ if __name__ == "__main__":
self.xy_step = xy_step
self.z_central = z_central
nwvl = 3
nx = 15
ny = 15
nwvl = 3 # number of wavelengths in the sub-cube (whatever, not used here)
nx = 15 # number of pixels on X
ny = 15 # number of pixels on Y
xy_step = muse_wfm.resolution_mas*1e-3 # comply with Muse resolution for this dummy cube
z_central = 0.8 # [um]
for_cube = MuseCube((nwvl,nx,ny), xy_step, z_central)
......
Markdown is supported
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