• API
• FAQ
• Tools
• Archive
SHARE
TWEET # Untitled a guest Jul 5th, 2017 87 Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
1. pro h5p1
2.
3. r=2.5 ;microns
4. uv = 4 ;m/s
5. T = 273 ;K
6. supersat = 1.002 ;2%
7. LWC = 1 ;g/m^3
8. Fk = 2.4E-2 ;J/msK
9. Fd = 2.21E-5 ;m^2/s
10. mu = 1.717E-5 ;kg/ms
11. E = 1
12. wd = 1.0 ;kg/m^3
14. g = 9.8 ;m/s^2
15. tu = 0 ;m/s
16. t = 0 ;s
17.
18. while tu lt uv do begin
19.
20. ;calculating Reynold's number
22. print, 're calculated'
23.
24. ;calculate diffusion coefficients
25. if re le 0.1 do begin
26. c = 24/re
27. print, 'first c while loop calculated'
28. endif
29.
30. if (re gt 0.1) and (re le 2) do begin
31. c = (24/re)*(1+(3/16)*re+(9/160)*re^2*ln(2*re))
32. print, 'second c while loop calculated'
33. endif
34.
35. if (re gt 2) and (re le 500) do begin
36. c = (24/re)(1+0.15*re^0.687)
37. print, 'third c while loop calculated'
38. endif
39.
40. if (re gt 500) and (re le 2E5) do begin
41. c = 0.44
42. print, 'fourth c while loop calculated'
43. endif
44.
45. ;calculating terminal velocity
47. print, 'terminal velocity calculated'
48.