szevis

 V=szesrc(par,u,v,band,nu)

 Same as pntsrc.m but for model of SZE cluster
 Calculate the complex visibilities expected for a source
 location in par and a set of baselines u,v

 par(1) = zero spacing flux at band center
 par(2) = Exp width
 par(3:4) = l,m position on field

0001 function V=szevis(par,u,v,band,nu)
0002 % V=szesrc(par,u,v,band,nu)
0003 %
0004 % Same as pntsrc.m but for model of SZE cluster
0005 % Calculate the complex visibilities expected for a source
0006 % location in par and a set of baselines u,v
0007 %
0008 % par(1) = zero spacing flux at band center
0009 % par(2) = Exp width
0010 % par(3:4) = l,m position on field
0011 
0012 % Calc relative SZ flux at each freq
0013 % Normalize to unity at band center
0014 % Note that this is not right if cluster is not at field center
0015 d=sze(nu)./sze(mean(nu(:)));
0016 % Combine with fit amp - the zero spacing flux
0017 s=cvec(par(1).*d);
0018 
0019 % Expand out to each vis
0020 s=s(band);
0021 
0022 % Compute vis amp rel to zero spacing for this baseline length
0023 r=exp(-2*pi*par(2)*pyth(u,v));
0024 
0025 % Dot product of src location and baseline vectors
0026 a=2*pi*(par(3)*u+par(4)*v);
0027 
0028 % Expected visibility
0029 V=s.*r.*complex(cos(a),sin(a));
0030 
0031 return
0032 
0033 function delIsz=sze(nu)
0034 Tcmb=2.726;
0035 fx=sz_fx(nu); delTsz=Tcmb.*fx;
0036 dIdT=planck_dIdT(nu,Tcmb);
0037 delIsz=delTsz.*dIdT;
0038 return