placalc.c

占星术4.0源码

	 1,	-1,	.711,	35.47,	491,	305.28,	4,
	 2,	-3,	.552,	161.15,	712,	71.12,	4,
	 2,	-2,	2.100,	161.15,	2370,	71.19,	4,
	 2,	-1,	3.724,	160.69,	899,	70.49,	4,
	 2,	0,	.729,	159.76,	763,	250.00,	4,
	 3,	-3,	.431,	105.37,	541,	15.53,	4,
	 3,	-2,	1.329,	104.78,	1157,	14.84,	4,
	 3,	-1,	.539,	278.95,	14,	282.00,	4,
	 4,	-2,	.484,	226.40,	234,	136.02,	4,
	 5,	-4,	.685,	-10.43,	849,	259.51,	4,
	 5,	-3,	2.810,	-10.14,	2954,	259.92,	4,
	 5,	-2,	7.356,	-12.22,	282,	255.43,	4,
	 5,	-1,	1.471,	-12.30,	1550,	77.75,	4,
	 5,	0,	.375,	-12.29,	472,	77.70,	4,
	 2,	-1,	.443,	218.48,	256,	128.36,	5,
	 4,	-2,	.374,	151.81,	397,	61.63,	5,
	 4,	-1,	.808,	145.93,	13,	35.00,	5,
	 1,	-1,	.697,	181.07,	708,	91.38,	6,
	 1,	0,	.574,	236.72,	75,	265.40,	6,
	 2,	-2,	.938,	36.98,	1185,	306.97,	6,
	 2,	-1,	3.275,	37.00,	3268,	306.99,	6,
	 2,	0,	.499,	31.91,	371,	126.90,	6,
	 3,	-1,	.353,	25.84,	347,	295.76,	6,
	 2,	-1,	.380,	239.87,	0,	0,	7,
	 ENDMARK
};

struct kor venuskor[] = {
	 -1,	2,	.264,	-19.20,	175,	251.10,	8,
	 -2,	5,	.361,	167.68,	55,	77.20,	8,
	 1,	-1,	4.889,	119.11,	2246,	29.11,	9,
	 2,	-2,	11.261,	148.23,	9772,	58.21,	9,
	 3,	-3,	7.128,	-2.57,	8271,	267.42,	9,
	 3,	-2,	3.446,	135.91,	737,	47.37,	9,
	 4,	-4,	1.034,	26.54,	1426,	296.49,	9,
	 4,	-3,	.677,	165.32,	445,	75.70,	9,
	 5,	-5,	.330,	56.88,	510,	-33.36,	9,
	 5,	-4,	1.575,	193.93,	1572,	104.21,	9,
	 5,	-3,	1.439,	138.08,	162,	229.90,	9,
	 6,	-6,	.143,	84.40,	236,	-5.80,	9,
	 6,	-5,	.205,	44.20,	256,	314.20,	9,
	 6,	-4,	.176,	164.30,	70,	75.70,	9,
	 8,	-5,	.231,	180.00,	25,	75.00,	9,
	 3,	-2,	.673,	221.62,	717,	131.60,	10,
	 3,	-1,	1.208,	237.57,	29,	149.00,	10,
	 1,	-1,	2.966,	208.09,	2991,	118.09,	11,
	 1,	0,	1.563,	268.31,	91,	-7.60,	11,
	 2,	-2,	.889,	145.16,	1335,	55.17,	11,
	 2,	-1,	.480,	171.01,	464,	80.95,	11,
	 3,	-2,	.169,	144.20,	250,	54.00,	11,
	 ENDMARK
};

struct kor marskor[] = {
	 -1,	1,	.115,	65.84,	684,	156.14,	12,
	 -1,	2,	.623,	246.03,	812,	155.77,	12,
	 -1,	3,	6.368,	57.60,	556,	-32.06,	12,
	 -1,	4,	.588,	57.24,	616,	147.28,	12,
	 -2,	5,	.138,	39.18,	157,	309.39,	12,
	 -2,	6,	.459,	217.58,	82,	128.10,	12,
	 -1,	-1,	.106,	33.60,	141,	303.45,	13,
	 -1,	0,	.873,	34.34,	1112,	304.05,	13,
	 -1,	1,	8.559,	35.10,	6947,	304.45,	13,
	 -1,	2,	13.966,	20.50,	2875,	113.20,	13,
	 -1,	3,	1.487,	22.18,	1619,	112.38,	13,
	 -1,	4,	.175,	22.46,	225,	112.15,	13,
	 -2,	2,	.150,	18.96,	484,	266.42,	13,
	 -2,	3,	7.355,	158.64,	6412,	68.62,	13,
	 -2,	4,	4.905,	154.09,	1985,	244.70,	13,
	 -2,	5,	.489,	154.39,	543,	244.50,	13,
	 -3,	3,	.216,	111.06,	389,	21.10,	13,
	 -3,	4,	.355,	110.64,	587,	19.17,	13,
	 -3,	5,	2.641,	280.58,	2038,	190.60,	13,
	 -3,	6,	.970,	276.06,	587,	6.75,	13,
	 -3,	7,	.100,	276.20,	116,	6.40,	13,
	 -4,	5,	.152,	232.48,	259,	142.60,	13,
	 -4,	6,	.264,	230.47,	387,	139.75,	13,
	 -4,	7,	1.156,	41.64,	749,	312.67,	13,
	 -4,	8,	.259,	37.92,	205,	128.80,	13,
	 -5,	8,	.172,	-8.99,	234,	260.70,	13,
	 -5,	9,	.575,	164.48,	308,	74.60,	13,
	 -6,	10,	.115,	113.70,	145,	23.53,	13,
	 -6,	11,	.363,	285.69,	144,	196.00,	13,
	 -7,	13,	.353,	48.83,	85,	319.10,	13,
	 -8,	15,	1.553,	170.14,	110,	81.00,	13,
	 -8,	16,	.148,	170.74,	154,	259.94,	13,
	 -9,	17,	.193,	293.70,	23,	22.80,	13,
	 1,	-3,	.382,	46.48,	521,	316.25,	14,
	 1,	-2,	3.144,	46.78,	3894,	316.39,	14,
	 1,	-1,	25.384,	48.96,	23116,	318.87,	14,
	 1,	0,	3.732,	-17.62,	1525,	117.81,	14,
	 1,	1,	.474,	-34.60,	531,	59.67,	14,
	 2,	-4,	.265,	192.88,	396,	103.12,	14,
	 2,	-3,	2.108,	192.72,	3042,	102.89,	14,
	 2,	-2,	16.035,	191.90,	22144,	101.99,	14,
	 2,	-1,	21.869,	188.35,	16624,	98.33,	14,
	 2,	0,	1.461,	189.66,	1478,	279.04,	14,
	 2,	1,	.167,	191.04,	224,	280.81,	14,
	 3,	-4,	.206,	167.11,	338,	76.13,	14,
	 3,	-3,	1.309,	168.27,	2141,	76.24,	14,
	 3,	-2,	2.607,	228.41,	3437,	139.74,	14,
	 3,	-1,	3.174,	207.20,	1915,	115.83,	14,
	 3,	0,	.232,	207.78,	240,	298.06,	14,
	 4,	-4,	.178,	127.25,	322,	36.16,	14,
	 4,	-3,	.241,	200.69,	389,	110.02,	14,
	 4,	-2,	.330,	267.57,	413,	179.86,	14,
	 4,	-1,	.416,	221.88,	184,	128.17,	14,
	 1,	-2,	.155,	-38.20,	191,	231.58,	15,
	 1,	-1,	1.351,	-34.10,	1345,	235.85,	15,
	 1,	0,	.884,	288.05,	111,	39.90,	15,
	 1,	1,	.132,	284.88,	144,	15.67,	15,
	 2,	-2,	.620,	35.15,	869,	305.30,	15,
	 2,	-1,	1.768,	32.50,	1661,	302.51,	15,
	 2,	0,	.125,	18.73,	103,	119.90,	15,
	 3,	-2,	.141,	47.59,	199,	318.06,	15,
	 3,	-1,	.281,	40.95,	248,	310.75,	15,
	 ENDMARK
};

#define NUM_MOON_CORR 93
/* moon correction data; revised 30-jul-88: all long. to 0.3" */
struct m45dat {		
	int  i0,i1,i2,i3;
	REAL8 lng,lat,par;
	} m45 [NUM_MOON_CORR] = {  
    /*	  l,     l',	F,	D,      Long,	   Lat,	        Par),*/
	{ 0,	 0,	0,	 4, 	 13.902,     14.06,	 0.2607},
	{ 0,	 0,	0,	 2,    2369.912,   2373.36,     28.2333},
	{ 1,	 0,	0,	 4, 	  1.979,      6.98,	 0.0433},
	{ 1,	 0,	0,	 2,     191.953,    192.72,	 3.0861},
	{ 1,	 0,	0,	 0,   22639.500,  22609.1,     186.5398},
	{ 1,	 0,	0,	-2,   -4586.465,  -4578.13,     34.3117},
	{ 1,	 0,	0,	-4,	-38.428,    -38.64,	 0.6008},
	{ 1,	 0,	0,	-6,	 -0.393,     -1.43,	 0.0086},
	{ 0,	 1,	0,	 4,	 -0.289,     -1.59,	-0.0053},
	{ 0,	 1,	0,	 2,	-24.420,    -25.10,	-0.3000},
	{ 0,	 1,	0,	 0,    -668.146,   -126.98,     -0.3997},
	{ 0,	 1,	0,	-2,    -165.145,   -165.06,      1.9178},
	{ 0,	 1,	0,	-4,	 -1.877,     -6.46,	 0.0339},
	{ 0,	 0,	0,	 3,	  0.403,     -4.01,	 0.0023},
	{ 0,	 0,	0,	 1,    -125.154,   -112.79,     -0.9781},
	{ 2,	 0,	0,	 4,	  0.213,      1.02,	 0.0054},
	{ 2,	 0,	0,	 2,	 14.387,     14.78,	 0.2833},
	{ 2,	 0,	0,	 0,	769.016,    767.96,	10.1657},
	{ 2,	 0,	0,	-2,    -211.656,   -152.53,     -0.3039},
	{ 2,	 0,	0,	-4,	-30.773,    -34.07,	 0.3722},
	{ 2,	 0,	0,	-6,	 -0.570,     -1.40,	 0.0109},
	{ 1,	 1,	0,	 2,	 -2.921,    -11.75,	-0.0484},
	{ 1,	 1,	0,	 0,    -109.673,   -115.18,     -0.9490},
	{ 1,	 1,	0,	-2,    -205.962,   -182.36,      1.4437},
	{ 1,	1,	0,	-4,	 -4.391,     -9.66,	 0.0673},
	{ 1,	-1,	0,	4,	  0.283,      1.53,	 0.0060},
	{ 1,	-1,	0,	2,	 14.577,     31.70,	 0.2302},
	{ 1,	-1,	0,	0,	147.687,    138.76,	 1.1528},
	{ 1,	-1,	0,	-2,	 28.475,     23.59,	-0.2257},
	{ 1,	-1,	0,	-4,	  0.636,      2.27,	-0.0102},
	{ 0,	2,	0,	2,	 -0.189,     -1.68,	-0.0028},
	{ 0,	2,	0,	0,	 -7.486,     -0.66,	-0.0086},
	{ 0,	2,	0,	-2,	 -8.096,    -16.35,	 0.0918},
	{ 0,	0,	2,	2,	 -5.741,     -0.04,	-0.0009},
	{ 0,	0,	2,	0,	-411.608,    -0.2,      -0.0124},
	{ 0,	0,	2,	-2,	-55.173,    -52.14,	-0.1052},
	{ 0,	0,	2,	-4,	  0.025,     -1.67,	 0.0031},
	{ 1,	0,	0,	1,	 -8.466,    -13.51,	-0.1093},
	{ 1,	0,	0,	-1,	 18.609,      3.59,	 0.0118},
	{ 1,	0,	0,	-3,	  3.215,      5.44,	-0.0386},
	{ 0,	1,	0,	1,	 18.023,     17.93,	 0.1494},
	{ 0,	1,	0,	-1,	  0.560,      0.32,	-0.0037},
	{ 3,	0,	0,	2,	  1.060,      2.96,	 0.0243},
	{ 3,	0,	0,	0,	 36.124,     50.64,	 0.6215},
	{ 3,	0,	0,	-2,	-13.193,    -16.40,	-0.1187},
	{ 3,	0,	0,	-4,	 -1.187,     -0.74,	 0.0074},
	{ 3,	0,	0,	-6,	 -0.293,     -0.31,	 0.0046},
	{ 2,	1,	0,	2,	 -0.290,     -1.45,	-0.0051},
	{ 2,	1,	0,	0,	 -7.649,    -10.56,	-0.1038},
	{ 2,	1,	0,	-2,	 -8.627,     -7.59,	-0.0192},
	{ 2,	1,	0,	-4,	 -2.740,     -2.54,	 0.0324},
	{ 2,	-1,	0,	2,	  1.181,      3.32,	 0.0213},
	{ 2,	-1,	0,	0,	  9.703,     11.67,	 0.1268},
	{ 2,	-1,	0,	-2,	 -2.494,     -1.17,	-0.0017},
	{ 2,	-1,	0,	-4,	  0.360,      0.20,	-0.0043},
	{ 1,	2,	0,	0,	 -1.167,     -1.25,	-0.0106},
	{ 1,	2,	0,	-2,	 -7.412,     -6.12,	 0.0484},
	{ 1,	2,	0,	-4,	 -0.311,     -0.65,	 0.0044},
	{ 1,	-2,	0,	2,	  0.757,      1.82,	 0.0112},
	{ 1,	-2,	0,	0,	  2.580,      2.32,	 0.0196},
	{ 1,	-2,	0,	-2,	  2.533,      2.40,	-0.0212},
	{ 0,	3,	0,	-2,	 -0.344,     -0.57,	 0.0036},
	{ 1,	0,	2,	2,	 -0.992,     -0.02,	 0},
	{ 1,	0,	2,	0,	-45.099,     -0.02,	-0.0010},
	{ 1,	0,	2,	-2,	 -0.179,     -9.52,	-0.0833},
	{ 1,	0,	-2,	2,	 -6.382,     -3.37,	-0.0481},
	{ 1,	0,	-2,	0,	 39.528,     85.13,	-0.7136},
	{ 1,	0,	-2,	-2,	  9.366,      0.71,     -0.0112},
	{ 0,	1,	2,	0,	  0.415,      0.10, 	 0.0013},
	{ 0,	1,	2,	-2,	 -2.152,     -2.26,	-0.0066},
	{ 0,	1,	-2,	2,	 -1.440,     -1.30,	 0.0014},
	{ 0,	1,	-2,	-2,	  0.384,      0.0,	 0.0},
	{ 2,	0,	0,	1,	 -0.586,     -1.20,	-0.0100},
	{ 2,	0,	0,	-1,	  1.750,      2.01,	 0.0155},
	{ 2,	0,	0,	-3,	  1.225,      0.91,     -0.0088},
	{ 1,	1,	0,	1,	  1.267,      1.52,	 0.0164},
	{ 1,	-1,	0,	-1,	 -1.089,      0.55,      0},
	{ 0,	0,	2,	-1,	  0.584,      8.84,	 0.0071},
	{ 4,	0,	0,	0,	  1.938,      3.60,	 0.0401},
	{ 4,	0,	0,	-2,	 -0.952,     -1.58,	-0.0130},
	{ 3,	1,	0,	0,	 -0.551,      0.94,	-0.0097},
	{ 3,	1,	0,	-2,	 -0.482,     -0.57,	-0.0045},
	{ 3,	-1,	0,	0,	  0.681,      0.96,      0.0115},
	{ 2,	0,	2,	0,	 -3.996,      0,	 0.0004},
	{ 2,	0,	2,	-2,	  0.557,     -0.75,     -0.0090},
	{ 2,	0,	-2,	2,	 -0.459,     -0.38,	-0.0053},
	{ 2,	0,	-2,	0,	 -1.298,      0.74,      0.0004},
	{ 2,	0,	-2,	-2,	  0.538,      1.14,	-0.0141},
	{ 1,	1,	-2,	-2,	  0.426,      0.07,    	-0.0006},
	{ 1,	-1,	2,	0,	 -0.304,      0.03,	 0.0003},
	{ 1,	-1,	-2,	2,	 -0.372,     -0.19,	-0.0027},
	{ 0,	0,	4,	0,	  0.418,      0,	 0},
	{ 2,	-1,	0,	-1,	 -0.352,     -0.37,     -0.0028}
};

# if MOON_TEST_CORR
/* moon additional correction terms */
struct m5dat {		
	REAL8 lng;
	int  i0,i1,i2,i3;
	} m5 [] = {  
    /*	  lng,	l,     l',	F,	D,      */
	0.127,	0,	0,	0,	6,
	-0.151,	0,	2,	0,	-4,
	-0.085,	0,	0,	2,	4,
	0.150,	0,	1,	0,	3,
	-0.091,	2,	1,	0,	-6,
	-0.103,	0,	3,	0,	0,
	-0.301,	1,	0,	2,	-4,
	0.202,	1,	0,	-2,	-4,
	0.137,	1,	1,	0,	-1,
	0.233,	1,	1,	0,	-3,
	-0.122,	1,	-1,	0,	1,
	-0.276,	1,	-1,	0,	-3,
	0.255,	0,	0,	2,	1,
	0.254,	0,	0,	2,	-3,
	-0.100,	3,	1,	0,	-4,
	-0.183,	3,	-1,	0,	-2,
	-0.297,	2,	2,	0,	-2,
	-0.161,	2,	2,	0,	-4,
	0.197,	2,	-2,	0,	0,
	0.254,	2,	-2,	0,	-2,
	-0.250,	1,	3,	0,	-2,
	-0.123,	2,	0,	2,	2,
	0.173,	2,	0,	-2,	-4,
	0.263,	1,	1,	2,	0,
	0.130,	3,	0,	0,	-1,
	0.113,	5,	0,	0,	0,
	0.092,	3,	0,	2,	-2,
	0,	99,	0,	0,	0	/* end mark */
	};
# endif	/* MOON_TEST_CORR */
#ifdef ASTROLOG
/* End contents of helconst.c */
#endif
#ifdef ASTROLOG
/* Begin contents of deltat.c */
#endif
/*****************************************************
$Header: deltat.c,v 1.10 93/01/27 14:37:06 alois Exp $
deltat.c
deltat(t): returns delta t (in julian days) from universal time t
is included by users
ET = UT +  deltat

  ---------------------------------------------------------------
  | Copyright Astrodienst Zurich AG and Alois Treindl, 1989.	|
  | The use of this source code is subject to regulations made	|
  | by Astrodienst Zurich. The code is NOT in the public domain.|
  |								|
  | This copyright notice must not be changed or removed	|
  | by any user of this program.				|
  ---------------------------------------------------------------

******************************************************/

double deltat (double jd_ad) /* Astrodienst relative julian date */
{ 
  double floor();
  static short int dt[] = {	/* in centiseconds */
	/* dt from 1637 to 2000, as tabulated in A.E. 
	the values 1620 - 1636 are not taken, as they fit
	badly the parabola 25.5 t*t for the next range. The
	best crossing point to switch to the parabola is
	1637, where we have fitted the value for continuity */
	6780,	6500,	6300,	
	6200,	6000,	5800,	5700,	5500,
	5400,	5300,	5100,	5000,	4900,
	4800,	4700,	4600,	4500,	4400,
	4300,	4200,	4100,	4000,	3800,	/* 1655 - 59 */
	3700,	3600,	3500,	3400,	3300,
	3200,	3100,	3000,	2800,	2700,
	2600,	2500,	2400,	2300,	2200,
	2100,	2000,	1900,	1800,	1700,
	1600,	1500,	1400,	1400,	1300,
	1200,	1200,	1100,	1100,	1000,
	1000,	1000,	900,	900,	900,
	900,	900,	900,	900,	900,
	900,	900,	900,	900,	900,	/* 1700 - 1704 */
	900,	900,	900,	1000,	1000,
	1000,	1000,	1000,	1000,	1000,
	1000,	1000,	1100,	1100,	1100,
	1100,	1100,	1100,	1100,	1100,
	1100,	1100,	1100,	1100,	1100,
	1100,	1100,	1100,	1100,	1200,	/* 1730	- 1734 */
	1200,	1200,	1200,	1200,	1200,
	1200,	1200,	1200,	1200,	1300,
	1300,	1300,	1300,	1300,	1300,
	1300,	1400,	1400,	1400,	1400,
	1400,	1400,	1400,	1500,	1500,
	1500,	1500,	1500,	1500,	1500,	/* 1760 - 1764 */
	1600,	1600,	1600,	1600,	1600,
	1600,	1600,	1600,	1600,	1600,
	1700,	1700,	1700,	1700,	1700,
	1700,	1700,	1700,	1700,	1700,
	1700,	1700,	1700,	1700,	1700,
	1700,	1700,	1600,	1600,	1600,	/* 1790 - 1794 */
	1600,	1500,	1500,	1400,	1400,
	1370,	1340,	1310,	1290,	1270,	/* 1800 - 1804 */
	1260,	1250,	1250,	1250,	1250,
	1250,	1250,	1250,	1250,	1250,
	1250,	1250,	1240,	1230,	1220,
	1200,	1170,	1140,	1110,	1060,
	1020,	960,	910,	860,	800,
	750,	700,	660,	630,	600,	/* 1830 - 1834 */
	580,	570,	560,	560,	560,
	570,	580,	590,	610,	620,
	630,	650,	660,	680,	690,
	710,	720,	730,	740,	750,
	760,	770,	770,	780,	780,
	788,	782,	754,	697,	640,	/* 1860 - 1864 */
	602,	541,	410,	292,	182,
	161,	10,	-102,	-128,	-269,
	-324,	-364,	-454,	-471,	-511,
	-540,	-542,	-520,	-546,	-546,
	-579,	-563,	-564,	-580,	-566,
	-587,	-601,	-619,	-664,	-644,	/* 1890 - 1894 */
	-647,	-609,	-576,	-466,	-374,
	-272,	-154,	-2,	124,	264,
	386,	537,	614,	775,	913,
	1046,	1153,	1336,	1465,	1601,
	1720,	1824,	1906,	2025,	2095,
	2116,	2225,	2241,	2303,	2349,	/* 1920 - 1924 */