气象家园投稿第四篇：编写程序
说起编程，平时我们都是写什么样的程序呢？算个平均、求个趋势，再检验下显著性？或者进行小波分析、EOF 分析？不知道大家的选择如何，反正我平时大多是靠 Matlab 来完成的，感觉二爷教的东西更有用啊喂！当然也有需要用到 Fortran 进行计算的，但谁会跑到 Linux 上写代码啊，在 Windows 下方便多了嘛！
不过呢，气象上用到的很多模式都需要在 Linux 的环境下运行，了解一些基本的 Linux 编程还是很有用的。模式的使用并不需要有很深的计算机知识，通常按照说明敲几行命令，说不定就 run 起来了。不过很多情况下还需要根据你的系统环境进行配置，有时甚至出现莫名其妙的问题而运行不了，这时候若是能对 Linux 编程略知一二，可能就帮了你的大忙了。这一期我会介绍一些有关 Linux 编程的基本知识，也稍微提一下 Makefile，但更多的内容还需要自己深入学习，相信大家不用多久就能成为模式小能手呢！那么第一步就从 hello world 开始吧！
开始编译吧！
在 Linux 下，我们用 vi 照样写代码，但写好的代码如何编译呢？Windwos 下只需要一个按钮就行，Linux 下可以用一行命令来实现。hello world 的代码我就不贴了，但编译命令是这样的：
1[N@Dell hello]$ ifort -o hello hello.f90
这里我的 Linux 装的是 Intel 的 Fortran 编译器，如果你想在你的 Linux 上编译 hello.f90，得看看你那的 Fortran 编译器是啥，然后把这里的 ifort 改成你的编译器的名称就行了。-o hello 会指定可执行文件名为 hello，你也可以改成 -o myhello”，而后面跟的hello.f90则是将要编译的源码文件。若是不指定-o，写成“ifort hello.f90，就会得到名为 a.out 的可执行文件。或许你的编译器是 gfortran？那就输入：
1[N@Dell hello]$ gfortran -o hello hello.f90
除了 -o 选项，还有 -O、-O2、-O3，这些“杠大欧”们会对你的程序进行优化，貌似很多数值模式都考虑把它给加上，以提高运算效率。还有更多的参数，若是需要用到就 man 一下吧。如果你没有 Fortran 编译器(很有可能是这样)，那么你可以去 Intel 官网下一个，他们有提供教育版的，安装过程也很方便。对了，我记得气象家园上有分享过来着。
加入函数库
为什么需要加入函数库呢？简而言之，外部函数库的功能就是帮助你完成一些 Fortran 本身无法完成的工作，比如读写 nc 文件、进行并行计算等等。通常下载安装完一个函数库以后，会得到一些库文件和头文件分别存放在 lib 和 include 文件夹下，以 netcdf 库(读写 nc 文件的库)为例：
1
2
3
4
5
6
7
8
9
10[N@Dell ~]$ ls /opt/netcdf/netcdf-latest/lib/
libnetcdf.a libnetcdff.so libnetcdf.la libnetcdf.so.11
libnetcdff.a libnetcdff.so.6 libnetcdf.settings libnetcdf.so.11.0.0
libnetcdff.la libnetcdff.so.6.0.1 libnetcdf.so pkgconfig
[N@Dell ~]$ ls /opt/netcdf/netcdf-latest/include/
netcdf4_f03.mod netcdf.h netcdf_nc_data.mod
netcdf4_nc_interfaces.mod netcdf.inc netcdf_nc_interfaces.mod
netcdf4_nf_interfaces.mod netcdf_mem.h netcdf_nf_data.mod
netcdf_f03.mod netcdf_meta.h netcdf_nf_interfaces.mod
netcdf_fortv2_c_interfaces.mod netcdf.mod typesizes.mod
那么如何使用 netcdf 库呢？首先需要在源代码里面 use netcdf，然后调用 netcdf 库函数进行读取。而在编译的时候则需要把库的名称和路径告诉编译器，即编译时加上几个参数。这里提供一个示例程序：
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
29program main
use netcdf
implicit none
character(50) :: ncfile
integer :: fid, varid, stat, year, recnum, start_num
integer :: i
real(kind=8), dimension(200) :: ice
real(kind=8), dimension(180) :: lat
real(kind=8), dimension(360) :: lon
!初始化
ncfile=’HadISST_ice.nc’
year=1979
recnum=12
start_num=(year-1870)*12+1
!打开文件，读取数据
stat=nf90_open(ncfile, NF90_NOWRITE, fid)
stat=nf90_inq_varid(fid, ‘sic’, varid)
stat=nf90_get_var(fid, varid, ice, start=(/180,10,start_num/), count=(/1,1,recnum/))
stat=nf90_inq_varid(fid, ‘latitude’, varid)
stat=nf90_get_var(fid, varid, lat)
stat=nf90_inq_varid(fid, ‘longitude’, varid)
stat=nf90_get_var(fid, varid, lon)
stat=nf90_close(fid)
!输出到屏幕
write(*,*)’ lat= ‘, lat(10), ‘ lon= ‘, lon(180)
do i=1,recnum
write(*,*)ice(i)
enddo
end program main
源码就不作过多说明了，我们来看看编译命令及运行结果：
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15[N@Dell readnc]$ ifort -o readnc readnc.f90 -lnetcdff -I/opt/netcdf/netcdf-latest/include -L/opt/netcdf/netcdf-latest/lib
[N@Dell readnc]$ ./readnc
lat= 80.5000000000000 lon= -0.500000000000000
0.920000016689301
0.959999978542328
0.850000023841858
0.959999978542328
0.870000004768372
0.910000026226044
0.680000007152557
0.589999973773956
0.259999990463257
0.560000002384186
0.720000028610229
0.750000000000000
编译的命令和之前相比多了三项。-lnetcdff 可以看成是 -l 和 netcdff 连在一起，其中的 -l 是 lib 的意思，netcdff 为库名，表示你将用到的库文件名为 libnetcdff.so 或 libnetcdff.a。库文件的位置由 -L 后的路径给出，头文件则从 -I 后的路径里去找。使用其他库时，也是类似地加上 -lxxx -Ixxx -Lxxx，然后就能正确调用外部函数了。当然，你可以加上多个 -l、-I 以及 -L。也许你也想读个 nc 文件试试，却发现没有 netcdf 库？等下载好了却装不成功？那就百度一下吧，或者参考下一期的内容。
.so / .a / .o文件
有人可能会问了，xxx.so 文件和 xxx.a 文件有啥区别啊，其实前者是动态函数库，而后者是静态函数库。两者的区别之一在于，是否需要从库文件中读取函数信息。对于静态函数库，它在编译的时候就直接整合到可执行文件中了，程序可以独立运行的，不需要再读 xxx.a 的内容。若是使用了动态函数库，可执行文件运行时还得去 xxx.so 里面读取函数。形象点描述的话，静态函数库是揣兜里的，随时可以使用，而动态函数是人家的，要用的时候还得找人家拿来看看。
也许你就要问了，前面使用 netcdf 库编译时加的是 -lnetcdff，可我们并不知道它到底用的 libnetcdff.so 还是 libnetcdff.a 啊？要知道，在 /opt/netcdf/netcdf-latest/lib 目录下，这两个文件可都是存在的！实际上我们可以不使用 -l 和 -L 的组合，而是这样写：
1
2[N@Dell readnc]$ ifort -o readnc readnc.f90 /opt/netcdf/netcdf-latest/lib/libnetcdff.so -I/opt/netcdf/netcdf-latest/include
[N@Dell readnc]$ ifort -o readnc readnc.f90 /opt/netcdf/netcdf-latest/lib/libnetcdff.a -I/opt/netcdf/netcdf-latest/include > ifort.log 2>&1
可以看到使用 libnetcdff.so 时，编译没有问题。而使用 libnetcdff.a 时则会出错，这里我将标准输出和标准错误输出都重定向了，不了解重定向没有关系，总之就是我把编译错误的信息放到 ifort.log 文件里了。当然我们不用查看 ifort.log，只需知道使用 libnetcdff.a 时确实出错了，也说明之前采用 -l 与 -L 组合时，实际是用了 libnetcdff.so 文件。至于为什么不能用静态库编译呢？大概是还需要其他函数库，而 libnetcdff.so 里有它们的信息但 libnetcdff.a 里面却没有吧。
那么我们来看另一个有关动态库的问题。还是以 netcdf 库为例，之前我们能够正常编译运行，但现在我把 .bashrc 文件里的某行注释掉(图中 netcdf 部分的 export LD_LIBRARY_PATH=… 那行)：
那么当我重新登录的时候，就会显示这样的结果：
1
2
3[N@Dell readnc]$ ifort -o readnc readnc.f90 -lnetcdff -I/opt/netcdf/netcdf-latest/include -L/opt/netcdf/netcdf-latest/lib
[N@Dell readnc]$ ./readnc
./readnc: error while loading shared libraries: libnetcdff.so.6: cannot open shared object file: No such file or directory
为什么呢？编译没错，却不能运行？那是因为我们编译时用的动态库，运行时就得再次访问 libnetcdff.so 文件，而程序需要由 LD_LIBRARY_PATH 这个环境变量提供动态库文件所在的位置。这样就懂了吧？注释掉那行就意味着程序再也找不着 libnetcdff.so 了。
动态库和静态库有各自的优点，动态库用起来方便一点，静态库独立一些。前面说过用 libnetcdff.a 编译时，因为缺少其他函数库的信息失败了，而 libnetcdff.so 则含有那些函数库的信息，ldd libnetcdff.so 就能看到。另一方面，动态库升级后不需要重新编译可执行文件，而静态库升级后需要重新编译才能使用新函数。不过，使用静态库的好处在于，编译得到的程序可以独立执行，不再需要访问库文件了。
这里我再稍微介绍下怎样生成 .a 及 .so 文件。不过先要提一下 .o 文件，即目标文件，相信大家都见过。以 hello world 为例，如果编译命令这样写，就会生成 .o 文件：
1
2
3
4[N@Dell hello]$ ifort -c hello.f90
[N@Dell hello]$ ifort -o hello hello.o
[N@Dell hello]$ ls
hello hello.f90 hello.o
反而多敲了一行，为什么一定要生成 .o 文件呢？其实静态库就是将很多 .o 和在一起形成的，动态库也差不多，有了库文件就能多次使用而不用每次都编译了。那如何形成库文件呢？这里我用另一个程序来演示一下，和 hello world 不同，它有三个源码文件。先不管我这个程序是做什么用的，只看看我是如何编译的吧：
1
2
3
4
5
6
7
8
9
10[N@Dell test]$ ls
main.c main.h solve_x.c solve_x.h solve_y.c solve_y.h
[N@Dell test]$ mpicc -c solve_x.c
[N@Dell test]$ mpicc -c solve_y.c
[N@Dell test]$ ar -r libsolve.a solve_x.o solve_y.o
ar: creating libsolve.a
[N@Dell test]$ mpicc -o solve main.c libsolve.a
[N@Dell test]$ mpirun -n 2 ./solve
Proc 0: x=1.994085,ex=0.000035,step=11
Proc 1: y=1.999241,ey=0.000002,step=11
这个程序是用 c 语言写的，还用到了并行，因此编译器是 mpicc，且运行时需要写成 mpirun -n 2 ./solve，但它跟普通的串行程序看起来区别不大。可以看到我先编译了 solve_x.c 和 solve_y.c，得到目标文件后把它们合成一个 libsolve.a，建立了自己的一个函数库。然后编译 main.c 时直接就用这个静态库了，最后执行 solve。如果考虑制作动态函数库，可以这样：
1
2
3
4
5
6
7
8[N@Dell test]$ ls
main.c main.h solve_x.c solve_x.h solve_y.c solve_y.h
[N@Dell test]$ mpicc -fPIC -shared -o libsolve.so solve_x.c solve_y.c
[N@Dell test]$ mpicc -o solve main.c libsolve.so
[N@Dell test]$ export LD_LIBRARY_PATH=$LD_LIBRARY_PATH:.
[N@Dell test]$ mpirun -n 2 ./solve
Proc 0: x=1.994085,ex=0.000035,step=11
Proc 1: y=1.999241,ey=0.000002,step=11
其过程和静态库生成差不太多，只是出现了 -fPIC 和 -shared 选项，同时还需要设置 LD_LIBRARY_PATH。相信大家百度一下 man 一下大致就懂，这里不再赘述。
如何调试？
Windows 下程序调试大家都会，编译器自带该功能。而 Linux 下可以用 GDB 来进行调试：
1
2[N@Dell readnc]$ ifort -g -o readnc readnc.f90 -lnetcdff -I/opt/netcdf/netcdf-latest/include -L/opt/netcdf/netcdf-latest/lib
[N@Dell readnc]$ gdb readnc
你会发现编译的时候多加了一个参数 -g，然后你才能用 gdb readnc 进行调试。说实话我平时在 Linux 下几乎不用 GDB，但也许你会需要它，所以这里还是列一下 GDB 的基本操作：
file [filename]：打开可执行文件 filename
info：查看相关信息，如 info breakpoints 查看所有断点
list：列出可执行文件对应的代码
break：设置断点，一般按源文件行号来设置，如 break 5
delete：按断点号删除断点
clear：按行号删除断点
run：运行可执行文件
continue：从断点处继续执行
next：单步执行
print：显示变量或表达式的值
whatis：查看变量类型
kill：终止正在调试的程序
quit：退出 GDB
大致就是这些，我就不演示了，用一用就熟悉啦！
Makefile
Makefile 是啥？为什么要使用它呢？想象一下你有个程序是由100个源文件编译链接而成，那我岂不是要敲很长的命令？若要重新编译，我还得再敲一遍？当然不是，我们可以把编译信息写到一个文件里，即 Makefile，然后敲入命令 make 就行了！而且当你改动了一些源码文件后再次编译时，make 还能够自动识别修改了哪些文件并只编译那些文件，这样可以节省不少时间。还是以 solve 程序为例来演示下 Makefile 的基本规则吧，我写的 Makefile 文件内容如下：
1
2
3
4
5
6
7
8
9
10
11
12CC = mpicc
solve: main.o solve_x.o solve_y.o
${CC} -o solve main.o solve_x.o solve_y.o
main.o: main.c
${CC} -c main.c
solve_x.o: solve_x.c
${CC} -c solve_x.c
solve_y.o: solve_y.c
${CC} -c solve_y.c
clean:
rm -f solve main.o solve_x.o solve_y.o
其中第一行为定义变量，这里我把编译器定为 mpicc，如果想换另一个编译器，只需要修改第一行。后面的内容则以这样的形式书写：
1
2target: object files
[tab] command
其中 target 为目标文件或标签，后面的 object files 为完成该 target 所需要用到的文件，两者以冒号隔开。第二行以 tab 键开头，后面紧跟完成该 target 的具体命令。这里共有5个 target，前四个分别生成 solve、main.o、solve_x.o、solve_y.o，最后一个是把这四个文件删除。如果你要执行某一条 target，可以敲入 make target，比如这样的操作：
1
2
3
4
5
6
7
8
9
10
11[N@Dell test]$ ls
libsolve.a main.c Makefile solve_x.h solve_y.h
libsolve.so main.h solve_x.c solve_y.c
[N@Dell test]$ make solve_x.o
mpicc -c solve_x.c
[N@Dell test]$ make solve
mpicc -c main.c
mpicc -c solve_y.c
mpicc -o solve main.o solve_x.o solve_y.o
[N@Dell test]$ make clean
rm -f solve main.o solve_x.o solve_y.o
当我们需要实现 solve_x.o 这个 target 的时候，它将 solve_x.c 编译成了 solve_x.o。实现 solve 这个 target 时，需要用到 main.o、solve_x.o、solve_y.o 三个文件，而 main.o 和 solve_y.o 并不存在，此时 make 就会先编译得到 main.o 与 solve_y.o，然后再得到 solve。如果你直接敲 make，那就会自动完成第一个目标。如果你要 make 的文件不叫 Makefile，而叫 MyMakefile，那就敲 make –f MyMakefile。
其实随便找一个模式看看，就发现它们基本都是采用 make(或者 gmake？)进行编译，因为源文件实在是太多了。不过模式的 Makefile 是很复杂的，我也看不太明白，但我知道它大致在干嘛，然后改一下关键的地方，嗯，编译通过啦！
结语
Linux 下写程序对于习惯了 Windows 的我们来说确实需要一段时间适应，适应了就一切都手到擒来。不过 Linux 编程的知识多了去了，我也有很多不懂的地方，真的是需要不断地学习啊！好了，这期就这样啦，我们下期见！