Commit a26fe20d authored by Pascal Parraud's avatar Pascal Parraud

Added NRLMSISE-00 Atmosphere model.

The NRLMSISE-00 model, developed by Mike Picone, Alan Hedin, and Doug
Drob (Naval Research Laboratory), calculates the neutral atmosphere
empirical model from the surface to lower exosphere (0 to 1000 km).

This a Java translation of both original FORTRAN distribution package,
provided by M. Picone & al, and C version, provided by Dominik
Brodowski.
parent 1fe3ec05
......@@ -33,3 +33,9 @@ The Apache Software Foundation (http://www.apache.org/)
This product depends on software developed by
The Hipparchus project (https://hipparchus.org/)
This product includes software translated from original work developed by
M. Picone, A.E. Hedin, D. Drob (Naval Research Laboratory)
This product includes software translated from original work developed by
Dominik Brodowski
......@@ -23,6 +23,7 @@ import org.hipparchus.geometry.euclidean.threed.Vector3D;
import org.orekit.errors.OrekitException;
import org.orekit.errors.OrekitMessages;
import org.orekit.forces.AbstractForceModel;
import org.orekit.forces.drag.atmosphere.Atmosphere;
import org.orekit.frames.Frame;
import org.orekit.frames.Transform;
import org.orekit.propagation.SpacecraftState;
......
......@@ -23,9 +23,7 @@ import org.orekit.bodies.OneAxisEllipsoid;
import org.orekit.errors.OrekitException;
import org.orekit.errors.OrekitMessages;
import org.orekit.frames.Frame;
import org.orekit.frames.Transform;
import org.orekit.time.AbsoluteDate;
import org.orekit.utils.PVCoordinates;
import org.orekit.utils.PVCoordinatesProvider;
......@@ -342,26 +340,6 @@ public class HarrisPriester implements Atmosphere {
return getDensity(sunInEarth, posInEarth);
}
/** Get the inertial velocity of atmosphere molecules.
* <p>
* Here the case is simplified : atmosphere is supposed to have a null velocity
* in earth frame.
* </p>
* @param date current date
* @param position current position
* @param frame the frame in which is defined the position
* @return velocity (m/s) (defined in the same frame as the position)
* @exception OrekitException if some frame conversion cannot be performed
*/
public Vector3D getVelocity(final AbsoluteDate date, final Vector3D position, final Frame frame)
throws OrekitException {
final Transform bodyToFrame = earth.getBodyFrame().getTransformTo(frame, date);
final Vector3D posInBody = bodyToFrame.getInverse().transformPosition(position);
final PVCoordinates pvBody = new PVCoordinates(posInBody, Vector3D.ZERO);
final PVCoordinates pvFrame = bodyToFrame.transformPVCoordinates(pvBody);
return pvFrame.getVelocity();
}
/** Get the height above the Earth for the given position.
* <p>
* The height computation is an approximation valid for the considered atmosphere.
......
This source diff could not be displayed because it is too large. You can view the blob instead.
/* Copyright 2002-2016 CS Systèmes d'Information
* Licensed to CS Systèmes d'Information (CS) under one or more
* contributor license agreements. See the NOTICE file distributed with
* this work for additional information regarding copyright ownership.
* CS licenses this file to You under the Apache License, Version 2.0
* (the "License"); you may not use this file except in compliance with
* the License. You may obtain a copy of the License at
*
* http://www.apache.org/licenses/LICENSE-2.0
*
* Unless required by applicable law or agreed to in writing, software
* distributed under the License is distributed on an "AS IS" BASIS,
* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
* See the License for the specific language governing permissions and
* limitations under the License.
*/
package org.orekit.forces.drag.atmosphere;
import java.io.Serializable;
import org.orekit.errors.OrekitException;
import org.orekit.time.AbsoluteDate;
/** Container for solar activity data, compatible with NRLMSISE-00 atmosphere model.
* <p>
* This model needs daily and average F10.7 solar fluxes and
* A<sub>p</sub> geomagnetic indices to compute the local density.
* </p>
*
* @author Pascal Parraud
*/
public interface NRLMSISE00InputParameters extends Serializable {
/** Gets the available data range minimum date.
* @return the minimum date.
* @exception OrekitException if data cannot be loaded
*/
AbsoluteDate getMinDate() throws OrekitException;
/** Gets the available data range maximum date.
* @return the maximum date.
* @exception OrekitException if data cannot be loaded
*/
AbsoluteDate getMaxDate() throws OrekitException;
/** Get the value of the daily F10.7 solar flux for previous day.
* @param date the current date
* @return the daily F10.7 flux for previous day
* @exception OrekitException if the date is out of range of available data
*/
double getDailyFlux(AbsoluteDate date) throws OrekitException;
/** Get the value of the 81 day average of F10.7 solar flux centered on current day.
* @param date the current date
* @return the 81 day average of F10.7 solar flux centered on current day
* @exception OrekitException if the date is out of range of available data
*/
double getAverageFlux(AbsoluteDate date) throws OrekitException;
/** Get the A<sub>p</sub> geomagnetic indices.
* <p>
* A<sub>p</sub> indices are provided as an array such as:
* <ul>
* <li>0 -> daily A<sub>p</sub></li>
* <li>1 -> 3 hr A<sub>p</sub> index for current time</li>
* <li>2 -> 3 hr A<sub>p</sub> index for 3 hrs before current time</li>
* <li>3 -> 3 hr A<sub>p</sub> index for 6 hrs before current time</li>
* <li>4 -> 3 hr A<sub>p</sub> index for 9 hrs before current time</li>
* <li>5 -> Average of eight 3 hr A<sub>p</sub> indices from 12 to 33 hrs
* prior to current time</li>
* <li>6 -> Average of eight 3 hr A<sub>p</sub> indices from 36 to 57 hrs
* prior to current time</li>
* </ul>
* </p>
* @param date the current date
* @return the array of A<sub>p</sub> indices
* @exception OrekitException if the date is out of range of available data
*/
double[] getAp(AbsoluteDate date) throws OrekitException;
}
/* Copyright 2002-2016 CS Systèmes d'Information
* Licensed to CS Systèmes d'Information (CS) under one or more
* contributor license agreements. See the NOTICE file distributed with
* this work for additional information regarding copyright ownership.
* CS licenses this file to You under the Apache License, Version 2.0
* (the "License"); you may not use this file except in compliance with
* the License. You may obtain a copy of the License at
*
* http://www.apache.org/licenses/LICENSE-2.0
*
* Unless required by applicable law or agreed to in writing, software
* distributed under the License is distributed on an "AS IS" BASIS,
* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
* See the License for the specific language governing permissions and
* limitations under the License.
*/
/**
*
* This package provides the atmosphere model interface and several implementations.
*
*/
package org.orekit.forces.drag.atmosphere;
/* Copyright 2002-2016 CS Systèmes d'Information
* Licensed to CS Systèmes d'Information (CS) under one or more
* contributor license agreements. See the NOTICE file distributed with
* this work for additional information regarding copyright ownership.
* CS licenses this file to You under the Apache License, Version 2.0
* (the "License"); you may not use this file except in compliance with
* the License. You may obtain a copy of the License at
*
* http://www.apache.org/licenses/LICENSE-2.0
*
* Unless required by applicable law or agreed to in writing, software
* distributed under the License is distributed on an "AS IS" BASIS,
* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
* See the License for the specific language governing permissions and
* limitations under the License.
*/
/**
*
* This package provides classes to get solar activity data,
* including solar flux and geomagnetic indices.
*
*/
package org.orekit.forces.drag.atmosphere.solaractivity;
......@@ -18,9 +18,5 @@
*
* This package provides all drag-related forces.
*
* <p>
* This includes both the atmosphere model interfaces and several implementations,
* as well as the force models that uses the atmosphere to compute drag.
* </p>
*/
package org.orekit.forces.drag;
......@@ -19,10 +19,10 @@ package org.orekit.propagation.semianalytical.dsst.forces;
import org.hipparchus.util.FastMath;
import org.hipparchus.util.MathUtils;
import org.orekit.errors.OrekitException;
import org.orekit.forces.drag.Atmosphere;
import org.orekit.forces.drag.DragForce;
import org.orekit.forces.drag.DragSensitive;
import org.orekit.forces.drag.IsotropicDrag;
import org.orekit.forces.drag.atmosphere.Atmosphere;
import org.orekit.propagation.SpacecraftState;
import org.orekit.propagation.events.EventDetector;
import org.orekit.utils.Constants;
......
......@@ -30,8 +30,8 @@ import java.util.TreeSet;
import org.hipparchus.exception.DummyLocalizable;
import org.orekit.errors.OrekitException;
import org.orekit.errors.OrekitMessages;
import org.orekit.forces.drag.DTM2000InputParameters;
import org.orekit.forces.drag.JB2006InputParameters;
import org.orekit.forces.drag.atmosphere.DTM2000InputParameters;
import org.orekit.forces.drag.atmosphere.JB2006InputParameters;
import org.orekit.time.AbsoluteDate;
import org.orekit.time.ChronologicalComparator;
import org.orekit.time.TimeScalesFactory;
......
......@@ -19,7 +19,7 @@ package org.orekit.estimation;
import org.orekit.errors.OrekitException;
import org.orekit.forces.ForceModel;
import org.orekit.forces.drag.DragForce;
import org.orekit.forces.drag.HarrisPriester;
import org.orekit.forces.drag.atmosphere.HarrisPriester;
import org.orekit.forces.gravity.HolmesFeatherstoneAttractionModel;
import org.orekit.forces.gravity.OceanTides;
import org.orekit.forces.gravity.Relativity;
......
......@@ -65,12 +65,12 @@ import org.orekit.estimation.measurements.Range;
import org.orekit.estimation.measurements.RangeRate;
import org.orekit.estimation.measurements.modifiers.AngularRadioRefractionModifier;
import org.orekit.estimation.measurements.modifiers.RangeTroposphericDelayModifier;
import org.orekit.forces.drag.Atmosphere;
import org.orekit.forces.drag.DTM2000;
import org.orekit.forces.drag.DragForce;
import org.orekit.forces.drag.DragSensitive;
import org.orekit.forces.drag.IsotropicDrag;
import org.orekit.forces.drag.MarshallSolarActivityFutureEstimation;
import org.orekit.forces.drag.atmosphere.Atmosphere;
import org.orekit.forces.drag.atmosphere.DTM2000;
import org.orekit.forces.drag.atmosphere.solaractivity.MarshallSolarActivityFutureEstimation;
import org.orekit.forces.gravity.HolmesFeatherstoneAttractionModel;
import org.orekit.forces.gravity.OceanTides;
import org.orekit.forces.gravity.Relativity;
......
......@@ -33,6 +33,9 @@ import org.orekit.bodies.OneAxisEllipsoid;
import org.orekit.errors.OrekitException;
import org.orekit.forces.AbstractForceModelTest;
import org.orekit.forces.BoxAndSolarArraySpacecraft;
import org.orekit.forces.drag.atmosphere.Atmosphere;
import org.orekit.forces.drag.atmosphere.HarrisPriester;
import org.orekit.forces.drag.atmosphere.SimpleExponentialAtmosphere;
import org.orekit.frames.Frame;
import org.orekit.frames.FramesFactory;
import org.orekit.orbits.CartesianOrbit;
......
......@@ -24,6 +24,7 @@ import org.junit.Test;
import org.orekit.Utils;
import org.orekit.bodies.OneAxisEllipsoid;
import org.orekit.errors.OrekitException;
import org.orekit.forces.drag.atmosphere.SimpleExponentialAtmosphere;
import org.orekit.frames.Frame;
import org.orekit.frames.FramesFactory;
import org.orekit.frames.Transform;
......
......@@ -41,11 +41,11 @@ import org.orekit.Utils;
import org.orekit.bodies.OneAxisEllipsoid;
import org.orekit.errors.OrekitException;
import org.orekit.forces.ForceModel;
import org.orekit.forces.drag.Atmosphere;
import org.orekit.forces.drag.DragForce;
import org.orekit.forces.drag.DragSensitive;
import org.orekit.forces.drag.IsotropicDrag;
import org.orekit.forces.drag.SimpleExponentialAtmosphere;
import org.orekit.forces.drag.atmosphere.Atmosphere;
import org.orekit.forces.drag.atmosphere.SimpleExponentialAtmosphere;
import org.orekit.forces.gravity.HolmesFeatherstoneAttractionModel;
import org.orekit.forces.gravity.NewtonianAttraction;
import org.orekit.forces.gravity.potential.GravityFieldFactory;
......
......@@ -31,11 +31,11 @@ import org.orekit.bodies.OneAxisEllipsoid;
import org.orekit.errors.OrekitException;
import org.orekit.errors.OrekitMessages;
import org.orekit.forces.ForceModel;
import org.orekit.forces.drag.Atmosphere;
import org.orekit.forces.drag.DragForce;
import org.orekit.forces.drag.DragSensitive;
import org.orekit.forces.drag.IsotropicDrag;
import org.orekit.forces.drag.SimpleExponentialAtmosphere;
import org.orekit.forces.drag.atmosphere.Atmosphere;
import org.orekit.forces.drag.atmosphere.SimpleExponentialAtmosphere;
import org.orekit.forces.gravity.HolmesFeatherstoneAttractionModel;
import org.orekit.forces.gravity.NewtonianAttraction;
import org.orekit.forces.gravity.potential.GravityFieldFactory;
......
......@@ -38,8 +38,8 @@ import org.orekit.errors.OrekitMessages;
import org.orekit.forces.ForceModel;
import org.orekit.forces.drag.DragForce;
import org.orekit.forces.drag.DragSensitive;
import org.orekit.forces.drag.HarrisPriester;
import org.orekit.forces.drag.IsotropicDrag;
import org.orekit.forces.drag.atmosphere.HarrisPriester;
import org.orekit.forces.gravity.HolmesFeatherstoneAttractionModel;
import org.orekit.forces.gravity.NewtonianAttraction;
import org.orekit.forces.gravity.ThirdBodyAttraction;
......
......@@ -42,8 +42,8 @@ import org.orekit.bodies.CelestialBody;
import org.orekit.bodies.CelestialBodyFactory;
import org.orekit.bodies.OneAxisEllipsoid;
import org.orekit.errors.OrekitException;
import org.orekit.forces.drag.Atmosphere;
import org.orekit.forces.drag.HarrisPriester;
import org.orekit.forces.drag.atmosphere.Atmosphere;
import org.orekit.forces.drag.atmosphere.HarrisPriester;
import org.orekit.forces.gravity.potential.GRGSFormatReader;
import org.orekit.forces.gravity.potential.GravityFieldFactory;
import org.orekit.forces.gravity.potential.ICGEMFormatReader;
......
......@@ -69,12 +69,12 @@ import org.orekit.estimation.measurements.Range;
import org.orekit.estimation.measurements.RangeRate;
import org.orekit.estimation.measurements.modifiers.AngularRadioRefractionModifier;
import org.orekit.estimation.measurements.modifiers.RangeTroposphericDelayModifier;
import org.orekit.forces.drag.Atmosphere;
import org.orekit.forces.drag.DTM2000;
import org.orekit.forces.drag.DragForce;
import org.orekit.forces.drag.DragSensitive;
import org.orekit.forces.drag.IsotropicDrag;
import org.orekit.forces.drag.MarshallSolarActivityFutureEstimation;
import org.orekit.forces.drag.atmosphere.Atmosphere;
import org.orekit.forces.drag.atmosphere.DTM2000;
import org.orekit.forces.drag.atmosphere.solaractivity.MarshallSolarActivityFutureEstimation;
import org.orekit.forces.gravity.HolmesFeatherstoneAttractionModel;
import org.orekit.forces.gravity.OceanTides;
import org.orekit.forces.gravity.Relativity;
......
......@@ -43,11 +43,11 @@ import org.orekit.bodies.OneAxisEllipsoid;
import org.orekit.data.DataProvidersManager;
import org.orekit.data.DirectoryCrawler;
import org.orekit.errors.OrekitException;
import org.orekit.forces.drag.Atmosphere;
import org.orekit.forces.drag.DragForce;
import org.orekit.forces.drag.DragSensitive;
import org.orekit.forces.drag.HarrisPriester;
import org.orekit.forces.drag.IsotropicDrag;
import org.orekit.forces.drag.atmosphere.Atmosphere;
import org.orekit.forces.drag.atmosphere.HarrisPriester;
import org.orekit.forces.gravity.HolmesFeatherstoneAttractionModel;
import org.orekit.forces.gravity.ThirdBodyAttraction;
import org.orekit.forces.gravity.potential.GravityFieldFactory;
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment