Fortran NetCDF - добавил новое измерение нужно заполнить нулями

Question

я добавил новое измерение в существующий файл NetCDF в Fortran, используя следующий код -

retval = nf_open(cfn,NF_WRITE,ncid) 
if (retval .ne. nf_noerr) call handle_err(retval) 
retval = nf_redef(ncid) 
if (retval .ne. nf_noerr) call handle_err(retval) 
retval = nf_def_dim(ncid,"xyz",len,dimid_xyz) 
if (retval .ne. nf_noerr) call handle_err(retval)   
retval = nf_enddef(ncid) 

Теперь я хочу, чтобы иметь возможность заполнить это измерение с значения нуля. Мощность этого множества равна мощности переменной в моей геополитической высоте. Кроме того, у меня есть три других измерения - время (неограниченное), широта, долгота и уровень.

Я посмотрел API NetCDF в Фортран, и я не уверен, что это API для call.When Я использую следующие API

retval = nf_put_var_real(ncid,dimid_xyz,xyzArray) 
    if (retval .ne. nf_noerr) call handle_err(retval) 

он заканчивает перезапись геопотенциала значения высоты с 0.0 (который единственная переменная в моем файле netCDF)

Как мне это сделать?

Answer 1

Как я понимаю, измерение отличается от переменной, размеры не могут иметь значения, но переменные могут - я думаю, довольно распространенной практикой может быть создание измерения, а также создание переменной с тем же именем. Затем вы можете дать переменной любые значения, которые вы хотите.

Ваш код может выглядеть

retval = nf_open(cfn,NF_WRITE,ncid) 
if (retval .ne. nf_noerr) call handle_err(retval) 
retval = nf_redef(ncid) 
if (retval .ne. nf_noerr) call handle_err(retval) 
retval = nf_def_dim(ncid,"xyz",len,dimid_xyz) 
if (retval .ne. nf_noerr) call handle_err(retval) 

retval = nf_def_var(ncid,"xyz",netcdf_real,1,[dimid_xyz], varid_xyz) 
if (retval .ne. nf_noerr) call handle_err(retval) 

retval = nf_enddef(ncid) 

retval = nf_put_vara(ncid,varid_xyz,[1],[len],arrayOfZero) 
if (retval .ne. nf_noerr) call handle_err(retval) 

Примечание Я рекомендую использовать переменную с именем len внутри fortran кода - это будет столкновение с внутренней одного и того же имени.