fortran使用netcdf 读写NC文件



先安装VS2010
再安装intel visual fortran2011

1、设置netcdf Path环境变量

去下载别人已编译好的netcdf库,http://bbs.lasg.ac.cn/bbs/viewthread.php?tid=955
NetCDF3.6.1Beta1.Win32.for.CVF.IVF.rar
为了简单操作,我将NetCDF3.6.1Beta1.Win32.for.CVF.IVF.rar里面文件结构根据说明重新调整了一下。

解压到C盘
C:\NetCDF3.6.1Beta1.Win32.for.IVF
得到3个文件夹:
bin,lib,include,开始设置visual studio。

将bin目录加入path
计算机——》右键,属性——》高级系统设置——》环境变量——》系统变量——》
在Path里面追加:C:\NetCDF3.6.1Beta1.Win32.for.IVF\bin;

注意修改path后,需要重启计算机或者先注销后登陆才会起作用。

2、设置Intel Visual Fortran Compiler环境

工具——》选项——》Intel Visual Fortran——》compilers
Libraries:
C:\NetCDF3.6.1Beta1.Win32.for.IVF\lib
Includes:
C:\NetCDF3.6.1Beta1.Win32.for.IVF\include

如图:



3、设置每个项目的依赖lib

项目上右键——》属性——》配置属性——》Linker——》General
Addtional Library Directories:
C:\NetCDF3.6.1Beta1.Win32.for.IVF\lib



4、读写文件Demo:

读NC文件:

program Console11
  use typeSizes
  use netcdf90
  implicit none
  include 'netcdf.inc'
    ! Variables
    character(len=50) filename
    integer ncid,ierr
    ! Body of Console3
    filename ='G:/si_monthly_1850-2012_r2.5x2.5.nc'
    print *, 'Hello World'
    
    ierr = nf90_open(filename,nf_nowrite,ncid)
    print *, ncid
end program Console11

写NC文件

! This program provides an elementary check of some of the parts of the 
!   Fortran 90 interface to netCDF 3.5. It is a Fortran 90 implementation
!   of the nctst.cpp program provided with the C++ interface to netcdf
!   (in the src/cxx directory of the netcdf distribution). 
!
program netcdfTest
  use typeSizes
  use netcdf90
  implicit none
  include 'netcdf.inc'
  
  ! netcdf related variables
  integer :: ncFileID,                                   &
             latDimID, lonDimID, frTimeDimID, timeDimID, &
             pressVarID, latVarID, lonVarID, frTimeVarID, refTimeVarID, scalarVarID
             
  ! Local variables
  integerparameter :: numLats = 4, numLons = 3, &
                        numFrTimes = 2, timeStringLen = 20
  character (len = *), parameter :: fileName = "example.nc"
  integer :: counter                      
  realdimension(numLons, numLats, numFrTimes) :: pressure
  integer (kind = FourByteInt), dimension(numFrTimes) :: frTimeVals
  real (kind = FourByteReal) fillVal;
  real (kind = FourByteReal), dimension(2) :: validRange;
  
  ! --------------------
  ! Code begins
  ! --------------------
  if(.not. byteSizesOK()) then
    print *, "Compiler does not appear to support required kinds of variables."
    stop
  end if
    
  ! Create the file
  call check(nf90_create(path = trim(fileName), cmode = nf90_clobber, ncid = ncFileID))
  
  ! Define the dimensions
  call check(nf90_def_dim(ncid = ncFileID, name = "lat",     len = numLats,        dimid = latDimID))
  call check(nf90_def_dim(ncid = ncFileID, name = "lon",     len = numLons,        dimid = lonDimID))
  call check(nf90_def_dim(ncid = ncFileID, name = "frtime",  len = nf90_unlimited, dimid = frTimeDimID))
  call check(nf90_def_dim(ncid = ncFileID, name = "timelen", len = timeStringLen,  dimid = timeDimID))
  ! Create variables and attributes
  call check(nf90_def_var(ncid = ncFileID, name = "P", xtype = nf90_float,     &
                     dimids = (/ lonDimID, latDimID, frTimeDimID /), varID = pressVarID) )
  call check(nf90_put_att(ncFileID, pressVarID, "long_name",   "pressure at maximum wind"))
  call check(nf90_put_att(ncFileID, pressVarID, "units",       "hectopascals") )
  ! Use 4-byte reals explicitly, to match 4-byte attribute type in test file
  validRange(1) = 0.
  validRange(2) = 1500
  call check(nf90_put_att(ncFileID, pressVarID, "valid_range", validRange))
  ! Use a 4-byte float constant, to match variable type
  fillVal = -9999.0
  call check(nf90_put_att(ncFileID, pressVarID,  "_FillValue", fillVal ) )
                      
  call check(nf90_def_var(ncFileID, "lat", nf90_float, dimids = latDimID, varID = latVarID) )
  call check(nf90_put_att(ncFileID, latVarID, "long_name""latitude"))
  call check(nf90_put_att(ncFileID, latVarID, "units""degrees_north"))
  call check(nf90_def_var(ncFileID, "lon", nf90_float, lonDimID, lonVarID) )
  call check(nf90_put_att(ncFileID, lonVarID, "long_name""longitude"))
  call check(nf90_put_att(ncFileID, lonVarID, "units",     "degrees_east"))
  call check(nf90_def_var(ncFileID, "frtime", nf90_int, frTimeDimID, frTimeVarID) )
  call check(nf90_put_att(ncFileID, frTimeVarID, "long_name""forecast time"))
  call check(nf90_put_att(ncFileID, frTimeVarID, "units",     "hours"))
  call check(nf90_def_var(ncFileID, "reftime", nf90_char, timeDimID, refTimeVarID) )
  call check(nf90_put_att(ncFileID, refTimeVarID, "long_name""reference time"))
  call check(nf90_put_att(ncFileID, refTimeVarID, "units",     "text_time"))
                     
  ! In the C++ interface the define a scalar variable - do we know how to do this? 
  call check(nf90_def_var(ncFileID, "ScalarVariable", nf90_real, scalarVarID))
  
  ! Global attributes
  call check(nf90_put_att(ncFileID, nf90_global, "history", &
                     "created by Unidata LDM from NPS broadcast"))
  call check(nf90_put_att(ncFileID, nf90_global, "title", &
                     "NMC Global Product Set: Pressure at Maximum Wind"))
  
  ! Leave define mode
  call check(nf90_enddef(ncfileID))
  
  ! Write the dimension variables
  call check(nf90_put_var(ncFileID, latVarId,     (/ -90., -87.5, -85., -82.5 /)) )
  call check(nf90_put_var(ncFileID, lonVarId,     (/ -180, -175, -170 /)      ) )
  ! Don't use anonymous array here, in case platform has 8-byte integers
  frTimeVals(1) = 12
  frTimeVals(2) = 18
  call check(nf90_put_var(ncFileID, frTimeVarId,  frTimeVals                  ) )
  call check(nf90_put_var(ncFileID, reftimeVarID, "1992-3-21 12:00"           ) )
  
  ! Write the pressure variable. Write a slab at a time to check incrementing.
  pressure = 949. + real(reshape( (/ (counter, counter = 1, numLats * numLons * numFrTimes) /),  &
                                    (/ numLons, numLats, numFrTimes /) ) )
  call check(nf90_put_var(ncFileID, pressVarID, pressure(:, :, 1:1)) )
  call check(nf90_put_var(ncFileID, pressVarID, pressure(:, :, 2:2), start = (/ 1, 1, 2 /)) )
  
  call check(nf90_put_var(ncFileID, scalarVarID, 10))
  call check(nf90_close(ncFileID))
contains
  ! Internal subroutine - checks error status after each netcdf, prints out text message each time
  !   an error code is returned. 
  subroutine check(status)
    integerintent ( in) :: status
    
    if(status /= nf90_noerr) then 
      print *, trim(nf90_strerror(status))
    end if
  end subroutine check  
end program netcdfTest


配置中出现问题

问题1:

1>LIBCMT.lib(invarg.obj) : error LNK2005: __initp_misc_invarg 已ò?经-在ú LIBCMTD.lib(invarg.obj) 中D定¨义ò?
1>LIBCMT.lib(invarg.obj) : error LNK2005: __call_reportfault 已ò?经-在ú LIBCMTD.lib(invarg.obj) 中D定¨义ò?
1>LIBCMT.lib(invarg.obj) : error LNK2005: __set_invalid_parameter_handler 已ò?经-在ú LIBCMTD.lib(invarg.obj) 中D定¨义ò?
1>LIBCMT.lib(invarg.obj) : error LNK2005: __get_invalid_parameter_handler 已ò?经-在ú LIBCMTD.lib(invarg.obj) 中D定¨义ò?
1>LIBCMT.lib(invarg.obj) : error LNK2005: __invoke_watson 已ò?经-在ú LIBCMTD.lib(invarg.obj) 中D定¨义ò?
1>LIBCMT.lib(invarg.obj) : error LNK2005: "void __cdecl _invoke_watson(unsigned short const *,unsigned short const *,unsigned short const *,unsigned int,unsigned int)" (?_invoke_watson@@YAXPBG00II@Z) 已ò?经-在ú LIBCMTD.lib(invarg.obj) 中D定¨义ò?
1>LIBCMT.lib(invarg.obj) : error LNK2005: __invalid_parameter 已ò?经-在ú LIBCMTD.lib(invarg.obj) 中D定¨义ò?
1>LIBCMT.lib(invarg.obj) : error LNK2005: "void __cdecl _invalid_parameter(unsigned short const *,unsigned short const *,unsigned short const *,unsigned int,unsigned int)" (?_invalid_parameter@@YAXPBG00II@Z) 已ò?经-在ú LIBCMTD.lib(invarg.obj) 中D定¨义ò?
1>LIBCMT.lib(invarg.obj) : error LNK2005: ___pInvalidArgHandler 已ò?经-在ú LIBCMTD.lib(invarg.obj) 中D定¨义ò?
1>LINK : warning LNK4098: 默?认è?库a“°LIBCMT”±与ó?其?他?库a的μ?使ê1用ó?冲3?突í?;£?请?使ê1用ó? /NODEFAULTLIB:library
1>Debug\Console1.exe : fatal error LNK1169: 找ò到μ?一ò?个?或ò多à个?多à重?定¨义ò?的μ?符·?号o?

解决办法:
项目上右键——》属性——》配置属性——》fortran——》General——》Input
Ignore specific library:LIBCMT.lib


问题2:

出现 "无法定位程序输入点_NF_ABORT@4于动态链接库netcdf.dll"
或 "无法启动此程序,因为计算机中丢失netcdf.dll ..."




解决办法:
将C:\NetCDF3.6.1Beta1.Win32.for.IVF\bin加入系统环境变量Path中
或者C:\NetCDF3.6.1Beta1.Win32.for.IVF\bin\netcdf.dll拷贝到工程的Debug目录下

注意将bin路径加入path后,需要重启计算机或者先注销后登陆才会起作用。

参考:

阅读更多
个人分类: fortran
想对作者说点什么? 我来说一句

netcdf4.3.3.1安装

2016年02月28日 2.94MB 下载

没有更多推荐了,返回首页

加入CSDN,享受更精准的内容推荐,与500万程序员共同成长!
关闭
关闭