Netcdf文件格式了解及如何创建、读取

godenflame135

转模式或做资料处理，只要涉及到将资料转换为模式所需要的格式或将模式输出文件进行处理，就会遇到文件格式的读写。而NETCDF格式是WRF所需要的，这里简单介绍下：
1.   一步生成读写源代码。在知道样本netcdf文件的情况,使用
ncdump –h   file1.nc  > file1.cdl  ！生成描述文件
ncgen   –f  file1.f  file1.cdl  ！生成fortran77源代码

则可以参考 file1.f  修改写出自己需要的代码。

2.     随意读取wrfout文件。WRF-ARW的UserGuide.pdf中在Tool一节介绍read_wrf_nc.f，想到这个文件的名字，相信我不多说，大家就知道里面应该有值得参考的句子。（当然，如果想看简单点的程序，请看netcdf软件自带的man中pdf文件介绍的examples（附件中））。

3.    附上f90方便快捷地创建netcdf文件和grads转netcdf格式文件的相关文献。

netcdf_mans.zip (1.36 MB, 下载次数: 741)

2011-12-6 00:20 上传

点击文件名下载附件
f90,f77,cpp,c,install,tutorial

read_wrf_nc.f (65.37 KB, 下载次数: 361)

2011-12-6 00:20 上传

点击文件名下载附件
WRF自带的读取nc格式wrfout的代码

实例分析如何用GRADS实现NetCDF格式转换.pdf (204.67 KB, 下载次数: 623)

2011-12-6 00:20 上传

点击文件名下载附件
文献

NetCDF数据的编程式访问.pdf (284.5 KB, 下载次数: 559)

2011-12-6 00:20 上传

点击文件名下载附件
文献

NetCDF格式的创建与应用.pdf (849.12 KB, 下载次数: 588)

2011-12-6 00:20 上传

点击文件名下载附件
文献

iapdzc

竟然要收费？？大家就别下这个了！拒绝收费！

寒江雪（王训

大气杨娃娃 发表于 2011-12-6 15:01

                               
登录/注册后可看大图

按你说的试了试，出现下列问题。。。
[root@localhost WRF2PBL]# pgf90 readwrf.f90 -I${NETCDF_INC} -L ...

如果是netcdf 4.0 以上，应该是lnercdff哈

冬日情愫+

先抢个沙发！版主分享了好东西了！

善人/jw

这么好的东西，怎么没有来顶啊

言深深

帖子很好，值得好好学习一下！
非常感谢@godenflame135的辛苦劳动，深深表示敬意

无径之林

好东西！感谢分享！

mofangbao

谢谢楼主的奉献  有大家才精彩

godenflame135

言深深 发表于 2011-12-6 08:20

                               
登录/注册后可看大图

帖子很好，值得好好学习一下！
非常感谢@godenflame135的辛苦劳动，深深表示敬意

坛主都如此认真负责，我们也来尽一点微薄之力。

大气杨娃娃

godenflame135 发表于 2011-12-6 12:47

                               
登录/注册后可看大图

坛主都如此认真负责，我们也来尽一点微薄之力。

通过ncdump和ncgen生成相应f文件后，pgf90该文件，发现找不到netcdf.inc。。。后来看网上说加上-I$NETCDF/include 语句，我试了却报错出现很多undefined reference。。。。
我netcdf安装好了，而且WRF编译运行也没有问题。
不解。。。

godenflame135

大气杨娃娃 发表于 2011-12-6 14:23

                               
登录/注册后可看大图

通过ncdump和ncgen生成相应f文件后，pgf90该文件，发现找不到netcdf.inc。。。后来看网上说加上-I$NETCDF ...

只加include是不足的。在环境变量中要设置一下：
export   NETCDF=/home/user/program_path/netcdf
export   NETCDF_INC=${NETCDF}/include
export   NETCDF_LIB=${NETCDF}/lib
export  PATH=${NETCDF}/bin:$PATH
正如大气洋娃娃说的，只加库是不够的，还要把目标文件（包含*.mod，*.o）的位置加上。
如： pgf90   file1.f   -I${NETCDF_INC}  -L${NETCDF_LIB}  -lnetcdf
这样编译就应该可以生成a.out 执行文件，只要你的代码正确无语法错误。


编外：现在fortran 90 的语句也有好多可以实现简捷编写，如数组赋初值可以  a（1：1000）=0.0  这样的语句来描述，风速质量控制where（a>50）a=-9999.0,这样的语句比起77来简便多了，另外netcdf函数fortran90版也比较简洁了，两句话用同一个函数可以写出一个变量，不像77那样还要分变量类型用不同的函数来写，所以修改生成的fortran 77语句时同时把use netcdf和include "netcdf.inc"加上就可以使用fortran90的新特征简捷高效写出自己想要的程序来。

大气杨娃娃

godenflame135 发表于 2011-12-6 14:39

                               
登录/注册后可看大图

只加include是不足的。在环境变量中要设置一下：
export   NETCDF=/home/user/program_path/netcdf
exp ...

按你说的试了试，出现下列问题。。。
[root@localhost WRF2PBL]# pgf90 readwrf.f90 -I${NETCDF_INC} -L${NETCDF_LIB} -Inetcdf
/tmp/pgf90SZ-ecfAxvHU.o: In function `get_wrf_var_':
readwrf.f90:(.text+0xfecc): undefined reference to `nf_open_'
readwrf.f90:(.text+0xffb2): undefined reference to `nf_inq_dimid_'
readwrf.f90:(.text+0xfff7): undefined reference to `nf_inq_dimlen_'
readwrf.f90:(.text+0x1029c): undefined reference to `nf_inq_varid_'
readwrf.f90:(.text+0x102e1): undefined reference to `nf_inq_varndims_'
readwrf.f90:(.text+0x10479): undefined reference to `nf_inq_vardimid_'
readwrf.f90:(.text+0x10502): undefined reference to `nf_inq_dimlen_'
readwrf.f90:(.text+0x106cf): undefined reference to `nf_get_var_text_'
readwrf.f90:(.text+0x109c6): undefined reference to `nf_get_var_real_'
readwrf.f90:(.text+0x10c8c): undefined reference to `nf_get_var_real_'
readwrf.f90:(.text+0x1103e): undefined reference to `nf_get_var_real_'
readwrf.f90:(.text+0x113f3): undefined reference to `nf_get_var_real_'
readwrf.f90:(.text+0x117a5): undefined reference to `nf_get_var_real_'
/tmp/pgf90SZ-ecfAxvHU.o:readwrf.f90:(.text+0x11d42): more undefined references to `nf_get_var_real_' follow
/tmp/pgf90SZ-ecfAxvHU.o: In function `get_wrf_var_':
readwrf.f90:(.text+0x1ab29): undefined reference to `nf_close_'
/tmp/pgf90SZ-ecfAxvHU.o: In function `get_wrf_var_check':
readwrf.f90:(.text+0x1cb8c): undefined reference to `nf_strerror_'


现在明白了 是-lnetcdf 不是-Inetcdf。。。。。。。。现在没有问题了