【已解决】用fortran读取grd文件时打不开。太长了记得要续行！！！

CKDDOUBI

如图，原数据是nc文件，楼主想用openGrADS的fwrite功能将数据转成grd文件进行处理。


楼主使用以下gs代码转成grd

1. 'reinit'
   
2. 'sdfopen C:\Users\CKDDOUBI\Desktop\Fortran\Meteorological_Statistics\ncepII.uwnd.mon.mean.nc'
   
3. 'set lev 200'
   
4. 'set time jan1981 jan2010'
   
5. 'set x 1 144'
   
6. 'set y 1 73'
   
7. 'set fwrite C:\Users\CKDDOUBI\Desktop\Fortran\Meteorological_Statistics\2\h200.grd'
   
8. 'set gxout fwrite'
   
9. 'd uwnd'
   
10. 'disable fwrite'

复制代码

之后使用以下fortran代码读取

1. program test2
   
2.   implicit none
   
3.   integer,parameter::nx=144,ny=73,nt=349   
   
4.   real u(nx,ny,nt)
   
5. open(1,file='C:\Users\CKDDOUBI\Desktop\Fortran\Meteorological_Statistics\2\h200.grd',form='unformatted',access="direct",recl=nx*ny*nt*4)
   
6.    read(1,rec=1) u
   
7.   close(1)
   
8.   open(2,file='C:\Users\CKDDOUBI\Desktop\Fortran\Meteorological_Statistics\2\test.txt')
   
9.    write(2,*) u
   
10.   close(2)
    
11. end program test2

复制代码

但此时一直报错

1. C:\Users\CKDDOUBI\Desktop\Fortran\Meteorological_Statistics\2\2.f90:5:132:
   
2. 
   
3.      open(1,file='C:\Users\CKDDOUBI\Desktop\Fortran\Meteorological_Statistics\2\h200.grd',form='unformatted', access="direct",recl=nx*ny*nt*4)
   
4. 
   
5.                                1
   
6. Error: Line truncated at (1) [-Werror=line-truncation]
   
7. C:\Users\CKDDOUBI\Desktop\Fortran\Meteorological_Statistics\2\2.f90:5:132:
   
8. 
   
9.      open(1,file='C:\Users\CKDDOUBI\Desktop\Fortran\Meteorological_Statistics\2\h200.grd',form='unformatted', access="direct",recl=nx*ny*nt*4)
   
10. 
    
11.                                 1
    
12. Error: Syntax error in OPEN statement at (1)
    
13. f951.exe: some warnings being treated as errors

复制代码

注，报错信息中的两个“1” 在recl下面
我不知道哪里出了问题，求各位帮帮忙，被困扰两天了

任性小王子

为什么不直接用NCL处理了，要转成grd文件格式勒

CKDDOUBI

任性小王子 发表于 2023-10-22 19:47
为什么不直接用NCL处理了，要转成grd文件格式勒

还不会NCL，学校里只教了Fortran和GrADS。QAQ
周围人都在自学python，唉。

CKDDOUBI

已解决，Fortran第5行太长了，要续行

1. program test2
   
2.   implicit none
   
3.   integer, parameter::nx = 144, ny = 73, nt = 349
   
4.   real u(nx, ny, nt)
   
5.   open (1, file='C:\Users\CKDDOUBI\Desktop\Fortran\Meteorological_Statistics&
   
6. &\2\h200.grd', form='unformatted', access='direct', recl=nx*ny*nt*4)
   
7.   read (1, rec=1) u
   
8.   close (1)
   
9. 
   
10. end program test2
    

program test2<br />   implicit none<br />   integer, parameter::nx = 144, ny = 73, nt = 349<br />   real u(nx, ny, nt)<br />   open (1, file='C:\Users\CKDDOUBI\Desktop\Fortran\Meteorological_Statistics&<br /> &\2\h200.grd', form='unformatted', access='direct', recl=nx*ny*nt*4)<br />   read (1, rec=1) u<br />   close (1)<br /> <br /> end program test2<br />

任性小王子

CKDDOUBI 发表于 2023-10-22 19:59
还不会NCL，学校里只教了Fortran和GrADS。QAQ
周围人都在自学python，唉。

解决就好，现在确实Python比较强大,z在处理数据方面无可挑剔