定量遙感：計算地方時和太陽高度角（C++代碼）

在定量遙感中，通常需要計算地方時和太陽高度角，本文采用C++語言實現。

#include <cmath>
#include <iostream>
#include <fstream>
using namespace std;
void main()
{
	 int JD,NF,Y,R,s[5],F[5];
	 float JF,WD;
	 float N0;
	 ifstream data1("d:\\result\\data-input.txt");
     ofstream out1("d:\\result\\6gaodujiao.txt");
	 data1>>JD>>JF;
	 data1>>WD;
     data1>>NF>>Y>>R;
	 out1<<JD<<"   "<<JF<<endl;
	 out1<<WD<<endl;
	 out1<<NF<<"   "<<Y<<"   "<<R<<endl;
	 out1<<"太陽方位角在正南爲0，由南向西爲0——+180，由南向東爲0——-180。午前取負值，午後取正值。"<<endl;
   N0=79.6764+0.2422*(NF-1985)-int((NF-1985)/4);
   float A,B,G,L;
	  A=NF/4;
	 double C;
	 C=32.8;
	 B=A-int(A);
  	 if(Y<=2)
	 C=30.6;
	 if(B=0 && Y>2 )
	 C=31.8;
	 G=int(30.6*Y-C+0.5)+R;
	 L=(JD+JF/60)/15;
	 float H[5], N[5], theta[5],Er[5],Ed[5],Et[5];
	 double Sd[5];
     float   St[5],  shj