generate_mjd

--------------------------------------------------------------
 dump to netcdf file
--------------------------------------------------------------

0001 %--------------------------------------------------------------
0002 % dump to netcdf file
0003 %--------------------------------------------------------------
0004 
0005 start_day = greg2mjulian(2007,4,1,1,0,0);
0006 end_day   = greg2mjulian(2007,7,1,0,0,0);
0007 
0008 time = start_day:(1./24):end_day;
0009 
0010 % open boundary forcing
0011 nc = netcdf('gom1v10_decelles_2007_time.nc', 'clobber');       
0012 
0013 
0014 % dimensions
0015 nc('time') = 0; 
0016 
0017 % variables
0018 nc{'time'} = ncfloat('time');
0019 nc{'time'}.long_name = 'time';  
0020 nc{'time'}.units     = 'days since 0.0';  
0021 nc{'time'}.time_zone = 'none';  
0022 
0023 nc{'Itime'} = ncint('time');
0024 nc{'Itime'}.units     = 'days since 0.0';  
0025 nc{'Itime'}.time_zone = 'none';  
0026 
0027 nc{'Itime2'} = ncint('time');
0028 nc{'Itime2'}.units     = 'msec since 00:00:00';
0029 nc{'Itime2'}.time_zone = 'none';  
0030 
0031 
0032 % dump dynamic data
0033 ntimes = numel(time);
0034 nc{'time'}(1:ntimes) = time(1:ntimes);
0035 nc{'Itime'}(1:ntimes) = floor(time(1:ntimes));
0036 nc{'Itime2'}(1:ntimes) = mod(time(1:ntimes),1)*24*3600*1000.;
0037 
0038 nc = close(nc);    
0039 
0040