【发布时间】：2015-10-22 12:31:30
【问题描述】：
您好，关于 stackoverflow 的第一个问题。
我已经坚持了 5 天。我想用 Fortran 写一个 netcdf 文件。
我正在使用 netcdf/3.6.3 我正在尝试输出一个名为 frech 的大型 43000x 18000 数组和一些较小的数组（一维数组）。
下面是我的代码示例，它是一个非常大的文件，问题不在于为变量输入值，问题在于结束变量的定义：
print*,"nunks is",nunks
print*,"neqns is",neqns
ok=nf90_create('michalek.nc', NF90_CLOBBER, ncid)
print *,"create ok=",ok
ok= nf90_def_dim(ncid,"nunks", nunks, nunks_dimid)
print *,"def nunks dimension ",ok
ok= nf90_def_dim(ncid,"neqns", neqns, neqns_dimid)
print *,"def neqns dimension ",ok
dimids=(/neqns_dimid, nunks_dimid/)
print *,dimids
ok= nf90_def_var(ncid,"frech", NF90_REAL, dimids, frech_varid)
print *,"def frech",ok
ok= nf90_def_var(ncid,"src", NF90_REAL, nunks_dimid, src_varid)
print *,"def src",ok
ok= nf90_def_var(ncid,"csrc", NF90_REAL, nunks_dimid, csrc_varid)
print *,"define csrc",ok
ok= nf90_def_var(ncid,"dat", NF90_REAL, neqns_dimid, dat_varid)
print *,"define dat",ok
ok= nf90_def_var(ncid,"cdat", NF90_REAL, neqns_dimid, cdat_varid)
print *,"define cdat",ok
ok= nf90_enddef(ncid)
print *,"end dif ", ok
ok= nf90_put_var(ncid, frech_varid, frech)
print *, 'frech put in ok=',ok
ok= nf90_put_var(ncid, src_varid, src)
print *, 'src put in ok=',ok
ok= nf90_put_var(ncid, csrc_varid, csrc)
print *, 'csrc put in ok=',ok
ok= nf90_put_var(ncid, dat_varid, dat)
print *, 'dat put in ok=',ok
ok= nf90_put_var(ncid, cdat_varid, cdat)
print *, 'cdat put in ok=',ok
ok= nf90_close(ncid)
print *, 'close?',ok
我知道正确读取文件时 ok=0 但是，当我到达结束文件定义的阶段时 (nf90_enddif) ok 返回为 =-62 并且未创建 netcdf 文件。我想这是一个太大的数组的问题，但我无法解决这个问题
上述代码的相关输出为：
nunks is 43894
neqns is 18144
create ok= 0
def nunks dimension 0
def neqns dimension 0
2 1
def frech 0
def src 0
define csrc 0
define dat 0
define cdat 0
end dif -62
frech put in ok= -39
src put in ok= -39
csrc put in ok= -39
dat put in ok= -39
cdat put in ok= -39
close? -62
感谢您的帮助！
问候 彼得:)
【问题讨论】：
-
您是否尝试过使用可变尺寸来查看较小尺寸是否成功？例如，选择 1000x1000 并逐渐增加以查看问题所在
-
这是个好主意。但我想避免对代码进行过多的试验，因为它需要 10 分钟才能执行。也许我应该尝试缩短代码的各个方面。
标签： arrays fortran fortran90 netcdf dataoutputstream