-
Notifications
You must be signed in to change notification settings - Fork 0
/
test_run.go
130 lines (118 loc) · 3.33 KB
/
test_run.go
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
package main
import (
"fmt"
"github.com/maltegrosse/go-spa"
"os"
"text/tabwriter"
"time"
)
func main() {
deltaUt1 := 0.
deltaT := 67.
longitude := -105.1786
latitude := 39.742476
elevation := 1830.14
pressure := 820.
temperature := 11.
slope := 30.
azmRotation := -10.
atmosRefract := 0.5667
loc, err := time.LoadLocation("America/Los_Angeles")
if err != nil {
fmt.Println(err)
return
}
dt := time.Date(2003, 10, 17, 12, 30, 30, 0, loc)
spa, err := spa.NewSpa(dt, latitude, longitude, elevation, pressure, temperature, deltaT, deltaUt1, slope, azmRotation, atmosRefract)
if err != nil {
fmt.Println(err)
return
}
writer := tabwriter.NewWriter(os.Stdout, 0, 8, 1, '\t', tabwriter.AlignRight)
_, err = fmt.Fprintln(writer, "- \tNREL spa_tester.c\tGO SPA")
if err != nil {
fmt.Println(err)
return
}
_, err = fmt.Fprintln(writer, "Julian Day", "\t", "2452930.312847", "\t", fmt.Sprintf("%.6f", spa.GetJd()), "\t")
if err != nil {
fmt.Println(err)
return
}
_, err = fmt.Fprintln(writer, "L", "\t", "24.01826", "\t", fmt.Sprintf("%.6f", spa.GetL()), "\t")
if err != nil {
fmt.Println(err)
return
}
_, err = fmt.Fprintln(writer, "B", "\t", "-0.0001011219", "\t", fmt.Sprintf("%.12f", spa.GetB()), "\t")
if err != nil {
fmt.Println(err)
return
}
_, err = fmt.Fprintln(writer, "R", "\t", "0.996542", "\t", fmt.Sprintf("%.6f", spa.GetR()), "\t")
if err != nil {
fmt.Println(err)
return
}
_, err = fmt.Fprintln(writer, "H", "\t", "11.105902", "\t", fmt.Sprintf("%.6f", spa.GetH()), "\t")
if err != nil {
fmt.Println(err)
return
}
_, err = fmt.Fprintln(writer, "Delta Psi", "\t", "-0.003998404", "\t", fmt.Sprintf("%.12f", spa.GetDelPsi()), "\t")
if err != nil {
fmt.Println(err)
return
}
_, err = fmt.Fprintln(writer, "Delta Epsilon", "\t", "0.001666568", "\t", fmt.Sprintf("%.12f", spa.GetDelEpsilon()), "\t")
if err != nil {
fmt.Println(err)
return
}
_, err = fmt.Fprintln(writer, "Epsilon", "\t", "23.440465", "\t", fmt.Sprintf("%.6f", spa.GetEpsilon()), "\t")
if err != nil {
fmt.Println(err)
return
}
_, err = fmt.Fprintln(writer, "Zenith", "\t", "50.111622", "\t", fmt.Sprintf("%.6f", spa.GetZenith()), "\t")
if err != nil {
fmt.Println(err)
return
}
_, err = fmt.Fprintln(writer, "Azimuth", "\t", "194.340241", "\t", fmt.Sprintf("%.6f", spa.GetAzimuth()), "\t")
if err != nil {
fmt.Println(err)
return
}
_, err = fmt.Fprintln(writer, "Incidence", "\t", "25.187000", "\t", fmt.Sprintf("%.6f", spa.GetIncidence()), "\t")
if err != nil {
fmt.Println(err)
return
}
err = writer.Flush()
if err != nil {
fmt.Println(err)
return
}
// use spa.Calculate() after the initial input values changed...
fmt.Println(spa.GetSunrise())
fmt.Println(spa.GetSunset())
/////////////////////////////////////////////
// The output of this program should be:
//
//Julian Day: 2452930.312847
//L: 2.401826e+01 degrees
//B: -1.011219e-04 degrees
//R: 0.996542 AU
//H: 11.105902 degrees
//Delta Psi: -3.998404e-03 degrees
//Delta Epsilon: 1.666568e-03 degrees
//Epsilon: 23.440465 degrees
//Zenith: 50.111622 degrees
//Azimuth: 194.340241 degrees
//Incidence: 25.187000 degrees
//Sunrise: 06:12:43 Local Time
//Sunset: 17:20:19 Local Time
//
/////////////////////////////////////////////
}