Diagnostics/DiagUtility/test_nc_add.m

xi=-180:4:180;
yi=-90:4:90;
[X,Y]=ndgrid(xi,yi);
q=sin(X*pi/180).*cos(Y*pi/180);

nc=netcdf('test.nc','clobber');

nc_add(nc,'lon','Longitude','degrees east',{'lon'},xi)
nc_add(nc,'lat','Latitude','degrees north',{'lat'},yi)
nc_add(nc,'time','Time','Years',{'time'},[])
nc_add(nc,'q','Random','nothing',{'lat' 'lon'},q)

nc_add(nc,'time','Time','Years',{'time'},1977,1)
nc_add(nc,'time','Time','Years',{'time'},1979,2)
nc_add(nc,'w','Random','nothing',{'time' 'lat' 'lon'},q+1,1)
nc_add(nc,'w','Random','nothing',{'time' 'lat' 'lon'},q+2,2)

nc=close(nc);
1	xi=-180:4:180;
2	yi=-90:4:90;
3	[X,Y]=ndgrid(xi,yi);
4	q=sin(Xpi/180).cos(Y*pi/180);
5
6	nc=netcdf('test.nc','clobber');
7
8	nc_add(nc,'lon','Longitude','degrees east',{'lon'},xi)
9	nc_add(nc,'lat','Latitude','degrees north',{'lat'},yi)
10	nc_add(nc,'time','Time','Years',{'time'},[])
11	nc_add(nc,'q','Random','nothing',{'lat' 'lon'},q)
12
13	nc_add(nc,'time','Time','Years',{'time'},1977,1)
14	nc_add(nc,'time','Time','Years',{'time'},1979,2)
15	nc_add(nc,'w','Random','nothing',{'time' 'lat' 'lon'},q+1,1)
16	nc_add(nc,'w','Random','nothing',{'time' 'lat' 'lon'},q+2,2)
17
18	nc=close(nc);