****************************************************** ****** IGRF GEOMAGNETIC FIELD MODEL 1945 - 2000 ****** *********** SHELLG L-VALUE CALCULATION ************* ************************************************************ This program allows you to produce B and L profiles in latitude, longitude, year or altitude. In each of the following windows you will be asked to enter one or more values, defining the conditions for your tables. In each window the current value(s) is (are) shown in the right upper corner (#...#). You can choose the current values by entering / at the prompt. If you enter a wrong character or a value outside the allowed parameter range, the program will ask you for a new entry. After your tables are displayed, you can change any parameter you wish to change and create a new profile. You can leave the program at any point by entering Ctrl Z. ************************* GOOD LUCK ************************* DO YOU WANT YOUR PROFILES #2# DISPLAYED ON YOUR MONITOR: ENTER 0 AT PROMPT STORED IN FILE OUTPUT.IGR: ENTER 1 AT PROMPT DISPLAYED AND STORED: ENTER 2 AT PROMPT ------------------------------------------------------------ Enter / to use previous value(s) (see # .. #); Ctrl Z to exit / SELECT YOUR VARIABLE: #LAST:3# 1 LATITUDE 3 ALTITUDE 2 LONGITUDE 4 YEAR ------------------------------------------------------------ Enter / to use previous value(s) (see # .. #); Ctrl Z to exit / CHOOSE YOUR VARIABLE RANGE: BEGIN, END, STEPWIDTH ? # 100.0, 1000.0, 100.0# ------------------------------------------------------------ Enter / to use previous value(s) (see # .. #); Ctrl Z to exit / GEOD LATITUDE ? !NORTH! [DEGREE,DECIMAL] # 45.1# ------------------------------------------------------------ Enter / to use previous value(s) (see # .. #); Ctrl Z to exit / GEOD LONGITUDE ? !EAST! [DEGREE,DECIMAL] # 293.1# ------------------------------------------------------------ Enter / to use previous value(s) (see # .. #); Ctrl Z to exit / YEAR(EPOCH) ? *decimal* #1985.5# ------------------------------------------------------------ Enter / to use previous value(s) (see # .. #); Ctrl Z to exit / OUTPUT OPTION: B OR B/B0 ? #0# if you enter 0, the absolute magnetic field strength will be listed, otherwise the field strength normalized to the field strength at the magnetic equator is listed ------------------------------------------------------------ Enter / to use previous value(s) (see # .. #); Ctrl Z to exit / H/km DIMO B/Gauss B-NORTH B-EAST B-DOWN DIP DEC L-VALUE C 100.00 0.3042 0.52942 0.16107 -.05615 0.50118 71.2 -19.2 3.216 2 200.00 0.3042 0.50404 0.15406 -.05222 0.47707 71.2 -18.7 3.264 2 300.00 0.3042 0.48030 0.14745 -.04865 0.45451 71.1 -18.3 3.313 1 400.00 0.3042 0.45808 0.14122 -.04539 0.43340 71.1 -17.8 3.361 1 500.00 0.3042 0.43724 0.13534 -.04241 0.41360 71.1 -17.4 3.409 1 600.00 0.3042 0.41769 0.12978 -.03968 0.39503 71.0 -17.0 3.457 1 700.00 0.3042 0.39931 0.12453 -.03717 0.37757 71.0 -16.6 3.505 1 800.00 0.3042 0.38202 0.11955 -.03487 0.36115 71.0 -16.3 3.553 1 900.00 0.3042 0.36574 0.11483 -.03275 0.34570 70.9 -15.9 3.601 1 1000.00 0.3042 0.35039 0.11036 -.03080 0.33113 70.9 -15.6 3.649 1 ------- International Geomagnetic Reference Field --- Edition 1995 ------- LATI= 45.1 LONGI= 293.1 I DIMO is Dipol I C=1 L and B0 correct ALT= 1000.0 YEAR=1985.5 I Moment in Gauss I =2 wrong, =3 approx. -------------------------------------------------------------------------- **** DO YOU WANT TO CONTINUE? ------------------------------------------------------------ "0" QUIT AND EXIT "1" NEW PARAMETERS ------------------------------------------------------------ 1 **** WHICH PARAMETER DO YOU WANT TO CHANGE? ------------------------------------------------------------ 0 NO FURTHER CHANGES, CALCULATE PROFILE 1 LATITUDE # 45.1# 5 DISPLAY OR STORE 2 LONGITUDE # 293.1# 6 SELECTION OF VARIABLE #H/km# 3 ALTITUDE # 1000.0# 7 VARIABLE RANGE 4 YEAR #1985.5# # 100.0, 1000.0, 100.0# 8 B OR B/B0 ------------------------------------------------------------ ENTER NUMBER 7 CHOOSE YOUR VARIABLE RANGE: BEGIN, END, STEPWIDTH ? # 100.0, 1000.0, 100.0# ------------------------------------------------------------ Enter / to use previous value(s) (see # .. #); Ctrl Z to exit 1000,10000,1000 **** WHICH PARAMETER DO YOU WANT TO CHANGE? ------------------------------------------------------------ 0 NO FURTHER CHANGES, CALCULATE PROFILE 1 LATITUDE # 45.1# 5 DISPLAY OR STORE 2 LONGITUDE # 293.1# 6 SELECTION OF VARIABLE #H/km# 3 ALTITUDE # 1000.0# 7 VARIABLE RANGE 4 YEAR #1985.5# # 1000.0, 10000.0, 1000.0# 8 B OR B/B0 ------------------------------------------------------------ ENTER NUMBER 0 H/km DIMO B/Gauss B-NORTH B-EAST B-DOWN DIP DEC L-VALUE C 1000.00 0.3042 0.35039 0.11036 -.03080 0.33113 70.9 -15.6 3.649 1 2000.00 0.3042 0.23588 0.07614 -.01761 0.22256 70.7 -13.0 4.134 1 3000.00 0.3042 0.16675 0.05465 -.01091 0.15716 70.5 -11.3 4.622 1 4000.00 0.3042 0.12239 0.04049 -.00718 0.11528 70.4 -10.1 5.107 1 5000.00 0.3042 0.09256 0.03080 -.00494 0.08715 70.3 -9.1 5.596 1 6000.00 0.3042 0.07174 0.02396 -.00353 0.06753 70.3 -8.4 6.092 1 7000.00 0.3042 0.05675 0.01899 -.00260 0.05341 70.3 -7.8 6.593 1 8000.00 0.3042 0.04567 0.01530 -.00196 0.04299 70.3 -7.3 7.093 1 9000.00 0.3042 0.03731 0.01250 -.00151 0.03512 70.3 -6.9 7.590 1 10000.00 0.3042 0.03088 0.01034 -.00118 0.02907 70.3 -6.5 8.084 1 ------- International Geomagnetic Reference Field --- Edition 1995 ------- LATI= 45.1 LONGI= 293.1 I DIMO is Dipol I C=1 L and B0 correct ALT=10000.0 YEAR=1985.5 I Moment in Gauss I =2 wrong, =3 approx. -------------------------------------------------------------------------- !! REMINDER: this field model does not include external sources !! **** DO YOU WANT TO CONTINUE? ------------------------------------------------------------ "0" QUIT AND EXIT "1" NEW PARAMETERS ------------------------------------------------------------ 1 **** WHICH PARAMETER DO YOU WANT TO CHANGE? ------------------------------------------------------------ 0 NO FURTHER CHANGES, CALCULATE PROFILE 1 LATITUDE # 45.1# 5 DISPLAY OR STORE 2 LONGITUDE # 293.1# 6 SELECTION OF VARIABLE #H/km# 3 ALTITUDE # 10000.0# 7 VARIABLE RANGE 4 YEAR #1985.5# # 1000.0, 10000.0, 1000.0# 8 B OR B/B0 ------------------------------------------------------------ ENTER NUMBER 3 ALTITUDE ? [KM] #10000.0# ------------------------------------------------------------ Enter / to use previous value(s) (see # .. #); Ctrl Z to exit 1000 **** WHICH PARAMETER DO YOU WANT TO CHANGE? ------------------------------------------------------------ 0 NO FURTHER CHANGES, CALCULATE PROFILE 1 LATITUDE # 45.1# 5 DISPLAY OR STORE 2 LONGITUDE # 293.1# 6 SELECTION OF VARIABLE #H/km# 3 ALTITUDE # 1000.0# 7 VARIABLE RANGE 4 YEAR #1985.5# # 1000.0, 10000.0, 1000.0# 8 B OR B/B0 ------------------------------------------------------------ ENTER NUMBER 1 GEOD LATITUDE ? !NORTH! [DEGREE,DECIMAL] # 45.1# ------------------------------------------------------------ Enter / to use previous value(s) (see # .. #); Ctrl Z to exit 30 **** WHICH PARAMETER DO YOU WANT TO CHANGE? ------------------------------------------------------------ 0 NO FURTHER CHANGES, CALCULATE PROFILE 1 LATITUDE # 30.0# 5 DISPLAY OR STORE 2 LONGITUDE # 293.1# 6 SELECTION OF VARIABLE #H/km# 3 ALTITUDE # 1000.0# 7 VARIABLE RANGE 4 YEAR #1985.5# # 1000.0, 10000.0, 1000.0# 8 B OR B/B0 ------------------------------------------------------------ ENTER NUMBER 6 SELECT YOUR VARIABLE: #LAST:3# 1 LATITUDE 3 ALTITUDE 2 LONGITUDE 4 YEAR ------------------------------------------------------------ Enter / to use previous value(s) (see # .. #); Ctrl Z to exit 4 CHOOSE YOUR VARIABLE RANGE: BEGIN, END, STEPWIDTH ? # 1000.0, 10000.0, 1000.0# ------------------------------------------------------------ Enter / to use previous value(s) (see # .. #); Ctrl Z to exit 1945,2000,5 **** WHICH PARAMETER DO YOU WANT TO CHANGE? ------------------------------------------------------------ 0 NO FURTHER CHANGES, CALCULATE PROFILE 1 LATITUDE # 30.0# 5 DISPLAY OR STORE 2 LONGITUDE # 293.1# 6 SELECTION OF VARIABLE #YEAR# 3 ALTITUDE # 1000.0# 7 VARIABLE RANGE 4 YEAR #1985.5# # 1945.0, 2000.0, 5.0# 8 B OR B/B0 ------------------------------------------------------------ ENTER NUMBER 8 OUTPUT OPTION: B OR B/B0 ? #0# if you enter 0, the absolute magnetic field strength will be listed, otherwise the field strength normalized to the field strength at the magnetic equator is listed ------------------------------------------------------------ Enter / to use previous value(s) (see # .. #); Ctrl Z to exit 1 **** WHICH PARAMETER DO YOU WANT TO CHANGE? ------------------------------------------------------------ 0 NO FURTHER CHANGES, CALCULATE PROFILE 1 LATITUDE # 30.0# 5 DISPLAY OR STORE 2 LONGITUDE # 293.1# 6 SELECTION OF VARIABLE #YEAR# 3 ALTITUDE # 1000.0# 7 VARIABLE RANGE 4 YEAR #1985.5# # 1945.0, 2000.0, 5.0# 8 B OR B/B0 ------------------------------------------------------------ ENTER NUMBER 0 YEAR DIMO B/B0 B-NORTH B-EAST B-DOWN DIP DEC L-VALUE C 1945.00 0.3122 11.336 0.14421 -.02172 0.28536 62.9 -8.6 2.220 1 1950.00 0.3118 11.153 0.14491 -.02244 0.28389 62.7 -8.8 2.209 1 1955.00 0.3113 10.966 0.14586 -.02342 0.28222 62.4 -9.1 2.199 1 1960.00 0.3104 10.696 0.14713 -.02436 0.27977 61.9 -9.4 2.183 1 1965.00 0.3095 10.430 0.14811 -.02543 0.27705 61.5 -9.7 2.167 1 1970.00 0.3083 10.109 0.14947 -.02637 0.27367 61.0 -10.0 2.148 1 1975.00 0.3070 9.753 0.15071 -.02735 0.26975 60.4 -10.3 2.127 1 1980.00 0.3057 9.407 0.15096 -.02849 0.26521 59.9 -10.7 2.108 1 1985.00 0.3043 9.059 0.15128 -.02957 0.26104 59.4 -11.1 2.087 1 1990.00 0.3032 8.736 0.15164 -.03045 0.25712 59.0 -11.4 2.068 1 1995.00 0.3021 8.438 0.15201 -.03129 0.25348 58.5 -11.6 2.049 1 2000.00 0.3010 8.136 0.15245 -.03213 0.24982 58.1 -11.9 2.029 1 ------- International Geomagnetic Reference Field --- Edition 1995 ------- LATI= 30.0 LONGI= 293.1 I DIMO is Dipol I C=1 L and B0 correct ALT= 1000.0 YEAR=2000.0 I Moment in Gauss I =2 wrong, =3 approx. -------------------------------------------------------------------------- **** DO YOU WANT TO CONTINUE? ------------------------------------------------------------ "0" QUIT AND EXIT "1" NEW PARAMETERS ------------------------------------------------------------ 1 **** WHICH PARAMETER DO YOU WANT TO CHANGE? ------------------------------------------------------------ 0 NO FURTHER CHANGES, CALCULATE PROFILE 1 LATITUDE # 30.0# 5 DISPLAY OR STORE 2 LONGITUDE # 293.1# 6 SELECTION OF VARIABLE #YEAR# 3 ALTITUDE # 1000.0# 7 VARIABLE RANGE 4 YEAR #2000.0# # 1945.0, 2000.0, 5.0# 8 B OR B/B0 ------------------------------------------------------------ ENTER NUMBER 4 YEAR(EPOCH) ? *decimal* #2000.0# ------------------------------------------------------------ Enter / to use previous value(s) (see # .. #); Ctrl Z to exit 1995 **** WHICH PARAMETER DO YOU WANT TO CHANGE? ------------------------------------------------------------ 0 NO FURTHER CHANGES, CALCULATE PROFILE 1 LATITUDE # 30.0# 5 DISPLAY OR STORE 2 LONGITUDE # 293.1# 6 SELECTION OF VARIABLE #YEAR# 3 ALTITUDE # 1000.0# 7 VARIABLE RANGE 4 YEAR #1995.0# # 1945.0, 2000.0, 5.0# 8 B OR B/B0 ------------------------------------------------------------ ENTER NUMBER 6 SELECT YOUR VARIABLE: #LAST:4# 1 LATITUDE 3 ALTITUDE 2 LONGITUDE 4 YEAR ------------------------------------------------------------ Enter / to use previous value(s) (see # .. #); Ctrl Z to exit 1 CHOOSE YOUR VARIABLE RANGE: BEGIN, END, STEPWIDTH ? # 1945.0, 2000.0, 5.0# ------------------------------------------------------------ Enter / to use previous value(s) (see # .. #); Ctrl Z to exit -80,80,20 **** WHICH PARAMETER DO YOU WANT TO CHANGE? ------------------------------------------------------------ 0 NO FURTHER CHANGES, CALCULATE PROFILE 1 LATITUDE # 30.0# 5 DISPLAY OR STORE 2 LONGITUDE # 293.1# 6 SELECTION OF VARIABLE #LATI# 3 ALTITUDE # 1000.0# 7 VARIABLE RANGE 4 YEAR #1995.0# # -80.0, 80.0, 20.0# 8 B OR B/B0 ------------------------------------------------------------ ENTER NUMBER 0 LATI DIMO B/B0 B-NORTH B-EAST B-DOWN DIP DEC L-VALUE C -80.00 0.3021 310.936 0.10062 0.04627 -.30864 -70.3 24.7 6.596 1 -60.00 0.3021 11.854 0.13016 0.02925 -.21172 -57.8 12.7 2.425 1 -40.00 0.3021 2.105 0.13885 0.00865 -.12325 -41.5 3.6 1.509 1 -20.00 0.3021 1.022 0.15771 -.00938 -.03933 -14.0 -3.4 1.241 1 0.00 0.3021 1.300 0.17571 -.02240 0.07077 21.8 -7.3 1.274 1 20.00 0.3021 3.697 0.16803 -.02961 0.19495 48.8 -10.0 1.633 1 40.00 0.3021 24.795 0.12767 -.03161 0.30420 66.6 -13.9 2.843 1 60.00 0.3021 966.671 0.06388 -.02775 0.36494 79.2 -23.5 9.204 1 80.00 0.30219999.999 0.01140 -.01756 0.37588 86.8 -57.0 509.552 3 ------- International Geomagnetic Reference Field --- Edition 1995 ------- LATI= 80.0 LONGI= 293.1 I DIMO is Dipol I C=1 L and B0 correct ALT= 1000.0 YEAR=1995.0 I Moment in Gauss I =2 wrong, =3 approx. -------------------------------------------------------------------------- **** DO YOU WANT TO CONTINUE? ------------------------------------------------------------ "0" QUIT AND EXIT "1" NEW PARAMETERS ------------------------------------------------------------ 0 FORTRAN STOP