7 |
|
* |
8 |
|
*} |
9 |
|
|
10 |
+ |
DEG : PI/180; |
11 |
+ |
|
12 |
|
RLAT = 0.66; {* Site latitude (radians) *} |
13 |
|
|
14 |
|
RLON = 2.13; {* Site longitude (radians) *} |
15 |
|
|
16 |
< |
RSM = 120 * PI/180; {* Standard meridian (radians) *} |
16 |
> |
RSM = 120 * DEG; {* Standard meridian (radians) *} |
17 |
|
|
18 |
|
|
19 |
|
stadj(jd) {* solar time adjustment (hours) *} |
36 |
|
{ sd is solar declination, st is solar time (hours) } |
37 |
|
= |
38 |
|
asin( sin(RLAT) * sin(sd) - cos(RLAT) * cos(sd) * |
39 |
< |
cos(st*PI/12) ) |
39 |
> |
cos(PI/12*st) ) |
40 |
|
; |
41 |
|
|
42 |
|
|
43 |
|
sazi(sd, st) {* the solar azimuth (radians) *} |
44 |
|
{ sd is solar declination, st is solar time (hours) } |
45 |
|
= |
46 |
< |
-atan2( cos(sd) * sin(st*PI/12), |
46 |
> |
-atan2( cos(sd) * sin(PI/12*st), |
47 |
|
-cos(RLAT)*sin(sd) - |
48 |
< |
sin(RLAT)*cos(sd)*cos(st*PI/12) ) |
48 |
> |
sin(RLAT)*cos(sd)*cos(PI/12*st) ) |
49 |
|
; |
50 |
|
|
51 |
|
|
52 |
< |
SAZI = sazi(SDEC, STIME) * 180/PI; |
52 |
> |
SAZI = sazi(SDEC, STIME) / DEG; |
53 |
|
|
54 |
< |
SALT = salt(SDEC, STIME) * 180/PI; |
54 |
> |
SALT = salt(SDEC, STIME) / DEG; |
55 |
|
|
56 |
|
SDEC = sdec(JDATE); |
57 |
|
|
58 |
|
STIME = TIME + stadj(JDATE); |
59 |
+ |
|
60 |
+ |
HALFDAY = acos(-tan(RLAT)*tan(SDEC)) * (12/PI); |
61 |
|
|
62 |
|
JDATE = DAY + select(MONTH,0,31,59,90,120,151,181,212,243,273,304,334); |