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)) { |
756 |
|
} |
757 |
|
cdlast = NULL; /* check for it in cache list */ |
758 |
|
for (cd = (SDMatCDst *)sdc->cdList; cd != NULL; |
759 |
< |
cdlast = cd, cd = (SDMatCDst *)cd->next) |
759 |
> |
cdlast = cd, cd = cd->next) |
760 |
|
if (cd->indx == myCD.indx && (cd->calen == myCD.calen) & |
761 |
|
(cd->ob_priv == myCD.ob_priv) & |
762 |
|
(cd->ob_vec == myCD.ob_vec)) |
775 |
|
} |
776 |
|
if (cdlast != NULL) { /* move entry to head of cache list */ |
777 |
|
cdlast->next = cd->next; |
778 |
< |
cd->next = sdc->cdList; |
778 |
> |
cd->next = (SDMatCDst *)sdc->cdList; |
779 |
|
sdc->cdList = (SDCDst *)cd; |
780 |
|
} |
781 |
|
return (SDCDst *)cd; /* ready to go */ |