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
  13. ad = 1.29 ;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
  21. re = (2*ad*r)/(mu*100000)
  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
  46. tu = sqrt((8*r*wd*g)/(3*ad*c))
  47. print, 'terminal velocity calculated'
  48.  
  49. ;calculating new radius
  50. r = ((E*uv*LWC*t)/(4*wd))+sqrt((2*(S-1)*t)/(Fk+Fd))
  51. print, 'radius calculated'
  52. t = t+1
  53. print, 'time incremented'
  54. endwhile
  55.  
  56. stop
  57. end
RAW Paste Data
We use cookies for various purposes including analytics. By continuing to use Pastebin, you agree to our use of cookies as described in the Cookies Policy. OK, I Understand
 
Top