49 |
|
double sundir[3]; |
50 |
|
double groundbr; |
51 |
|
double F2; |
52 |
< |
double solarbr; |
52 |
> |
double solarbr = -1.0; |
53 |
|
|
54 |
|
char *progname; |
55 |
|
char errmsg[128]; |
91 |
|
cloudy = 0; |
92 |
|
dosun = argv[i][0] == '+'; |
93 |
|
break; |
94 |
+ |
case 'r': |
95 |
+ |
solarbr = atof(argv[++i]); |
96 |
+ |
break; |
97 |
|
case 'c': |
98 |
|
cloudy = argv[i][0] == '+' ? 2 : 1; |
99 |
|
dosun = 0; |
178 |
|
0.45*sundir[2]*sundir[2]); |
179 |
|
groundbr = zenithbr*normsc(altitude)/F2/PI; |
180 |
|
printf("# Ground ambient level: %f\n", groundbr); |
181 |
< |
if (sundir[2] > 0.0) { |
182 |
< |
if (sundir[2] > .16) |
183 |
< |
solarbr = (1.5e9/SUNEFFICACY) * |
184 |
< |
(1.147 - .147/sundir[2]); |
182 |
< |
else |
183 |
< |
solarbr = 1.5e9/SUNEFFICACY*(1.147-.147/.16); |
181 |
> |
if (sundir[2] > 0.0 && solarbr != 0.0) { |
182 |
> |
if (solarbr < 0.0) |
183 |
> |
solarbr = 1.5e9/SUNEFFICACY * |
184 |
> |
(1.147 - .147/(sundir[2]>.16?sundir[2]:.16)); |
185 |
|
groundbr += solarbr*6e-5*sundir[2]/PI; |
186 |
|
} else |
187 |
|
dosun = 0; |