GlobalPressureTemperatureModel.java

  1. /* Copyright 2002-2019 CS Systèmes d'Information
  2.  * Licensed to CS Systèmes d'Information (CS) under one or more
  3.  * contributor license agreements.  See the NOTICE file distributed with
  4.  * this work for additional information regarding copyright ownership.
  5.  * CS licenses this file to You under the Apache License, Version 2.0
  6.  * (the "License"); you may not use this file except in compliance with
  7.  * the License.  You may obtain a copy of the License at
  8.  *
  9.  *   http://www.apache.org/licenses/LICENSE-2.0
  10.  *
  11.  * Unless required by applicable law or agreed to in writing, software
  12.  * distributed under the License is distributed on an "AS IS" BASIS,
  13.  * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
  14.  * See the License for the specific language governing permissions and
  15.  * limitations under the License.
  16.  */
  17. package org.orekit.models.earth;

  19. import org.hipparchus.util.CombinatoricsUtils;
  20. import org.hipparchus.util.FastMath;
  21. import org.orekit.forces.gravity.potential.GravityFieldFactory;
  22. import org.orekit.frames.Frame;
  23. import org.orekit.time.AbsoluteDate;
  24. import org.orekit.time.DateTimeComponents;
  25. import org.orekit.time.TimeScalesFactory;

  27. /** The Global Pressure and Temperature model.
  28.  * This model is an empirical model that provides the temperature and the pressure depending
  29.  * the latitude and the longitude of the station.
  30.  * <p>
  31.  * The Global Pressure and Temperature model is based on spherical harmonics up
  32.  * to degree and order of 9. The residual values ​​of this model can reach 20 hPa
  33.  * for pressure and 10 ° C for temperature. They are significant for higher latitudes and
  34.  * small near the equator (Böhm, 2007)
  35.  * </p>
  36.  *
  37.  * @see J. Böhm, R. Heinkelmann, and H. Schuh (2007),
  38.  * Short Note: A global model of pressure and temperature for geodetic applications. J Geod,
  39.  * doi:10.1007/s00190-007-0135-3.
  40.  *
  41.  * @author Bryan Cazabonne
  42.  *
  43.  */
  44. public class GlobalPressureTemperatureModel implements WeatherModel {

  46.     /** Temperature gradient (°C/m). */
  47.     private static final double TEMPERATURE_GRADIENT = -6.5e-3;

  49.     /** Geodetic latitude, in radians. */
  50.     private final double latitude;

  52.     /** Geodetic longitude, in radians. */
  53.     private final double longitude;

  55.     /** Temperature site, in kelvins. */
  56.     private double temperature;

  58.     /** Pressure site, in hPa. */
  59.     private double pressure;

  61.     /** Body frame related to body shape. */
  62.     private final Frame bodyFrame;

  64.     /** Build a new instance.
  65.      * <p>
  66.      * At the initialization the values of the pressure and the temperature are set to NaN.
  67.      * The user has to call {@link #weatherParameters(double, AbsoluteDate)} method before using
  68.      * the values of the pressure and the temperature.
  69.      * </p>
  70.      * @param latitude geodetic latitude, in radians
  71.      * @param longitude geodetic longitude, in radians
  72.      * @param bodyFrame the frame to attach to the ellipsoid. The origin is at
  73.      *                  the center of mass, the z axis is the minor axis.
  74.      */
  75.     public GlobalPressureTemperatureModel(final double latitude, final double longitude, final Frame bodyFrame) {
  76.         this.bodyFrame   = bodyFrame;
  77.         this.latitude    = latitude;
  78.         this.longitude   = longitude;
  79.         this.temperature = Double.NaN;
  80.         this.pressure    = Double.NaN;
  81.     }

  83.     /** Get the atmospheric temperature of the station depending its position.
  84.      * @return the temperature in kelvins
  85.      */
  86.     public double getTemperature() {
  87.         return temperature;
  88.     }

  90.     /** Get the atmospheric pressure of the station depending its position.
  91.      * @return the pressure in hPa
  92.      */
  93.     public double getPressure() {
  94.         return pressure;
  95.     }

  97.     @Override
  98.     public void weatherParameters(final double height, final AbsoluteDate date) {

  100.         // Day of year computation
  101.         final DateTimeComponents dtc = date.getComponents(TimeScalesFactory.getUTC());
  102.         final int dofyear = dtc.getDate().getDayOfYear();

  104.         // Reference day: 28 January 1980 (Niell, 1996)
  105.         final int t0 = 28;
  106.         final double coef = ((dofyear + 1 - t0) / 365.25) * 2 * FastMath.PI;
  107.         final double cosCoef = FastMath.cos(coef);

  109.         // Compute Legendre Polynomials Pnm(sin(phi))
  110.         final int degree = 9;
  111.         final int order  = 9;
  112.         final LegendrePolynomials p = new LegendrePolynomials(degree, order);

  114.         // Geoid for height computation
  115.         final Geoid geoid = new Geoid(GravityFieldFactory.getNormalizedProvider(degree, order),
  116.                                       ReferenceEllipsoid.getWgs84(bodyFrame));

  118.         // Corrected height
  119.         final double correctedheight = FastMath.max(0.0, height - geoid.getUndulation(latitude, longitude, date));

  121.         // Eq. 4 (Ref)
  122.         double meanT0      = 0.0;
  123.         double amplitudeT0 = 0.0;
  124.         double meanP0      = 0.0;
  125.         double amplitudeP0 = 0.0;
  126.         final ABCoefficients abCoef = new ABCoefficients();
  127.         int j = 0;
  128.         for (int n = 0; n <= 9; n++) {
  129.             for (int m = 0; m <= n; m++) {
  130.                 final double pCosmLambda = p.getPnm(n, m) * FastMath.cos(m * longitude);
  131.                 final double pSinmLambda = p.getPnm(n, m) * FastMath.sin(m * longitude);

  133.                 meanT0      = meanT0 +
  134.                                 (abCoef.getAnmTemperatureMean(j) * pCosmLambda + abCoef.getBnmTemperatureMean(j) * pSinmLambda);
  135.                 amplitudeT0 = amplitudeT0 +
  136.                                 (abCoef.getAnmTemperatureAmpl(j) * pCosmLambda + abCoef.getBnmTemperatureAmpl(j) * pSinmLambda);
  137.                 meanP0      = meanP0 +
  138.                                 (abCoef.getAnmPressureMean(j) * pCosmLambda + abCoef.getBnmPressureMean(j) * pSinmLambda);
  139.                 amplitudeP0 = amplitudeP0 +
  140.                                 (abCoef.getAnmPressureAmpl(j) * pCosmLambda + abCoef.getBnmPressureAmpl(j) * pSinmLambda);

  142.                 j = j + 1;
  143.             }
  144.         }

  146.         // Eq. 3 (Ref)
  147.         final double temp0 = meanT0 + amplitudeT0 * cosCoef;
  148.         final double pres0 = meanP0 + amplitudeP0 * cosCoef;

  150.         // Compute pressure and temperature Eq. 1 and 2 (Ref)
  151.         final double degrees = temp0 + TEMPERATURE_GRADIENT * correctedheight;
  152.         this.temperature = degrees + 273.15;
  153.         this.pressure    = pres0 * FastMath.pow(1.0 - correctedheight * 0.0000226, 5.225);
  154.     }

  156.     /** Computes the P<sub>nm</sub>(sin(&#934)) coefficients of Eq. 4 (Ref).
  157.      *  The computation of the Legendre polynomials is performed following:
  158.      *  Heiskanen and Moritz, Physical Geodesy, 1967, eq. 1-62
  159.      */
  160.     private class LegendrePolynomials {

  162.         /** Array for the Legendre polynomials. */
  163.         private double[][] pCoef;

  165.         /** Create Legendre polynomials for the given degree and order.
  166.          * @param degree degree of the spherical harmonics
  167.          * @param order order of the spherical harmonics
  168.          */
  169.         LegendrePolynomials(final int degree, final int order) {

  171.             this.pCoef = new double[degree + 1][order + 1];

  173.             final double t  = FastMath.sin(latitude);
  174.             final double t2 = t * t;

  176.             for (int n = 0; n <= degree; n++) {

  178.                 // m shall be <= n (Heiskanen and Moritz, 1967, pp 21)
  179.                 for (int m = 0; m <= FastMath.min(n, order); m++) {

  181.                     // r = int((n - m) / 2)
  182.                     final int r = (int) (n - m) / 2;
  183.                     double sum = 0.;
  184.                     for (int k = 0; k <= r; k++) {
  185.                         final double term = FastMath.pow(-1.0, k) * CombinatoricsUtils.factorialDouble(2 * n - 2 * k) /
  186.                                         (CombinatoricsUtils.factorialDouble(k) * CombinatoricsUtils.factorialDouble(n - k) *
  187.                                          CombinatoricsUtils.factorialDouble(n - m - 2 * k)) *
  188.                                          FastMath.pow(t, n - m - 2 * k);
  189.                         sum = sum + term;
  190.                     }

  192.                     pCoef[n][m] = FastMath.pow(2, -n) * FastMath.pow(1.0 - t2, 0.5 * m) * sum;

  194.                 }

  196.             }

  198.         }

  200.         /** Return the coefficient P<sub>nm</sub>.
  201.          * @param n index
  202.          * @param m index
  203.          * @return The coefficient P<sub>nm</sub>
  204.          */
  205.         public double getPnm(final int n, final int m) {
  206.             return pCoef[n][m];
  207.         }

  209.     }

  211.     private static class ABCoefficients {

  213.         /** Mean A<sub>nm</sub> coefficients for the pressure. */
  214.         private static final double[] A_PRESSURE_MEAN = {
  215.             1.0108e+03,
  216.             8.4886e+00,
  217.             1.4799e+00,
  218.             -1.3897e+01,
  219.             3.7516e-03,
  220.             -1.4936e-01,
  221.             1.2232e+01,
  222.             -7.6615e-01,
  223.             -6.7699e-02,
  224.             8.1002e-03,
  225.             -1.5874e+01,
  226.             3.6614e-01,
  227.             -6.7807e-02,
  228.             -3.6309e-03,
  229.             5.9966e-04,
  230.             4.8163e+00,
  231.             -3.7363e-01,
  232.             -7.2071e-02,
  233.             1.9998e-03,
  234.             -6.2385e-04,
  235.             -3.7916e-04,
  236.             4.7609e+00,
  237.             -3.9534e-01,
  238.             8.6667e-03,
  239.             1.1569e-02,
  240.             1.1441e-03,
  241.             -1.4193e-04,
  242.             -8.5723e-05,
  243.             6.5008e-01,
  244.             -5.0889e-01,
  245.             -1.5754e-02,
  246.             -2.8305e-03,
  247.             5.7458e-04,
  248.             3.2577e-05,
  249.             -9.6052e-06,
  250.             -2.7974e-06,
  251.             1.3530e+00,
  252.             -2.7271e-01,
  253.             -3.0276e-04,
  254.             3.6286e-03,
  255.             -2.0398e-04,
  256.             1.5846e-05,
  257.             -7.7787e-06,
  258.             1.1210e-06,
  259.             9.9020e-08,
  260.             5.5046e-01,
  261.             -2.7312e-01,
  262.             3.2532e-03,
  263.             -2.4277e-03,
  264.             1.1596e-04,
  265.             2.6421e-07,
  266.             -1.3263e-06,
  267.             2.7322e-07,
  268.             1.4058e-07,
  269.             4.9414e-09
  270.         };

  272.         /** Mean B<sub>nm</sub> coefficients for the pressure. */
  273.         private static final double[] B_PRESSURE_MEAN = {
  274.             0.0000e+00,
  275.             0.0000e+00,
  276.             -1.2878e+00,
  277.             0.0000e+00,
  278.             7.0444e-01,
  279.             3.3222e-01,
  280.             0.0000e+00,
  281.             -2.9636e-01,
  282.             7.2248e-03,
  283.             7.9655e-03,
  284.             0.0000e+00,
  285.             1.0854e+00,
  286.             1.1145e-02,
  287.             -3.6513e-02,
  288.             3.1527e-03,
  289.             0.0000e+00,
  290.             -4.8434e-01,
  291.             5.2023e-02,
  292.             -1.3091e-02,
  293.             1.8515e-03,
  294.             1.5422e-04,
  295.             0.0000e+00,
  296.             6.8298e-01,
  297.             2.5261e-03,
  298.             -9.9703e-04,
  299.             -1.0829e-03,
  300.             +1.7688e-04,
  301.             -3.1418e-05,
  302.             +0.0000e+00,
  303.             -3.7018e-01,
  304.             4.3234e-02,
  305.             7.2559e-03,
  306.             3.1516e-04,
  307.             2.0024e-05,
  308.             -8.0581e-06,
  309.             -2.3653e-06,
  310.             0.0000e+00,
  311.             1.0298e-01,
  312.             -1.5086e-02,
  313.             5.6186e-03,
  314.             3.2613e-05,
  315.             4.0567e-05,
  316.             -1.3925e-06,
  317.             -3.6219e-07,
  318.             -2.0176e-08,
  319.             0.0000e+00,
  320.             -1.8364e-01,
  321.             1.8508e-02,
  322.             7.5016e-04,
  323.             -9.6139e-05,
  324.             -3.1995e-06,
  325.             1.3868e-07,
  326.             -1.9486e-07,
  327.             3.0165e-10,
  328.             -6.4376e-10
  329.         };

  331.         /** Amplitude A<sub>nm</sub> coefficients for the pressure. */
  332.         private static final double[] A_PRESSURE_AMPLITUDE = {
  333.             -1.0444e-01,
  334.             1.6618e-01,
  335.             -6.3974e-02,
  336.             1.0922e+00,
  337.             5.7472e-01,
  338.             -3.0277e-01,
  339.             -3.5087e+00,
  340.             7.1264e-03,
  341.             -1.4030e-01,
  342.             3.7050e-02,
  343.             4.0208e-01,
  344.             -3.0431e-01,
  345.             -1.3292e-01,
  346.             4.6746e-03,
  347.             -1.5902e-04,
  348.             2.8624e+00,
  349.             -3.9315e-01,
  350.             -6.4371e-02,
  351.             1.6444e-02,
  352.             -2.3403e-03,
  353.             4.2127e-05,
  354.             1.9945e+00,
  355.             -6.0907e-01,
  356.             -3.5386e-02,
  357.             -1.0910e-03,
  358.             -1.2799e-04,
  359.             4.0970e-05,
  360.             2.2131e-05,
  361.             -5.3292e-01,
  362.             -2.9765e-01,
  363.             -3.2877e-02,
  364.             1.7691e-03,
  365.             5.9692e-05,
  366.             3.1725e-05,
  367.             2.0741e-05,
  368.             -3.7622e-07,
  369.             2.6372e+00,
  370.             -3.1165e-01,
  371.             1.6439e-02,
  372.             2.1633e-04,
  373.             1.7485e-04,
  374.             2.1587e-05,
  375.             6.1064e-06,
  376.             -1.3755e-08,
  377.             -7.8748e-08,
  378.             -5.9152e-01,
  379.             -1.7676e-01,
  380.             8.1807e-03,
  381.             1.0445e-03,
  382.             2.3432e-04,
  383.             9.3421e-06,
  384.             2.8104e-06,
  385.             -1.5788e-07,
  386.             -3.0648e-08,
  387.             2.6421e-10
  388.         };

  390.         /** Amplitude B<sub>nm</sub> coefficients for the pressure. */
  391.         private static final double[] B_PRESSURE_AMPLITUDE = {
  392.             0.0000e+00,
  393.             0.0000e+00,
  394.             9.3340e-01,
  395.             0.0000e+00,
  396.             8.2346e-01,
  397.             2.2082e-01,
  398.             0.0000e+00,
  399.             9.6177e-01,
  400.             -1.5650e-02,
  401.             1.2708e-03,
  402.             0.0000e+00,
  403.             -3.9913e-01,
  404.             2.8020e-02,
  405.             2.8334e-02,
  406.             8.5980e-04,
  407.             0.0000e+00,
  408.             3.0545e-01,
  409.             -2.1691e-02,
  410.             6.4067e-04,
  411.             -3.6528e-05,
  412.             -1.1166e-04,
  413.             0.0000e+00,
  414.             -7.6974e-02,
  415.             -1.8986e-02,
  416.             +5.6896e-03,
  417.             -2.4159e-04,
  418.             -2.3033e-04,
  419.             -9.6783e-06,
  420.             0.0000e+00,
  421.             -1.0218e-01,
  422.             -1.3916e-02,
  423.             -4.1025e-03,
  424.             -5.1340e-05,
  425.             -7.0114e-05,
  426.             -3.3152e-07,
  427.             1.6901e-06,
  428.             0.0000e+00,
  429.             -1.2422e-02,
  430.             +2.5072e-03,
  431.             +1.1205e-03,
  432.             -1.3034e-04,
  433.             -2.3971e-05,
  434.             -2.6622e-06,
  435.             5.7852e-07,
  436.             4.5847e-08,
  437.             0.0000e+00,
  438.             4.4777e-02,
  439.             -3.0421e-03,
  440.             2.6062e-05,
  441.             -7.2421e-05,
  442.             1.9119e-06,
  443.             3.9236e-07,
  444.             2.2390e-07,
  445.             2.9765e-09,
  446.             -4.6452e-09
  447.         };

  449.         /** Mean A<sub>nm</sub> coefficients for the temperature. */
  450.         private static final double[] A_TEMPERATURE_MEAN = {
  451.             1.6257e+01,
  452.             2.1224e+00,
  453.             9.2569e-01,
  454.             -2.5974e+01,
  455.             1.4510e+00,
  456.             9.2468e-02,
  457.             -5.3192e-01,
  458.             2.1094e-01,
  459.             -6.9210e-02,
  460.             -3.4060e-02,
  461.             -4.6569e+00,
  462.             2.6385e-01,
  463.             -3.6093e-02,
  464.             1.0198e-02,
  465.             -1.8783e-03,
  466.             7.4983e-01,
  467.             1.1741e-01,
  468.             3.9940e-02,
  469.             5.1348e-03,
  470.             5.9111e-03,
  471.             8.6133e-06,
  472.             6.3057e-01,
  473.             1.5203e-01,
  474.             3.9702e-02,
  475.             4.6334e-03,
  476.             2.4406e-04,
  477.             1.5189e-04,
  478.             1.9581e-07,
  479.             5.4414e-01,
  480.             3.5722e-01,
  481.             5.2763e-02,
  482.             4.1147e-03,
  483.             -2.7239e-04,
  484.             -5.9957e-05,
  485.             1.6394e-06,
  486.             -7.3045e-07,
  487.             -2.9394e+00,
  488.             5.5579e-02,
  489.             1.8852e-02,
  490.             3.4272e-03,
  491.             -2.3193e-05,
  492.             -2.9349e-05,
  493.             3.6397e-07,
  494.             2.0490e-06,
  495.             -6.4719e-08,
  496.             -5.2225e-01,
  497.             2.0799e-01,
  498.             1.3477e-03,
  499.             3.1613e-04,
  500.             -2.2285e-04,
  501.             -1.8137e-05,
  502.             -1.5177e-07,
  503.             6.1343e-07,
  504.             7.8566e-08,
  505.             1.0749e-09
  506.         };

  508.         /** Mean B<sub>nm</sub> coefficients for the temperature. */
  509.         private static final double[] B_TEMPERATURE_MEAN = {
  510.             0.0000e+00,
  511.             0.0000e+00,
  512.             1.0210e+00,
  513.             0.0000e+00,
  514.             6.0194e-01,
  515.             1.2292e-01,
  516.             0.0000e+00,
  517.             -4.2184e-01,
  518.             1.8230e-01,
  519.             4.2329e-02,
  520.             0.0000e+00,
  521.             9.3312e-02,
  522.             9.5346e-02,
  523.             -1.9724e-03,
  524.             5.8776e-03,
  525.             0.0000e+00,
  526.             -2.0940e-01,
  527.             3.4199e-02,
  528.             -5.7672e-03,
  529.             -2.1590e-03,
  530.             5.6815e-04,
  531.             0.0000e+00,
  532.             2.2858e-01,
  533.             1.2283e-02,
  534.             -9.3679e-03,
  535.             -1.4233e-03,
  536.             -1.5962e-04,
  537.             4.0160e-05,
  538.             0.0000e+00,
  539.             3.6353e-02,
  540.             -9.4263e-04,
  541.             -3.6762e-03,
  542.             5.8608e-05,
  543.             -2.6391e-05,
  544.             3.2095e-06,
  545.             -1.1605e-06,
  546.             0.0000e+00,
  547.             1.6306e-01,
  548.             1.3293e-02,
  549.             -1.1395e-03,
  550.             5.1097e-05,
  551.             3.3977e-05,
  552.             7.6449e-06,
  553.             -1.7602e-07,
  554.             -7.6558e-08,
  555.             0.0000e+00,
  556.             -4.5415e-02,
  557.             -1.8027e-02,
  558.             3.6561e-04,
  559.             -1.1274e-04,
  560.             1.3047e-05,
  561.             2.0001e-06,
  562.             -1.5152e-07,
  563.             -2.7807e-08,
  564.             7.7491e-09
  565.         };

  567.         /** Amplitude A<sub>nm</sub> coefficients for the temperature. */
  568.         private static final double[] A_TEMPERATURE_AMPLITUDE = {
  569.             -1.8654e+00,
  570.             -9.0041e+00,
  571.             -1.2974e-01,
  572.             -3.6053e+00,
  573.             2.0284e-02,
  574.             2.1872e-01,
  575.             -1.3015e+00,
  576.             4.0355e-01,
  577.             2.2216e-01,
  578.             -4.0605e-03,
  579.             1.9623e+00,
  580.             4.2887e-01,
  581.             2.1437e-01,
  582.             -1.0061e-02,
  583.             -1.1368e-03,
  584.             -6.9235e-02,
  585.             5.6758e-01,
  586.             1.1917e-01,
  587.             -7.0765e-03,
  588.             3.0017e-04,
  589.             3.0601e-04,
  590.             1.6559e+00,
  591.             2.0722e-01,
  592.             6.0013e-02,
  593.             1.7023e-04,
  594.             -9.2424e-04,
  595.             1.1269e-05,
  596.             -6.9911e-06,
  597.             -2.0886e+00,
  598.             -6.7879e-02,
  599.             -8.5922e-04,
  600.             -1.6087e-03,
  601.             -4.5549e-05,
  602.             3.3178e-05,
  603.             -6.1715e-06,
  604.             -1.4446e-06,
  605.             -3.7210e-01,
  606.             1.5775e-01,
  607.             -1.7827e-03,
  608.             -4.4396e-04,
  609.             2.2844e-04,
  610.             -1.1215e-05,
  611.             -2.1120e-06,
  612.             -9.6421e-07,
  613.             -1.4170e-08,
  614.             7.8720e-01,
  615.             -4.4238e-02,
  616.             -1.5120e-03,
  617.             -9.4119e-04,
  618.             4.0645e-06,
  619.             -4.9253e-06,
  620.             -1.8656e-06,
  621.             -4.0736e-07,
  622.             -4.9594e-08,
  623.             1.6134e-09
  624.         };

  626.         /** Amplitude B<sub>nm</sub> coefficients for the temperature. */
  627.         private static final double[] B_TEMPERATURE_AMPLITUDE = {
  628.             0.0000e+00,
  629.             0.0000e+00,
  630.             -8.9895e-01,
  631.             0.0000e+00,
  632.             -1.0790e+00,
  633.             -1.2699e-01,
  634.             0.0000e+00,
  635.             -5.9033e-01,
  636.             3.4865e-02,
  637.             -3.2614e-02,
  638.             0.0000e+00,
  639.             -2.4310e-02,
  640.             1.5607e-02,
  641.             -2.9833e-02,
  642.             -5.9048e-03,
  643.             0.0000e+00,
  644.             2.8383e-01,
  645.             4.0509e-02,
  646.             -1.8834e-02,
  647.             -1.2654e-03,
  648.             -1.3794e-04,
  649.             0.0000e+00,
  650.             1.3306e-01,
  651.             3.4960e-02,
  652.             -3.6799e-03,
  653.             -3.5626e-04,
  654.             1.4814e-04,
  655.             3.7932e-06,
  656.             0.0000e+00,
  657.             2.0801e-01,
  658.             6.5640e-03,
  659.             -3.4893e-03,
  660.             -2.7395e-04,
  661.             7.4296e-05,
  662.             -7.9927e-06,
  663.             -1.0277e-06,
  664.             0.0000e+00,
  665.             3.6515e-02,
  666.             -7.4319e-03,
  667.             -6.2873e-04,
  668.             8.2461e-05,
  669.             3.1095e-05,
  670.             -5.3860e-07,
  671.             -1.2055e-07,
  672.             -1.1517e-07,
  673.             0.0000e+00,
  674.             3.1404e-02,
  675.             1.5580e-02,
  676.             -1.1428e-03,
  677.             3.3529e-05,
  678.             1.0387e-05,
  679.             -1.9378e-06,
  680.             -2.7327e-07,
  681.             7.5833e-09,
  682.             -9.2323e-09
  683.         };

  685.         /** Build a new instance. */
  686.         ABCoefficients() {

  688.         }

  690.         /** Get the value of the mean A<sub>nm</sub> pressure coefficient for the given index.
  691.          * @param index index
  692.          * @return the mean A<sub>nm</sub> pressure coefficient for the given index
  693.          */
  694.         public double getAnmPressureMean(final int index) {
  695.             return A_PRESSURE_MEAN[index];
  696.         }

  698.         /** Get the value of the mean B<sub>nm</sub> pressure coefficient for the given index.
  699.          * @param index index
  700.          * @return the mean B<sub>nm</sub> pressure coefficient for the given index
  701.          */
  702.         public double getBnmPressureMean(final int index) {
  703.             return B_PRESSURE_MEAN[index];
  704.         }

  706.         /** Get the value of the amplitude A<sub>nm</sub> pressure coefficient for the given index.
  707.          * @param index index
  708.          * @return the amplitude A<sub>nm</sub> pressure coefficient for the given index.
  709.          */
  710.         public double getAnmPressureAmpl(final int index) {
  711.             return A_PRESSURE_AMPLITUDE[index];
  712.         }

  714.         /** Get the value of the amplitude B<sub>nm</sub> pressure coefficient for the given index.
  715.          * @param index index
  716.          * @return the amplitude B<sub>nm</sub> pressure coefficient for the given index
  717.          */
  718.         public double getBnmPressureAmpl(final int index) {
  719.             return B_PRESSURE_AMPLITUDE[index];
  720.         }

  722.         /** Get the value of the mean A<sub>nm</sub> temperature coefficient for the given index.
  723.          * @param index index
  724.          * @return the mean A<sub>nm</sub> temperature coefficient for the given index
  725.          */
  726.         public double getAnmTemperatureMean(final int index) {
  727.             return A_TEMPERATURE_MEAN[index];
  728.         }

  730.         /** Get the value of the mean B<sub>nm</sub> temperature coefficient for the given index.
  731.          * @param index index
  732.          * @return the mean B<sub>nm</sub> temperature coefficient for the given index
  733.          */
  734.         public double getBnmTemperatureMean(final int index) {
  735.             return B_TEMPERATURE_MEAN[index];
  736.         }

  738.         /** Get the value of the amplitude A<sub>nm</sub> temperature coefficient for the given index.
  739.          * @param index index
  740.          * @return the amplitude A<sub>nm</sub> temperature coefficient for the given index.
  741.          */
  742.         public double getAnmTemperatureAmpl(final int index) {
  743.             return A_TEMPERATURE_AMPLITUDE[index];
  744.         }

  746.         /** Get the value of the amplitude B<sub>nm</sub> temperature coefficient for the given index.
  747.          * @param index index
  748.          * @return the amplitude B<sub>nm</sub> temperature coefficient for the given index
  749.          */
  750.         public double getBnmTemperatureAmpl(final int index) {
  751.             return B_TEMPERATURE_AMPLITUDE[index];
  752.         }
  753.     }

  755. }
Created with JaCoCo 0.8.2.201808211720