q.txt 2.0 KB

12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152535455565758596061626364656667686970717273747576777879808182838485868788899091
  1. Water cycle in PlaSim gq gqdt dq dqdt
  2. =========================================
  3. GRIDPOINTA
  4. INVLEGA
  5. gq <- sq (sp2fc)
  6. -----
  7. gq <- gq (fc2gp) gq = ps * q gq
  8. CALCGP
  9. zq <- gq / gp
  10. gqm <- (gq(lev+1) - gq(lev)) * zsdotp
  11. gqn <- (gqm(lev) - gqm(lev-1)) * rdsig
  12. -----
  13. guq <- gu * gq
  14. gvq <- gv * gq
  15. gqn <- gqn (gp2fc) gqn
  16. guq <- guq (gp2fc)
  17. gvq <- gvq (gp2fc)
  18. QTEND
  19. sqf <- gqn & guq & gvq (Legendre)
  20. -----
  21. sqt <- sqf (mpsumsc)
  22. dq <- gq * (psurf/dp) dq
  23. =========================================
  24. SPECTRALA
  25. aqm <- sqm
  26. sqm <- sqm & aqm ??? sqm == aqm ???
  27. sqp <- aqm + delt2 * sqt
  28. sq <- sqp (mpgallsp)
  29. =========================================
  30. GRIDPOINTD
  31. INVLEGD
  32. gq <- sq (sp2fc)
  33. -----
  34. gq <- gq (fc2gp) gq = gqn ?
  35. dq <- gq * (psurf/dp)
  36. TRACER_MAIN
  37. gq <- gq (SL advection)
  38. -----
  39. gqdt <- gq (SL term)
  40. gqdt <- dqdt * (dp/ww/psurf)
  41. MISCSTEP
  42. FIXER
  43. gqdt
  44. dq <- fix negative values
  45. .....
  46. -----
  47. FLUXSTEP
  48. MKEVAP
  49. zqn <- dq
  50. dqdt <- zqn & dqdt & dq
  51. .....
  52. VDIFF
  53. dqdt <- zqn
  54. .....
  55. -----
  56. RADSTEP
  57. RAINSTEP
  58. dqt <- sqt (mkdqtgp)
  59. KUO
  60. MKDCA
  61. MKLSP
  62. MKRAIN
  63. -----
  64. SURFSTEP
  65. -----
  66. gqdt <- gqdt & dqdt ??? gqdt vorher 0 ???
  67. gqdt <- gqdt (gp2fc)
  68. sqf <- gqdt (fc2sp)
  69. sqt <- sqf (mpsumsc)
  70. =========================================
  71. OUTSP
  72. write sqout
  73. =========================================
  74. SPECTRALD
  75. zsq <- sqp (mpgallsp)
  76. sqp <- sqp + delt2 * sqt
  77. MKDQTGP
  78. dqt <- sqt
  79. -----
  80. sqt <- sqp & tdiss
  81. zqgp <- zsq (sp2fc)
  82. zqgp <- zqgp (fc2gp)
  83. zqgp <- zqgp * (psurf/dp)
  84. sqp <- sqp + delt2 * sqt
  85. sqm <- sqm + pnu * sqp
  86. sq <- sqp (mügallsp)
  87. MKDQTGP
  88. dqt <- sqt
  89. -----
  90. OUTGP