FACT++  1.0
double palDt ( double  epoch)

Definition at line 97 of file palDt.c.

References t.

Referenced by palPlantu(), palRdplan(), and t_dat().

97  {
98 
99  double t,w,s;
100 
101  /* Centuries since 1800 */
102  t = (epoch - 1800.0) / 100.0;
103 
104  /* Select model */
105  if ( epoch >= 1708.185161980887 ) {
106 
107  /* Post-1708: use McCarthy & Babcock */
108  w = t - 0.19;
109  s = 5.156 + 13.3066 * w * w;
110 
111  } else if ( epoch >= 979.0258204760233 ) {
112 
113  /* 978-1708: use Stephenson & Morrison's 948-1600 model */
114  s = 25.5 * t * t;
115 
116  } else {
117 
118  /* Pre-979: use Stephenson & Morrison's 390 BC to AD 948 model */
119  s = 1360.0 + (320.0 + 44.3*t) * t;
120 
121  }
122 
123  return s;
124 
125 }
double w
Definition: palObs.c:168
TT t
Definition: test_client.c:26

+ Here is the caller graph for this function: