root/DPWavesProc/trunk/adcp/diwasp_1_1GD/private/velz.m
Revision 168 (checked in by cbc, 16 years ago) |
---|
Line | |
---|---|
1 | function trm=velz(ffreqs,dirs,wns,z,depth,xpos,ypos) |
2 | |
3 | Kz=sinh(z*wns)./sinh(depth*wns); |
4 | %include a maximum cuttoff for the velocity response function |
5 | Kz(find(Kz<0.1))=0.1; |
6 | Kz(find(isnan(Kz)))=1; |
7 | trm=-i*(ffreqs.*Kz)*ones(size(dirs)); |
Note: See TracBrowser for help on using the browser.