view freqCalcCyl4No_D4Mach.pd @ 0:ba6878e40772

Demo of Aeolian tone sound synthesis model
author Rod Selfridge <r.selfridge@qmul.ac.uk>
date Tue, 17 May 2016 09:50:06 +0100
parents
children 971da3e24318
line wrap: on
line source
#N canvas 159 62 1370 784 10;
#X text 582 246 Mass Density of Air;
#X obj 704 247 *~ 1.225;
#X obj 704 322 /~ 1.81e-05;
#X obj 704 277 *~ 1;
#X obj 704 440 /~ 1;
#X obj 704 416 *~;
#X obj 223 63 inlet;
#X text 737 366 Reynolds Number;
#X text 584 433 Strouhal Number;
#X obj 388 61 inlet;
#X text 734 302 Viscosity;
#X obj 223 84 s \$0-diameter;
#X obj 765 246 r \$0-diameter;
#X obj 772 421 r \$0-diameter;
#X obj 120 85 s~ \$0-windSpeed;
#X obj 771 398 r~ \$0-windSpeed;
#X obj 704 219 r~ \$0-windSpeed;
#X obj 704 464 s~ \$0-myFreq;
#X obj 1052 375 *~;
#X obj 717 346 s~ \$0-myRey;
#X obj 636 57 inlet;
#X text 668 502 Calculates the frequency;
#X text 669 515 based on the physics of;
#X obj 317 425 r~ \$0-windSpeed;
#X obj 357 487 r \$0-diameter;
#X obj 460 60 inlet;
#X obj 534 59 inlet;
#X text 293 367 Calculates the gain of the dipole source;
#X text 459 37 Azimuth;
#X text 641 39 Pan;
#X text 533 41 Elevation;
#X text 384 42 Distance;
#X text 218 42 Diameter;
#X text 133 41 Wind Speed;
#X obj 315 60 inlet;
#X text 311 41 Length;
#X obj 704 387 newStrou;
#X obj 1052 444 throw~ l;
#X obj 1079 422 throw~ r;
#X obj 1052 400 fcpan;
#X obj 315 81 s \$0-leng;
#X obj 388 82 s \$0-dist;
#X obj 400 331 s~ \$0-dragDiGain;
#X obj 313 352 s~ \$0-liftDiGain;
#X obj 304 406 r~ \$0-str;
#X obj 343 465 r \$0-leng;
#X obj 370 508 r \$0-dist;
#X obj 304 601 s~ \$0-liftNoiseGain;
#X obj 604 418 s~ \$0-str;
#X obj 636 126 s \$0-thisPan;
#X obj 1052 178 r~ \$0-liftNoiseGain;
#X obj 1079 375 r \$0-thisPan;
#X obj 731 66 inlet;
#X obj 803 65 inlet;
#X obj 877 64 inlet;
#X obj 731 89 s \$0-rotaX;
#X obj 803 87 s \$0-rotaY;
#X obj 877 87 s \$0-rotaZ;
#X text 727 47 X Rotation;
#X text 797 46 Y Rotation;
#X text 872 45 Z Rotation;
#X text 1033 152 Noise spectrum calculation;
#X obj 1067 323 hip~;
#X obj 1067 237 rpole~ 0.99;
#X obj 321 163 r~ \$0-windSpeed;
#X obj 339 202 r \$0-diameter;
#X obj 313 142 r~ \$0-str;
#X obj 331 183 r \$0-leng;
#X obj 356 241 r \$0-dist;
#N canvas 83 79 1238 721 GoldfreqGain 0;
#N canvas 759 125 286 294 pow6 0;
#X obj 47 67 *~;
#X obj 76 67 *~;
#X obj 76 35 inlet~;
#X obj 47 163 outlet~;
#X obj 136 184 snapshot~;
#X obj 136 159 metro 100;
#X obj 136 137 tgl 15 0 empty empty empty 17 7 0 10 -262144 -1 -1 1
1;
#X floatatom 136 211 10 0 0 0 - - -;
#X obj 216 210 outlet;
#X obj 136 111 loadbang;
#X obj 106 67 *~;
#X obj 47 93 *~;
#X obj 47 118 *~;
#X obj 372 202 *;
#X obj 372 228 *;
#X obj 382 22 hsl 128 15 0 10 0 0 empty empty empty -2 -8 0 10 -262144
-1 -1 0 1;
#X floatatom 374 288 10 0 0 0 - - -;
#X obj 303 20 t f b;
#X floatatom 455 55 5 0 0 0 - - -;
#X obj 372 260 *;
#X obj 400 202 *;
#X msg 303 -4 6;
#X obj 434 201 *;
#X connect 0 0 11 0;
#X connect 1 0 11 1;
#X connect 2 0 0 0;
#X connect 2 0 0 1;
#X connect 2 0 1 1;
#X connect 2 0 10 0;
#X connect 2 0 10 1;
#X connect 2 0 1 0;
#X connect 4 0 7 0;
#X connect 4 0 8 0;
#X connect 5 0 4 0;
#X connect 6 0 5 0;
#X connect 9 0 6 0;
#X connect 10 0 12 1;
#X connect 11 0 12 0;
#X connect 12 0 4 0;
#X connect 12 0 3 0;
#X connect 13 0 14 0;
#X connect 14 0 19 0;
#X connect 15 0 13 0;
#X connect 15 0 13 1;
#X connect 15 0 18 0;
#X connect 15 0 20 0;
#X connect 15 0 20 1;
#X connect 17 0 13 0;
#X connect 17 0 13 1;
#X connect 17 0 20 0;
#X connect 17 0 20 1;
#X connect 17 0 22 0;
#X connect 17 0 22 1;
#X connect 17 1 13 0;
#X connect 17 1 20 0;
#X connect 17 1 22 0;
#X connect 19 0 16 0;
#X connect 20 0 14 1;
#X connect 21 0 17 0;
#X connect 22 0 19 1;
#X restore 219 150 pd pow6;
#X obj 173 140 *~;
#X obj 219 176 *~;
#X obj 411 161 *;
#X text 587 13 Distance;
#X floatatom 567 73 5 0 0 0 - - -;
#X text 724 595 Calculates Drag Dipole gain;
#X text 48 531 Calculates Lift Dipole gain;
#X obj 173 100 inlet~;
#X obj 219 13 inlet~;
#X text 328 170 Diameter;
#X obj 291 169 inlet;
#X obj 567 29 inlet;
#X obj 219 590 outlet~;
#X obj 758 570 outlet~;
#X obj 240 263 inlet;
#X text 278 258 Length;
#X obj 219 289 *~ 1;
#X obj 234 469 /~ 360;
#X obj 234 518 *~;
#X obj 219 541 *~;
#X obj 868 30 inlet;
#X obj 967 32 inlet;
#X obj 1063 35 inlet;
#X floatatom 868 77 5 0 0 0 - - -;
#X floatatom 967 78 5 0 0 0 - - -;
#X floatatom 1064 79 5 0 0 0 - - -;
#X text 864 11 X Rotation;
#X text 961 13 Y Rotation;
#X text 1058 16 Z Rotation;
#X obj 873 53 s \$0-this_rotaX;
#X obj 974 54 s \$0-thisrotaY;
#X obj 1071 58 s \$0-thisrotaZ;
#X text 74 100 Strouhal number;
#X text 149 37 Wind speed;
#X text 266 211 Corr Length;
#X obj 219 348 *~ 2.50663;
#X obj 219 238 *~ 1.225;
#X text 93 234 Mass Density of Air;
#X text 145 343 Sqrt(2xpi);
#X obj 306 495 /~ 360;
#X obj 306 541 *~;
#X obj 306 517 cos~;
#X obj 219 567 *~;
#X obj 411 191 * 4.05513e+07;
#X text 362 210 Speed of Sound to power 3;
#X obj 411 244 * 32;
#X obj 525 213 cos~;
#X obj 525 236 *~;
#X obj 473 237 sig~ 1;
#X obj 510 267 -~;
#N canvas 1150 227 286 289 pow4 0;
#X obj 47 67 *~;
#X obj 76 67 *~;
#X obj 76 35 inlet~;
#X obj 47 163 outlet~;
#X obj 136 184 snapshot~;
#X obj 136 159 metro 100;
#X obj 136 137 tgl 15 0 empty empty empty 17 7 0 10 -262144 -1 -1 1
1;
#X floatatom 136 211 10 0 0 0 - - -;
#X obj 216 210 outlet;
#X obj 136 111 loadbang;
#X obj 47 93 *~;
#X connect 0 0 10 0;
#X connect 1 0 10 1;
#X connect 2 0 0 0;
#X connect 2 0 0 1;
#X connect 2 0 1 1;
#X connect 2 0 1 0;
#X connect 4 0 7 0;
#X connect 4 0 8 0;
#X connect 5 0 4 0;
#X connect 6 0 5 0;
#X connect 9 0 6 0;
#X connect 10 0 3 0;
#X connect 10 0 4 0;
#X restore 510 290 pd pow4;
#X obj 411 272 sig~;
#X obj 411 322 *~;
#X obj 219 385 /~;
#X obj 133 375 s~ \$0-aGainG;
#X obj 872 246 r~ \$0-aGainG;
#X obj 758 512 *~ 0.1;
#X obj 723 340 +~ 90;
#X obj 723 363 /~ 360;
#X obj 723 409 *~;
#X obj 723 385 sin~;
#X obj 773 359 /~ 360;
#X obj 773 405 *~;
#X obj 773 381 cos~;
#X obj 758 436 *~;
#X text 763 12 Azimuth;
#X text 662 14 Elevation;
#X obj 665 32 inlet;
#X obj 767 30 inlet;
#X floatatom 665 79 5 0 0 0 - - -;
#X floatatom 767 79 5 0 0 0 - - -;
#X obj 773 294 r \$0-this_azim;
#X obj 670 55 s \$0-this_elev;
#X obj 774 55 s \$0-this_azim;
#X obj 525 129 r \$0-this_elev;
#X obj 574 51 s \$0-thisDist;
#X obj 887 268 r \$0-thisDist;
#X obj 411 135 r \$0-thisDist;
#X obj 887 292 *;
#X obj 887 317 * 4.05513e+07;
#X text 892 335 Speed of Sound to power 3;
#X obj 887 370 * 32;
#X obj 887 398 sig~;
#X obj 872 423 /~;
#X text 595 191 Remember - Goldstein moves cylinder;
#X obj 525 194 /~ 360;
#X obj 758 475 *~;
#X obj 478 51 s~ \$0-thisColL;
#X obj 270 150 r~ \$0-thisColL;
#X obj 270 190 *~ 1;
#X obj 219 210 *~;
#X obj 478 28 inlet~;
#X floatatom 268 86 5 0 0 0 - - -;
#X obj 234 495 sin~;
#X floatatom 278 34 5 0 0 0 - - -;
#N canvas 921 315 426 401 aveTimeCal 0;
#X text 224 38 Elevation;
#X obj 226 56 r \$0-elev;
#X floatatom 240 81 5 0 0 0 - - -;
#X obj 205 17 r~ MY_MACH;
#X obj 226 135 cos;
#X obj 226 107 * 0.01745;
#X obj 205 158 *~ 1;
#X obj 161 131 sig~ 1;
#X obj 161 185 -~;
#X obj 121 104 r~ \$0-myFreq;
#X obj 106 82 sig~ 6.28318;
#X obj 106 128 /~;
#X obj 106 219 *~;
#X obj 36 204 r myMet;
#X obj 106 242 snapshot~;
#X floatatom 36 293 0 0 0 0 - - -;
#X obj 106 263 * 44100;
#X obj 129 287 outlet;
#X obj 235 200 r myMet;
#X obj 235 221 snapshot~;
#X floatatom 235 245 0 0 0 0 - - -;
#X connect 1 0 2 0;
#X connect 1 0 5 0;
#X connect 3 0 6 0;
#X connect 4 0 6 1;
#X connect 5 0 4 0;
#X connect 6 0 8 1;
#X connect 7 0 8 0;
#X connect 8 0 12 1;
#X connect 9 0 11 1;
#X connect 10 0 11 0;
#X connect 11 0 12 0;
#X connect 12 0 14 0;
#X connect 13 0 14 0;
#X connect 14 0 16 0;
#X connect 16 0 15 0;
#X connect 16 0 17 0;
#X connect 18 0 19 0;
#X connect 19 0 20 0;
#X restore 268 15 pd aveTimeCal;
#X obj 219 59 average~ 1000;
#X obj 219 122 line~;
#X obj 525 173 line~;
#X obj 723 273 r \$0-this_elev;
#X obj 723 317 line~;
#X msg 723 295 \$1 5;
#X msg 525 151 \$1 5;
#X obj 773 338 line~;
#X msg 773 316 \$1 5;
#X obj 306 427 r \$0-this_azim;
#X obj 234 401 r \$0-this_elev;
#X obj 234 445 line~;
#X msg 234 423 \$1 5;
#X obj 306 471 line~;
#X msg 306 449 \$1 5;
#X obj 219 317 /~ 400;
#X text -53 316 Scale factor to bring in line with CFD calcs;
#X obj 685 572 s~ plotOPD;
#X msg 219 100 \$1 15;
#X obj 565 217 r~ MY_MACH;
#X connect 0 0 2 0;
#X connect 1 0 2 1;
#X connect 2 0 91 0;
#X connect 3 0 44 0;
#X connect 8 0 1 0;
#X connect 8 0 1 1;
#X connect 9 0 97 0;
#X connect 11 0 90 1;
#X connect 12 0 5 0;
#X connect 12 0 76 0;
#X connect 15 0 17 1;
#X connect 17 0 112 0;
#X connect 18 0 94 0;
#X connect 19 0 20 1;
#X connect 20 0 43 0;
#X connect 21 0 24 0;
#X connect 21 0 30 0;
#X connect 22 0 25 0;
#X connect 22 0 31 0;
#X connect 23 0 26 0;
#X connect 23 0 32 0;
#X connect 36 0 54 0;
#X connect 36 0 55 0;
#X connect 37 0 17 0;
#X connect 40 0 42 0;
#X connect 41 0 43 1;
#X connect 42 0 41 0;
#X connect 42 0 41 1;
#X connect 43 0 13 0;
#X connect 44 0 46 0;
#X connect 46 0 52 0;
#X connect 47 0 48 0;
#X connect 48 0 50 1;
#X connect 49 0 50 0;
#X connect 50 0 51 0;
#X connect 51 0 53 1;
#X connect 52 0 53 0;
#X connect 53 0 54 1;
#X connect 54 0 20 0;
#X connect 56 0 84 0;
#X connect 57 0 14 0;
#X connect 57 0 114 0;
#X connect 58 0 59 0;
#X connect 59 0 61 0;
#X connect 60 0 65 0;
#X connect 61 0 60 0;
#X connect 61 0 60 1;
#X connect 62 0 64 0;
#X connect 63 0 65 1;
#X connect 64 0 63 0;
#X connect 64 0 63 1;
#X connect 65 0 87 0;
#X connect 68 0 70 0;
#X connect 68 0 73 0;
#X connect 69 0 71 0;
#X connect 69 0 74 0;
#X connect 72 0 105 0;
#X connect 75 0 103 0;
#X connect 77 0 79 0;
#X connect 77 0 79 1;
#X connect 78 0 3 0;
#X connect 78 0 3 1;
#X connect 79 0 80 0;
#X connect 80 0 82 0;
#X connect 82 0 83 0;
#X connect 83 0 84 1;
#X connect 84 0 87 1;
#X connect 86 0 47 0;
#X connect 87 0 57 0;
#X connect 89 0 90 0;
#X connect 90 0 91 1;
#X connect 91 0 37 0;
#X connect 92 0 88 0;
#X connect 94 0 19 0;
#X connect 94 0 19 1;
#X connect 96 0 95 0;
#X connect 96 0 97 0;
#X connect 97 0 93 0;
#X connect 97 0 115 0;
#X connect 98 0 0 0;
#X connect 99 0 86 0;
#X connect 100 0 102 0;
#X connect 101 0 58 0;
#X connect 102 0 101 0;
#X connect 103 0 99 0;
#X connect 104 0 62 0;
#X connect 105 0 104 0;
#X connect 106 0 111 0;
#X connect 107 0 109 0;
#X connect 108 0 18 0;
#X connect 109 0 108 0;
#X connect 110 0 40 0;
#X connect 111 0 110 0;
#X connect 112 0 36 0;
#X connect 115 0 98 0;
#X connect 116 0 48 1;
#X restore 313 307 pd GoldfreqGain;
#N canvas 300 30 1006 797 GoldnoiseGain 0;
#X text -21037 364 Mass Density of Air \, squared;
#X obj -20777 459 *;
#X text -20777 412 Distance;
#N canvas 759 125 235 275 pow8 0;
#X obj 47 67 *~;
#X obj 76 67 *~;
#X obj 76 35 inlet~;
#X obj 47 147 outlet~;
#X obj 128 146 snapshot~;
#X obj 128 121 metro 100;
#X obj 128 99 tgl 15 0 empty empty empty 17 7 0 10 -262144 -1 -1 1
1;
#X floatatom 128 173 10 0 0 0 - - -;
#X obj 128 73 loadbang;
#X obj 47 93 *~;
#X obj 47 121 *~;
#X connect 0 0 9 0;
#X connect 1 0 9 1;
#X connect 2 0 0 0;
#X connect 2 0 0 1;
#X connect 2 0 1 1;
#X connect 2 0 1 0;
#X connect 4 0 7 0;
#X connect 5 0 4 0;
#X connect 6 0 5 0;
#X connect 8 0 6 0;
#X connect 9 0 10 0;
#X connect 9 0 10 1;
#X connect 10 0 3 0;
#X connect 10 0 4 0;
#X restore -20846 224 pd pow8;
#X obj -20846 137 inlet~;
#X obj -20807 286 inlet;
#X obj -20777 433 inlet;
#X obj -20861 617 outlet~;
#X text -20858 648 Gain proportional to windspeed to power 8 for spectrum
;
#X text -20771 287 Diameter;
#X text -20847 118 Wind speed;
#X text -21027 627 Calculates wake gain;
#X text -20399 16 Azimuth;
#X obj -20398 33 inlet;
#X obj -20342 33 inlet;
#X text -20343 16 Elevation;
#X floatatom -20304 57 5 0 0 0 - - -;
#X floatatom -20470 58 5 0 0 0 - - -;
#X obj -20398 117 s~ \$0-this_azim;
#X obj -20342 95 s~ \$0-this_elev;
#X obj -20861 371 *~ 1.225;
#X obj -20587 239 r~ \$0-this_azim;
#N canvas 759 125 250 309 pow5 0;
#X obj 47 67 *~;
#X obj 76 35 inlet~;
#X obj 47 170 outlet~;
#X obj 136 173 snapshot~;
#X obj 136 148 metro 100;
#X obj 136 126 tgl 15 0 empty empty empty 17 7 0 10 -262144 -1 -1 1
1;
#X floatatom 136 200 10 0 0 0 - - -;
#X obj 136 100 loadbang;
#X obj 47 87 *~;
#X obj 47 118 *~;
#X connect 0 0 8 0;
#X connect 0 0 8 1;
#X connect 1 0 0 0;
#X connect 1 0 0 1;
#X connect 1 0 9 1;
#X connect 3 0 6 0;
#X connect 4 0 3 0;
#X connect 5 0 4 0;
#X connect 7 0 5 0;
#X connect 8 0 9 0;
#X connect 9 0 2 0;
#X connect 9 0 3 0;
#X restore -20705 475 pd pow5;
#X obj -20705 451 sig~ 343.56;
#N canvas 875 261 250 309 pow5 0;
#X obj 47 67 *~;
#X obj 76 35 inlet~;
#X obj 47 170 outlet~;
#X obj 136 173 snapshot~;
#X obj 136 148 metro 100;
#X obj 136 126 tgl 15 0 empty empty empty 17 7 0 10 -262144 -1 -1 1
1;
#X floatatom 136 200 10 0 0 0 - - -;
#X obj 136 100 loadbang;
#X obj 47 87 *~;
#X obj 47 118 *~;
#X connect 0 0 8 0;
#X connect 0 0 8 1;
#X connect 1 0 0 0;
#X connect 1 0 0 1;
#X connect 1 0 9 1;
#X connect 3 0 6 0;
#X connect 4 0 3 0;
#X connect 5 0 4 0;
#X connect 7 0 5 0;
#X connect 8 0 9 0;
#X connect 9 0 2 0;
#X connect 9 0 3 0;
#X restore -20670 401 pd pow5;
#X obj -20624 216 /~ 360;
#X obj -20624 243 cos~;
#X obj -20624 349 *~;
#X obj -20670 348 sig~ 1;
#X obj -20670 375 -~;
#X obj -20705 322 sig~ 3.1416;
#X obj -20705 347 *~;
#X obj -20720 299 sig~ 16;
#X obj -20720 371 *~;
#X obj -20720 429 *~;
#X obj -20720 502 *~;
#X obj -20777 530 *~;
#X obj -20861 574 /~;
#X obj -20624 194 -~;
#X text -20543 190 Angle goes with the flow in Goldstein hence 90 -
theta so that theta is same angle as dipole angles.;
#X obj -20624 149 sig~ 180;
#X obj -20609 172 r~ \$0-this_elev;
#X obj -20916 137 inlet~;
#X obj -20916 164 *~;
#X obj -20828 253 inlet~;
#X obj -20828 312 *~ 1;
#X obj -20861 251 *~;
#X obj -20861 347 *~;
#X obj -20809 392 inlet;
#X obj -20861 444 *~ 1;
#X obj -20861 528 *~ 2.50663;
#X text -20782 254 Corrolation Length;
#X text -20806 371 Length;
#X obj -20846 159 average~ 64;
#X floatatom -20757 165 5 0 0 0 - - -;
#N canvas 801 264 426 401 aveTimeCal 0;
#X text 224 38 Elevation;
#X obj 226 56 r \$0-elev;
#X floatatom 240 81 5 0 0 0 - - -;
#X obj 205 17 r~ MY_MACH;
#X obj 226 135 cos;
#X obj 226 107 * 0.01745;
#X obj 205 158 *~ 1;
#X obj 161 131 sig~ 1;
#X obj 161 185 -~;
#X obj 121 104 r~ \$0-myFreq;
#X obj 106 82 sig~ 6.28318;
#X obj 106 128 /~;
#X obj 106 219 *~;
#X obj 36 204 r myMet;
#X obj 106 242 snapshot~;
#X floatatom 44 283 0 0 0 0 - - -;
#X obj 106 263 * 44100;
#X obj 106 292 outlet;
#X connect 1 0 2 0;
#X connect 1 0 5 0;
#X connect 3 0 6 0;
#X connect 4 0 6 1;
#X connect 5 0 4 0;
#X connect 6 0 8 1;
#X connect 7 0 8 0;
#X connect 8 0 12 1;
#X connect 9 0 11 1;
#X connect 10 0 11 0;
#X connect 11 0 12 0;
#X connect 12 0 14 0;
#X connect 13 0 14 0;
#X connect 14 0 16 0;
#X connect 16 0 15 0;
#X connect 16 0 17 0;
#X restore -20766 139 pd aveTimeCal;
#X obj -20846 203 line~;
#X msg -20846 181 \$1 5;
#X obj -20342 75 line~;
#X msg -20342 53 \$1 5;
#X obj -20777 502 line~;
#X msg -20777 480 \$1 5;
#X obj -20398 86 line~;
#X msg -20398 64 \$1 5;
#X obj -20587 285 cos~;
#X obj -20861 595 *~ 1;
#X obj -20587 262 /~ 360;
#X obj -20624 306 *~;
#X text -21131 486 Scale factor to bring in line with CFD calcs;
#X text -20929 530 Sqrt(2*pi);
#X obj -20861 488 /~ 10000;
#X obj -20609 326 r~ MY_MACH;
#X connect 1 0 61 0;
#X connect 3 0 46 1;
#X connect 4 0 53 0;
#X connect 5 0 45 1;
#X connect 6 0 1 0;
#X connect 6 0 1 1;
#X connect 13 0 17 0;
#X connect 13 0 63 0;
#X connect 14 0 16 0;
#X connect 14 0 59 0;
#X connect 20 0 49 0;
#X connect 21 0 66 0;
#X connect 22 0 35 1;
#X connect 23 0 22 0;
#X connect 24 0 34 1;
#X connect 25 0 26 0;
#X connect 26 0 67 0;
#X connect 27 0 29 1;
#X connect 28 0 29 0;
#X connect 29 0 24 0;
#X connect 30 0 31 0;
#X connect 30 0 31 1;
#X connect 31 0 33 1;
#X connect 32 0 33 0;
#X connect 33 0 34 0;
#X connect 34 0 35 0;
#X connect 35 0 36 1;
#X connect 36 0 37 1;
#X connect 37 0 65 0;
#X connect 38 0 25 0;
#X connect 40 0 38 0;
#X connect 41 0 38 1;
#X connect 42 0 43 0;
#X connect 42 0 43 1;
#X connect 43 0 46 0;
#X connect 44 0 45 0;
#X connect 45 0 47 1;
#X connect 46 0 47 0;
#X connect 47 0 20 0;
#X connect 48 0 49 1;
#X connect 49 0 70 0;
#X connect 50 0 37 0;
#X connect 53 0 57 0;
#X connect 55 0 54 0;
#X connect 55 0 53 0;
#X connect 56 0 3 0;
#X connect 57 0 56 0;
#X connect 58 0 19 0;
#X connect 59 0 58 0;
#X connect 60 0 36 0;
#X connect 61 0 60 0;
#X connect 62 0 18 0;
#X connect 63 0 62 0;
#X connect 64 0 67 1;
#X connect 65 0 7 0;
#X connect 66 0 64 0;
#X connect 67 0 27 0;
#X connect 70 0 50 0;
#X connect 71 0 27 1;
#X restore 304 573 pd GoldnoiseGain;
#X text 295 620 Calculates the gain of the wake;
#X obj 460 80 s \$0-azim;
#X obj 534 80 s \$0-elev;
#X obj 365 262 r \$0-elev;
#X obj 373 283 r \$0-azim;
#X obj 397 552 r \$0-elev;
#X obj 383 529 r \$0-azim;
#N canvas 743 232 444 428 corrLen 0;
#X obj 103 71 inlet~;
#X obj 88 200 outlet~;
#X obj 88 143 *~;
#X obj 124 93 sig~ -0.245;
#X obj 103 119 pow~;
#X text 234 22 Y = 10^b . X^m;
#X text 235 40 b = 1.536;
#X text 235 56 m = -0.245;
#X obj 88 169 max~ 1;
#X obj 88 46 sig~ 44.67;
#X obj 229 190 r myMet;
#X obj 229 211 snapshot~;
#X floatatom 230 236 5 0 0 0 - - -;
#X text 56 298 This patch calculates the correlation length associated
with the dipole source. Equation of line derived from Norberg 2001
;
#X connect 0 0 4 0;
#X connect 2 0 8 0;
#X connect 3 0 4 1;
#X connect 4 0 2 1;
#X connect 8 0 1 0;
#X connect 8 0 11 0;
#X connect 9 0 2 0;
#X connect 10 0 11 0;
#X connect 11 0 12 0;
#X restore 956 85 pd corrLen Calc;
#X obj 956 62 r~ \$0-myRey;
#X obj 956 109 s~ \$0-myCor;
#X obj 347 222 r~ \$0-myCor;
#X obj 1067 350 hip~;
#X obj 330 445 r~ \$0-myCor;
#X obj 636 81 max 0.01;
#X obj 636 104 min 0.99;
#N canvas 93 119 976 679 Dipoles 0;
#X obj 118 66 r~ \$0-myFreq;
#X obj 397 51 r~ \$0-myFreq;
#X obj 668 44 r~ \$0-myFreq;
#X obj 816 126 r~ \$0-liftDiGain;
#X obj 668 68 *~ 5;
#X obj 397 75 *~ 3;
#X obj 139 407 r~ \$0-myFreq;
#X obj 253 492 r~ \$0-dragDiGain;
#X obj 139 435 *~ 2;
#X obj 409 406 r~ \$0-myFreq;
#X obj 550 464 r~ \$0-dragDiGain;
#X obj 409 434 *~ 4;
#X obj 550 487 *~ 0.125;
#X obj 210 183 *~;
#X obj 210 248 throw~ l;
#X obj 237 225 throw~ r;
#X obj 210 204 fcpan;
#X obj 210 108 vcf~ 3;
#X obj 237 182 r \$0-thisPan;
#X obj 210 67 r~ \$0-thisNoise;
#X obj 480 192 *~;
#X obj 480 257 throw~ l;
#X obj 507 234 throw~ r;
#X obj 480 213 fcpan;
#X obj 480 116 vcf~ 3;
#X obj 507 191 r \$0-thisPan;
#X obj 480 48 r~ \$0-thisNoise;
#X obj 750 183 *~;
#X obj 750 248 throw~ l;
#X obj 777 225 throw~ r;
#X obj 750 204 fcpan;
#X obj 750 107 vcf~ 3;
#X obj 777 182 r \$0-thisPan;
#X obj 750 38 r~ \$0-thisNoise;
#X obj 218 525 *~;
#X obj 218 590 throw~ l;
#X obj 245 567 throw~ r;
#X obj 218 546 fcpan;
#X obj 218 449 vcf~ 3;
#X obj 245 524 r \$0-thisPan;
#X obj 218 405 r~ \$0-thisNoise;
#X obj 492 523 *~;
#X obj 492 588 throw~ l;
#X obj 519 565 throw~ r;
#X obj 492 544 fcpan;
#X obj 492 447 vcf~ 3;
#X obj 519 522 r \$0-thisPan;
#X obj 492 373 r~ \$0-thisNoise;
#X text 273 27 Generates the sound associated with the lift dipole.
;
#X text 176 373 Generates the sound associated with the drag dipole.
;
#X obj 270 87 r \$0-Qval;
#X obj 523 72 r \$0-Qval;
#X obj 792 62 r \$0-Qval;
#X obj 538 394 r \$0-Qval;
#X obj 263 427 r \$0-Qval;
#X floatatom 277 122 5 0 0 0 - - -;
#X obj 210 134 vcf~ 3;
#X obj 480 138 vcf~ 3;
#X obj 750 131 vcf~ 3;
#X obj 492 468 vcf~ 3;
#X obj 218 470 vcf~ 3;
#X obj 527 143 r~ \$0-liftDiGain;
#X obj 225 156 r~ \$0-liftDiGain;
#X obj 816 148 *~ 0.1;
#X obj 523 95 * 2;
#X obj 792 84 * 4;
#X obj 533 421 * 2;
#X obj 527 165 *~ 0.6;
#X text 179 282 Fundamental;
#X text 444 288 3rd Harmonic;
#X text 729 288 5th Harmonic;
#X text 151 618 Fundamental Drag tone = 2nd Harmonic;
#X text 420 616 2 Harmonic Drag tone = 4th Harmonic;
#X connect 0 0 17 1;
#X connect 0 0 56 1;
#X connect 1 0 5 0;
#X connect 2 0 4 0;
#X connect 3 0 63 0;
#X connect 4 0 31 1;
#X connect 4 0 58 1;
#X connect 5 0 24 1;
#X connect 5 0 57 1;
#X connect 6 0 8 0;
#X connect 7 0 34 1;
#X connect 8 0 38 1;
#X connect 8 0 60 1;
#X connect 9 0 11 0;
#X connect 10 0 12 0;
#X connect 11 0 45 1;
#X connect 11 0 59 1;
#X connect 12 0 41 1;
#X connect 13 0 16 0;
#X connect 16 0 14 0;
#X connect 16 1 15 0;
#X connect 17 0 56 0;
#X connect 18 0 16 1;
#X connect 19 0 17 0;
#X connect 20 0 23 0;
#X connect 23 0 21 0;
#X connect 23 1 22 0;
#X connect 24 0 57 0;
#X connect 25 0 23 1;
#X connect 26 0 24 0;
#X connect 27 0 30 0;
#X connect 30 0 28 0;
#X connect 30 1 29 0;
#X connect 31 0 58 0;
#X connect 32 0 30 1;
#X connect 33 0 31 0;
#X connect 34 0 37 0;
#X connect 37 0 35 0;
#X connect 37 1 36 0;
#X connect 38 0 60 0;
#X connect 39 0 37 1;
#X connect 40 0 38 0;
#X connect 41 0 44 0;
#X connect 44 0 42 0;
#X connect 44 1 43 0;
#X connect 45 0 59 0;
#X connect 46 0 44 1;
#X connect 47 0 45 0;
#X connect 50 0 17 2;
#X connect 50 0 55 0;
#X connect 50 0 56 2;
#X connect 51 0 64 0;
#X connect 52 0 65 0;
#X connect 53 0 66 0;
#X connect 54 0 38 2;
#X connect 54 0 60 2;
#X connect 56 0 13 0;
#X connect 57 0 20 0;
#X connect 58 0 27 0;
#X connect 59 0 41 0;
#X connect 60 0 34 0;
#X connect 61 0 67 0;
#X connect 62 0 13 1;
#X connect 63 0 27 1;
#X connect 64 0 24 2;
#X connect 64 0 57 2;
#X connect 65 0 31 2;
#X connect 65 0 58 2;
#X connect 66 0 45 2;
#X connect 66 0 59 2;
#X connect 67 0 20 1;
#X restore 660 614 pd Dipoles;
#X obj 1088 69 noise~;
#X obj 1088 92 s~ \$0-thisNoise;
#X obj 1067 214 r~ \$0-thisNoise;
#X text 669 528 the object and environment;
#X obj 1107 280 r myMet;
#X obj 1098 302 snapshot~;
#N canvas 308 58 1027 635 StBandwidth 0;
#X obj 23 70 inlet;
#X obj 460 365 swap 1;
#X obj 460 391 /;
#X floatatom 534 431 0 0 0 0 - - -;
#X obj 460 342 / 100;
#X obj 461 434 outlet;
#X text 437 453 Q value;
#X obj 531 -36 r~ \$0-myRey;
#X obj 477 -36 r myMet;
#X obj 477 -14 snapshot~;
#X obj 203 164 spigot;
#X obj 219 187 spigot;
#X obj 236 140 > 0;
#X obj 219 233 + 0.9797;
#X floatatom 194 270 0 0 0 0 - - -;
#X obj 219 208 * 4.624e-05;
#N canvas 106 95 810 724 linearCalc 0;
#X obj 240 36 inlet;
#X obj 248 579 outlet;
#X floatatom 263 126 0 0 0 0 - - -;
#X obj 240 59 * 38.7462;
#X obj 240 81 + 0.6254;
#X text 223 3 Peak dependant on TU value;
#X obj 339 146 inlet;
#X obj 215 435 < 193260;
#X obj 182 458 spigot;
#X obj 308 457 spigot;
#X obj 341 433 > 193260;
#X obj 125 164 - 7.685;
#X obj 125 187 / 48293;
#X obj 140 210 - 145021;
#X obj 125 255 *;
#X obj 125 278 + 7.685;
#X floatatom 81 327 0 0 0 0 - - -;
#X floatatom 123 486 0 0 0 0 - - -;
#X obj 426 221 - 193260;
#X obj 411 150 swap 2.95;
#X obj 411 174 -;
#X obj 411 197 / 50060;
#X obj 411 267 *;
#X obj 140 232 t f b;
#X obj 411 302 +;
#X obj 370 273 t f b;
#X floatatom 409 328 0 0 0 0 - - -;
#X obj 426 242 t f b;
#X floatatom 405 530 0 0 0 0 - - -;
#X obj 240 102 min 200;
#X text 548 50 Calculates the peak value of DeltaF / F depending on
the TU values. The higher the TU the higher the DeltaF/F value. A limit
of 200 has been made as the value is a % and bandwidth greater than
200% would invlove -ve frequencies;
#X text -9 102 Linear up slope of the second section;
#X text 478 174 Linear down slope of the second section;
#X text 391 391 Peak of the DeltaF / F was found to be at a Reynolds
number of 193260;
#X obj 559 256 hsl 128 15 0 1e+06 0 0 empty empty empty -2 -8 0 10
-262144 -1 -1 0 1;
#X obj 65 15 hsl 128 15 0 20 0 0 empty empty empty -2 -8 0 10 -262144
-1 -1 0 1;
#X floatatom 566 289 0 0 0 0 - - -;
#X floatatom 66 55 0 0 0 0 - - -;
#X obj 379 460 < 243320;
#X obj 346 484 spigot;
#X obj 200 483 spigot;
#X obj 233 460 > 145021;
#X obj 392 60 *;
#X obj 403 88 +;
#X floatatom 463 120 0 0 0 0 - - -;
#X obj 495 26 * 0.05;
#X obj 424 61 * 7.89;
#X obj 452 95 + 21.64;
#X connect 3 0 4 0;
#X connect 4 0 29 0;
#X connect 7 0 8 1;
#X connect 8 0 40 0;
#X connect 9 0 39 0;
#X connect 10 0 9 1;
#X connect 11 0 12 0;
#X connect 12 0 14 0;
#X connect 13 0 23 0;
#X connect 14 0 15 0;
#X connect 15 0 16 0;
#X connect 15 0 8 0;
#X connect 18 0 27 0;
#X connect 19 0 20 0;
#X connect 19 1 20 1;
#X connect 20 0 21 0;
#X connect 21 0 22 0;
#X connect 22 0 24 0;
#X connect 23 0 14 1;
#X connect 23 1 14 0;
#X connect 24 0 26 0;
#X connect 24 0 9 0;
#X connect 25 0 24 1;
#X connect 25 1 24 0;
#X connect 27 0 22 1;
#X connect 27 1 22 0;
#X connect 29 0 2 0;
#X connect 29 0 19 0;
#X connect 29 0 25 0;
#X connect 29 0 11 0;
#X connect 34 0 10 0;
#X connect 34 0 18 0;
#X connect 34 0 13 0;
#X connect 34 0 7 0;
#X connect 34 0 36 0;
#X connect 34 0 38 0;
#X connect 34 0 41 0;
#X connect 35 0 3 0;
#X connect 35 0 37 0;
#X connect 35 0 42 1;
#X connect 35 0 42 0;
#X connect 35 0 46 0;
#X connect 38 0 39 1;
#X connect 39 0 1 0;
#X connect 39 0 28 0;
#X connect 40 0 1 0;
#X connect 40 0 17 0;
#X connect 41 0 40 1;
#X connect 42 0 45 0;
#X connect 43 0 47 0;
#X connect 45 0 43 0;
#X connect 46 0 43 1;
#X connect 47 0 44 0;
#X restore 23 92 pd linearCalc;
#X floatatom -37 134 0 0 0 0 - - -;
#X obj 56 116 > 145021;
#X obj 72 139 < 243320;
#X obj 23 140 spigot;
#X obj 39 163 spigot;
#X obj 602 150 spigot;
#X obj 618 173 spigot;
#X obj 651 149 < 1e+06;
#X obj 618 216 * 1.227e-10;
#X obj 618 196 *;
#X obj 690 215 * -8.553e-05;
#X obj 618 261 + 16.5;
#X obj 618 238 +;
#X floatatom 640 291 0 0 0 0 - - -;
#X floatatom 565 18 0 0 0 0 - - -;
#X floatatom 54 188 0 0 0 0 - - -;
#X text -84 298 First section of the Reynolds number - DeltaF/F relationship
is linear and did not appear to be reliant on current TU values;
#X text 565 329 Third section of the Reynolds number - DeltaF/F relationship
is linear and did not appear to be reliant on current TU values;
#X text 114 -16 Middle section of Reynolds number - DeltaF/F relationship
is also linear but depends on TU Value;
#X obj 461 412 max 0.005;
#X floatatom 402 435 0 0 0 0 - - -;
#X obj 707 42 > 243320;
#X obj 712 97 > 145021;
#X obj 635 126 > 193260;
#X obj 252 163 < 193260;
#X text 184 511 This patch calculates the Q values to be applied to
the bandpass filters. Future development may link this to the Turbulence
Intensity but not enough data available at time of writing to do this.
;
#X connect 0 0 16 0;
#X connect 1 0 2 0;
#X connect 1 1 2 1;
#X connect 2 0 3 0;
#X connect 2 0 36 0;
#X connect 4 0 1 0;
#X connect 7 0 9 0;
#X connect 8 0 9 0;
#X connect 9 0 24 0;
#X connect 9 0 22 0;
#X connect 9 0 12 0;
#X connect 9 0 10 0;
#X connect 9 0 31 0;
#X connect 9 0 40 0;
#X connect 9 0 41 0;
#X connect 10 0 11 0;
#X connect 11 0 15 0;
#X connect 12 0 10 1;
#X connect 13 0 14 0;
#X connect 13 0 4 0;
#X connect 15 0 13 0;
#X connect 16 0 17 0;
#X connect 16 0 20 0;
#X connect 18 0 20 1;
#X connect 19 0 21 1;
#X connect 20 0 21 0;
#X connect 21 0 32 0;
#X connect 22 0 23 0;
#X connect 23 0 27 0;
#X connect 23 0 26 1;
#X connect 23 0 26 0;
#X connect 24 0 23 1;
#X connect 25 0 29 0;
#X connect 26 0 25 0;
#X connect 27 0 29 1;
#X connect 28 0 30 0;
#X connect 28 0 4 0;
#X connect 29 0 28 0;
#X connect 36 0 5 0;
#X connect 36 0 37 0;
#X connect 40 0 22 1;
#X connect 41 0 11 1;
#X restore 1189 91 pd StBandwidth;
#X obj 1189 67 inlet;
#X obj 1189 114 s \$0-Qval;
#X obj 1098 260 r~ \$0-myFreq;
#X obj 120 65 vd~ varwind;
#X obj 60 106 s~ MY_MACH;
#X obj 817 469 outlet~;
#X obj 872 468 outlet~;
#X obj 931 470 outlet~;
#X obj 60 84 /~ 343.56;
#X connect 1 0 3 0;
#X connect 2 0 19 0;
#X connect 2 0 36 0;
#X connect 2 0 100 0;
#X connect 3 0 2 0;
#X connect 4 0 17 0;
#X connect 4 0 99 0;
#X connect 5 0 4 0;
#X connect 6 0 11 0;
#X connect 9 0 41 0;
#X connect 12 0 3 1;
#X connect 13 0 4 1;
#X connect 15 0 5 1;
#X connect 16 0 1 0;
#X connect 18 0 39 0;
#X connect 20 0 84 0;
#X connect 23 0 70 1;
#X connect 24 0 70 4;
#X connect 25 0 72 0;
#X connect 26 0 73 0;
#X connect 34 0 40 0;
#X connect 36 0 48 0;
#X connect 36 0 5 0;
#X connect 36 0 101 0;
#X connect 39 0 37 0;
#X connect 39 1 38 0;
#X connect 44 0 70 0;
#X connect 45 0 70 3;
#X connect 46 0 70 5;
#X connect 50 0 18 0;
#X connect 51 0 39 1;
#X connect 52 0 55 0;
#X connect 53 0 56 0;
#X connect 54 0 57 0;
#X connect 62 0 82 0;
#X connect 63 0 62 0;
#X connect 64 0 69 1;
#X connect 65 0 69 3;
#X connect 66 0 69 0;
#X connect 67 0 69 2;
#X connect 68 0 69 5;
#X connect 69 0 43 0;
#X connect 69 1 42 0;
#X connect 70 0 47 0;
#X connect 74 0 69 6;
#X connect 75 0 69 7;
#X connect 76 0 70 7;
#X connect 77 0 70 6;
#X connect 78 0 80 0;
#X connect 79 0 78 0;
#X connect 81 0 69 4;
#X connect 82 0 18 1;
#X connect 83 0 70 2;
#X connect 84 0 85 0;
#X connect 85 0 49 0;
#X connect 87 0 88 0;
#X connect 89 0 63 0;
#X connect 91 0 92 0;
#X connect 92 0 62 1;
#X connect 92 0 82 1;
#X connect 93 0 95 0;
#X connect 94 0 93 0;
#X connect 96 0 92 0;
#X connect 97 0 14 0;
#X connect 97 0 102 0;
#X connect 102 0 98 0;