matlab在nc中加入新变量

%% 清楚变量
clear,clc;

%% 读nc
file2=['restart_0113.nc'];%%%%%
b=ncinfo(file2);
ss=struct(b.Variables);
name2=cellstr(char(ss.Name));
for i=33
    eval(['variable_',char(name2(i)),'=ncread(file2,name2{i});']);%%读取nc变量中第33个bv变量数据，赋值给variable_bv
end
%% 打开nc
copyfile(which('Bohai_restart_okaice0110-0113.nc'),'Bohai_restart_okaice0113.nc');
ncid = netcdf.open('Bohai_restart_okaice0113.nc', 'WRITE');

% 新变量是以原本nc中的维度进行设置的，所以此处需要先把维度ID提取出来。
nele_dim=netcdf.inqDimID(ncid,'nele');
node_dim=netcdf.inqDimID(ncid,'node');

% nc进入定义模式
netcdf.reDef(ncid);

%-----------------------给维度分配变量-------------------------------------
bv_var_id = netcdf.defVar(ncid,'bv','float',[node_dim,siglev_dim,time_dim]);
hs_var_id = netcdf.defVar(ncid,'hs','float',[node_dim,time_dim]);

%-----------------------赋值-----------------------------------------------
netcdf.putVar(ncid,bv_var_id,variable_bv);

%------------------------描述及单位信息------------------------------------
%bv
netcdf.putAtt(ncid,bv_var_id,'long_name','wave-induced vertical kinematic viscosity or diffusivity');
netcdf.putAtt(ncid,bv_var_id,'units','m 2 s-1');
netcdf.putAtt(ncid,bv_var_id,'grid','fvcom_grid');

%hs
netcdf.putAtt(ncid,hs_var_id,'long_name','Significant Wave Height');
netcdf.putAtt(ncid,hs_var_id,'units','meters'); 

%% 关闭nc
netcdf.close(ncid);
toc