| 1 |
PRO PWIN, file1, file2, startday, offdays, ndays, stepdays, tstep |
| 2 |
|
| 3 |
ppd=long(86400/tstep) |
| 4 |
|
| 5 |
g=readfits(file1) |
| 6 |
gf=readfits(file2) |
| 7 |
i0=offdays*ppd |
| 8 |
ndt=ndays*ppd |
| 9 |
|
| 10 |
g0=g(i0:i0+ndt-1) |
| 11 |
gx=gf(i0:i0+ndt-1) |
| 12 |
gfdc=total(rebin(gx,1)) |
| 13 |
print,startday+offdays,ndays,total(rebin(g0,1)),gfdc,1.0/gfdc |
| 14 |
px=abs(fft(gx,1))^2 |
| 15 |
df=1d6/tstep/ndt |
| 16 |
nx=long(200/df) |
| 17 |
openw,1,'pwin.'+strtrim(string(ndays),2)+'.'+strtrim(string(startday+offdays),2) |
| 18 |
for j=-nx,nx do printf,1,j*df,px(abs(j))/px(0) |
| 19 |
close,1 |
| 20 |
|
| 21 |
nperiods=ndays/stepdays |
| 22 |
ndt=stepdays*ppd |
| 23 |
df=1d6/tstep/ndt |
| 24 |
nx=long(200/df) |
| 25 |
|
| 26 |
for i3=0,nperiods-1 do begin |
| 27 |
i=i0+ndt*i3 |
| 28 |
g0=g(i:i+ndt-1) |
| 29 |
gx=gf(i:i+ndt-1) |
| 30 |
gfdc=total(rebin(gx,1)) |
| 31 |
print,startday+offdays+stepdays*i3,stepdays,total(rebin(g0,1)),gfdc,1.0/gfdc |
| 32 |
px=abs(fft(gx,1))^2 |
| 33 |
openw,1,'pwin.'+stringn(stepdays,2)+'.'+stringn(startday+offdays+stepdays*i3,4) |
| 34 |
for j=-nx,nx do printf,1,j*df,px(abs(j))/px(0) |
| 35 |
close,1 |
| 36 |
end |
| 37 |
|
| 38 |
end |
