netcdf_matlab_examples/modis-aqua/Commands.m

%
%  Ed Hill
%  $Id$
%

%  Data from:
%    http://mqabi.gsfc.nasa.gov/cgi-bin/MQABI/subset/map_04_week_D1.sst.cgi

matlab -nojvm

  !./pre_process
  
  clear all
  files = [ ...
      '2002065'; '2002073'; '2002089'; '2002097'; '2002105'; ...
      '2002113'; ...
      '2002145'; '2002153'; '2002161'; '2002169'; '2002177'; ...
      '2002185'; '2002193'; '2002201'; '2002209'; '2002217'; ...
          ];
  lats = linspace(47.0, 29.0, 410);
  lons = linspace(-78.0, -54.0, 546);
  i = 1;
  nt = size(files,1);
  days = [1:nt];
  all = zeros(nt,length(lats),length(lons));
  for i = 1:size(files,1)
    
    var = load([ 'sst-' files(i,:) ]);
    vt = var;
    for j = 1:prod(size(vt))
      if vt(j) == -99
        vt(j) = NaN;
      end
    end
    all(i,:,:) = vt;
    days(i) = str2double(files(i,5:7));

    % surf(lons,lats,vt), view(2), shading interp
    % pause(1)

  end

  id = 'SST';
  units = 'deg C (weekly average)';

  nc = netcdf('modis_aqua_gulfstr.nc', 'clobber');
  nc.reference = 'Various MODIS-Aqua data sets of the Gulf Stream';
  nc.author = 'Ed Hill <eh3@mit.edu>';
  nc.date = 'Jun 16, 2004';
  
  nc('T') = nt;
  nc('Y') = length(lats);
  nc('X') = length(lons);
  nc{'T'}  = 'T';
  nc{'Y'}  = 'Y';
  nc{'X'}  = 'X';
  nc{ id } = { 'T', 'Y', 'X' };
  nc{'X'}.uniquename = 'X';
  nc{'X'}.long_name = 'longitude';
  nc{'X'}.gridtype = ncint(1);
  nc{'X'}.units = 'degree_east';
  nc{'Y'}.uniquename = 'Y';
  nc{'Y'}.long_name = 'latitude';
  nc{'Y'}.gridtype = ncint(0);
  nc{'Y'}.units = 'degree_north';
  nc{'T'}.uniquename = 'T';
  nc{'T'}.long_name = 'Time';
  nc{'T'}.gridtype = ncint(0);
  nc{'T'}.units = 'days since Jan 1, 2002';
  nc{ id }.units = units;
  nc{ id }.long_name = id;
  nc{ id }.missing_value = ncdouble(NaN);
  nc{ id }.FillValue_ = ncdouble(0.);
  nc{'Y'}(:) = lats;
  nc{'X'}(:) = lons;
  nc{'T'}(:) = days;
  nc{ id }(:) = permute(all,[ 1 2 3 ]);
  nc = close(nc);

  ! scp ./modis_aqua_gulfstr.nc channel.mit.edu:/home/edhill/INGRID_PEOPLE/EH3/modis-aqua
1	edhill	1.1	%
2			% Ed Hill
3			% $Id$
4			%
5
6			% Data from:
7			% http://mqabi.gsfc.nasa.gov/cgi-bin/MQABI/subset/map_04_week_D1.sst.cgi
8
9			matlab -nojvm
10
11			!./pre_process
12
13			clear all
14			files = [ ...
15			'2002065'; '2002073'; '2002089'; '2002097'; '2002105'; ...
16			'2002113'; ...
17			'2002145'; '2002153'; '2002161'; '2002169'; '2002177'; ...
18			'2002185'; '2002193'; '2002201'; '2002209'; '2002217'; ...
19			];
20			lats = linspace(47.0, 29.0, 410);
21			lons = linspace(-78.0, -54.0, 546);
22			i = 1;
23			nt = size(files,1);
24			days = [1:nt];
25			all = zeros(nt,length(lats),length(lons));
26			for i = 1:size(files,1)
27
28			var = load([ 'sst-' files(i,:) ]);
29			vt = var;
30			for j = 1:prod(size(vt))
31			if vt(j) == -99
32			vt(j) = NaN;
33			end
34			end
35			all(i,:,:) = vt;
36			days(i) = str2double(files(i,5:7));
37
38			% surf(lons,lats,vt), view(2), shading interp
39			% pause(1)
40
41			end
42
43			id = 'SST';
44			units = 'deg C (weekly average)';
45
46			nc = netcdf('modis_aqua_gulfstr.nc', 'clobber');
47			nc.reference = 'Various MODIS-Aqua data sets of the Gulf Stream';
48			nc.author = 'Ed Hill <eh3@mit.edu>';
49			nc.date = 'Jun 16, 2004';
50
51			nc('T') = nt;
52			nc('Y') = length(lats);
53			nc('X') = length(lons);
54			nc{'T'} = 'T';
55			nc{'Y'} = 'Y';
56			nc{'X'} = 'X';
57			nc{ id } = { 'T', 'Y', 'X' };
58			nc{'X'}.uniquename = 'X';
59			nc{'X'}.long_name = 'longitude';
60			nc{'X'}.gridtype = ncint(1);
61			nc{'X'}.units = 'degree_east';
62			nc{'Y'}.uniquename = 'Y';
63			nc{'Y'}.long_name = 'latitude';
64			nc{'Y'}.gridtype = ncint(0);
65			nc{'Y'}.units = 'degree_north';
66			nc{'T'}.uniquename = 'T';
67			nc{'T'}.long_name = 'Time';
68			nc{'T'}.gridtype = ncint(0);
69			nc{'T'}.units = 'days since Jan 1, 2002';
70			nc{ id }.units = units;
71			nc{ id }.long_name = id;
72			nc{ id }.missing_value = ncdouble(NaN);
73			nc{ id }.FillValue_ = ncdouble(0.);
74			nc{'Y'}(:) = lats;
75			nc{'X'}(:) = lons;
76			nc{'T'}(:) = days;
77			nc{ id }(:) = permute(all,[ 1 2 3 ]);
78			nc = close(nc);
79
80			! scp ./modis_aqua_gulfstr.nc channel.mit.edu:/home/edhill/INGRID_PEOPLE/EH3/modis-aqua