极光和气辉辐射的全球气辉模型
glowaurora的Python项目详细描述
使用python、matlab或mpi新建Modern Fortran GLOW。
发光0.973
辉光0.973气辉/极光model——现在是Python!
安装
需要Fortran编译器:
- Linux系统:
apt install gfortran
- macos:
brew install gcc
- Windows
pip install -e .
示例
自我测试f2py
这个自测应该没有错误。 这将测试来自python的fortran代码。:
pytest -v
容积排放率图
生成github站点上看到的图:
python RunGlow.py
选项包括:
-t --simtime
模拟时间-c --latlon
地理坐标(纬度,经度)[度]-q --flux
总通量
以体积排放率和中间过程为模型 给定一次电子沉淀输入。 你可以做得更多 通常在下一节中用作特征剖面。
生产/损耗率特征曲线
这需要两个步骤:
- 生成单位输入微分数通量与能量的关系
- 计算电离层能量沉积,从而产生/损失 模拟动力学化学的速率(总共12)
这由脚本gridaurora/MakeIonoEigenprofile.py
Matlab访问Glow
Matlab可以通过Python接口使用Glow,如示例所示
glow.m
。
论文
(感谢Stephen Kaeppler指出这些内容)
附录(一般用户不需要)
从stan solomon下载glow v0.973源代码
斯坦的团队已经使用现代Fortran发布了Glow v0.98,但是 原版的。
wget -r -np -nc -nH --cut-dirs=4 --random-wait --wait 1 -R "index.html*" http://download.hao.ucar.edu/pub/stans/glow/v0.973/
下载斯坦的IRI文件副本
斯坦稍微调整了一下IRI90,这是他使用的拷贝。
wget -r -np -nc -nH --cut-dirs=3 --random-wait --wait 1 -R "index.html*" http://download.hao.ucar.edu/pub/stans/iri/
自行编译Fortran代码
fortran程序自己使用时会吐出大量文本作为它的 输出:
cd bin
cmake ../src
cmake --build .
FORTRAN自检
极光示例
./auroraexample < aurexample.in > aurtest.dat
高能示例
./hexexample < hexexample.in > hextest.dat
注释
许可
目录fortran/
中的原始Fortran代码,从
http://download.hao.ucar.edu/pub/stans/glow/:“这个软件是
辉光模型。使用受开源学术研究的制约
包含在文件glowlicense.txt中的许可协议。“
python代码和对原始fortran代码的修改:gnu affero gplv3+