Comments (5)
from geopack.
from geopack.
Hi Sheng,
I am running my code as follows:
from geopack import geopack
import numpy as np
import pandas as pd
glats = [45.871,48.027,50.199,52.035,53.856,54.929,56.376,57.707,58.763,62.824]
glons = [264.917,263.635,263.960,264.537,265.340,264.713,265.300,265.794,265.920,267.890]
data = pd.Series([2.008e+03, 1.000e+00, 2.100e+01, 3.500e+01, -2.140e+00,
1.560e+00, 2.200e-01, -3.514e+02, -2.200e+00, -7.900e+00,
7.750e+00, 3.131e+04, -4.000e+00, 1.000e+00, 1.000e+00,
-3.380e-01, 1.850e+00, 0.000e+00, 0.000e+00, 0.000e+00,
0.000e+00, 0.000e+00, 0.000e+00])
data.columns = ['IYEAR', 'IDAY', 'IHOUR', 'MIN', 'BXGSM', 'BYGSM', 'BZGSM', 'VXGSE',
'VYGSE', 'VZGSE', 'DEN', 'TEMP', 'SYMH', 'IMFFLAG', 'ISWFLAG', 'TILT',
'Pdyn', 'W1', 'W2', 'W3', 'W4', 'W5', 'W6']
ut = data['IHOUR']*3600 + data['MIN']*60
ps = geopack.recalc(ut)
Re = 6371.2
x, y,= geopack.geodgeo(100,np.deg2rad(glats[0]),1)
xgsm, ygsm, zgsm = geopack.geogsm(x,y,100,1)
xgsm /= Re
ygsm /= Re
zgsm /= Re
parmod = [data['Pdyn'][0],data['SYMH'][0],data['BYGSM'][0],data['BZGSM'][0]]
geopack.trace(xgsm,ygsm,zgsm,1,parmod=parmod,rlim=100,exname='t96')
Thanks,
Rhys
from geopack.
from geopack.
from geopack.
Related Issues (18)
- Geopack.trace not providing vector of output values HOT 2
- Releasing new version to PyPI? HOT 1
- How to reference? HOT 1
- scipy and M1 arm64
- Feature Request: Evaluate models at multiple points HOT 1
- np.float is deprecated since numpy v1.24.0 HOT 1
- Getting latest bug fixes into PyPI? HOT 2
- geopack on conda-forge? HOT 1
- Bug report on t96.py HOT 6
- T01 IDL vs Python differences HOT 4
- FAILED geopack/test_geopack1.py::KnownValues::test_to_known_values - AssertionError: -0.53356312486214 != -0.5335631237417825
- Geopack says it has no attribute called geogsm, igrf_geo, etc. HOT 2
- Update to IGRF-13 HOT 2
- Consider adding project to list at https://heliopython.org/projects/ HOT 2
- inconsistent comparison results HOT 2
- some errors in t96.py HOT 2
- Difference in output of t96 and t01 model between IDL and Python geopacks! HOT 4
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from geopack.