shift_solid_thermal.cpp File Reference

#include "support/configCosmos.h"
#include "support/jsonclass.h"
#include "support/datalib.h"
#include "support/timelib.h"
#include "support/convertlib.h"

Include dependency graph for shift_solid_thermal.cpp:

Functions
int	main (int argc, char *argv[])

Function Documentation

int main	(	int	argc,
		char *	argv[]
	)

{
    string filename = argv[1];
    double deltat = atof(argv[2])/86400.;

    string input;

    FILE *fp = fopen(filename.c_str(), "r");
    int ic;
    int32_t id;
    int32_t im;
    int32_t iy;
    int32_t ihh;
    int32_t imm;
    int32_t iss;
    float x;
    float y;
    float z;
    float a;
    float b;
    float c;
    double initialutc = 0.;

    while ((ic=fscanf(fp, "%d/%d/%d %d:%d:%d %f %f %f %f %f %f\n", &im, &id, &iy, &ihh, &imm, &iss, &x, &y, &z, &a, &b, &c)) != EOF)
    {
        locstruc loc;
        loc.pos.eci.utc = cal2mjd(iy, im, id, ihh, imm, iss);
        if (initialutc == 0.)
        {
            initialutc = loc.pos.eci.utc;
        }
        loc.pos.eci.s.col[0] = x;
        loc.pos.eci.s.col[1] = y;
        loc.pos.eci.s.col[2] = z;
        loc.pos.eci.pass++;
        pos_eci(loc);
        avector abc(RADOF(a), RADOF(b), RADOF(c));
        loc.att.icrf.utc = loc.pos.eci.utc;
        loc.att.icrf.s = q_euler2quaternion(abc);
        loc.att.icrf.pass++;
        att_icrf(loc);

        loc.pos.geoc.utc += deltat;
        loc.pos.geoc.pass++;
        pos_geoc(loc);
        loc.att.geoc.utc += deltat;
        loc.att.geoc.pass++;
        att_geoc(loc);

        abc = a_quaternion2euler(loc.att.icrf.s);
        calstruc date = mjd2cal(loc.pos.eci.utc);

        printf("%d/%d/%d %d:%02d:%02d\t%.0f\t%f\t%f\t%f\t%f\t%f\t%f\n", date.month, date.dom, date.year, date.hour, date.minute, date.second, 86400.*(loc.pos.eci.utc - initialutc), loc.pos.eci.s.col[0], loc.pos.eci.s.col[1], loc.pos.eci.s.col[2], DEGOF(abc.h), DEGOF(abc.e), DEGOF(abc.b));
    };

    fclose(fp);

}