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; |