1 |
|
2 |
****************************************************** |
3 |
****** IGRF GEOMAGNETIC FIELD MODEL 1945 - 2000 ****** |
4 |
*********** SHELLG L-VALUE CALCULATION ************* |
5 |
************************************************************ |
6 |
This program allows you to produce B and L profiles in |
7 |
latitude, longitude, year or altitude. |
8 |
In each of the following windows you will be asked to enter |
9 |
one or more values, defining the conditions for your tables. |
10 |
In each window the current value(s) is (are) shown in the right |
11 |
upper corner (#...#). You can choose the current values by |
12 |
entering / at the prompt. |
13 |
If you enter a wrong character or a value outside the allowed |
14 |
parameter range, the program will ask you for a new entry. |
15 |
After your tables are displayed, you can change any parameter |
16 |
you wish to change and create a new profile. |
17 |
You can leave the program at any point by entering Ctrl Z. |
18 |
************************* GOOD LUCK ************************* |
19 |
|
20 |
DO YOU WANT YOUR PROFILES #2# |
21 |
DISPLAYED ON YOUR MONITOR: ENTER 0 AT PROMPT |
22 |
STORED IN FILE OUTPUT.IGR: ENTER 1 AT PROMPT |
23 |
DISPLAYED AND STORED: ENTER 2 AT PROMPT |
24 |
------------------------------------------------------------ |
25 |
Enter / to use previous value(s) (see # .. #); Ctrl Z to exit |
26 |
/ |
27 |
|
28 |
|
29 |
SELECT YOUR VARIABLE: #LAST:3# |
30 |
|
31 |
1 LATITUDE 3 ALTITUDE |
32 |
2 LONGITUDE 4 YEAR |
33 |
------------------------------------------------------------ |
34 |
Enter / to use previous value(s) (see # .. #); Ctrl Z to exit |
35 |
/ |
36 |
|
37 |
|
38 |
CHOOSE YOUR VARIABLE RANGE: BEGIN, END, STEPWIDTH ? |
39 |
# 100.0, 1000.0, 100.0# |
40 |
------------------------------------------------------------ |
41 |
Enter / to use previous value(s) (see # .. #); Ctrl Z to exit |
42 |
/ |
43 |
|
44 |
|
45 |
GEOD LATITUDE ? !NORTH! [DEGREE,DECIMAL] # 45.1# |
46 |
------------------------------------------------------------ |
47 |
Enter / to use previous value(s) (see # .. #); Ctrl Z to exit |
48 |
/ |
49 |
|
50 |
|
51 |
GEOD LONGITUDE ? !EAST! [DEGREE,DECIMAL] # 293.1# |
52 |
------------------------------------------------------------ |
53 |
Enter / to use previous value(s) (see # .. #); Ctrl Z to exit |
54 |
/ |
55 |
|
56 |
|
57 |
YEAR(EPOCH) ? *decimal* #1985.5# |
58 |
------------------------------------------------------------ |
59 |
Enter / to use previous value(s) (see # .. #); Ctrl Z to exit |
60 |
/ |
61 |
|
62 |
|
63 |
OUTPUT OPTION: B OR B/B0 ? #0# |
64 |
|
65 |
if you enter 0, the absolute magnetic field strength |
66 |
will be listed, otherwise the field strength normalized |
67 |
to the field strength at the magnetic equator is listed |
68 |
------------------------------------------------------------ |
69 |
Enter / to use previous value(s) (see # .. #); Ctrl Z to exit |
70 |
/ |
71 |
|
72 |
|
73 |
|
74 |
|
75 |
|
76 |
|
77 |
|
78 |
|
79 |
|
80 |
|
81 |
|
82 |
|
83 |
H/km DIMO B/Gauss B-NORTH B-EAST B-DOWN DIP DEC L-VALUE C |
84 |
100.00 0.3042 0.52942 0.16107 -.05615 0.50118 71.2 -19.2 3.216 2 |
85 |
200.00 0.3042 0.50404 0.15406 -.05222 0.47707 71.2 -18.7 3.264 2 |
86 |
300.00 0.3042 0.48030 0.14745 -.04865 0.45451 71.1 -18.3 3.313 1 |
87 |
400.00 0.3042 0.45808 0.14122 -.04539 0.43340 71.1 -17.8 3.361 1 |
88 |
500.00 0.3042 0.43724 0.13534 -.04241 0.41360 71.1 -17.4 3.409 1 |
89 |
600.00 0.3042 0.41769 0.12978 -.03968 0.39503 71.0 -17.0 3.457 1 |
90 |
700.00 0.3042 0.39931 0.12453 -.03717 0.37757 71.0 -16.6 3.505 1 |
91 |
800.00 0.3042 0.38202 0.11955 -.03487 0.36115 71.0 -16.3 3.553 1 |
92 |
900.00 0.3042 0.36574 0.11483 -.03275 0.34570 70.9 -15.9 3.601 1 |
93 |
1000.00 0.3042 0.35039 0.11036 -.03080 0.33113 70.9 -15.6 3.649 1 |
94 |
------- International Geomagnetic Reference Field --- Edition 1995 ------- |
95 |
LATI= 45.1 LONGI= 293.1 I DIMO is Dipol I C=1 L and B0 correct |
96 |
ALT= 1000.0 YEAR=1985.5 I Moment in Gauss I =2 wrong, =3 approx. |
97 |
-------------------------------------------------------------------------- |
98 |
|
99 |
**** DO YOU WANT TO CONTINUE? |
100 |
------------------------------------------------------------ |
101 |
"0" QUIT AND EXIT "1" NEW PARAMETERS |
102 |
------------------------------------------------------------ |
103 |
1 |
104 |
|
105 |
|
106 |
**** WHICH PARAMETER DO YOU WANT TO CHANGE? |
107 |
------------------------------------------------------------ |
108 |
0 NO FURTHER CHANGES, CALCULATE PROFILE |
109 |
1 LATITUDE # 45.1# 5 DISPLAY OR STORE |
110 |
2 LONGITUDE # 293.1# 6 SELECTION OF VARIABLE #H/km# |
111 |
3 ALTITUDE # 1000.0# 7 VARIABLE RANGE |
112 |
4 YEAR #1985.5# # 100.0, 1000.0, 100.0# |
113 |
8 B OR B/B0 |
114 |
------------------------------------------------------------ |
115 |
ENTER NUMBER |
116 |
7 |
117 |
|
118 |
|
119 |
CHOOSE YOUR VARIABLE RANGE: BEGIN, END, STEPWIDTH ? |
120 |
# 100.0, 1000.0, 100.0# |
121 |
------------------------------------------------------------ |
122 |
Enter / to use previous value(s) (see # .. #); Ctrl Z to exit |
123 |
1000,10000,1000 |
124 |
|
125 |
|
126 |
**** WHICH PARAMETER DO YOU WANT TO CHANGE? |
127 |
------------------------------------------------------------ |
128 |
0 NO FURTHER CHANGES, CALCULATE PROFILE |
129 |
1 LATITUDE # 45.1# 5 DISPLAY OR STORE |
130 |
2 LONGITUDE # 293.1# 6 SELECTION OF VARIABLE #H/km# |
131 |
3 ALTITUDE # 1000.0# 7 VARIABLE RANGE |
132 |
4 YEAR #1985.5# # 1000.0, 10000.0, 1000.0# |
133 |
8 B OR B/B0 |
134 |
------------------------------------------------------------ |
135 |
ENTER NUMBER |
136 |
0 |
137 |
|
138 |
|
139 |
|
140 |
|
141 |
|
142 |
|
143 |
|
144 |
|
145 |
|
146 |
|
147 |
|
148 |
|
149 |
H/km DIMO B/Gauss B-NORTH B-EAST B-DOWN DIP DEC L-VALUE C |
150 |
1000.00 0.3042 0.35039 0.11036 -.03080 0.33113 70.9 -15.6 3.649 1 |
151 |
2000.00 0.3042 0.23588 0.07614 -.01761 0.22256 70.7 -13.0 4.134 1 |
152 |
3000.00 0.3042 0.16675 0.05465 -.01091 0.15716 70.5 -11.3 4.622 1 |
153 |
4000.00 0.3042 0.12239 0.04049 -.00718 0.11528 70.4 -10.1 5.107 1 |
154 |
5000.00 0.3042 0.09256 0.03080 -.00494 0.08715 70.3 -9.1 5.596 1 |
155 |
6000.00 0.3042 0.07174 0.02396 -.00353 0.06753 70.3 -8.4 6.092 1 |
156 |
7000.00 0.3042 0.05675 0.01899 -.00260 0.05341 70.3 -7.8 6.593 1 |
157 |
8000.00 0.3042 0.04567 0.01530 -.00196 0.04299 70.3 -7.3 7.093 1 |
158 |
9000.00 0.3042 0.03731 0.01250 -.00151 0.03512 70.3 -6.9 7.590 1 |
159 |
10000.00 0.3042 0.03088 0.01034 -.00118 0.02907 70.3 -6.5 8.084 1 |
160 |
------- International Geomagnetic Reference Field --- Edition 1995 ------- |
161 |
LATI= 45.1 LONGI= 293.1 I DIMO is Dipol I C=1 L and B0 correct |
162 |
ALT=10000.0 YEAR=1985.5 I Moment in Gauss I =2 wrong, =3 approx. |
163 |
-------------------------------------------------------------------------- |
164 |
!! REMINDER: this field model does not include external sources !! |
165 |
|
166 |
**** DO YOU WANT TO CONTINUE? |
167 |
------------------------------------------------------------ |
168 |
"0" QUIT AND EXIT "1" NEW PARAMETERS |
169 |
------------------------------------------------------------ |
170 |
1 |
171 |
|
172 |
|
173 |
**** WHICH PARAMETER DO YOU WANT TO CHANGE? |
174 |
------------------------------------------------------------ |
175 |
0 NO FURTHER CHANGES, CALCULATE PROFILE |
176 |
1 LATITUDE # 45.1# 5 DISPLAY OR STORE |
177 |
2 LONGITUDE # 293.1# 6 SELECTION OF VARIABLE #H/km# |
178 |
3 ALTITUDE # 10000.0# 7 VARIABLE RANGE |
179 |
4 YEAR #1985.5# # 1000.0, 10000.0, 1000.0# |
180 |
8 B OR B/B0 |
181 |
------------------------------------------------------------ |
182 |
ENTER NUMBER |
183 |
3 |
184 |
|
185 |
|
186 |
ALTITUDE ? [KM] #10000.0# |
187 |
------------------------------------------------------------ |
188 |
Enter / to use previous value(s) (see # .. #); Ctrl Z to exit |
189 |
1000 |
190 |
|
191 |
|
192 |
**** WHICH PARAMETER DO YOU WANT TO CHANGE? |
193 |
------------------------------------------------------------ |
194 |
0 NO FURTHER CHANGES, CALCULATE PROFILE |
195 |
1 LATITUDE # 45.1# 5 DISPLAY OR STORE |
196 |
2 LONGITUDE # 293.1# 6 SELECTION OF VARIABLE #H/km# |
197 |
3 ALTITUDE # 1000.0# 7 VARIABLE RANGE |
198 |
4 YEAR #1985.5# # 1000.0, 10000.0, 1000.0# |
199 |
8 B OR B/B0 |
200 |
------------------------------------------------------------ |
201 |
ENTER NUMBER |
202 |
1 |
203 |
|
204 |
|
205 |
GEOD LATITUDE ? !NORTH! [DEGREE,DECIMAL] # 45.1# |
206 |
------------------------------------------------------------ |
207 |
Enter / to use previous value(s) (see # .. #); Ctrl Z to exit |
208 |
30 |
209 |
|
210 |
|
211 |
**** WHICH PARAMETER DO YOU WANT TO CHANGE? |
212 |
------------------------------------------------------------ |
213 |
0 NO FURTHER CHANGES, CALCULATE PROFILE |
214 |
1 LATITUDE # 30.0# 5 DISPLAY OR STORE |
215 |
2 LONGITUDE # 293.1# 6 SELECTION OF VARIABLE #H/km# |
216 |
3 ALTITUDE # 1000.0# 7 VARIABLE RANGE |
217 |
4 YEAR #1985.5# # 1000.0, 10000.0, 1000.0# |
218 |
8 B OR B/B0 |
219 |
------------------------------------------------------------ |
220 |
ENTER NUMBER |
221 |
6 |
222 |
|
223 |
|
224 |
SELECT YOUR VARIABLE: #LAST:3# |
225 |
|
226 |
1 LATITUDE 3 ALTITUDE |
227 |
2 LONGITUDE 4 YEAR |
228 |
------------------------------------------------------------ |
229 |
Enter / to use previous value(s) (see # .. #); Ctrl Z to exit |
230 |
4 |
231 |
|
232 |
|
233 |
CHOOSE YOUR VARIABLE RANGE: BEGIN, END, STEPWIDTH ? |
234 |
# 1000.0, 10000.0, 1000.0# |
235 |
------------------------------------------------------------ |
236 |
Enter / to use previous value(s) (see # .. #); Ctrl Z to exit |
237 |
1945,2000,5 |
238 |
|
239 |
|
240 |
**** WHICH PARAMETER DO YOU WANT TO CHANGE? |
241 |
------------------------------------------------------------ |
242 |
0 NO FURTHER CHANGES, CALCULATE PROFILE |
243 |
1 LATITUDE # 30.0# 5 DISPLAY OR STORE |
244 |
2 LONGITUDE # 293.1# 6 SELECTION OF VARIABLE #YEAR# |
245 |
3 ALTITUDE # 1000.0# 7 VARIABLE RANGE |
246 |
4 YEAR #1985.5# # 1945.0, 2000.0, 5.0# |
247 |
8 B OR B/B0 |
248 |
------------------------------------------------------------ |
249 |
ENTER NUMBER |
250 |
8 |
251 |
|
252 |
|
253 |
OUTPUT OPTION: B OR B/B0 ? #0# |
254 |
|
255 |
if you enter 0, the absolute magnetic field strength |
256 |
will be listed, otherwise the field strength normalized |
257 |
to the field strength at the magnetic equator is listed |
258 |
------------------------------------------------------------ |
259 |
Enter / to use previous value(s) (see # .. #); Ctrl Z to exit |
260 |
1 |
261 |
|
262 |
|
263 |
**** WHICH PARAMETER DO YOU WANT TO CHANGE? |
264 |
------------------------------------------------------------ |
265 |
0 NO FURTHER CHANGES, CALCULATE PROFILE |
266 |
1 LATITUDE # 30.0# 5 DISPLAY OR STORE |
267 |
2 LONGITUDE # 293.1# 6 SELECTION OF VARIABLE #YEAR# |
268 |
3 ALTITUDE # 1000.0# 7 VARIABLE RANGE |
269 |
4 YEAR #1985.5# # 1945.0, 2000.0, 5.0# |
270 |
8 B OR B/B0 |
271 |
------------------------------------------------------------ |
272 |
ENTER NUMBER |
273 |
0 |
274 |
|
275 |
|
276 |
|
277 |
|
278 |
|
279 |
|
280 |
|
281 |
|
282 |
|
283 |
|
284 |
|
285 |
|
286 |
YEAR DIMO B/B0 B-NORTH B-EAST B-DOWN DIP DEC L-VALUE C |
287 |
1945.00 0.3122 11.336 0.14421 -.02172 0.28536 62.9 -8.6 2.220 1 |
288 |
1950.00 0.3118 11.153 0.14491 -.02244 0.28389 62.7 -8.8 2.209 1 |
289 |
1955.00 0.3113 10.966 0.14586 -.02342 0.28222 62.4 -9.1 2.199 1 |
290 |
1960.00 0.3104 10.696 0.14713 -.02436 0.27977 61.9 -9.4 2.183 1 |
291 |
1965.00 0.3095 10.430 0.14811 -.02543 0.27705 61.5 -9.7 2.167 1 |
292 |
1970.00 0.3083 10.109 0.14947 -.02637 0.27367 61.0 -10.0 2.148 1 |
293 |
1975.00 0.3070 9.753 0.15071 -.02735 0.26975 60.4 -10.3 2.127 1 |
294 |
1980.00 0.3057 9.407 0.15096 -.02849 0.26521 59.9 -10.7 2.108 1 |
295 |
1985.00 0.3043 9.059 0.15128 -.02957 0.26104 59.4 -11.1 2.087 1 |
296 |
1990.00 0.3032 8.736 0.15164 -.03045 0.25712 59.0 -11.4 2.068 1 |
297 |
1995.00 0.3021 8.438 0.15201 -.03129 0.25348 58.5 -11.6 2.049 1 |
298 |
2000.00 0.3010 8.136 0.15245 -.03213 0.24982 58.1 -11.9 2.029 1 |
299 |
------- International Geomagnetic Reference Field --- Edition 1995 ------- |
300 |
LATI= 30.0 LONGI= 293.1 I DIMO is Dipol I C=1 L and B0 correct |
301 |
ALT= 1000.0 YEAR=2000.0 I Moment in Gauss I =2 wrong, =3 approx. |
302 |
-------------------------------------------------------------------------- |
303 |
|
304 |
**** DO YOU WANT TO CONTINUE? |
305 |
------------------------------------------------------------ |
306 |
"0" QUIT AND EXIT "1" NEW PARAMETERS |
307 |
------------------------------------------------------------ |
308 |
1 |
309 |
|
310 |
|
311 |
**** WHICH PARAMETER DO YOU WANT TO CHANGE? |
312 |
------------------------------------------------------------ |
313 |
0 NO FURTHER CHANGES, CALCULATE PROFILE |
314 |
1 LATITUDE # 30.0# 5 DISPLAY OR STORE |
315 |
2 LONGITUDE # 293.1# 6 SELECTION OF VARIABLE #YEAR# |
316 |
3 ALTITUDE # 1000.0# 7 VARIABLE RANGE |
317 |
4 YEAR #2000.0# # 1945.0, 2000.0, 5.0# |
318 |
8 B OR B/B0 |
319 |
------------------------------------------------------------ |
320 |
ENTER NUMBER |
321 |
4 |
322 |
|
323 |
|
324 |
YEAR(EPOCH) ? *decimal* #2000.0# |
325 |
------------------------------------------------------------ |
326 |
Enter / to use previous value(s) (see # .. #); Ctrl Z to exit |
327 |
1995 |
328 |
|
329 |
|
330 |
**** WHICH PARAMETER DO YOU WANT TO CHANGE? |
331 |
------------------------------------------------------------ |
332 |
0 NO FURTHER CHANGES, CALCULATE PROFILE |
333 |
1 LATITUDE # 30.0# 5 DISPLAY OR STORE |
334 |
2 LONGITUDE # 293.1# 6 SELECTION OF VARIABLE #YEAR# |
335 |
3 ALTITUDE # 1000.0# 7 VARIABLE RANGE |
336 |
4 YEAR #1995.0# # 1945.0, 2000.0, 5.0# |
337 |
8 B OR B/B0 |
338 |
------------------------------------------------------------ |
339 |
ENTER NUMBER |
340 |
6 |
341 |
|
342 |
|
343 |
SELECT YOUR VARIABLE: #LAST:4# |
344 |
|
345 |
1 LATITUDE 3 ALTITUDE |
346 |
2 LONGITUDE 4 YEAR |
347 |
------------------------------------------------------------ |
348 |
Enter / to use previous value(s) (see # .. #); Ctrl Z to exit |
349 |
1 |
350 |
|
351 |
|
352 |
CHOOSE YOUR VARIABLE RANGE: BEGIN, END, STEPWIDTH ? |
353 |
# 1945.0, 2000.0, 5.0# |
354 |
------------------------------------------------------------ |
355 |
Enter / to use previous value(s) (see # .. #); Ctrl Z to exit |
356 |
-80,80,20 |
357 |
|
358 |
|
359 |
**** WHICH PARAMETER DO YOU WANT TO CHANGE? |
360 |
------------------------------------------------------------ |
361 |
0 NO FURTHER CHANGES, CALCULATE PROFILE |
362 |
1 LATITUDE # 30.0# 5 DISPLAY OR STORE |
363 |
2 LONGITUDE # 293.1# 6 SELECTION OF VARIABLE #LATI# |
364 |
3 ALTITUDE # 1000.0# 7 VARIABLE RANGE |
365 |
4 YEAR #1995.0# # -80.0, 80.0, 20.0# |
366 |
8 B OR B/B0 |
367 |
------------------------------------------------------------ |
368 |
ENTER NUMBER |
369 |
0 |
370 |
|
371 |
|
372 |
|
373 |
|
374 |
|
375 |
|
376 |
|
377 |
|
378 |
|
379 |
|
380 |
|
381 |
|
382 |
LATI DIMO B/B0 B-NORTH B-EAST B-DOWN DIP DEC L-VALUE C |
383 |
-80.00 0.3021 310.936 0.10062 0.04627 -.30864 -70.3 24.7 6.596 1 |
384 |
-60.00 0.3021 11.854 0.13016 0.02925 -.21172 -57.8 12.7 2.425 1 |
385 |
-40.00 0.3021 2.105 0.13885 0.00865 -.12325 -41.5 3.6 1.509 1 |
386 |
-20.00 0.3021 1.022 0.15771 -.00938 -.03933 -14.0 -3.4 1.241 1 |
387 |
0.00 0.3021 1.300 0.17571 -.02240 0.07077 21.8 -7.3 1.274 1 |
388 |
20.00 0.3021 3.697 0.16803 -.02961 0.19495 48.8 -10.0 1.633 1 |
389 |
40.00 0.3021 24.795 0.12767 -.03161 0.30420 66.6 -13.9 2.843 1 |
390 |
60.00 0.3021 966.671 0.06388 -.02775 0.36494 79.2 -23.5 9.204 1 |
391 |
80.00 0.30219999.999 0.01140 -.01756 0.37588 86.8 -57.0 509.552 3 |
392 |
------- International Geomagnetic Reference Field --- Edition 1995 ------- |
393 |
LATI= 80.0 LONGI= 293.1 I DIMO is Dipol I C=1 L and B0 correct |
394 |
ALT= 1000.0 YEAR=1995.0 I Moment in Gauss I =2 wrong, =3 approx. |
395 |
-------------------------------------------------------------------------- |
396 |
|
397 |
**** DO YOU WANT TO CONTINUE? |
398 |
------------------------------------------------------------ |
399 |
"0" QUIT AND EXIT "1" NEW PARAMETERS |
400 |
------------------------------------------------------------ |
401 |
0 |
402 |
FORTRAN STOP |
403 |
|