218 |
|
if (li == last_li) /* cached latitude? */ |
219 |
|
return last_ohm; |
220 |
|
last_li = li; |
221 |
– |
theta1 = M_PI/180. * ab->lat[li+1].tmin; |
222 |
– |
if (ab->lat[li].nphis == 1) /* special case */ |
223 |
– |
return last_ohm = M_PI*(1. - sq(cos(theta1))); |
221 |
|
theta = M_PI/180. * ab->lat[li].tmin; |
222 |
+ |
theta1 = M_PI/180. * ab->lat[li+1].tmin; |
223 |
|
return last_ohm = M_PI*(sq(cos(theta)) - sq(cos(theta1))) / |
224 |
|
(double)ab->lat[li].nphis; |
225 |
|
} |
333 |
|
ezxml_child(ezxml_child(wbb, |
334 |
|
"ThetaBounds"), "UpperTheta"))); |
335 |
|
if (!i) |
336 |
< |
abase_list[nabases].lat[i].tmin = |
339 |
< |
-abase_list[nabases].lat[i+1].tmin; |
336 |
> |
abase_list[nabases].lat[0].tmin = 0; |
337 |
|
else if (!fequal(atof(ezxml_txt(ezxml_child(ezxml_child(wbb, |
338 |
|
"ThetaBounds"), "LowerTheta"))), |
339 |
|
abase_list[nabases].lat[i].tmin)) { |