single_slice

0001 function [x,z,c,rho,K,C_zero]=single_slice(p,S,tw);
0002 
0003 tic,
0004 %fname='input1';
0005 
0006 
0007 [xg,IIza,IIzb]=hl_int_IIz(p,S,tw);
0008 [xg,IIha,IIhb,IIhay,IIhby]=hl_int_IIh(p,S,tw);
0009 [xg,Iza,Izb]=hl_int_Iz(p,S,tw);
0010 [xg,Iha,Ihb,Ihay,Ihby]=hl_int_Ih(p,S,tw);
0011 
0012 
0013 [tmp,tmp]=size(Iha);
0014 K=zeros(tmp-1,tmp-1);
0015 
0016 
0017 %[ok,p]=read_input(fname);
0018 
0019 out=p.out_directory;
0020 grid=load([out 'grids']);
0021 z=grid.z(S);
0022 c=grid.c(S);
0023 rho=grid.rho(S);
0024 omega=p.w(tw);
0025 
0026 
0027 tmpa=Iza.*IIza + Iha.*IIha + Ihay.*IIhay;
0028 tmpb=Izb.*IIzb + Ihb.*IIhb + Ihby.*IIhby;
0029 
0030 
0031 tmpa=tmpa(2:end,2:end);
0032 tmpb=tmpb(2:end,2:end);
0033 
0034 
0035 
0036 K=2*omega*(2*pi)^10*calc_mzero(omega,p)*rho*(tmpa+tmpb);
0037 
0038 
0039 
0040 
0041 x=xg(2:end);
0042 
0043 tm=toc;
0044 disp(['finished slice ' num2str(S) ' in ' num2str(tm) ' seconds']);
0045 
0046 [C_zero]=compute_Czero(p,tw);
0047 
0048 return
0049 
0050 %
0051 % kernel for phase
0052 %
0053 phi=angle(C_zero);
0054 dp=cos(phi)*(imag(K)-tan(phi)*real(K))/abs(C_zero);
0055 
0056 
0057 dp=dp+fliplr(flipud(dp));
0058 
0059 
0060 dp=dp*1e24;  %% convert to 1/Mm^3;
0061 
0062 
0063 
0064 
0065