physicslib.h

Go to the documentation of this file.

/********************************************************************
* Copyright (C) 2015 by Interstel Technologies, Inc.
*   and Hawaii Space Flight Laboratory.
*
* This file is part of the COSMOS/core that is the central
* module for COSMOS. For more information on COSMOS go to
* <http://cosmos-project.com>
*
* The COSMOS/core software is licenced under the
* GNU Lesser General Public License (LGPL) version 3 licence.
*
* You should have received a copy of the
* GNU Lesser General Public License
* If not, go to <http://www.gnu.org/licenses/>
*
* COSMOS/core is free software: you can redistribute it and/or
* modify it under the terms of the GNU Lesser General Public License
* as published by the Free Software Foundation, either version 3 of
* the License, or (at your option) any later version.
*
* COSMOS/core is distributed in the hope that it will be useful, but
* WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
* Lesser General Public License for more details.
*
* Refer to the "licences" folder for further information on the
* condititons and terms to use this software.
********************************************************************/

#ifndef _ORBITLIB_H
#define _ORBITLIB_H 1


#include "support/configCosmos.h"

#include "support/cosmos-defs.h"
#include "math/mathlib.h"
#include "physics/nrlmsise-00.h"
#include "support/convertlib.h"
#include "support/demlib.h"
#include "physics/physicsdef.h"
#include "support/jsondef.h"

#include <fcntl.h>
#include <cmath>
#include <time.h>
#include <errno.h>



void pleph_(double [], long *, long *, double []);
void dpleph_(double [], long *, long *, double []);

rvector gravity_vector(svector pos,int model,uint32_t degree);
double gravity_potential(double lon, double lat, double r,int model,uint32_t degree);
rvector gravity_accel(posstruc pos, int model, uint32_t degree);
rvector gravity_accel2(posstruc pos, int model, uint32_t degree);
double gravity(double radius, double colat, double elon, int model, uint32_t degree);
int32_t gravity_params(int model);
double nplgndr(uint32_t l, uint32_t m, double x);
svector groundstation(locstruc &satellite, locstruc &groundstation);
void simulate_hardware(cosmosstruc *cinfo, locstruc &loc);
void simulate_hardware(cosmosstruc *cinfo, vector <locstruc> &locvec);
void initialize_imu(uint16_t index, devspecstruc &devspec, locstruc &loc);
void simulate_imu(int index, cosmosstruc *root, locstruc &loc);
int32_t pos_accel(physicsstruc &physics, locstruc &loc);
void att_accel(physicsstruc &physics, locstruc &loc);
void geod2icrf(posstruc *pos);
double msis86_density(posstruc pos,float f107avg,float f107,float magidx);
double msis00_density(posstruc pos,float f107avg,float f107,float magidx);
void orbit_init_tle(int32_t mode, double dt, double mjd, cosmosstruc *root);
void orbit_init_eci(int32_t mode, double dt, double mjd, cartpos ipos, cosmosstruc *root);
void orbit_init_shape(int32_t mode, double dt, double mjd, double altitude, double angle, double hour, cosmosstruc *root);
void propagate(cosmosstruc *root, double mjd);
double rearth(double lat);
int update_eci(cosmosstruc *root, double utc, cartpos pos);

void hardware_init_eci(cosmosstruc *cinfo, locstruc &loc);
void gauss_jackson_setup(gj_handle &gjh, uint32_t order, double utc, double &dt);
void gauss_jackson_init_tle(gj_handle &gjh, uint32_t order, int32_t mode, double dt, double mjd, cosmosstruc *cinfo);
void gauss_jackson_init_eci(gj_handle &gjh, uint32_t order, int32_t mode, double dt, double mjd, cartpos ipos, qatt iatt, physicsstruc &physics, locstruc &loc);
void gauss_jackson_init_stk(gj_handle &gjh, uint32_t order, int32_t mode, double dt, double mjd, stkstruc &stk, physicsstruc &physics, locstruc &loc);
void gauss_jackson_init(gj_handle &gjh, uint32_t order, int32_t mode, double dt, double mjd, double altitude, double angle, double hour, locstruc &iloc, physicsstruc &physics, locstruc &loc);
locstruc gauss_jackson_converge_orbit(gj_handle &gjh, physicsstruc &physics);
void gauss_jackson_converge_hardware(gj_handle &gjh, physicsstruc &physics);
vector<locstruc> gauss_jackson_propagate(gj_handle &gjh, physicsstruc &physics, locstruc &loc, double mjd);
int orbit_propagate(cosmosstruc *root, double mjd);
int orbit_init(int32_t mode, double dt, double mjd, string ofile, cosmosstruc *root);
//void SolidTide(posstruc pos, double dc[5][4], double ds[5][4]);


#endif