发一个简单的fortran读取nc文件的程序

半个瞎子

         上次要用fortran读取nc数据，好像在论坛里没有找到相关的帖子，所以就想着发一个简单的读取程序，希望能帮助到像我一样的菜鸟。我就直接贴上来了。

• Subroutine read_rain_nc(rfile,lat,lon,RAIN)
• implicit none
• include 'app/netcdf/include/netcdf.inc'            ！这里用到了netcdf库
• integer,parameter :: nx=720,ny=361,nt=4
• integer ncid,ierr,varidx,varidy,varid,err
• character(len=50) rfile
• real*4 lat(ny),lon(nx),RAIN(nx,ny,nt)
• ierr=NF_OPEN(trim(rfile),NF_NOWRITE,ncid)
• ierr=NF_INQ_VARID(ncid,'longitude',varidx)
• ierr=NF_INQ_VARID(ncid,'latitude',varidy)
• ierr=NF_GET_VAR_REAL(ncid,varidx,lon)
• ierr=NF_GET_VAR_REAL(ncid,varidy,lat)
• ierr=NF_INQ_VARID(ncid,'APCP_surface',varid)
• ierr=NF_GET_VAR_REAL(ncid,varid,RAIN)
• ! print *,RAIN(1,1,1),RAIN(720,361,4)
• err=NF_CLOSE(ncid)
• end subroutine read_rain_nc
  

来自群组: 成于大气 信达天下

数值分析

可不可有介绍这些函数汉语的呀

andrewsoong

VERY GOOD ！

lovechang1314

谢谢分享
订楼主

letsdo8

这个子程序有改进的地方，nx、ny、和nt应该放进子程序哑元中，在子程序外面（调用前）给值。这样通用性更强。

半个瞎子

letsdo8 发表于 2014-5-21 17:10
这个子程序有改进的地方，nx、ny、和nt应该放进子程序哑元中，在子程序外面（调用前）给值。这样通用性更强 ...

你是说的allocate数组么？初学者，还需要一步步改进，谢谢你的意见。

任我行

初学者，还需要一步步改进，谢谢你的意见。

问题少女

歇息诶分享啦~不过我一直使用IDL来读.nc和.hdf的数据的，fortran没有试过，嘿嘿，有空试一下咯

黑胡子777

学习学习。。

mymaple

感激不尽，非常感谢楼主的贡献。。。

limoxi

请问这是在cvf中编译的吗？ivf是不是不识别这些命令？