Hi Kalyani, The weights are computed as follows LANCZOS(A,F1,F2,N) nn = (N-1)/2 weight(0) = 2.*(f2-f1) do k = 1,nn pik = pi*k aa = pik/nn sigma = sin(aa)/aa aa = 2.*pik*f1 bb = 2.*pik*f2 weight( k) = ( sin(bb)/pik - sin(aa)/pik )*sigma end do And then they are applied as follows (arg_1 is the data, and there are other loops over the x, y, z dimensions of arg_1, which accounts for the i, j, k indices) DO 100 l=res_lo_ss(T_AXIS)+nn, res_hi_ss(T_AXIS)-nn do m = 1,nn result(i,j,k,l) = result(i,j,k,l) + weight(m)* . ( arg_1(i,j,k,l-m) + arg_1(i,j,k,l+m) ) end do 100 end do C Kalyani Devasena wrote: Dear ferret users, |