function Helio() { 
  with (Math) {
	var V = 172.74 + 0.00111588*d;
	sV = sin (V*PI/180);
	cV = cos (V*PI/180);
             N = (317.018 + 0.0334443 * d - 0.814 *  sV ) /360;
             N = (N - floor(N)) *2*PI;
             B = 6.366 * sin(N) + 0.211 * sin(2*N) + 0.011 * sin(3*N);
             r = 9.5695 - 0.5304*cos(N)  - 0.0147 *cos(2*N);
             lambda = 50.07 + 0.0334982 * d - 0.814 * sV + B;
			 
	}
}


function Geoz() { 
  with (Math) {

             var J = (50.391 + 0.95214946 * d + 0.814 * sV)/360;  // synodische Länge
             J = (J - floor(J)) *2*PI;
	K = J - B*PI/180 + Cs;
	delta = sqrt( r*r + Rerde*Rerde - 2*r*Rerde*cos(K));
             phi = asin(Rerde * sin(K)/delta);
	lgeoz = lambda*PI/180 - phi;  //  geoz

	beta = 2.49 * sin((lambda - 113.55)*PI/180); 
 
             bgeoz = beta*PI*r/(180*delta); 

	var x = cos(lgeoz)*cos(bgeoz);
	var y = cos(EKL) * cos(bgeoz) * sin(lgeoz) - sin(EKL) * sin(bgeoz);
	var z = sin(EKL) * cos(bgeoz) * sin(lgeoz) + cos(EKL) * sin(bgeoz); 
	alphasa = atan(y/x); if (x<0) alphasa = alphasa + PI
	deltasa = asin(z);
	

	var laenge = lambda*PI/180;
	x = cos(laenge);
	RAs = atan ((cos(EKL)*sin(laenge) - sin(EKL)*tan(bgeoz))/x);
	if (x<0) RAs = RAs + PI
	DEs = asin(sin(EKL)*cos(bgeoz)*sin(laenge) + cos(EKL)*sin(bgeoz));
	}

}


function Phys() { 
  with (Math) {
  
  ///  neu
  
  
  
  var i_ring=28.075216-0.012998*T;
  i_ring=i_ring*PI/180;
  var omega_ring=169.508470+1.394681*T;
  omega_ring=omega_ring*PI/180;
  
  var b_ring_erde=sin(i_ring)*cos(bgeoz)*sin(lgeoz-omega_ring)-cos(i_ring)*sin(bgeoz);
  B0=asin(b_ring_erde);
  
  var b_ring_sonne=sin(i_ring)*cos(beta*PI/180)*sin(lambda*PI/180-omega_ring)-cos(i_ring)*sin(beta*PI/180);
  B0s=asin(b_ring_sonne);
  
  
  
  
  
	A0 = (40.59 + 1.285E-4*d)*PI/180;
	D0 = (83.54 + 1.11E-5*d)*PI/180;
	A0RA = A0-alphasa;
	B0 = B0; // asin(-sin(D0)*sin(deltasa) - cos(D0)*cos(deltasa)*cos(A0RA));
	A0RAs = A0-RAs;
	B0s =B0s; // asin(- sin(D0)*sin(DEs) - cos(D0)*cos(DEs)*cos(A0RAs));

	x = sin(A0RA);
	K0 = atan((sin(D0)*cos(A0RA) - cos(D0)*tan(deltasa))/x); 
	if (x<0) K0 = K0 + PI
	DAequator = 165.46/delta;
	DPol = sqrt(1 - 0.2037*cos(B0)*cos(B0));
	DPol = DAequator * DPol;
	Rax = 375.4/delta;
	Ray = 375.4*abs(sin(B0))/delta;
	mv = -8.68 + 5 * 0.43429 * log(r*delta) + 0.044*phi*180/PI - 2.60 * sin(abs(B0)) + 1.25 * sin(B0) * sin(B0);
	}

}