function Convert2NC(grd,filename) % E. Di Lorenzo (edl@ucsd.edu) % nc = netcdf(filename, 'clobber'); if isempty(nc), return, end %% Global attributes: nc.type = ncchar(''); nc.title = ncchar(''); nc.out_file = ncchar(''); nc.grd_file = ncchar(''); nc.inp_file = ncchar(''); nc.history = ncchar(''); %% Dimensions: nc('lon') = grd.L; nc('lat') = grd.P; nc('time') = grd.T; %% Variables and attributes: nc{'time'} = ncdouble('time'); %% 1 element. nc{'time'}.long_name = ncchar('time since initialization'); nc{'time'}.units = ncchar('second'); nc{'time'}.field = ncchar('time, scalar, series'); nc{'u'} = ncfloat('time', 'lat', 'lon'); %% 189600 elements. nc{'u'}.long_name = ncchar('u-momentum component'); nc{'u'}.units = ncchar('cm second-1'); nc{'u'}.field = ncchar('u-velocity, scalar, series'); nc{'v'} = ncfloat('time', 'lat', 'lon'); %% 190400 elements. nc{'v'}.long_name = ncchar('v-momentum component'); nc{'v'}.units = ncchar('cm second-1'); nc{'v'}.field = ncchar('v-velocity, scalar, series'); endef(nc) close(nc) return