After complete tutorial I and tutorial II, you can follow the README file to run the simulation.

The Trelis software interface is shown below.

Panels that are frequently used::

1. Command Panel: For quick geometry creation, segmentation, combination and mesh division
2. Command Line: Complete the same function as command panel but by the way of the command lines.
3. Tool bar: The first five lines are create, open, save, import and export files. The middle part includes running the journal script file and adjusting the observation mode.
4. Model Tree: Each part of the model is displayed at the level of body-volume-surface-curve-vertex. And each element has a corresponding ID number for designation.
   More details can be found in Official Guide. Note that each step of Trelis has a corresponding script command, which is convenient for saving and sharing.

2.Create Geometry

2.1. Create a volume

First create a 600km x 600km x 400km volume. Since the default origin is in the center of the new volume, we will move the volume down by 200km. Thus the upside surface represents the surface.
The following script can be copied and pasted directly to the command line.

1
2
3
4
5
6
7
8
9
10
11
12
13

# Like shell script,Trelis uses '#' for comments# ----------------------------------------------------------------------
# Create block
# ----------------------------------------------------------------------

# Block is 600 km x 600 km x 400 km
# -300 km <= x <= 300 km
# -300 km <= y <= 300 km
# -400 km <= z <= 0 km
reset
brick x 600000 y 600000 z 400000

# Translate block so the top is at z=0
volume 1 move x 0 y 0 z -200000

In addition to the above method, Trelis can also create a volume by:

creating vertexes–> creating surfaces by specified vertex–>create a volume by specify surfaces
Although it is more cumbersome, an example also shows here for reference.

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54

# ----------------------------------------------------------------------
# Set units to SI.
# ----------------------------------------------------------------------
${Units('si')}
# ----------------------------------------------------------------------
# Create vertex
# ----------------------------------------------------------------------
# A cuboid has 8 corner points
reset
create vertex  { 300*km}  { 300*km} {   0*km}
create vertex  { 300*km}  {-300*km} {   0*km}
create vertex  {-300*km}  {-300*km} {   0*km}
create vertex  {-300*km}  { 300*km} {   0*km}
create vertex  { 300*km}  { 300*km} {-400*km}
create vertex  { 300*km}  {-300*km} {-400*km}
create vertex  {-300*km}  {-300*km} {-400*km}
create vertex  {-300*km}  { 300*km} {-400*km}
# ----------------------------------------------------------------------
# Create curve
# ----------------------------------------------------------------------
# 8 corner points connected to 12 side lines 
# top 4 lines
create curve vertex 1 2  
create curve vertex 2 3  
create curve vertex 3 4
create curve vertex 4 1 
# four bottom lines
create curve vertex 5 6  
create curve vertex 6 7  
create curve vertex 7 8
create curve vertex 8 5 
# four side lines (note that the id number of the element in the process of creating a new geometry may also change, so we can give names to important elements.)
create curve vertex 1 5 
create curve vertex 6 9 
create curve vertex 7 10 
create curve vertex 11 15
# ----------------------------------------------------------------------
# Create surface
# ----------------------------------------------------------------------
# 12 lines connected to 6 surfaces
# Top surface
create surface curve 1 2 3 4  
# Bottom surface
create surface curve 5 6 7 8
# Side surfaces
create surface curve 1 5 9 10
create surface curve 2 6 10 11
create surface curve 3 7 11 12
create surface curve 4 8 9 12
# ----------------------------------------------------------------------
# Create volume
# ----------------------------------------------------------------------
# Create the volume 
create volume surface 1 2 3 4 5 6

The results of the above two ways are exactly the same. For more complicated examples, you can choose your own way according to the situation. Note that if you use the second method to establish the geometry, the id number of most elements will change for the left part of this manual.

2.2. Split the cuboid into two parts

Here we create a plane to divide the cuboid into a upper part as the elastic layer and a lower part as the viscoelastic layer. And a fault plane will be created for simulation.

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32

# ----------------------------------------------------------------------
# Create interface surfaces
# ----------------------------------------------------------------------
create vertex 0 -300000 0
create vertex 0 300000 0
create vertex -200000  0 -400000
create planar surface with plane vertex 9 10 11   
surface 7 name "fault_surface"
delete vertex 9 10 11
create planar surface with plane zplane offset -20000
surface 8 name "material_interface"

# ----------------------------------------------------------------------
# Divide volumes using interface surfaces
# ----------------------------------------------------------------------

webcut volume 1 with plane surface fault_surface
webcut volume 1 with plane surface material_interface
webcut volume 4 with plane surface material_interface
volume 1 name "elastic_xpos"
volume 4 name "elastic_xneg"
volume 5 name "visco_xpos"
volume 6 name "visco_xneg"
# ----------------------------------------------------------------------
# create fault_plane
# ----------------------------------------------------------------------
create planar surface with plane yplane offset -20000 
create planar surface with plane yplane offset 20000 
split surface 22 with surface 39
split surface 42 with surface 40
surface 43 name "fault_plane"
delete body 2 3 7 8

A more detailed segmentation process can be found in the official tutorial . (A proper vpn tool is required)

2.3. imprint & merge

After creating the split the volume, the points and lines of each sub volumes are probably not the same, so they need to be merged, which takes a certain amount of time under complex geometry.

1
2
3
4
5

# ----------------------------------------------------------------------
# Imprint all volumes, merging surfaces
# ----------------------------------------------------------------------
imprint all with volume all
merge all

3.Meshing

We need to specify the size of each curve, surface or volume before meshing them. And the map meshing scheme is set by default. You can use the command panal or type help volume scheme in the command line to view all the possible meshing schemes. The tetmesh scheme is applied in this example.

1
2
3
4
5
6
7
8
9
10
11
12

volume all scheme tetmesh 
curve all scheme default
# ----------------------------------------------------------------------
# Set discretization size
# ----------------------------------------------------------------------
surface 43 size 1000
volume all size 100000
# ----------------------------------------------------------------------
# Generate the mesh
# ----------------------------------------------------------------------
volume all scheme tetmesh proximity layers off geometric sizing on
mesh volume all

Now the mesh part has completed, and we recommend you to do mesh QA. You can use the command panal or type

1

quality volume all shape global draw mesh

in the command line to do mesh QA.Generally the minimum value should be greater than 0.2 for a reliable Pylith calculation.

The mesh quality of this example is as follows

4.Declare each block and the boundary surfaces

Then we need to declare each block and the bounding surfaces for calculations and determining the boundary conditions of Pylith.

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81

# ----------------------------------------------------------------------
# Create blocks for materials
# ----------------------------------------------------------------------
block 1 volume 1 4
block 1 name "elastic"
block 2 volume 5 6
block 2 name "viscoelastic"

# ----------------------------------------------------------------------
# Create nodeset for fault
# ----------------------------------------------------------------------
group "fault" add node in fault_plane
nodeset 10 group fault
nodeset 10 name "fault"

# ----------------------------------------------------------------------
# Create sideset for faces in fault
# ----------------------------------------------------------------------
group "fault_faces" add face in fault_plane
sideset 10 group fault_faces
sideset 10 name "fault_faces"

# ----------------------------------------------------------------------
# Create nodeset for -x face
# ----------------------------------------------------------------------
group "face_xpos" add node in surface 20
group "face_xpos" add node in surface 28
nodeset 11 group face_xpos
nodeset 11 name "face_xpos"

# ----------------------------------------------------------------------
# Create nodeset for +x face
# ----------------------------------------------------------------------
group "face_xneg" add node in surface 30
group "face_xneg" add node in surface 38
nodeset 12 group face_xneg
nodeset 12 name "face_xneg"

# ----------------------------------------------------------------------
# Create nodeset for +y face
# ----------------------------------------------------------------------
group "face_ypos" add node in surface 21
group "face_ypos" add node in surface 27
group "face_ypos" add node in surface 33
group "face_ypos" add node in surface 35
nodeset 13 group face_ypos
nodeset 13 name "face_ypos"

# ----------------------------------------------------------------------
# Create nodeset for -y face
# ----------------------------------------------------------------------
group "face_yneg" add node in surface 23
group "face_yneg" add node in surface 25
group "face_yneg" add node in surface 31
group "face_yneg" add node in surface 37
nodeset 14 group face_yneg
nodeset 14 name "face_yneg"

# ----------------------------------------------------------------------
# Create nodeset for -z face
# ----------------------------------------------------------------------
group "face_zneg" add node in surface 12
group "face_zneg" add node in surface 16
nodeset 15 group face_zneg
nodeset 15 name "face_zneg"

# ----------------------------------------------------------------------
# Create nodeset for -z face w/o fault
# ----------------------------------------------------------------------
group "face_zneg_nofault" add node in face_zneg
group "face_zneg_nofault" remove node in fault
nodeset 16 group face_zneg_nofault
nodeset 16 name "face_zneg_nofault"

# ----------------------------------------------------------------------
# Create nodeset for +z face
# ----------------------------------------------------------------------
group "face_zpos" add node in surface 10
group "face_zpos" add node in surface 17
nodeset 17 group face_zpos
nodeset 17 name "face_zpos"

5. Export the mesh file

At this point, we only need to exported the mesh file in exo format.

1
2
3
4

# ----------------------------------------------------------------------
# Export exodus file
# ----------------------------------------------------------------------
export mesh "mesh_1km.exo" dimension 3 overwrite

From the output of the Command Line you can see how many elements and nodes you have in the generated mesh file. Generally speaking, the more nodes, the longer the calculation time.

Reference
[1] Aagaard, B.T., Kientz, S., Knepley, M.G., Strand, L., Williams, C., 2017. PyLith User Manual: version 2.2.1. Davis: California: Computational Infrastructure of Geodynamics
[2] Tutorial of Trelis

1.1 Download

Click here to the download page
You can request a free trial on the official website for three months after the registration.

1.2 Installation

• For MacOS:
  Just double click the package and complete the installation
• For Linux

  1 2	cd ~/Software/ sudo dpkg -i Trelis-16.5.0-Lin64.deb

  Use the code you get from the trial request during Trelis Activation.

2 Calculation: Pylith

2.1 Download

• For MacOS:
  https://geodynamics.org/cig/software/pylith/git/v2.2.2/pylith-2.2.2-darwin-10.13.6.tgz
• For linux:
  https://geodynamics.org/cig/software/pylith/git/v2.2.2/pylith-2.2.2-linux-x86_64.tgz

2.2 Configuation

1. Unzip the file to where you want, for example

   1 2	tar zxvf pylith-2.2.1-linux-x86_64.tgz -C ~/Software/ # for macos it's pylith-2.2.2-darwin-10.13.6.tgz

2. Configuation

   1 2 3

   cd ~/Software/pylith-2.2.1-linux-x86_64 # for macos it's cd ~/Software/pylith-2.2.2-darwin-10.13.6] source setup.sh

   If the output showes up as

   1	Ready to run Pylith

   Then the pylith congiguation has completed. And everytime before running pylith you need run source setup.sh under corresponding direcotory. If you want to take it a easier way. You can add the following script into your ~/.bashrc

   1 2	# for pylith alias setpylith='tp=`pwd`; cd ~/Software/pylith-2.2.1-linux-x86_64 ; source setup.sh; cd $tp'

   Thus you just need to type setpylith for every new terminal window.

3 Visualization: Paraview

• For MacOS:

  1	brew cask install paraview

• For Linux:

  1	sudo apt install paraview

这是因为没有对图片进行栅格化,只需执行

1

set(gcf,'renderer','painters')

栅格化后再进行保存,输出为eps格式,便能得到无限放大不失真的矢量图了.

另外,若只是想生成高分辨率的png图片,可以执行

1

print -dpng FILENAME -r400

其中-r400是可指定的ppi值,越高则越清晰,生成的图片文件越大.

参考链接

Mainly commands:

• psxy - to plot fault geometry with fault slip
• pscoast - to plot basemap
• psvelo - to plot GPS velocity

Material needed：

1. Fault_geo1.dat: Fault geometry 1 with coseismic slip
2. Fault_geo2.dat: Fault geometry 2 with coseismic slip
3. GPS_coseis.dat: GPS data corresponding to the coseismic slip.
4. Plot_fault.ksh: Plotting script.

Material Format

The format for Fault_geo*.dat is

1
2
3
4
5
6

-Z ${slip magnitute}
${lon1} ${lat1}
${lon2} ${lat2}
${lon3} ${lat3}
${lon4} ${lat4}
...

The format for GPS_coseis.dat is

1

${lon}  ${lat} ${east_component} ${north_component} ${east_significant} ${noth_significant} 

Plotting script

Next we’ll mainly focus on our plotting script.

1. Give name to some initial parameters

1
2
3
4
5
6
7
8

#!/bin/bash
ps=plot.ps
fault_geo1=fault_geo1.dat
fault_geo2=fault_geo2.dat
gps_name=GPS_coseis.dat
cpt=fault_slip.cpt
R=-157/-144/55/62
J=M5i

2. Plot basic picture

1
2

gmt makecpt -Chot.cpt -G0.4/1 -Ic -T0/22 > $cpt
gmt pscoast -R$R -J$J -Bx4 -By2 -Da -G229/195/148 -S156/233/253 -K -P > $ps

Using makecpt to define the colormap: Choose hot colormaphen, cut the range to 0.41 with -G, invert colormap with -Ic, choose 022m’s slip range with -G, finally output to $cpt.(see makecpt)
Using pscoast to plot the land and ocean with different color.(see pscoast)

3. Plot fault geometry

1
2
3
4
5
6

gmt psxy $fault_geo1 -R$R -J$J  -W0.1p -L -C$cpt -O -K >> $ps
gmt psxy $fault_geo2 -R$R -J$J  -W0.1p -L -C$cpt -O -K >> $ps
gmt psscale -D3.5i/0.5i+w2.5i/0.25i+h+jCT+m -C$cpt \
        -Bx5+l"Coseismic slip (m)" -O -K -P >> $ps
gmt pscoast -R$R -J$J -Bx4 -By2 -Da -W229/195/148 -O -K -P >> $ps

Here we use psxy to plot the fault geometry, giving different color to fault patches according to slip maganitude and assigned colormap.(see psxy)
And we plot the coast line again to avoid everything being covered by fault geometry.

4. Plot GPS vectors

1
2
3
4
5
6
7
8

    awk '{print $1,$2,$3,$4,$5,$6,0.0}' ${gps_name} |\
    gmt psvelo    -R${R}  -W0.25p,black  -Gblack  -L  -Se0.3c/0.95/25  -J${J} -A0.3c+p1.5p+e -P -V -O -K >> $ps

    gmt psvelo <<END  -h2  -R${R}  -W0.25p,black  -Gblack  -L  -Se0.3c/0.95/12 -J${J} -A0.3c+p1.5p+e -P -V -O  >> $ps
         Long.   Lat.   Evel   Nvel   Esig   Nsig  CorEN SITE
          (deg)  (deg)    (m)        (m)
      -155.5      61.5     2.0    0      0     0     0     2m
END

We use psvelo to plot GPS displacement with the first two lines.(see psvelo)
The left 5 lines are used to plot displacement’s legend.

5. Convert to png

1
2

gmt psconvert $ps -A -Tg
rm $ps gmt.* $cpt

Finally, we can convert the .ps file to .png file for better Visualization and remove files we don’t need anymore.

The output picture should looks like this.

这通常是由于网格间距过小,程序在调用PETSC计算时,无法在指定次数内迭代收敛到给定的收敛阈值,可以

1. 在建立网格文件时,扩大网格间距(mesh size),直至迭代计算能顺利完成.
2. 调大迭代次数: 在pylithapp.cfg中找到Convergence paameters对应模块,调整ksp_max_it的数值(但不宜超过500,否则即使顺利完成计算,结果仍十分不稳定)
3. 调整收敛阈值(不建议): 同调大迭代次数类似,调整对应部分的ksp_rtol与ksp_atol即可.

Could not find final slip at (* * *) using spatial database Final slip

这是由于修改后的断层滑移分布与指定断层的linear slip function不匹配,给定的滑移量在同一材料范围内无法满足线性关系,一般可以

1. 修改断层滑移分布至满足线性分布,示例文件中一般放在spatialdb/finalslip.spatialdb类似的文件中
2. 不采用linear slip function,而使用默认的滑移分布求解方案,直接将对应*.cfg文件中指定slip.query_type的一行注释掉即可.

Mesh missing group of vertices ‘***’ for subdomain output

在修改范例的geometry文档时,对应边界条件surface的ID号可能也会改变,需要在mesh_***.jou或create_bc.jou中,找到add group对应行,重新指定边界surface的ID号.
当然更好的办法是直接在创建surface的时候给其起好名字(如name 'surface_example'),这样之后创建边界条件时,就可以指定其名字而非之后会变动的ID号了

Material id ‘*‘ for cell ‘*‘ does not match the id of any available materials or interfaces

这是由于读取的网格文件对应材料的block ID与*.cfg中确定问题时的material ID不匹配造成的,调整为相吻合的ID号即可

Surface *** has no facet

具体原因仍不明白,但个人摸索的办法是建立geometry时,对应surface的控制点坐标要尽量是整数,如50km,100km等,如果有找到好办法的欢迎评论告诉我.

• 适用于暗影精灵的4k屏线(~RMB120,可淘宝定制屏线)
• 兼容性较好的4k分辨率15.6寸屏幕–本文所用型号为京东方BOE06D7(~RMB350)
• 十字螺丝刀,镊子,翘片或信用卡

二.换屏流程

1 拆下机身底面

• 首先卸下底面所有螺丝
  
• 回到机身C面,使用卡片或薄片四周环边翘开卡扣
  
• 回到机身D面,从开合笔记本一侧抬起背板,便能卸下
  

2 断开相关接口

• 断开电源线,以免拆机过程出现短路
  
• 断开屏线:先向上撬开卡扣再抽出屏线
  
• 断开无线网卡天线,由于屏幕这边连着天线,故需要一并拆下,可使用小镊子拆下两个金属头
  

3 拆下屏幕框

• 由于暗影精灵系屏幕与框之间有双面胶固定,需要使用金属薄片或卡片从内侧探入,刮开双面胶(勿探入过深一面伤及屏线和原装屏四周的绝缘胶布).
  
• 刮开后用手探入掰开框四边的卡扣,可轻度使劲反复多次
  
• 拆下后见下图
  

4 卸下屏线

• 卸下图中固定屏幕四角的螺丝,沿轴翻到另一侧,以方便继续卸屏线
  
• 卸下连接摄像头的线头
  
• 卸下固定屏线的金属框
  
• 至此便能完全连带原屏线取下原屏幕

5 换新

• 取出新屏线和新屏幕,依次倒序执行步骤4和2
  

6 测试

• 完成后,别急着直接装回去,开机看看能否亮机
  
• 若无法亮机,检查所有的接口特别是屏线接口是否插好.

三.完成

至此大功告成,可以好好使用你的新屏幕了!本人在黑苹果下开启hidpi,分辨率调为1920*1080,分辨率甚是丝滑哈哈

说明

• 若只是想更换成更高色域的1080p屏幕,则完全不需要更换屏线,在步骤2中只需断开电源线,步骤4中卸下固定屏幕四角的螺丝并断开屏线接口,替换屏幕即可.
• 对于其他型号的笔记本更换屏幕,也可作为一定的参考,祝换屏成功!

• 全球热流数据是根据钻井,温泉,直接采样相关岩石计算导出的热流值.
• 数据来源于NOAA的geothermal数据.
• 下载数据视不同的路径格式可能为.gmt,.xyz格式
• 数据下载:http://www.geophysik.rwth-aachen.de/IHFC/heatflow.html

  程序

1. 初始参数设定

   1 2 3 4 5 6 7

   #!/bin/bash PS=plot/global_geothermal.ps in_gmt=data/Global.gmt in_xyz=data/global_heatflow.xyz R=-180/180/-90/90 J=X18d/9d gmt psbasemap -R$R -J$J -B30 -X2c -Y2c -BwSEN -K > $PS #设置底图参数

2. 将数据转为正常的.xyz格式,以便之后绘制散点图(由于数据在大洋区域较为分散,不适合做画平滑热流图,读者可自行尝试)

   1	gmt convert $in_gmt -a2=Heat_Flow > $in_xyz

3. 绘制全球海岸线

   1	gmt pscoast -J$J -B30 -R$R -A4 -Gwhite -W0.1p -Dhigh -Swhite -K -O >> $PS

4. 以散点形式绘制全球热流值

   1 2 3 4 5

   cat $in_xyz | awk '{if($3<50){print $1,$2}}'| gmt psxy -R$R -J$J -Sc0.07 -G117/0/251 -K -O >> $PS cat $in_xyz | awk '{if($3<100 && $3 >=50){print $1,$2}}'| gmt psxy -R$R -J$J -Sc0.07 -G0/182/255 -K -O >> $PS cat $in_xyz | awk '{if($3 >= 100 && $3 < 150){print $1,$2}}'| gmt psxy -R$R -J$J -Sc0.07 -G0/255/58 -K -O >> $PS cat $in_xyz | awk '{if($3 >= 150 && $3 < 200){print $1,$2}}'| gmt psxy -R$R -J$J -Sc0.07 -G255/201/0 -K -O >> $PS cat $in_xyz | awk '{if($3 >= 200){print $1,$2}}'| gmt psxy -R$R -J$J -Sc0.07 -G255/0/0 -K -O >> $PS

5. 绘制图例

   1 2 3 4 5 6 7 8 9 10 11 12 13 14

   gmt pslegend -R$R -J$J -DjTL+w6i/0.3i+o1c/8c -O <<EOF >> $PS # Legend for $PS # G is vertical gap, V is vertical line, N sets # of columns, D draws horizontal line. # H is header, L is label, S is symbol, T is paragraph text, M is map scale. # G 0.02i N 5 S 0.1i c 0.15 117/0/251 0.25p 0.3i 0-50 S 0.1i c 0.15 0/182/255 0.25p 0.3i 50-100 S 0.1i c 0.15 0/255/58 0.25p 0.3i 100-150 S 0.1i c 0.15 255/201/0 0.01p 0.3i 150-200 S 0.2i c 0.15 255/0/0 0.01p 0.3i >200mW/m2 G 0.02i EOF

6. 将成图转为png格式,并删除不必要的文件

   1 2	gmt psconvert $PS -A -Tg -P rm *.history $PS

成图示例

关键报错如下

1
2
3
4
5
6
7
8
9
10

-- The C compiler identification is unknown
CMake Error at CMakeLists.txt:13 (project):
  No CMAKE_C_COMPILER could be found.

  Tell CMake where to find the compiler by setting either the environment
  variable "CC" or the CMake cache entry CMAKE_C_COMPILER to the full path to
  the compiler, or to the compiler name if it is in the PATH.
...
...
make: *** No targets specified and no makefile found.  Stop.

提示找不到C/CXX compiler路径,搜了一下似乎常见于Visual studio问题,并找到方法多为以下1,2两种

解决方案

1. 安装build-essential

   1	$ sudo apt-get install build-essential

2. 安装g++

   1	$ sudo apt-get install g++

3. 手动添加路径
   若二者均已安装仍报错,那么找到gcc位置,手动添加下路径即可

   1 2 3

   $ whereis gcc gcc: /usr/bin/gcc /usr/lib/gcc $ export CC=/usr/bin/gcc # 视上一行输出填写路径

再次编译应该就能顺利完成了.

参考链接

环境:

os: Debian 9.3 stretch
git: 2.11.0
node: 8.9.3
npm: 5.8.0
hexo-cli:1.1.0

1. 关联github

1.1 设置git

1
2

$ git config --global user.name "username"
$ git config --global user.email "username@example.com"

设置git私钥

安装openssh,生成SSHkey

1
2

$ sudo apt install openssh
$ ssh-keygen -t rsa -b 4096 -C "your_email@example.com"

确保ssh-agent可运行

1

$ ssh-agent -s

添加SSH key

1

$ ssh-add ~/.ssh/id_rsa

拷贝 ~/.ssh/id_rsa.pub 中的内容，然后，在 GitHub 右上方点击头像，选择 ”Settings”，在右边的 ”Personal settings” 侧边栏选择 ”SSH Keys”。点击 ”Add key” 按钮添加，接着粘贴 key。

测试:

1

$ ssh -T git@github.com

2.安装node.js和hexo

2.1安装node.js

法1:
安装nvm:

1

$ curl -o- https://raw.githubusercontent.com/creationix/nvm/v0.33.2/install.sh | bash

重启Terminal,安装Node.js

1
2

$ nvm ls-remote
$ nvm install v7.2.0

法2:
包管理器安装

1
2

sudo apt install npm
sudo apt install nodejs

2.2 安装Hexo

1

$ sudo npm install -g hexo-cli      # 注：这里不用 sudo 会报错

3. 恢复hexo

3.1 初始化hexo

1

hexo init

3.2 拷贝所需文件

1. _config.yml
2. theme/
3. source/

3.3 安装必要组件

安装hexo-deployer-git,以使用hexo d部署到git上

1

npm install hexo-deployer-git --save

安装hexo-generator-feed,以建立RSS订阅

1

npm install hexo-generator-feed --save

安装hexo-generator-sitemap,以建立站点地图

1

npm install hexo-generator-sitemap --save

4. 验证

进行发布部署

1

hexo d -g

若成功就说明大功告成了!

参考链接

archlinux

1

sudo pacman -S paraview #需要事先添加archlinuxcn源

debian/ubuntu

1

sudo apt install paraview

目前使用包管理器安装时的版本分别为

• archlinux：paraview5.2.0
• debian/ubuntu：paraview5.1.2
  而当前最新版本的paraview为5.4.1，若想要安装最新的paraview，archlinux用户可以选择使用yaourt安装第三方aur包，或者使用源码编译

二.源码编译安装

在此以debian环境为例

1. 下载源码并安装相关依赖

• paraview5.4.1下载链接（linux）
• 安装相关依赖

  1 2	sudo apt install cmake qt5-default libxt-dev g++ gcc cmake-curses-gui mesa-common-dev ... openmpi-common openmpi-bin libopenmpi-dev

2. 使用CMake进行配置

建议建立一个单独的build路径，将source code拷贝到此处再进行配置。如在源码路径paraview下

1
2
3
4
5

cd ..
mkdir paraview_build 
cd paraview_build
#进行拷贝操作
ccmake ../paraview

将会出现以下界面

接下来

1. 选择合适的参数（上下方向键移动，enter修改）
2. 按c运行configure
3. 根据反馈信息调整参数或安装缺少的依赖，回到步骤1
4. 一切正常后，按g应用

具体各个参数的含义可见参考来源，其中对于报错缺少的文件，可使用

1

apt-file search $FILE#(如Qt5CoreConfig.cmake) 

来查找需要安装哪个包

3.编译

CMake结束后将会产生一个makefile，只需在当前目录执行make即可开始编译，视具体情况会花费几~几十分钟。可使用多线程进行加速，以四核处理器为例。

1

make -j 4

编译完成后你就已经可以通过

1

./bin/paraview

运行paraview了

4. 安装

接下来我们将paraview安装到本机环境中，只需执行

1

sudo make install

paraview便会自动安装到之前configure步骤CMAKE_INSTALL_PREFIX的指定目录。

至此安装完毕，对于其他linux发行版，windows下或macOS的安装步骤大同小异，在此不再赘述，可参考官方教程。

1.参考来源
2.各平台最新版paraview下载
3.paraview快速教程

1

-Bxinfo[/yinfo][:."title string":][W|w][E|e][S|s][N|n]

• 其中 xinfo 格式为：
• [a]tick[m|c][ftick[m|c]][gtick[m|c]][l|p][:"axis label":][:,"unit label":]
• 其中a,f,g分别为刻度标记，刻度和网格间距

-C 指定绘图拷贝数

缺省为1

-H 告知头记录数

缺省为1，若数据文件头大于1，如4，必须写-H4

-J 选择地图投影

• 线性投影 -JX width/height
  • -JXwidth[d] 当x,y用地理坐标的度表示时。
  • -JXwidthl 缩放前数值的log10
  • -JXwidthp power 缩放前数值的指数
• 莫卡托投影 -JMwidth
• 正形投影 -JGwidth 指定投影中心
• 圆锥投影 -JBwidth 指定投影中心和两条标准平行线
• Eckert投影 -JK[f|s]lon0/width 指定中央子午线

-R 指定绘图范围

-Rxmin/xmax/ymin/ymax

-Rxlleft/ylleft/xuright/yuright[r] 加r则是左下角和右上角坐标

-U 标注绘图日期

-U"optinal string"

-K 和 -O

规则

• 如果一张图只要一个命令就能画出来，则不需要-K -O
• 如果一张图要两个命令画出，则第一个-K 第二个 -O
• 如果一张图要多个命令画出，则第一个用-K 中间用-K-O，最后一个用-O
• 第一个命令使用重定向符号 >创建新文件并写入PS命令
• 余下所有绘图命令使用重定向符号>>向PS文件中追加新的PS命令。

避免犯错技巧

例

1
2
3
4
5
6
7
8
9
10
11
12
13
14

#!/bin/bash
#
# GMT template
#
PS=example.ps
J=X6i #表明是6英寸
R=0/1/0/1
# Begin GMT Plot
psxy -J$J -R$R -T -K > $PS
# now begin your real plot commands
# every command should be with "-K" "-O" and ">>"
# End GMT Plot
psxy -J$J -R$R -T -O >> $PS
# Plot finished. Now you can convert PS to other format