#include "tm5.inc"
! ----------------------------------------------------------------------
!
! NAME
!   const_ec_v - vertical level definition of the driving model
!
! USAGE
!   use const_ec_v
!
! REFERENCES
!    Home > Products > Data Services > Technical information > Model Levels
!    http://www.ecmwf.int/products/data/technical/model_levels/index.html
!
! ----------------------------------------------------------------------

module const_ec_v

  implicit none
  
  ! --- in/out -----------------------------
  
  private 

  public :: lme, a_ec, b_ec
  
  
  ! --- const --------------------------------
  
  ! number of levels:
  integer, parameter :: lme = 31
      
  ! Hybrid coordinate coefficients at half level interfaces,
  ! specifying 31 vertical ECMWF levels.
  ! Coefficient a is given in [Pa] .
  ! Pressure at a certain hybrid level is defined as: p = a + b.SP

  real, parameter  :: a_ec(0:lme) = (/   &
           0.000000,   2000.000000,   4000.000000,   6000.000000,   8000.000000, &
        9976.135361,  11820.539617,  13431.393926,  14736.356909,  15689.207458, &
       16266.610500,  16465.005734,  16297.619332,  15791.598604,  14985.269630, &
       13925.517858,  12665.291662,  11261.228878,   9771.406290,   8253.212096, &
        6761.341326,   5345.914240,   4050.717678,   2911.569385,   1954.805296, &
        1195.889791,    638.148911,    271.626545,     72.063577,      0.000000, &
           0.000000,      0.000000 /)

  real, parameter  :: b_ec(0:lme) = (/   &
       0.0000000000,      0.0000000000,  0.0000000000,   0.0000000000, 0.0000000000, &
       0.0003908582,      0.0029197006,  0.0091941320,   0.0203191555, 0.0369748598, &
       0.0594876397,      0.0878949492,  0.1220035886,   0.1614415235, 0.2057032385, &
       0.2541886223,      0.3062353873,  0.3611450218,   0.4182022749, 0.4766881754, &
       0.5358865832,      0.5950842740,  0.6535645569,   0.7105944258, 0.7654052430, &
       0.8171669567,      0.8649558510,  0.9077158297,   0.9442132326, 0.9729851852, &
       0.9922814815,      1.0000000000 /)
          
   
end module const_ec_v