| 36 |
|
extern double s_meridian; |
| 37 |
|
/* required values */ |
| 38 |
|
int month, day; /* date */ |
| 39 |
< |
double hour; /* standard time */ |
| 39 |
> |
double hour; /* time */ |
| 40 |
> |
int tsolar; /* 0=standard, 1=solar */ |
| 41 |
|
double altitude, azimuth; /* or solar angles */ |
| 42 |
|
/* default values */ |
| 43 |
|
int cloudy = 0; /* 1=standard, 2=uniform */ |
| 77 |
|
month = atoi(argv[1]); |
| 78 |
|
day = atoi(argv[2]); |
| 79 |
|
hour = atof(argv[3]); |
| 80 |
+ |
tsolar = argv[3][0] == '+'; |
| 81 |
|
} |
| 82 |
|
for (i = 4; i < argc; i++) |
| 83 |
|
if (argv[i][0] == '-' || argv[i][0] == '+') |
| 136 |
|
|
| 137 |
|
jd = jdate(month, day); /* Julian date */ |
| 138 |
|
sd = sdec(jd); /* solar declination */ |
| 139 |
< |
st = hour + stadj(jd); /* solar time */ |
| 139 |
> |
if (tsolar) /* solar time */ |
| 140 |
> |
st = hour; |
| 141 |
> |
else |
| 142 |
> |
st = hour + stadj(jd); |
| 143 |
|
altitude = salt(sd, st); |
| 144 |
|
azimuth = sazi(sd, st); |
| 145 |
|
} |
