const_ec_v.F90 2.9 KB

12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152535455565758596061626364656667686970717273
  1. #include "tm5.inc"
  2. ! ----------------------------------------------------------------------
  3. !
  4. ! NAME
  5. ! const_ec_v - vertical level definition of the ecwmf model
  6. !
  7. ! USAGE
  8. ! use const_ec_v
  9. !
  10. ! CHANGES
  11. ! may 2001, Arjo Segers, KNMI
  12. ! Rebuilt from 'constants_ecmwf_##_##_##.f90'
  13. !
  14. ! ----------------------------------------------------------------------
  15. module const_ec_v
  16. implicit none
  17. ! --- in/out -----------------------------
  18. private
  19. public :: lme, a_ec, b_ec
  20. ! --- const --------------------------------
  21. ! number of levels:
  22. integer, parameter :: lme = 60
  23. ! Hybrid coordinate coefficients at half level interfaces,
  24. ! specifying 60 vertical ECMWF levels.
  25. ! Coefficient a is given in [Pa] .
  26. ! Pressure at a certain hybrid level is defined as: p = a + b.SP
  27. real, parameter :: a_ec(0:lme) = (/ &
  28. 0.000000, 20.000000, 38.425343, 63.647804, &
  29. 95.636963, 134.483307, 180.584351, 234.779053, &
  30. 298.495789, 373.971924, 464.618134, 575.651001, &
  31. 713.218079, 883.660522, 1094.834717, 1356.474609, &
  32. 1680.640259, 2082.273926, 2579.888672, 3196.421631, &
  33. 3960.291504, 4906.708496, 6018.019531, 7306.631348, &
  34. 8765.053711, 10376.126953, 12077.446289, 13775.325195, &
  35. 15379.805664, 16819.474609, 18045.183594, 19027.695313, &
  36. 19755.109375, 20222.205078, 20429.863281, 20384.480469, &
  37. 20097.402344, 19584.330078, 18864.750000, 17961.357422, &
  38. 16899.468750, 15706.447266, 14411.124023, 13043.218750, &
  39. 11632.758789, 10209.500977, 8802.356445, 7438.803223, &
  40. 6144.314941, 4941.778320, 3850.913330, 2887.696533, &
  41. 2063.779785, 1385.912598, 855.361755, 467.333588, &
  42. 210.393890, 65.889244, 7.367743, 0.000000, &
  43. 0.000000/)
  44. real, parameter :: b_ec(0:lme) = (/ &
  45. 0.00000000, 0.00000000, 0.00000000, 0.00000000, &
  46. 0.00000000, 0.00000000, 0.00000000, 0.00000000, &
  47. 0.00000000, 0.00000000, 0.00000000, 0.00000000, &
  48. 0.00000000, 0.00000000, 0.00000000, 0.00000000, &
  49. 0.00000000, 0.00000000, 0.00000000, 0.00000000, &
  50. 0.00000000, 0.00000000, 0.00000000, 0.00000000, &
  51. 0.00007582, 0.00046139, 0.00181516, 0.00508112, &
  52. 0.01114291, 0.02067788, 0.03412116, 0.05169041, &
  53. 0.07353383, 0.09967469, 0.13002251, 0.16438432, &
  54. 0.20247590, 0.24393314, 0.28832296, 0.33515489, &
  55. 0.38389215, 0.43396294, 0.48477158, 0.53570992, &
  56. 0.58616841, 0.63554746, 0.68326861, 0.72878581, &
  57. 0.77159661, 0.81125343, 0.84737492, 0.87965691, &
  58. 0.90788388, 0.93194032, 0.95182151, 0.96764523, &
  59. 0.97966272, 0.98827010, 0.99401945, 0.99763012, &
  60. 1.00000000 /)
  61. end module const_ec_v