loadANTS_simple.m
author A.M. Thurnherr <ant@ldeo.columbia.edu>
Tue, 21 Feb 2012 14:30:23 -0500
changeset 5 2e9ce6753135
parent 3 5e3caf9c0f2e
permissions -rw-r--r--
.
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
3
5e3caf9c0f2e first whoosher version
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
     1
% Usage: D = loadANTS(filename[,fieldname[,...]]);
5e3caf9c0f2e first whoosher version
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
     2
%
5e3caf9c0f2e first whoosher version
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
     3
% Examples:
5e3caf9c0f2e first whoosher version
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
     4
%	prof = loadANTS('070.2db');
5e3caf9c0f2e first whoosher version
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
     5
%		- read all fields from file 070.2db into matrix prof
5e3caf9c0f2e first whoosher version
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
     6
%		- each field becomes a column
5e3caf9c0f2e first whoosher version
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
     7
%	prof = loadANTS('070.2db','press','temp','salin');
5e3caf9c0f2e first whoosher version
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
     8
%		- read pressure, temperature, and salinity fields from file 070.2db
5e3caf9c0f2e first whoosher version
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
     9
%		  into matrix prof
5e3caf9c0f2e first whoosher version
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    10
%		- 1st column is pressure, 2nd column is temperature, 3rd column is
5e3caf9c0f2e first whoosher version
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    11
%		  salinity
5e3caf9c0f2e first whoosher version
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    12
5e3caf9c0f2e first whoosher version
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    13
%======================================================================
5e3caf9c0f2e first whoosher version
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    14
%                    L O A D A N T S . M 
5e3caf9c0f2e first whoosher version
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    15
%                    doc: Fri Sep 17 11:12:44 2010
5e3caf9c0f2e first whoosher version
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    16
%                    dlm: Fri Sep 17 11:41:01 2010
5e3caf9c0f2e first whoosher version
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    17
%                    (c) 2010 A.M. Thurnherr
5e3caf9c0f2e first whoosher version
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    18
%                    uE-Info: 25 29 NIL 0 0 72 2 2 4 NIL ofnI
5e3caf9c0f2e first whoosher version
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    19
%======================================================================
5e3caf9c0f2e first whoosher version
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    20
5e3caf9c0f2e first whoosher version
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    21
function D = loadANTS(file,varargin)
5e3caf9c0f2e first whoosher version
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    22
5e3caf9c0f2e first whoosher version
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    23
if nargin==1
5e3caf9c0f2e first whoosher version
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    24
	cmd = sprintf('data -Qw %s',file);
5e3caf9c0f2e first whoosher version
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    25
	[status,dta] = system(cmd);
5e3caf9c0f2e first whoosher version
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    26
	if status~=0
5e3caf9c0f2e first whoosher version
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    27
		error(sprintf('cmd "%s" failed',cmd));
5e3caf9c0f2e first whoosher version
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    28
	end
5e3caf9c0f2e first whoosher version
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    29
	nfields = sscanf(dta,'%d');
5e3caf9c0f2e first whoosher version
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    30
	cmd = sprintf('Cat -Q %s',file);
5e3caf9c0f2e first whoosher version
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    31
else
5e3caf9c0f2e first whoosher version
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    32
	nfields = nargin-1;
5e3caf9c0f2e first whoosher version
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    33
	cmd = 'list -Q';
5e3caf9c0f2e first whoosher version
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    34
	for i=1:nfields
5e3caf9c0f2e first whoosher version
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    35
	  cmd = sprintf('%s %s',cmd,varargin{i});
5e3caf9c0f2e first whoosher version
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    36
	end
5e3caf9c0f2e first whoosher version
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    37
	cmd = sprintf('%s %s',cmd,file);
5e3caf9c0f2e first whoosher version
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    38
end
5e3caf9c0f2e first whoosher version
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    39
5e3caf9c0f2e first whoosher version
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    40
[status,dta] = system(cmd);
5e3caf9c0f2e first whoosher version
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    41
if status~=0
5e3caf9c0f2e first whoosher version
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    42
	error(sprintf('cmd "%s" failed',cmd));
5e3caf9c0f2e first whoosher version
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    43
end
5e3caf9c0f2e first whoosher version
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    44
5e3caf9c0f2e first whoosher version
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    45
D = sscanf(dta,'%f',[nfields,inf])';
5e3caf9c0f2e first whoosher version
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    46