Added config parameters for estimated cd

src/main/java/fr/cs/examples/estimation/OrbitDetermination.java

@@ -619,6 +619,7 @@ public class OrbitDetermination { // Class 1
             final double  area        = parser.getDouble(ParameterKey.DRAG_AREA);
             final boolean cdEstimated = parser.getBoolean(ParameterKey.DRAG_CD_ESTIMATED);
             System.out.println(cd);
+            System.out.println(parser.getDouble(ParameterKey.DRAG_CD_MAX));
             MarshallSolarActivityFutureEstimation msafe =
                             new MarshallSolarActivityFutureEstimation(MarshallSolarActivityFutureEstimation.DEFAULT_SUPPORTED_NAMES,
                                                                       MarshallSolarActivityFutureEstimation.StrengthLevel.AVERAGE);
@@ -627,13 +628,23 @@ public class OrbitDetermination { // Class 1
             //Atmosphere atmosphere = new DTM2000(msafe, CelestialBodyFactory.getSun(), body);
             Atmosphere atmosphere = new NRLMSISE00(msafe, CelestialBodyFactory.getSun(), body);
             propagatorBuilder.addForceModel(new DragForce(atmosphere, new IsotropicDrag(area, cd)));
+
+
             if (cdEstimated) {
                 for (final ParameterDriver driver : propagatorBuilder.getPropagationParametersDrivers().getDrivers()) {
                     if (driver.getName().equals(DragSensitive.DRAG_COEFFICIENT)) {
                         driver.setSelected(true);
                     }
+                    if (parser.containsKey(ParameterKey.DRAG_CD_MIN)) {
+                        driver.setMinValue(parser.getDouble(ParameterKey.DRAG_CD_MIN));
+                    }
+
+                    if (parser.containsKey(ParameterKey.DRAG_CD_MAX)) {
+                        driver.setMinValue(parser.getDouble(ParameterKey.DRAG_CD_MAX));
+                    }
                 }
             }
+
         }
 
         // solar radiation pressure
@@ -2240,6 +2251,8 @@ public class OrbitDetermination { // Class 1
         DRAG_CD,
         DRAG_CD_ESTIMATED,
         DRAG_AREA,
+        DRAG_CD_MIN,
+        DRAG_CD_MAX,
         SOLAR_RADIATION_PRESSURE,
         SOLAR_RADIATION_PRESSURE_CR,
         SOLAR_RADIATION_PRESSURE_CR_ESTIMATED,

src/main/resources/maxvalier.in

@@ -133,6 +133,9 @@ drag.cd =  2
 drag.cd.estimated = true
 ## Drag area (m^2)
 drag.area = 0.25
+## Maximum and minimum values of the drag coefficient
+drag.cd.min = 0
+drag.cd.max = 20
 
 ## Solar Radiation Pressure (true/false) [false]
 solar.radiation.pressure = false
@@ -588,11 +591,11 @@ estimator.max.iterations                                = 600
 estimator.max.evaluations                               = 600
 
 # comma-separated list of measurements files (in the same directory as this file)
-measurements.files = generated-doppler-39-CGBSAT-VHF.dat
+measurements.files = 2019-06-11T19_58_49_145.961_4171_42778.dat
 
 
 # generated-doppler-F4KLD-UNIVERSITE-PAUL-SABATIER-Toulouse-III.dat
-# 
+# generated-doppler-39-CGBSAT-VHF.dat
 # generated-doppler-39-CGBSAT-VHF-fixeddrag.dat
 # generated-doppler-Marcs.dat 
 #  generated-doppler-ZL1WJQ.dat
@@ -619,7 +622,7 @@ measurements.files = generated-doppler-39-CGBSAT-VHF.dat
 
 #Measured Data
 # 2019-06-12T07_20_42_145.961_4171_42778.dat
-# 2019-06-11T19_58_49_145.961_4171_42778.dat , 2019-06-11T21_32_55_145.961_4171_42778.dat
+#  , 2019-06-11T21_32_55_145.961_4171_42778.dat
 #
 
 # base name of the output files (log and residuals), no files created if empty