JB2008 Model Instantiation

Orekit Python Wrapper usage
2

Apologies for this post, I have now solved the issue.

For anyone who might encounter the same problem this is my current implementation:

        from org.orekit.models.earth.atmosphere.data import JB2008SpaceEnvironmentData, SOLFSMYDataLoader, DtcDataLoader
        from org.orekit.models.earth.atmosphere import JB2008, JB2008InputParameters
        from org.orekit.data import DataSource
        import requests
        from java.io import File

        # Function to download file and return a java.io.File object
        def download_file(url, local_filename):
            with requests.get(url, stream=True) as r:
                r.raise_for_status()
                with open(local_filename, 'wb') as f:
                    for chunk in r.iter_content(chunk_size=8192):
                        f.write(chunk)
            return File(local_filename)

        # Download SOLFSMY and DTCFILE files
        solfsmy_file = download_file("https://sol.spacenvironment.net/JB2008/indices/SOLFSMY.TXT", "external/jb08_inputs/SOLFSMY.TXT")
        dtcfile_file = download_file("https://sol.spacenvironment.net/JB2008/indices/DTCFILE.TXT", "external/jb08_inputs/DTCFILE.TXT")

        # Create DataSource instances
        solfsmy_data_source = DataSource(solfsmy_file)
        dtcfile_data_source = DataSource(dtcfile_file)

        jb08_data = JB2008SpaceEnvironmentData(solfsmy_data_source,
                                            dtcfile_data_source)
        from org.orekit.time import TimeScalesFactory
        utc = TimeScalesFactory.getUTC()
        atmosphere = JB2008(jb08_data, sun, wgs84Ellipsoid, utc)
        isotropicDrag = IsotropicDrag(float(cross_section), float(cd))
        dragForce = DragForce(atmosphere, isotropicDrag)
        force_models.append(dragForce)

Best,
Charles

4 Likes