又过了半年,已经完全忘了当时8月份的python工作了。这半年,先是跑合同,然后张罗出海,完了又搞了两次验收,还做了两次汇报,开了无数的会,忙坏了。
现在,得重新张罗电磁这档子事儿了。下面记录下最近的操作。
1 打开编程环境。
都忘了怎么打开环境了,翻一下以前的博客(这是第三篇),终于想起来了。
打开anaconda,然后激活环境e1(总共装了2个,1个e1,1个e2,忘了哪个能用了,随便激活1个)。然后打开spyder。
下面是spyder界面:
这是4.2.5版本,python用的3.7
2 打开simpeg文件
装这个anaconda主要是为了跑simpeg,里面最感兴趣的就是频率域三维电磁了。下面就逐步分析一下这个文件。
2.1 引入依赖包
这里的代码是如下。首先来盘点一下:
- 从discretize引入treemesh包,实现树形网格;
- 从discretize.utils引入mkvc和refine_tree_xyz,实现向量组装和细化数网格;
- 从SimPEG.utils引入plot2Ddata和su