关于nc文件的读取.docVIP

精品文档 . 关于nc文件的读取 2012-11-02 21:47:01|??分类：? \o 信号处理DSP 信号处理DSP?|??标签： nc文件?? panoply??|举报|字号?订阅 最近在学习冰后回弹模型（http://www.atmosp.physics.utoronto.ca/~peltier/data.php），里面给出的数据为nc（netcdf）格式的，读取是个问题，在言深深的帮助下，成功读取，在这里将其总结出来，以备后用也希望可以方便有用之人。 有几种方式可以进行读取： 1、Excel的一个插件， NetCDF4Excel可以完成相关的查看工作，需要的可以进行下载。 安装之后，我们可以得到，一个快捷方式，如图： ?打开之后，最上面标题栏点击加载项，我们可以看到下面的情形： ?? 然后打开nc文件，我们就可以进行查看操作： ? 2、?通过Grads进行相关的读取，主要是在言深深同学的帮助下完成，深深的感谢！ 第一步，安装相关的 grads软件，这个这里不做介绍，相关可以参考 气象论坛网站； 第二步，打开grads软件，如图，回车之后的那个窗口千万不能关，那个窗口是图像显示的，关闭之后输入命令之后直接就退出了，刚开始我老犯这个低级错误（对菜鸟而言）。 ?? ??第三步，可以通过两个命令简单的看到数据的头文件包含信息： ①sdfopen d:\1.nc? ②q ctlinfo 1.nc文件 ? 然后编写gs文件，如下： reinit sdfopen D:\1.nc set gxout fwrite set fwrite D:\1.dat set x 1 360 set y 1 180 set z 1 set t 1 1 d dsea_250 disable fwrite ps1：路径自己可以按需更改，不要有空格。 ps2：这里的dsea_250可不是随便取的，看上面的ctlinfo心里里最后一行，必须要一致！ 然后就可以生成一个1.dat的文件。 第四步，通过fortran程序将1.dat读为txt文档 ? Program nc2txt ? ? parameter(m=180,n=360,o=1) ? ? dimension ? x(m,n,o) ? ? open(1,file=1.dat,access=direct,form=binary,recl=n*m*o*4) ? ? ? ? read(1,rec=1) (((x(i,j,k),j=1,n),i=1,m),k=1,o) ? ? close(1) ? ?? ? ? open(2,file=Ice-5gGeoid.txt) ? ? ? ?do i=1,m ? ? ? ? ? !print*,i ? ? ? ? ? do j=1,n ? ? ? ? ? ? ?write(2,(2f10.1,2x,f15.6)) -(i-1)*1.+180,(j)*1.,x(i,j,1) ? ? ? ? ? enddo ? ? ? ?enddo ? ? close(2) ? End program nc2txt ? ?? ps：这里的参数要与自己的数据一致，比如mno，以及最后的输出，有的是1°，有的是2.5°。 然后我们就可以得到相关的数据文件如下： ? 至此，数据读取操作完成，需要对有需要的人有帮助，再次对言深深表示感谢。 此外grads还可以通过那个窗口绘图，具体如下： （1）sdfopen d:\1.nc ? (2) ?d dsea_250 ? Nc2TxT（版权归作者所有，本着学习的精神） 3、一个国外学者对我的回复： In the meantime, may I suggest that you investigate the following ?free software: ncview and ncdump (for Linux) panoply?(for Windows)? /tools/panoply/ Personally,I use NCAR Graphics NCL programs to create plots from netcdf files. That is also free software available from NCAR (National ?Center for Atmospheric Research in Boulder, Colorado, USA) for a variety of platforms (Linux and Windows)。 /Download/install.shtml ?ncview and panoply are both viewers that d