correctmodel.h

计算摄动力模型中关于海潮固体潮极潮以及利用行星星历表计算太阳月亮位置的程序

// CorrectModel.h: interface for the CCorrectModel class.
//
//////////////////////////////////////////////////////////////////////

#if !defined(AFX_CORRECTMODEL_H__19E77690_EF5D_4BF3_B948_9961C516E5D3__INCLUDED_)
#define AFX_CORRECTMODEL_H__19E77690_EF5D_4BF3_B948_9961C516E5D3__INCLUDED_

#if _MSC_VER > 1000
#pragma once
#endif // _MSC_VER > 1000

#include "ReadJPL_EPH.h"
//////////////////////////////////////////////////////////////////////////////////////////////////////
typedef struct  
{
  int N,M;
  double C1,C2,S1,S2; // C1=C+ C2=C- S1=S+ S2=S-
  double Theths;  // 相位角
}
OTideType;
typedef CArray <OTideType,OTideType&> TideArray;
///////////////////////////////////////////////////////////////////////////////////////////////////////
class CCorrectModel  
{
public:
	double GetFnm(int N, int M);
	double GetKn(int N);
	void GetOceanTide(double DeltCoef[]);
	CCorrectModel();
	virtual ~CCorrectModel();
public:
	void CalSunPos();
	void CalMoonPos();
	void SetJUL(double JD);
	void ReadPole();
	void GetPoleTide(double DeltCoef[]);
	void GetSoildTide(double DeltCoef[]);
	double PI;
private:
	CReadJPL_EPH m_ReadJPL;
	double  m_X, m_Y, m_Z; // the position of Satellite
	double  m_Xm, m_Ym, m_Zm; // the position of moon
	double  m_Xs, m_Ys, m_Zs; // the position of sun
	double  Xp,Yp;  // the parameters of pole 
	int m_year,m_month,m_day,m_hour,m_minute; // calculation time of perturb
	double m_second;
	double DayRest;
	double ET; // JUL of calculation for perturb
	double CoefDot[5]; //保存由于潮汐引起的位系数摄动:C20,C21,C22,S21,S22
    TideArray m_OTide; // save the information of ocean tide including N,M,C+,C-,S+,S- 

};

#endif // !defined(AFX_CORRECTMODEL_H__19E77690_EF5D_4BF3_B948_9961C516E5D3__INCLUDED_)