! $Id: ropp_io_ecmwf.f90 4452 2015-01-29 14:42:02Z idculv $

MODULE ropp_io_ecmwf

!****m* ECMWF/ropp_io_ecmwf *
!
! NAME
!   ropp_io_ecmwf
!
! SYNOPSIS
!   USE ropp_io_ecmwf
!
! DESCRIPTION
!   This module contains the Ak and Bk values for the 60L, 91L and 137L 
!   ECMWF models.
!
! AUTHOR
!   Met Office, Exeter, UK, using data from ECMWF, Reading, UK.
!   Any comments on this software should be given via the ROM SAF
!   Helpdesk at http://www.romsaf.org
!
! COPYRIGHT
!   (c) EUMETSAT. All rights reserved.
!   For further details please refer to the file COPYRIGHT
!   which you should have received as part of this distribution.
!
!****

  USE typesizes,     wp => EightByteReal


! 60 level data
! -----------------

! Standard 60 level Ak (Pa), from ECMWF's namelist NAMVV1 in vtable_L60, 
! approximately the same as the numbers in from
! http://www.ecmwf.int/products/data/technical/model_levels/model_def_60.html
  REAL(wp), PARAMETER, DIMENSION(61)  :: Ak60=(/ &
    0.000000_wp, 0.000000_wp, 7.367743_wp, 65.889244_wp, 210.393890_wp, 467.333588_wp, &
    855.361755_wp, 1385.912598_wp, 2063.779785_wp, 2887.696533_wp, 3850.913330_wp, &
    4941.778320_wp, 6144.314941_wp, 7438.803223_wp, 8802.356445_wp, 10209.500977_wp, &
    11632.758789_wp, 13043.218750_wp, 14411.124023_wp, 15706.447266_wp, 16899.468750_wp, &
    17961.357422_wp, 18864.750000_wp, 19584.330078_wp, 20097.402344_wp, 20384.480469_wp, &
    20429.863281_wp, 20222.205078_wp, 19755.109375_wp, 19027.695313_wp, 18045.183594_wp, &
    16819.474609_wp, 15379.805664_wp, 13775.325195_wp, 12077.446289_wp, 10376.126953_wp, &
    8765.053711_wp, 7306.631348_wp, 6018.019531_wp, 4906.708496_wp, 3960.291504_wp, &
    3196.421631_wp, 2579.888672_wp, 2082.273926_wp, 1680.640259_wp, 1356.474609_wp, &
    1094.834717_wp, 883.660522_wp, 713.218079_wp, 575.651001_wp, 464.618134_wp, &
    373.971924_wp, 298.495789_wp, 234.779053_wp, 180.584351_wp, 134.483307_wp, &
    95.636963_wp, 63.647804_wp, 38.425343_wp, 20.000000_wp, 0.000000_wp /)

! Standard 60 level Bk (no unit), from ECMWF's namelist NAMVV1 in vtable_L60, 
! approximately the same as the numbers in from
! http://www.ecmwf.int/products/data/technical/model_levels/model_def_60.html
  REAL(wp), PARAMETER, DIMENSION(61)  :: Bk60=(/ &
    1.0000000000_wp, 0.9976301193_wp, 0.9940194488_wp, 0.9882701039_wp, 0.9796627164_wp, &
    0.9676452279_wp, 0.9518215060_wp, 0.9319403172_wp, 0.9078838825_wp, 0.8796569109_wp, &
    0.8473749161_wp, 0.8112534285_wp, 0.7715966105_wp, 0.7287858129_wp, 0.6832686067_wp, &
    0.6355474591_wp, 0.5861684084_wp, 0.5357099175_wp, 0.4847715795_wp, 0.4339629412_wp, &
    0.3838921487_wp, 0.3351548910_wp, 0.2883229554_wp, 0.2439331412_wp, 0.2024759352_wp, &
    0.1643843204_wp, 0.1300225109_wp, 0.0996746942_wp, 0.0735338330_wp, 0.0516904071_wp, &
    0.0341211632_wp, 0.0206778757_wp, 0.0111429105_wp, 0.0050811190_wp, 0.0018151561_wp, &
    0.0004613950_wp, 0.0000758235_wp, 0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, &
    0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, &
    0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, &
    0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, &
    0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, &
    0.0000000000_wp /)


! 91 level data
! -------------

! Standard 91 level Ak (Pa), from ECMWF's namelist NAMVV1 in vtable_L91, 
! approximately the same as the numbers in 
! http://www.ecmwf.int/products/data/technical/model_levels/model_def_91.html
  REAL(wp), PARAMETER, DIMENSION(92)  :: Ak91=(/ &
    0.000000_wp, 0.003160_wp, 6.575628_wp, 54.208336_wp, 162.043427_wp, 336.772369_wp, &
    576.314148_wp, 895.193542_wp, 1297.656128_wp, 1784.854614_wp, 2356.202637_wp, &
    3010.146973_wp, 3743.464355_wp, 4550.215820_wp, 5422.802734_wp, 6353.920898_wp, &
    7335.164551_wp, 8356.252930_wp, 9405.222656_wp, 10471.310547_wp, 11543.166992_wp, &
    12608.383789_wp, 13653.219727_wp, 14665.645508_wp, 15633.566406_wp, 16544.585938_wp, &
    17385.595703_wp, 18141.296875_wp, 18798.822266_wp, 19348.775391_wp, 19785.357422_wp, &
    20107.031250_wp, 20319.011719_wp, 20426.218750_wp, 20434.158203_wp, 20348.916016_wp, &
    20175.394531_wp, 19919.796875_wp, 19587.513672_wp, 19184.544922_wp, 18716.968750_wp, &
    18191.029297_wp, 17613.281250_wp, 16990.623047_wp, 16329.560547_wp, 15638.053711_wp, &
    14922.685547_wp, 14192.009766_wp, 13453.225586_wp, 12713.897461_wp, 11982.662109_wp, &
    11262.484375_wp, 10558.881836_wp, 9873.560547_wp, 9208.305664_wp, 8564.624023_wp, &
    7942.926270_wp, 7341.469727_wp, 6759.727051_wp, 6199.839355_wp, 5663.156250_wp, &
    5150.859863_wp, 4663.776367_wp, 4202.416504_wp, 3767.196045_wp, 3358.425781_wp, &
    2976.302246_wp, 2620.898438_wp, 2292.155518_wp, 1989.874390_wp, 1713.709595_wp, &
    1463.163940_wp, 1237.585449_wp, 1036.166504_wp, 857.945801_wp, 701.813354_wp, &
    566.519226_wp, 450.685791_wp, 352.824493_wp, 271.356506_wp, 204.637451_wp, &
    150.986023_wp, 108.715561_wp, 76.167656_wp, 51.746601_wp, 33.952858_wp, 21.413612_wp, &
    12.908319_wp, 7.387186_wp, 3.980832_wp, 2.000040_wp, 0.000000_wp /)

! Standard 91 level Bk (no unit), from ECMWF's namelist NAMVV1 in vtable_L91, 
! approximately the same as the numbers in 
! http://www.ecmwf.int/products/data/technical/model_levels/model_def_91.html
  REAL(wp), PARAMETER, DIMENSION(92)  :: Bk91=(/ &
    1.0000000000_wp, 0.9976301193_wp, 0.9942041636_wp, 0.9891529679_wp, 0.9822381139_wp, &
    0.9734660387_wp, 0.9630070925_wp, 0.9502738118_wp, 0.9351570606_wp, 0.9176509380_wp, &
    0.8977672458_wp, 0.8755183816_wp, 0.8509504199_wp, 0.8241854310_wp, 0.7953847647_wp, &
    0.7646794915_wp, 0.7322238088_wp, 0.6982235909_wp, 0.6629335880_wp, 0.6265588999_wp, &
    0.5893171430_wp, 0.5514584780_wp, 0.5132797360_wp, 0.4750164151_wp, 0.4369063377_wp, &
    0.3992047608_wp, 0.3622025847_wp, 0.3263293803_wp, 0.2919934094_wp, 0.2595544755_wp, &
    0.2293148786_wp, 0.2015201300_wp, 0.1760910898_wp, 0.1529335231_wp, 0.1319348216_wp, &
    0.1129790097_wp, 0.0959640592_wp, 0.0807773024_wp, 0.0673161894_wp, 0.0554742292_wp, &
    0.0451461338_wp, 0.0362269096_wp, 0.0286103487_wp, 0.0221886449_wp, 0.0168595780_wp, &
    0.0125082629_wp, 0.0090349903_wp, 0.0063221669_wp, 0.0042670486_wp, 0.0027647193_wp, &
    0.0017010754_wp, 0.0010001344_wp, 0.0005483841_wp, 0.0002788848_wp, 0.0001313641_wp, &
    0.0000546672_wp, 0.0000139116_wp, 0.0000002724_wp, 0.0000000000_wp, 0.0000000000_wp, &
    0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, &
    0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, &
    0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, &
    0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, &
    0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, &
    0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, &
    0.0000000000_wp, 0.0000000000_wp /)


! 137 level data
! ------------------

! Standard 137 level Ak (Pa), from ECMWF's namelist NAMVV1 in vtable_L137, 
! approximately the same as the numbers in 
! http://www.ecmwf.int/products/data/technical/model_levels/model_def_137.html
  REAL(wp), PARAMETER, DIMENSION(138)  :: Ak137=(/ &
    0.000000_wp, 0.000000_wp, 3.757813_wp, 22.835938_wp, 62.781250_wp, 122.101563_wp, &
    202.484375_wp, 302.476563_wp, 424.414063_wp, 568.062500_wp, 734.992188_wp, &
    926.507813_wp, 1143.250000_wp, 1387.546875_wp, 1659.476563_wp, 1961.500000_wp, &
    2294.242188_wp, 2659.140625_wp, 3057.265625_wp, 3489.234375_wp, 3955.960938_wp, &
    4457.375000_wp, 4993.796875_wp, 5564.382813_wp, 6168.531250_wp, 6804.421875_wp, &
    7470.343750_wp, 8163.375000_wp, 8880.453125_wp, 9617.515625_wp, 10370.175781_wp, &
    11133.304688_wp, 11901.339844_wp, 12668.257813_wp, 13427.769531_wp, 14173.324219_wp, &
    14898.453125_wp, 15596.695313_wp, 16262.046875_wp, 16888.687500_wp, 17471.839844_wp, &
    18006.925781_wp, 18489.707031_wp, 18917.460938_wp, 19290.226563_wp, 19608.572266_wp, &
    19874.025391_wp, 20087.085938_wp, 20249.511719_wp, 20361.816406_wp, 20425.718750_wp, &
    20442.078125_wp, 20412.308594_wp, 20337.863281_wp, 20219.664063_wp, 20059.931641_wp, &
    19859.390625_wp, 19620.042969_wp, 19343.511719_wp, 19031.289063_wp, 18685.718750_wp, &
    18308.433594_wp, 17901.621094_wp, 17467.613281_wp, 17008.789063_wp, 16527.322266_wp, &
    16026.115234_wp, 15508.256836_wp, 14975.615234_wp, 14432.139648_wp, 13881.331055_wp, &
    13324.668945_wp, 12766.873047_wp, 12211.547852_wp, 11660.067383_wp, 11116.662109_wp, &
    10584.631836_wp, 10065.978516_wp, 9562.682617_wp, 9076.400391_wp, 8608.525391_wp, &
    8159.354004_wp, 7727.412109_wp, 7311.869141_wp, 6911.870605_wp, 6526.946777_wp, &
    6156.074219_wp, 5798.344727_wp, 5452.990723_wp, 5119.895020_wp, 4799.149414_wp, &
    4490.817383_wp, 4194.930664_wp, 3911.490479_wp, 3640.468262_wp, 3381.743652_wp, &
    3135.119385_wp, 2900.391357_wp, 2677.348145_wp, 2465.770508_wp, 2265.431641_wp, &
    2076.095947_wp, 1897.519287_wp, 1729.448975_wp, 1571.622925_wp, 1423.770142_wp, &
    1285.610352_wp, 1156.853638_wp, 1037.201172_wp, 926.344910_wp, 823.967834_wp, &
    729.744141_wp, 643.339905_wp, 564.413452_wp, 492.616028_wp, 427.592499_wp, &
    368.982361_wp, 316.420746_wp, 269.539581_wp, 227.968948_wp, 191.338562_wp, &
    159.279404_wp, 131.425507_wp, 107.415741_wp, 86.895882_wp, 69.520576_wp, &
    54.955463_wp, 42.879242_wp, 32.985710_wp, 24.985718_wp, 18.608931_wp, 13.605424_wp, &
    9.746966_wp, 6.827977_wp, 4.666084_wp, 3.102241_wp, 2.000365_wp, 0.000000_wp /)

! Standard 137 level Bk (no unit), from ECMWF's namelist NAMVV1 in vtable_L137, 
! approximately the same as the numbers in 
! http://www.ecmwf.int/products/data/technical/model_levels/model_def_137.html
  REAL(wp), PARAMETER, DIMENSION(138)  :: Bk137=(/ &
    1.0000000000_wp, 0.9976301193_wp, 0.9950025082_wp, 0.9919840097_wp, 0.9884995222_wp, &
    0.9845418930_wp, 0.9800716043_wp, 0.9750784039_wp, 0.9695134163_wp, 0.9633517265_wp, &
    0.9565495253_wp, 0.9490644336_wp, 0.9408595562_wp, 0.9318807721_wp, 0.9220956564_wp, &
    0.9114481807_wp, 0.8999004960_wp, 0.8874075413_wp, 0.8739292622_wp, 0.8594318032_wp, &
    0.8438811302_wp, 0.8272560835_wp, 0.8095360398_wp, 0.7907167673_wp, 0.7707975507_wp, &
    0.7497965693_wp, 0.7277387381_wp, 0.7046689987_wp, 0.6806430221_wp, 0.6557359695_wp, &
    0.6300358176_wp, 0.6036480665_wp, 0.5766921639_wp, 0.5493011475_wp, 0.5216192007_wp, &
    0.4938003123_wp, 0.4660032988_wp, 0.4383912086_wp, 0.4111247659_wp, 0.3843633235_wp, &
    0.3582541943_wp, 0.3329390883_wp, 0.3085984588_wp, 0.2853540182_wp, 0.2632418871_wp, &
    0.2422440052_wp, 0.2223328650_wp, 0.2034912109_wp, 0.1856894493_wp, 0.1689104140_wp, &
    0.1531250328_wp, 0.1383128911_wp, 0.1244481280_wp, 0.1115046516_wp, 0.0994616672_wp, &
    0.0882857367_wp, 0.0779583082_wp, 0.0684482530_wp, 0.0597284138_wp, 0.0517730154_wp, &
    0.0445479602_wp, 0.0380263999_wp, 0.0321760960_wp, 0.0269635208_wp, 0.0223548450_wp, &
    0.0183184519_wp, 0.0148156285_wp, 0.0118060224_wp, 0.0092614600_wp, 0.0071333768_wp, &
    0.0053778146_wp, 0.0039709536_wp, 0.0028571242_wp, 0.0019918380_wp, 0.0013528055_wp, &
    0.0008896979_wp, 0.0005615553_wp, 0.0003403797_wp, 0.0001985774_wp, 0.0001119143_wp, &
    0.0000589220_wp, 0.0000243480_wp, 0.0000067607_wp, 0.0000000382_wp, 0.0000000000_wp, &
    0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, &
    0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, &
    0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, &
    0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, &
    0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, &
    0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, &
    0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, &
    0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, &
    0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, &
    0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp, &
    0.0000000000_wp, 0.0000000000_wp, 0.0000000000_wp /)


END MODULE ropp_io_ecmwf