root/adcp/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.
