////////////////////////////////////////////////////////////////
//
//	人工的にコード位相を生成
//
////////////////////////////////////////////////////////////////

#include <stdio.h>
#include <math.h>

const double PI       = 3.1415926535898;//円周率
const double CS       = 299792458.0;//光速
const double F1		  = 1575420000.0;
const double F2       = 1227600000.0;

const int SATMAX             = 36;//GPSの最大総数
const int GPS                = 32;//GPSの最大総数

#define DIV1     1//何個(縦横)にグリッドを区切るか
#define DIVMAX  4//グリッドを格納する配列の最大値を設定

extern FILE *dd_code;

extern int Enable[DIVMAX][DIVMAX],Enable_sat[DIVMAX][DIVMAX][SATMAX];//使用可能衛星
extern int Actual_satn;//実際のＧＰＳ，ＧＥＯ衛星の運用数を入れる
extern int Enable_sat_stock[DIVMAX][DIVMAX][SATMAX];//2000/11/04
//予備
extern int Enable2[DIVMAX][DIVMAX];
extern int Enable_sat2[DIVMAX][DIVMAX][SATMAX];//使用可能衛星
extern int Total_sat;//トータル衛星数
extern int Ref_sat;//最大仰角衛星のこと
extern int UNUSE_GPS[16];//アルマナック情報のない衛星をストックする
extern int Com_sat[SATMAX],Com_satn;
extern int Integer_r1[SATMAX],Integer_m1[SATMAX],DD_integer1[SATMAX];
extern int QZSS;
extern int CS_check2[SATMAX],CS_check[SATMAX];

extern double Pr1ref[SATMAX],Pr1mov[SATMAX],Pr2ref[SATMAX],Pr2mov[SATMAX];
extern double DD_code1[SATMAX];
extern double Sat_xpos[SATMAX],Sat_ypos[SATMAX],Sat_zpos[SATMAX];//地球中心座標系
extern double POSx[DIVMAX][DIVMAX],POSy[DIVMAX][DIVMAX],POSz[DIVMAX][DIVMAX];//解析位置
extern double Cn[SATMAX],GPSTIME;

void make_dd_code()
{
	int i,k=0,j;
	int prn;

	//再度基準局と移動局における共通衛星を抜き出す
	for(i=0;i<Enable[0][0];i++){
		for(j=0;j<Enable[DIV1][DIV1];j++){
			if(Enable_sat[0][0][i] == Enable_sat[DIV1][DIV1][j]){
				prn = Enable_sat[0][0][i];
				if(CS_check[prn] == 1 && CS_check2[prn] == 1)
					Com_sat[k++] = Enable_sat[0][0][i];
			}
		}
	}
	Com_satn = k;

	//二重位相差の計算
	for(i=0;i<Com_satn;i++){
		prn = Com_sat[i];
		if(Com_sat[i] != Ref_sat){
			DD_code1[prn] = (Pr1mov[prn] - Pr1ref[prn])
				-(Pr1mov[Ref_sat] - Pr1ref[Ref_sat]);
		}
	}
}