全局核心等离子体模型的Python包装器
PyGCPM的Python项目详细描述
PyGCPM公司
PyGCPM是用于全局核心等离子模型的python3包装器 (Gallagher等人,2000年)。此模块使用来自 https://plasmasphere.nasa.gov/models/
原始代码被稍微修改了一下,这样生成的代码就更少了 输出到终端,并传递使用的数据文件的路径 通过C包装器,以便Fortran代码能够找到它们。在
值得注意的是,该模型的IGRF组件需要 更新-2010年以后可能会发生奇怪的事情。我会的 解决这个问题。在
安装
使用pip
:
python3 -m pip install PyGCPM --user
或者
^{pr2}$本次回购:
git clone https://github.com/mattkjames7/PyGCPM.git
cd PyGCPM/
python3 setup.py install --user
使用
到目前为止,只有两个工作函数GCPM
和{PlotMLTSlice
函数正在做一些奇怪的事情。值得注意的是
这个模型需要一些时间来运行大量的点。在
GCPM
这一功能在太阳磁层(SM)中占有一席之地 坐标系以及日期和时间将生成一个模型 每个位置的输出。在
importPyGCPMne,nH,nHe,nO=PyGCPM.GCPM(x,y,z,Date,ut,Kp=Kp,Verbose=Verbose)
其中输入为:
Variable | Data Type | Description |
---|---|---|
^{ | ^{ | Scalar or array x position(s) in RE |
^{ | ^{ | Scalar or array y position(s) in RE |
^{ | ^{ | Scalar or array z position(s) in RE |
^{ | ^{ | Integer date in format yyyymmdd - if provided with a scalar then all positions will use the same date, otherwise an array can be provided for a time series |
^{ | ^{ | Floating point hours where ^{ |
^{ | ^{ | (optional) Kp index, as with ^{ |
^{ | ^{ | (optional) If ^{ |
功能输出:
^{tb2}$PlotEqSlice
此函数将绘制穿过SM赤道的切片(zSM=0)。在
importPyGCPMPyGCPM.PlotEqSlice(Date,ut,Parameter='ne',Rmax=10.0,dR=0.5,Kp=1.0,fig=None,maps=[1,1,0,0],zlog=True,cmap='gnuplot',scale=None,Verbose=False)^{3}$
参考文献
Gallagher,D.L.,Craven,p.D.,&Comfort,R.H.(2000年8月)。全球核心等离子体模型。 J、 地球物理学。自然科学,105(A8),18819-18833。doi:10.1029/1999JA000241
- 项目
标签: