|
Revision 334
(checked in by gdusek, 14 years ago)
|
Total update to the DPWP code. Significant changes made to direspec, specmultiplot, radialtouvw and IMLM code. 6/14/10
|
| Line | |
|---|
| 1 |
function trm=velz(ffreqs,dirs,wns,z,depth,xpos,ypos) |
|---|
| 2 |
|
|---|
| 3 |
%Edits were made on 4/20/10 to fix a bug in the cutoff limits in the KZ |
|---|
| 4 |
%parts of the transfer function. |
|---|
| 5 |
|
|---|
| 6 |
%set the cutoff value (for high frequencies) |
|---|
| 7 |
cutoff=0.1; |
|---|
| 8 |
|
|---|
| 9 |
|
|---|
| 10 |
KZ=sinh(z*wns)./sinh(depth*wns); |
|---|
| 11 |
FLKZ=ffreqs.*KZ; |
|---|
| 12 |
|
|---|
| 13 |
%find the cutoff and replace the values below it |
|---|
| 14 |
cutindex=find(FLKZ < cutoff); |
|---|
| 15 |
FLKZ(cutindex)=cutoff; |
|---|
| 16 |
|
|---|
| 17 |
|
|---|
| 18 |
trm=-1i*(FLKZ)*ones(size(dirs)); |
|---|
