1   /* Copyright 2002-2025 CS GROUP
2    * Licensed to CS GROUP (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.troposphere;
18  
19  import org.hipparchus.util.FastMath;
20  import org.junit.jupiter.api.Test;
21  import org.orekit.bodies.GeodeticPoint;
22  import org.orekit.time.AbsoluteDate;
23  import org.orekit.time.TimeScalesFactory;
24  import org.orekit.utils.TrackingCoordinates;
25  
26  public class ViennaOneMappingFunctionTest extends AbstractMappingFunctionTest<ViennaOne> {
27  
28      protected ViennaOne buildMappingFunction() {
29          return new ViennaOne(new ConstantViennaAProvider(new ViennaACoefficients(0.00127683, 0.00060955)),
30                               new ConstantAzimuthalGradientProvider(null),
31                               new ConstantTroposphericModel(new TroposphericDelay(2.0966, 0.2140, 0, 0)),
32                               TimeScalesFactory.getUTC());
33      }
34  
35      @Test
36      public void testMappingFactors() {
37          // Site (NRAO, Green Bank, WV): latitude:  38°
38          //                              longitude: 280°
39          //                              height:    824.17 m
40          //
41          // Date: MJD 55055 -> 12 August 2009 at 0h UT
42          //
43          // Ref for the inputs:    Petit, G. and Luzum, B. (eds.), IERS Conventions (2010),
44          //                        IERS Technical Note No. 36, BKG (2010)
45          //
46          // Values: ah  = 0.00127683
47          //         aw  = 0.00060955
48          //         zhd = 2.0966 m
49          //         zwd = 0.2140 m
50          //
51          // Values taken from: http://vmf.geo.tuwien.ac.at/trop_products/GRID/2.5x2/VMF1/VMF1_OP/2009/VMFG_20090812.H00
52          //
53          // Expected mapping factors : hydrostatic -> 3.425088
54          //                                    wet -> 3.448300
55          //
56          // Expected outputs are obtained by performing the Matlab script vmf1_ht.m provided by TU WIEN:
57          // http://vmf.geo.tuwien.ac.at/codes/
58          //
59          doTestMappingFactors(AbsoluteDate.createMJDDate(55055, 0, TimeScalesFactory.getUTC()),
60                               new GeodeticPoint(FastMath.toRadians(38.0), FastMath.toRadians(280.0), 824.17),
61                               new TrackingCoordinates(0.0, 0.5 * FastMath.PI - 1.278564131, 0.0),
62                               3.425088, 3.448300);
63      }
64  
65      @Test
66      @Override
67      public void testDerivatives() {
68          doTestDerivatives(5.0e-16, 2.0e-19, 1.0e-100, 3.0e-8, 1.0e-100);
69      }
70  
71  }