全局核心等离子体模型的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)

其中输入为:

VariableData TypeDescription
^{}^{} or ^{}Scalar or array x position(s) in RE
^{}^{} or ^{}Scalar or array y position(s) in RE
^{}^{} or ^{}Scalar or array z position(s) in RE
^{}^{} or ^{}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
^{}^{} or ^{}Floating point hours where ^{}. As with ^{} - if provided with a scalar then all positions will use the same time, otherwise an array can be provided for a time series
^{}^{} or ^{}(optional) Kp index, as with ^{} and ^{} - this can be a scalar or an array
^{}^{}(optional) If ^{} then the function will display its progress in the terminal

功能输出:

^{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

欢迎加入QQ群-->: 979659372 Python中文网_新手群

推荐PyPI第三方库


热门话题
java返回。p12文件到客户端,而不创建密钥库文件   java创建聚合表的最佳方法是什么?   html Java Thymeleaf:从一个表单获取数据到另一个表单,然后在提交时使用这两个表单的数据   对象声明和初始化级别的java Jackson@JsonIgnoreProperties   java如何从一串数字中获取数值?   java JavaFX Hibernate返回空指针异常   java我需要在Netbeans中更改什么,以便它构建一个混合的JavaFX应用程序?   如何保护我的Java AES加密密钥   即使我关机,Java线程调度程序仍然有效   java HTTP状态500找不到元素“bean”的声明   java如何将电子邮件另存为。带有XUnsent的eml:1正在使用消息。写   AmazonWeb服务如何使用JavaSDK使Cloudfront中的缓存失效   java将月份转换为毫秒的最佳方法   在一行中使用相同的Java SAXParser   java获取异常组织。openqa。硒。SessionNotCreatedException