Skip to content

New pyiriplas #47

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
wants to merge 2 commits into
base: master
Choose a base branch
from
Open

New pyiriplas #47

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
969b830
update pyiriplas with 2020 code and constants
maaijke Dec 18, 2023
263a5ba
removed old igrf code
maaijke Dec 18, 2023
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
6 changes: 3 additions & 3 deletions RMextract/PosTools.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -315,7 +315,7 @@ def getPP(h=ION_HEIGHT,mPosition=[0.,0.,0.],direction=[0.,0.,0.]):

a = dx*dx + dy*dy + dz*dz;
b = x*dx + y*dy + z*dz;
alpha = (-b + sqrt(b*b - a*c))/a;
alpha = (-b + np.sqrt(b*b - a*c))/a;
pp_x = stationX + alpha*direction[0];
pp_y = stationY + alpha*direction[1]
pp_z = stationZ + alpha*direction[2];
Expand All @@ -324,7 +324,7 @@ def getPP(h=ION_HEIGHT,mPosition=[0.,0.,0.],direction=[0.,0.,0.]):
normal_y = pp_y * ion_ellipsoid_a2_inv;
normal_z = pp_z * ion_ellipsoid_b2_inv;
norm_normal2 = normal_x*normal_x + normal_y*normal_y + normal_z*normal_z;
norm_normal = sqrt(norm_normal2);
norm_normal = np.sqrt(norm_normal2);
sin_lat2 = normal_z*normal_z / norm_normal2;


Expand Down Expand Up @@ -567,7 +567,7 @@ def getProfile(source_pos, stat_pos, time, h = np.arange(60,20200,10)*u.km):
Rpp = fRloc(hdist)
#am = 1./((R_earthkm**2-hdist**2-(R_earthkm+h.to(u.km).value)**2)/(-2.*hdist*(R_earthkm+h.to(u.km).value)))
am = 1./((R_local**2-hdist**2-Rpp**2)/(-2.*hdist*Rpp))
return latpp, lonpp, latdir, londir, am
return latpp, lonpp, latdir, londir, am, azel



Expand Down
3,376 changes: 1,701 additions & 1,675 deletions RMextract/pyiriplas/Iris2017.for → RMextract/pyiriplas/Iris2020.for
View file
Open in desktop

Large diffs are not rendered by default.

4,561 changes: 2,862 additions & 1,699 deletions RMextract/pyiriplas/apf107.dat
View file
Open in desktop

Large diffs are not rendered by default.

42,368 changes: 21,184 additions & 21,184 deletions RMextract/pyiriplas/aprz.dat
View file
Open in desktop

Large diffs are not rendered by default.

159 changes: 81 additions & 78 deletions RMextract/pyiriplas/f107_12.dat
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,78 +1,81 @@
14,03,2020,

1,1948,12,2021

1948 191.4 189.5 183.3 173.9 168.3 171.7 172.8 179.6 186.5 191.3 186.2 178.3
1949 176.6 172.1 173.0 175.5 179.6 178.6 175.0 171.2 165.1 161.8 159.9 159.7
1950 158.2 155.3 148.9 139.9 133.8 129.6 124.7 121.5 120.3 119.5 119.7 120.4
1951 119.2 117.6 117.1 115.8 112.2 112.3 115.4 114.7 112.1 111.8 108.8 102.0
1952 96.3 94.6 92.7 88.8 86.9 85.2 83.1 81.0 80.9 81.3 80.4 79.7
1953 79.2 78.5 76.9 75.0 74.3 73.2 71.8 70.8 71.1 71.6 70.8 70.3
1954 70.2 69.8 69.3 69.2 69.4 69.7 70.4 71.7 72.0 72.4 73.9 75.8
1955 77.0 78.8 81.6 84.4 88.4 92.7 97.8 104.6 110.9 117.2 123.3 129.2
1956 134.9 141.4 150.9 159.4 170.1 180.6 187.6 188.7 191.1 194.8 195.4 203.1
1957 211.1 217.7 222.5 228.5 232.0 230.8 233.5 234.2 234.6 237.9 242.1 241.3
1958 241.7 246.2 245.9 242.4 238.1 235.9 234.6 230.2 231.6 230.1 229.5 227.4
1959 223.0 228.9 223.1 217.5 213.2 210.3 208.1 197.3 196.1 194.0 188.8 184.9
1960 179.3 178.6 174.2 167.7 166.3 162.4 157.3 149.4 145.2 143.9 138.0 134.6
1961 132.0 128.1 122.7 116.7 112.6 107.2 102.8 102.3 101.5 102.4 102.6 101.6
1962 99.1 94.8 93.3 91.3 89.9 89.7 89.3 88.5 85.7 84.3 84.1 83.2
1963 82.4 82.3 82.9 82.2 80.8 80.5 80.4 80.1 79.6 79.4 79.0 77.8
1964 76.8 75.9 74.9 73.9 73.1 72.8 72.9 73.0 72.9 72.5 73.1 73.5
1965 74.0 74.8 75.3 76.1 75.8 76.1 77.0 77.6 78.8 80.7 83.3 85.0
1966 86.8 89.7 92.1 93.9 96.3 100.1 104.9 109.8 116.0 119.8 123.5 124.7
1967 126.6 130.5 131.2 135.5 137.1 142.9 146.3 147.4 149.2 146.7 148.4 147.2
1968 146.6 147.8 148.8 151.2 148.9 150.0 149.6 146.6 147.1 147.6 150.4 150.5
1969 149.7 149.8 147.6 149.6 150.5 149.2 152.2 153.7 154.8 153.0 154.2 156.3
1970 154.3 152.8 153.9 154.1 153.8 154.4 156.2 155.9 150.6 147.1 141.9 137.2
1971 133.2 131.8 130.7 126.9 124.4 121.9 117.4 116.2 116.6 116.7 117.5 119.1
1972 121.6 120.0 120.7 123.1 121.0 119.4 118.7 117.3 115.2 112.9 113.4 110.2
1973 105.9 102.7 100.6 97.8 94.9 93.9 92.8 91.4 88.1 88.6 88.2 87.2
1974 87.4 87.3 87.7 86.2 85.8 86.0 84.9 85.1 84.2 82.9 82.4 81.4
1975 81.0 81.0 80.6 79.6 78.0 76.7 75.5 73.9 75.2 76.1 75.8 76.5
1976 76.3 76.4 74.2 73.7 73.8 73.3 73.8 73.9 74.7 74.7 75.5 77.1
1977 77.9 78.9 81.3 83.0 84.0 85.3 88.8 94.2 97.9 104.6 110.9 115.7
1978 118.6 118.3 124.6 128.4 132.8 139.6 148.6 157.9 157.9 161.0 163.2 164.4
1979 164.3 167.5 175.9 181.4 187.7 190.4 194.1 193.5 188.6 189.6 193.2 198.2
1980 198.4 198.7 201.9 198.3 196.4 198.4 196.7 198.4 197.8 201.3 199.5 193.3
1981 197.5 197.9 201.9 203.2 205.4 207.4 200.5 201.5 202.6 199.1 194.4 192.7
1982 195.3 190.2 185.3 184.0 180.3 177.0 172.8 167.7 160.5 153.6 151.5 148.6
1983 149.2 147.6 142.7 137.8 130.3 123.9 118.4 118.5 120.0 120.9 123.6 120.6
1984 114.5 108.8 106.1 102.3 99.1 100.0 99.4 96.6 91.7 88.4 84.3 80.1
1985 78.3 76.8 75.4 74.8 74.2 74.3 73.9 75.3 75.9 75.9 75.9 74.8
1986 74.5 72.7 72.8 73.8 74.1 74.3 73.8 74.0 73.3 73.0 74.8 75.7
1987 77.2 78.8 80.0 81.3 82.3 84.4 86.7 89.0 93.1 96.4 99.0 104.1
1988 109.4 115.5 119.8 125.8 130.8 136.5 146.8 156.1 164.8 170.0 177.2 184.6
1989 185.4 191.7 198.8 204.2 209.0 209.1 214.9 207.5 204.5 204.9 203.3 203.7
1990 198.9 204.7 200.0 196.1 193.2 192.1 192.6 191.9 200.4 201.4 202.1 202.2
1991 203.8 204.9 201.9 205.9 205.4 206.9 207.4 208.0 205.2 198.1 194.7 184.5
1992 178.6 169.5 166.1 162.6 160.4 158.2 148.6 143.6 136.5 131.1 126.4 125.7
1993 125.7 120.4 117.7 117.7 114.8 111.6 110.5 109.3 105.9 100.3 97.9 94.3
1994 91.0 90.2 89.5 90.2 88.3 87.5 85.4 83.5 82.3 80.7 81.2 80.7
1995 80.5 79.6 79.1 79.4 78.3 77.5 76.8 76.1 75.2 73.6 73.5 73.3
1996 72.7 72.7 72.3 72.0 71.8 71.8 71.8 71.7 72.2 72.8 72.6 73.2
1997 73.4 74.0 75.8 76.8 79.0 80.6 82.0 83.1 84.7 89.3 91.2 94.4
1998 98.3 104.0 107.5 108.6 113.4 117.2 119.7 121.9 124.6 127.4 129.3 134.1
1999 139.6 144.4 145.3 147.9 151.0 151.3 152.3 154.6 163.3 166.2 170.1 174.8
2000 178.8 177.1 174.8 178.9 179.6 178.3 175.8 178.7 177.9 177.1 173.5 173.3
2001 169.7 167.5 170.5 172.2 178.0 180.1 188.0 188.6 185.8 189.6 189.7 188.2
2002 188.9 194.9 197.6 190.7 188.0 185.7 176.2 167.2 162.3 159.0 153.1 150.0
2003 149.3 146.3 139.3 137.3 135.7 130.0 126.6 124.2 124.3 118.9 117.8 118.5
2004 117.7 116.6 114.3 114.3 111.1 107.4 104.5 104.1 103.5 101.4 100.8 102.0
2005 101.3 99.2 98.0 95.8 92.9 91.3 90.2 88.7 87.1 88.5 87.4 85.3
2006 84.3 83.3 82.1 79.4 80.6 80.7 80.2 79.5 79.5 79.7 77.7 77.4
2007 77.1 76.8 76.1 75.0 74.3 73.6 72.9 72.2 72.2 72.2 71.8 70.5
2008 70.2 69.9 69.4 69.5 69.6 69.6 68.9 68.6 68.6 68.4 68.6 68.7
2009 68.8 68.9 69.1 69.3 69.5 70.2 71.2 72.6 73.6 74.2 74.6 74.6
2010 75.3 75.8 76.9 78.3 79.1 79.9 79.8 81.1 83.7 86.3 88.2 90.0
2011 92.0 92.8 95.5 100.2 107.2 111.8 115.8 117.2 116.4 116.3 116.1 119.4
2012 123.5 127.4 128.2 127.1 124.7 121.0 118.5 118.1 118.1 119.2 121.3 121.5
2013 119.5 119.5 117.5 116.8 118.0 121.8 125.5 129.2 134.3 134.9 134.3 133.8
2014 135.1 135.0 139.1 143.5 146.2 147.5 145.7 144.4 139.2 137.3 135.6 134.9
2015 135.0 133.2 131.1 126.9 123.6 119.8 116.6 113.8 111.4 108.7 105.5 102.2
2016 99.8 97.7 96.9 95.2 92.9 90.9 87.7 85.7 83.3 82.8 81.9 80.7
2017 80.2 79.4 80.2 78.4 77.7 76.9 76.2 75.9 74.0 75.3 75.0 75.0
2018 74.7 73.8 73.1 70.8 69.9 69.5 69.7 69.9 70.1 70.5 70.9 70.7
2019 70.1 69.9 69.5 69.2 69.2 69.5 69.7 69.8 71.5 71.9 72.4 73.0
2020 73.6 74.2 74.9 75.6 76.3 77.1 77.9 78.8 79.7 80.6 81.6 82.5
2021 83.6 84.6 85.7 86.7 87.9 89.0 90.1 91.3 92.5 93.7 95.0 96.2
08,03,2023,

1,1948,12,2024

1948 191.4 189.5 183.3 173.9 168.3 171.7 172.8 179.6 186.5 191.3 186.2 178.3
1949 176.6 172.1 173.0 175.5 179.6 178.6 175.0 171.2 165.1 161.8 159.9 159.7
1950 158.2 155.3 148.9 139.9 133.8 129.6 124.7 121.5 120.3 119.5 119.7 120.4
1951 119.2 117.6 117.1 115.8 112.2 112.3 115.4 114.7 112.1 111.8 108.8 102.0
1952 96.3 94.6 92.7 88.8 86.9 85.2 83.1 81.0 80.9 81.3 80.4 79.7
1953 79.2 78.5 76.9 75.0 74.3 73.2 71.8 70.8 71.1 71.6 70.8 70.3
1954 70.2 69.8 69.3 69.2 69.4 69.7 70.4 71.7 72.0 72.4 73.9 75.8
1955 77.0 78.8 81.6 84.4 88.4 92.7 97.8 104.6 110.9 117.2 123.3 129.2
1956 134.9 141.4 150.9 159.4 170.1 180.6 187.6 188.7 191.1 194.8 195.4 203.1
1957 211.1 217.7 222.5 228.5 232.0 230.8 233.5 234.2 234.6 237.9 242.1 241.3
1958 241.7 246.2 245.9 242.4 238.1 235.9 234.6 230.2 231.6 230.1 229.5 227.4
1959 223.0 228.9 223.1 217.5 213.2 210.3 208.1 197.3 196.1 194.0 188.8 184.9
1960 179.3 178.6 174.2 167.7 166.3 162.4 157.3 149.4 145.2 143.9 138.0 134.6
1961 132.0 128.1 122.7 116.7 112.6 107.2 102.8 102.3 101.5 102.4 102.6 101.6
1962 99.1 94.8 93.3 91.3 89.9 89.7 89.3 88.5 85.7 84.3 84.1 83.2
1963 82.4 82.3 82.9 82.2 80.8 80.5 80.4 80.1 79.6 79.4 79.0 77.8
1964 76.8 75.9 74.9 73.9 73.1 72.8 72.9 73.0 72.9 72.5 73.1 73.5
1965 74.0 74.8 75.3 76.1 75.8 76.1 77.0 77.6 78.8 80.7 83.3 85.0
1966 86.8 89.7 92.1 93.9 96.3 100.1 104.9 109.8 116.0 119.8 123.5 124.7
1967 126.6 130.5 131.2 135.5 137.1 142.9 146.3 147.4 149.2 146.7 148.4 147.2
1968 146.6 147.8 148.8 151.2 148.9 150.0 149.6 146.6 147.1 147.6 150.4 150.5
1969 149.7 149.8 147.6 149.6 150.5 149.2 152.2 153.7 154.8 153.0 154.2 156.3
1970 154.3 152.8 153.9 154.1 153.8 154.4 156.2 155.9 150.6 147.1 141.9 137.2
1971 133.2 131.8 130.7 126.9 124.4 121.9 117.4 116.2 116.6 116.7 117.5 119.1
1972 121.6 120.0 120.7 123.1 121.0 119.4 118.7 117.3 115.2 112.9 113.4 110.2
1973 105.9 102.7 100.6 97.8 94.9 93.9 92.8 91.4 88.1 88.6 88.2 87.2
1974 87.4 87.3 87.7 86.2 85.8 86.0 84.9 85.1 84.2 82.9 82.4 81.4
1975 81.0 81.0 80.6 79.6 78.0 76.7 75.5 73.9 75.2 76.1 75.8 76.5
1976 76.3 76.4 74.2 73.7 73.8 73.3 73.8 73.9 74.7 74.7 75.5 77.1
1977 77.9 78.9 81.3 83.0 84.0 85.3 88.8 94.2 97.9 104.6 110.9 115.7
1978 118.6 118.3 124.6 128.4 132.8 139.6 148.6 157.9 157.9 161.0 163.2 164.4
1979 164.3 167.5 175.9 181.4 187.7 190.4 194.1 193.5 188.6 189.6 193.2 198.2
1980 198.4 198.7 201.9 198.3 196.4 198.4 196.7 198.4 197.8 201.3 199.5 193.3
1981 197.5 197.9 201.9 203.2 205.4 207.4 200.5 201.5 202.6 199.1 194.4 192.7
1982 195.3 190.2 185.3 184.0 180.3 177.0 172.8 167.7 160.5 153.6 151.5 148.6
1983 149.2 147.6 142.7 137.8 130.3 123.9 118.4 118.5 120.0 120.9 123.6 120.6
1984 114.5 108.8 106.1 102.3 99.1 100.0 99.4 96.6 91.7 88.4 84.3 80.1
1985 78.3 76.8 75.4 74.8 74.2 74.3 73.9 75.3 75.9 75.9 75.9 74.8
1986 74.5 72.7 72.8 73.8 74.1 74.3 73.8 74.0 73.3 73.0 74.8 75.7
1987 77.2 78.8 80.0 81.3 82.3 84.4 86.7 89.0 93.1 96.4 99.0 104.1
1988 109.4 115.5 119.8 125.8 130.8 136.5 146.8 156.1 164.8 170.0 177.2 184.6
1989 185.4 191.7 198.8 204.2 209.0 209.1 214.9 207.5 204.5 204.9 203.3 203.7
1990 198.9 204.7 200.0 196.1 193.2 192.1 192.6 191.9 200.4 201.4 202.1 202.2
1991 203.8 204.9 201.9 205.9 205.4 206.9 207.4 208.0 205.2 198.1 194.7 184.5
1992 178.6 169.5 166.1 162.6 160.4 158.2 148.6 143.6 136.5 131.1 126.4 125.7
1993 125.7 120.4 117.7 117.7 114.8 111.6 110.5 109.3 105.9 100.3 97.9 94.3
1994 91.0 90.2 89.5 90.2 88.3 87.5 85.4 83.5 82.3 80.7 81.2 80.7
1995 80.5 79.6 79.1 79.4 78.3 77.5 76.8 76.1 75.2 73.6 73.5 73.3
1996 72.7 72.7 72.3 72.0 71.8 71.8 71.8 71.7 72.2 72.8 72.6 73.2
1997 73.4 74.0 75.8 76.8 79.0 80.6 82.0 83.1 84.7 89.3 91.2 94.4
1998 98.3 104.0 107.5 108.6 113.4 117.2 119.7 121.9 124.6 127.4 129.3 134.1
1999 139.6 144.4 145.3 147.9 151.0 151.3 152.3 154.6 163.3 166.2 170.1 174.8
2000 178.8 177.1 174.8 178.9 179.6 178.3 175.8 178.7 177.9 177.1 173.5 173.3
2001 169.7 167.5 170.5 172.2 178.0 180.1 188.0 188.6 185.8 189.6 189.7 188.2
2002 188.9 194.9 197.6 190.7 188.0 185.7 176.2 167.2 162.3 159.0 153.1 150.0
2003 149.3 146.3 139.3 137.3 135.7 130.0 126.6 124.2 124.3 118.9 117.8 118.5
2004 117.7 116.6 114.3 114.3 111.1 107.4 104.5 104.1 103.5 101.4 100.8 102.0
2005 101.3 99.2 98.0 95.8 92.9 91.3 90.2 88.7 87.1 88.5 87.4 85.3
2006 84.3 83.3 82.1 79.4 80.6 80.7 80.2 79.5 79.5 79.7 77.7 77.4
2007 77.1 76.8 76.1 75.0 74.3 73.6 72.9 72.2 72.2 72.2 71.8 70.5
2008 70.2 69.9 69.4 69.5 69.6 69.6 68.9 68.6 68.6 68.4 68.6 68.7
2009 68.8 68.9 69.1 69.3 69.5 70.2 71.2 72.6 73.6 74.2 74.6 74.6
2010 75.3 75.8 76.9 78.3 79.1 79.9 79.8 81.1 83.7 86.3 88.2 90.0
2011 92.0 92.8 95.5 100.2 107.2 111.8 115.8 117.2 116.4 116.3 116.1 119.4
2012 123.5 127.4 128.2 127.1 124.7 121.0 118.5 118.1 118.1 119.2 121.3 121.5
2013 119.5 119.5 117.5 116.8 118.0 121.8 125.5 129.2 134.3 134.9 134.3 133.8
2014 135.1 135.0 139.1 143.5 146.2 147.5 145.7 144.4 139.2 137.3 135.6 134.9
2015 135.0 133.2 131.1 126.9 123.6 119.8 116.6 113.8 111.4 108.7 105.5 102.2
2016 99.8 97.7 96.9 95.2 92.9 90.9 87.7 85.7 83.3 82.8 81.9 80.7
2017 80.2 79.4 80.2 78.4 77.7 76.9 76.2 75.9 74.0 75.3 75.0 75.0
2018 74.7 73.8 73.1 70.8 69.9 69.5 69.7 69.9 70.1 70.5 70.9 70.7
2019 70.1 69.9 69.5 69.2 69.2 69.5 69.7 69.8 69.8 69.9 69.6 69.5
2020 69.5 69.9 70.1 70.5 72.2 73.4 73.8 73.9 74.7 75.1 74.5 75.8
2021 77.7 78.4 79.3 80.7 81.5 82.1 83.3 86.4 89.5 94.1 99.4 101.7
2022 105.6 108.0 111.7 114.3 116.5 122.9 128.5 135.6 133.8 136.0 140.7 142.9
2023 147.5 149.7 154.2 156.3 160.6 162.7 166.7 168.7 172.5 174.3 177.8 179.4
2024 182.5 183.9 186.6 187.9 190.1 191.1 192.9 193.7 195.0 195.5 196.3 196.6
Loading