fortran读netcdf文件求解[已经解决]

infocert

有个nc文件，头格式如下：

1. netcdf VEGET {
   
2. dimensions:
   
3.         lon = 64 ;
   
4.         lat = 32 ;
   
5.         time = UNLIMITED ; // (1992 currently)
   
6. variables:
   
7.         float lon(lon) ;
   
8.                 lon:long_name = "longitude coordinate" ;
   
9.                 lon:standard_name = "longitude" ;
   
10.                 lon:units = "degrees_east" ;
    
11.                 lon:axis = "X" ;
    
12.         float lat(lat) ;
    
13.                 lat:long_name = "latitude coordinate" ;
    
14.                 lat:standard_name = "latitude" ;
    
15.                 lat:units = "degrees_north" ;
    
16.                 lat:axis = "Y" ;
    
17.         float time(time) ;
    
18.                 time:long_name = "time" ;
    
19.                 time:standard_name = "time" ;
    
20.                 time:units = "years since 0001-01-01 00:00:00" ;
    
21.                 time:time_origin = "01-JAN-0001 00:00:00" ;
    
22.         float vegfrac(time, lat, lon) ;
    
23.                 vegfrac:long_name = "fraction de culture par maille" ;
    
24.                 vegfrac:standard_name = "fraction_culture" ;
    
25.                 vegfrac:units = "%" ;
    
26.                 vegfrac:missing_value = -99.99f ;
    
27. }
    

复制代码

我读取这个nc所用的fortran如下[已经修改为正确代码]：

1. PROGRAM crop_grille_ecbilt
   
2. 
   
3.   USE NETCDF
   
4.   implicit none
   
5. 
   
6.   !variable pour l'ouverture ecriture du netcdf
   
7.   character(256) :: filenameIn, filenameOut
   
8.   integer :: intputID, unlimDimID, varID, ID, nbTime, iyr
   
9. 
   
10.   !variable pour le calcul des moyennes, etc...
    
11.   real*4, dimension(:), allocatable :: lon, lat, time, ValueVar
    
12.   real*8, dimension(:,:,:), allocatable :: vegfrac
    
13. 
    
14.   !variables pour le changement de grille
    
15.   !real*8 farea(64,32,1)
    
16.   integer nlat,nlon
    
17. 
    
18.   integer :: outputID
    
19.   integer :: outdimid, outvarid
    
20. 
    
21.   if(iargc().eq.2) then
    
22.     call getarg(1,filenameIn)
    
23.     call getarg(2,filenameOut)
    
24.   else
    
25.     write(*,'(A)') "Usage: ./ConvertVegetYo [input.nc] [output.nc]"
    
26.     stop
    
27.   endif
    
28. 
    
29.   write(*,*) "Open ", trim(filenameIn)
    
30.   call check(nf90_open(filenameIn, nf90_nowrite, intputID))
    
31.   call check(nf90_inquire(intputID, unlimitedDimId = unlimDimID))
    
32. 
    
33.   write(*,*) "Read longitude"
    
34.   call check(nf90_inq_varid(intputID, "lon", ID))
    
35.   allocate(lon(64))
    
36.   call check(nf90_get_var(intputID, ID, lon, (/ 1 /), (/ 64 /)))
    
37. 
    
38.   write(*,*) "Read latitude"
    
39.   call check(nf90_inq_varid(intputID, "lat", ID))
    
40.   allocate(lat(32))
    
41.   call check(nf90_get_var(intputID, ID, lat, (/ 1 /), (/ 32 /)))
    
42. 
    
43.   call check(nf90_inq_varid(intputID, "time", ID))
    
44.   allocate(time(1992))
    
45.   call check(nf90_get_var(intputID, ID, time, (/ 1 /), (/ 1992 /)))
    
46. 
    
47.   call check(nf90_inq_varid(intputID, "vegfrac", varID))
    
48.   allocate(vegfrac(64,32,1992))
    
49.   call check(nf90_get_var(intputID, varID, vegfrac, (/ 1, 1, 1 /), (/ 64,32,1992 /)))
    
50.   
    
51.   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    
52.   !!  Create netcdf START
    
53.   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    
54.   !write target file
    
55.   write(*,*) "Create ", trim(filenameOut)
    
56.   call check(nf90_create(filenameOut, nf90_clobber, outputID))
    
57.   
    
58.   !define dimension
    
59.   call check(nf90_def_dim(outputID, "lon", 64, outdimid))
    
60.   call check(nf90_def_var(outputID, "lon", NF90_FLOAT, (/ 1 /), outvarid))
    
61.   call check(nf90_put_att(outputID, outvarid, "long_name", "longitude coordinate"))
    
62.   call check(nf90_put_att(outputID, outvarid, "standard_name", "longitude"))
    
63.   call check(nf90_put_att(outputID, outvarid, "units", "degrees_east"))
    
64.   call check(nf90_put_att(outputID, outvarid, "axis", "X"))
    
65. 
    
66.   call check(nf90_def_dim(outputID, "lat", 32, outdimid))
    
67.   call check(nf90_def_var(outputID, "lat", NF90_FLOAT, (/ 2 /), outvarid))
    
68.   call check(nf90_put_att(outputID, outvarid, "long_name", "latitude coordinate"))
    
69.   call check(nf90_put_att(outputID, outvarid, "standard_name", "latitude"))
    
70.   call check(nf90_put_att(outputID, outvarid, "units", "degrees_north"))
    
71.   call check(nf90_put_att(outputID, outvarid, "axis", "Y"))
    
72.   
    
73.   !NF90_UNLIMITED
    
74.   call check(nf90_def_dim(outputID, "time", NF90_UNLIMITED, outdimid))
    
75.   call check(nf90_def_var(outputID, "time", NF90_FLOAT, (/ 3 /), outvarid))
    
76.   call check(nf90_put_att(outputID, outvarid, "long_name", "time"))
    
77.   call check(nf90_put_att(outputID, outvarid, "standard_name", "time"))
    
78.   call check(nf90_put_att(outputID, outvarid, "units", "years since 0001-01-01 00:00:00"))
    
79.   call check(nf90_put_att(outputID, outvarid, "time_origin", "01-JAN-0001 00:00:00"))
    
80. 
    
81.   !define variable
    
82.   call check(nf90_def_var(outputID, "vegfrac", NF90_FLOAT, (/ 1, 2, 3 /), outvarid))
    
83.   call check(nf90_put_att(outputID, outvarid, "long_name", "fraction de culture par maille"))
    
84.   call check(nf90_put_att(outputID, outvarid, "standard_name", "fraction_culture"))
    
85.   call check(nf90_put_att(outputID, outvarid, "units", "%"))
    
86.   call check(nf90_put_att(outputID, outvarid, "missing_value", real(-99.99)))
    
87. 
    
88.   !finish the configuration of the output file
    
89.   call check(nf90_enddef(outputID))
    
90.   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    
91.   !! Create END
    
92.   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    
93. 
    
94.   do iyr=1,1992
    
95.     do nlat=1,32
    
96.       do nlon=1,64
    
97.         !farea(nlon,nlat,iyr)=vegfrac(nlon,nlat,1)
    
98.         write(*,*) nlon,",",nlat,",",iyr," => ", farea(nlon,nlat,iyr)
    
99.       end do
    
100.     end do
     
101.   end do
     
102.   
     
103.   call check(nf90_close(intputID))
     
104. 
     
105. END PROGRAM crop_grille_ecbilt
     
106. 
     
107. 
     
108. SUBROUTINE check(status)
     
109.   USE NETCDF
     
110.   IMPLICIT NONE
     
111. 
     
112.   INTEGER, INTENT (IN) :: status
     
113.   if(status /= nf90_noerr) then
     
114.     write(*,*)"Error  : ", trim(nf90_strerror(status))
     
115.     stop
     
116.   end if
     
117. END SUBROUTINE check
     

复制代码

我想问，这段代码到底对不对？编译成功，运行时，只打印出第一年的数据，第二年的数为空，然后程序就出错
出错信息为：

1. forrtl: severe (174): SIGSEGV, segmentation fault occurred
   
2. Image              PC                Routine            Line        Source
   
3. a.out              0000000000401740  Unknown               Unknown  Unknown
   
4. a.out              0000000000400E1C  Unknown               Unknown  Unknown
   
5. libc.so.6          00002B1417944CDD  Unknown               Unknown  Unknown
   
6. a.out              0000000000400D19  Unknown               Unknown  Unknown
   

复制代码

请论坛中的前辈帮忙解解，多谢了！！！

改变读数据方式后成功，读vegfrac时不在循环里读，直接一次性读出来

kongfeng0824

我一直就没安装上netcdf，不知道你有安装的软件可以传我一份吗？多谢了！

qxtlyf

第一年的数据处理的正确吗？？？

infocert

kongfeng0824 发表于 2013-3-2 01:46

                               
登录/注册后可看大图

我一直就没安装上netcdf，不知道你有安装的软件可以传我一份吗？多谢了！

我就是直接在大型机上安装的netcdf4.2/intel 12.1.4，你可以试着自己安装看看，不会的地方可以联系我

infocert

qxtlyf 发表于 2013-3-2 09:39

                               
登录/注册后可看大图

第一年的数据处理的正确吗？？？

没太仔细看，全数据确实是nc里的数据，不像第二年全是0.000000000

infocert

修改第53,54行为：

1. farea(nlon,nlat,iyr)=vegfrac(nlon,nlat,1)
   
2.         write(*,*) nlon,",",nlat,",",iyr," => ", farea(nlon,nlat,1)

复制代码

现在第二年的结果不再是0.000000了，但第二年只打印几十行，就出现1楼原贴的错误信息，求解，请各位前辈支招，小弟初学，非常感谢！！！
还有个问题，nlon,nlat,iyr这三个变量只是lon,lat,time数组的索引，farea(nlon,nlat,iyr)是不是应该变成farea(lon(nlon),lat(nlat),time(iyr))，同样vegfrac(nlon,nlat,1)是不是也应该改变？

火_旦

你好 请问你这个只是读取netcdf数据的文件吗？  有没有做过修改某些数值后 重新替换掉原来的netcdf文件？

infocert

火_旦 发表于 2013-11-26 14:20

                               
登录/注册后可看大图

你好 请问你这个只是读取netcdf数据的文件吗？  有没有做过修改某些数值后 重新替换掉原来的netcdf文件？

别用fortran做这些，太麻烦，用ncl比较快

火_旦

infocert 发表于 2013-11-27 12:41

                               
登录/注册后可看大图

别用fortran做这些，太麻烦，用ncl比较快

好的  谢谢您的建议  试一下NCL

凝露灵动

请问你的问题后来是怎么解决的？