| 647 |
|
const int mo = erec.date.month+1; |
| 648 |
|
const int da = erec.date.day; |
| 649 |
|
const double hr = erec.date.hour; |
| 650 |
< |
double aod = erec.optdepth; |
| 650 |
> |
double aod = erec.optdepth * 1e3; |
| 651 |
> |
if (aod >= 999.0) { |
| 652 |
> |
fprintf(stderr, "aod not set, using default value %.3f\n", AOD0_CA); |
| 653 |
> |
aod = AOD0_CA; |
| 654 |
> |
} |
| 655 |
|
double cc = erec.skycover; |
| 656 |
+ |
if (cc >= 99.0) { |
| 657 |
+ |
fprintf(stderr, "skycover not set, using default value 0.0\n"); |
| 658 |
+ |
cc = 0.0; |
| 659 |
+ |
} |
| 660 |
|
double sda, sta, st; |
| 661 |
|
int sun_in_sky; |
| 662 |
|
|
| 655 |
– |
if (aod == 0.0) { |
| 656 |
– |
aod = AOD0_CA; |
| 657 |
– |
fprintf(stderr, "aod is zero, using default value %.3f\n", AOD0_CA); |
| 658 |
– |
} |
| 663 |
|
/* compute solar position */ |
| 664 |
|
if ((mo == 2) & (da == 29)) { |
| 665 |
|
julian_date = 60; |