用Fortran读nc数据的程序

瞌睡猫喵

好东西，收藏先~

river

很好用的东西，谢谢

hyingniu

我怎么才能下载啊，这个文件对我很有用啊！

hyingniu

楼主，能否破例给发一个，急需啊！

hyingniu

我怎样才能知道自己的权限呢？还有，我什么条件才能下载啊？

stopfire

顶一个~~~~~~~~~~

米妮儿恩

我无法下载（自动回复：请不要使用迅雷等下载工具，使用IE右键另存为就可以），能不能给我发一下谢谢，qzsyanni@163.com

dolcevitak

同问：
按照楼主所说，第一步生成read_netcdf.f90后，第2步中每个变量会有以下一段程序代码（以相对湿度为例）
status = nf_inq_var(ncid,  22,dummy,xtype,ndim,dimids,natts)
           if ( status /= nf_noerr ) write (*,*) nf_strerror(status)
      do j = 1,ndim
      status = nf_inq_dim(ncid,dimids(j),dummy,len)
           if ( status /= nf_noerr ) write (*,*) nf_strerror(status)
      start(j) = 1 ; count(j) = len
      enddo
      status = nf_get_vara_real(ncid,  22,start,count,RH)
请教楼主这段的作用是什么，注解说是提取数据，但是我没看明白怎么个提取
还有如果我想把RH这个量单独提取出来，并转换成二进制文件，最后应该如何写程序

hyingniu

我的信箱是hyingniu@dlou.edu.cn，麻烦楼主，能不能给发一个呀？还有，我试着写了一个读取nc
f90程序，不知道为什么出错，请大侠给看看呗！

        !---Longitude(经度),Latitude(纬度)
        program read_ncdata
                implicit none
                !---读取某一层高度场数据
                include                'netcdf.inc'
               
!                real*4, parameter                ::        add_offset=117965.0, scale_factor=1.0                                !偏移量,放大系数
!                integer*4, parameter        ::        missing_value=32766, lats=73, lons=144, levs=17                        !缺测值,纬度数,经度数,高度场层数           
                integer*4                                ::        ncid,status,uwnd
                character*4                                ::        year_str
                character*200                        ::        in_file_str
!                character*200                        ::        out_file_str
                integer*4                                ::        year,days
                integer*4                                ::        i,j,l                !---
                real*4,allocatable                ::        uwind(:,:,:)        !动态数组，尚未分配控件,实传递到读数据子程序空间，读一年数据
!                integer*4                                ::        get_days_of_year
                external handle_err       
!                external nf_open       
!                external nf_inq_varid
!                external nf_get_var_real       
!                external nf_close       
!                days = get_days_of_year(year)
!                allocate(uwind(lons,lats,days))  !给动态数组分配控件，lon:0~357.5;lat:90~-90;lev:1~17;day:1~365
               
       
                !生成所要读取的nc文件的文件名
!                write(year_str,'(i4)')year
!                in_file_str = 'e:\work\data\slp\slp.'//year_str//'.nc'       
                in_file_str='D:\doctor\holland\data\uwnd.sig995.1949.nc'                       
                       
                !----------读netcdf数据------------------------------------------------
                status=nf_open(trim(in_file_str),nf_nowrite,ncid)   !打开nc
                if(status/=nf_noerr) call handle_err(status)
               
                status=nf_inq_varid(ncid,'uwnd',uwnd)        !获取变量信息
                if(status/=nf_noerr) call handle_err(status)
               
                status = nf_get_var_real(ncid,uwnd,uwind) !---读取变量
                if(status/=nf_noerr) call handle_err(status)
               
                status=nf_close(ncid)   ! 关闭nc
                if(status/=nf_noerr) call handle_err(status)

                !-------------------------------------------------------------------------
                !生成数据到数组
                do l = 1,1460 !---对天数循环
                        do j = 0,25 !---对纬度格数循环
                                do i = 100,125 !---对经度格数循环
                                        uwind(i,j,l)=uwind(i,j,l)                               
                                enddo
                        enddo
                enddo
               
        end

错误提示如下：

Deleting intermediate files and output files for project 'readnc - Win32 Debug'.
--------------------Configuration: readnc - Win32 Debug--------------------
Compiling Fortran...
D:\doctor\holland\readnc\read_nc_slp.f90
Linking...
read_nc_slp.obj : error LNK2001: unresolved external symbol _HANDLE_ERR@4
Debug/readnc.exe : fatal error LNK1120: 1 unresolved externals
Error executing link.exe.

readnc.exe - 2 error(s), 0 warning(s)
拜托哦！

wlzhongouc

hyingniu 发表于 2012-11-2 17:13

                               
登录/注册后可看大图

我的信箱是，麻烦楼主，能不能给发一个呀？还有，我试着写了一个读取nc
f90程序，不知道为什么出错，请大侠 ...

我只是用过这个读取nc数据的程序  自己没有写过读取nc数据的程序  要写还不容易啊