restore, '/home/aosup/adopt/measures/adsec_calib/if_functions/20110704_173413/adsec_save.sav', /ver
slaves = adsec.mir_act2act[indgen(total(adsec_shell.n_act_ring[0:2]))]
print, intersection(slaves, adsec.act_w_cl, slcl)
cmd = float(slave_commands(slcl, adsec_save.ff_p_svec, adsec_save=adsec_save))
writefits, '~/fernando_m2c.fits', cmd

qui c'e' un problema: il tracking number: 20110704_175709 e' preso/ridotto male dall'interferometro

ora procediamo con gli slave nel fagiolo che non si riesce a riprendere in modo da oscurare le pupille del W da aggiungere: [437,438,439,425,424,430,429,428,427,426]

m2cf = '/towerdata/adsec_calib/M2C/KL_MAG585_20110706/m2c.fits'
m2c = readfits(m2cf)
offmat = readfits('/home/aosup/Source/Magellan/calib/adsec/585m/data/ff_matrix.fits')
ocl = readfits('/home/aosup/Source/Magellan/calib/adsec/585m/data/ff_list.fits')
slaves = [adsec.mir_act2act[indgen(total(adsec_shell.n_act_ring[0:3]))] , 437,438,439,425,424,430,429,428,427,426]
print, intersection(slaves, adsec.act_w_cl, slcl)
new_m2c = slave_commands(slcl,m2c)
new = adsec.ff_matrix ## m2c
old=new*0 & old[*,ocl] = offmat ## m2c[*,ocl] 
new_slaved = adsec.ff_matrix ## new_m2c
print, minmax(old[*, slcl])
print, minmax(new[*, slcl])
print, minmax(new_slaved[*, slcl])
writefits, '/towerdata/M2C/KL_SLAVE_v1/m2c.fits',new_m2c



f1='/home/aosup/adopt/measures/adsec_calib/flat/20121124_200912/20121124_100modes2.sav'
f2='/home/aosup/adopt/measures/adsec_calib/flat/20111125_164852/20111125_400modes.sav'
restore, f1, /ver & ff1 = flattened_status.current
restore, f2, /ver & ff2 = flattened_status.current
cl = adsec.act_w_cl
df = (ff2-ff1) & df -= mean(df[cl])
ffred = adsec.ff_matrix[adsec.act_w_cl, *]
ffred = ffred[*, adsec.act_w_cl]
dpos = pseudo_invert(ffred, /dou, n_mod=7) ## (df[cl] - mean(df[cl]))
dp = fltarr(672) & dp[cl]=dpos & print, minmax(adsec.ff_matrix ## dp)
x = adsec.act_coordinates[0,*]
y = adsec.act_coordinates[1,*]
kk =  surf_fit(x[cl]/adsec_shell.out_radius, y[cl]/adsec_shell.out_radius, dpos, /zern, indgen(22)+1, coeff=cc)
cc[0:5]=0
cc[8:*]=0
 uu= dpos-2*(pm_offload.pmz ## cc)[cl]
dp = fltarr(672) & dp[cl]=uu & print, minmax(adsec.ff_matrix ## dp)
writefits, '~/comandi.fits', dp

-- MarcoXompero - 24 Nov 2012
Topic revision: r3 - 25 Nov 2012, MarcoXompero
This site is powered by FoswikiCopyright © by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding Foswiki? Send feedback