小编jkp*_*jkp的帖子

我在 Fortran 中编写了一个代码来读取具有 4 维数据 [时间、水平、经度、纬度] 的 NetCDF 文件。但是，我的代码产生错误

NetCDF: Start+count exceeds dimension bound

在我使用的任何 4-d NetCDF 文件上。例如，http://people.sc.fsu.edu/~jburkardt/f_src/netcdf/pres_temp_4D.nc上的文件有压力和温度。我在下面粘贴我的代码。请建议出了什么问题。

程序 rw_nc4d_main

  仅使用 rw_nc4d：read_nc4

  隐式无

  字符（LEN = 50）:: ncfn
  字符（LEN = 15）:: vname

  ncfn = 'pres_temp_4D.nc'
  vname = '压力'

  调用 read_nc4(ncfn, vname)    

结束程序 rw_nc4d_main
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!!!!!!!!!!!!!!!!
模块 rw_nc4d

使用 netcdf

隐式无

包含

  子程序 read_nc4(fname,vin_name)   

    隐式无

    字符（LEN = *），意图（输入）:: fname
    字符（LEN = *），意图（输入）:: vin_name

    ！局部变量
    整数 :: ncid、var_id、ndim、nvar、nattr、unlim_id
    字符(LEN=15) :: dname
    整数 :: dlength
    整数 :: ii、状态、lx、ly、lz、lt、lzp1 
    真实的:: SF, ofs
    真实，维度（：，：，：，：），可分配：：vin …

我通过在各个地方包含印刷语句来测试一些fortran代码.我还想在print语句中包含一个行号,例如:

...
write(*,*)'[current line #] I am here'
.....
.....
write(*,*)'[current line #] I am here too'

目前,我通过查看状态栏手动插入行号.有没有相同的宏？

我正在尝试使用read.table读取R中的文件(ascii)该文件如下所示:

  DAILY MAXIMUM TEMPARATURE  
  YEAR DAY MT DT   LAT. 66.5   67.5   68.5   69.5   70.5
  1969 001 01 01   6.5  99.90  99.90  31.90  99.90  99.90 
  1969 001 01 01   7.5  99.90  20.90  99.90  99.90  23.90
  1969 001 01 01   8.5  99.90  99.90  30.90  99.90  18.90
  .....
  ..... 
  YEAR DAY MT DT   LAT. 66.5   67.5   68.5   69.5   70.5
  1969 001 01 02   6.5  21.90  99.90  99.90  99.90  99.90 
  1969 001 01 02   7.5  99.90  33.90  99.90  99.90  99.90
  1969 001 01 02 …