mlosch/interp_llc/extract_levitus_monthly.m

function [ts,x,y,z]=extract_levitus_monthly(suff,months,lon_lo,lon_hi,lat_lo,lat_hi)
%
% Extract monthly Levitus data from netcdf files 
%
% e.g.
%   [T,x,y,z]=extract_levitus_monthly('temp',1,-180,180,-90,90);
%   [S,x,y,z]=extract_levitus_monthly('salt',1:12,0,360,-90,90);

% change a name convention:
  if strcmp(suff,'temp'); suff = 'theta'; end
  
  
  ncquiet;
  load DATAPATH datapath
  fname = fullfile(datapath,['levitus_month_' suff '.cdf']);
  theResult=ncload(fname);
  if isempty(theResult);
    error([fname ' not found'])
  end
  
  if size(X,1) > 1; X = X'; end
  if size(Y,1) > 1; Y = Y'; end
  if size(Z,1) > 1; Z = Z'; end
  if size(T,1) > 1; T = T'; end

  i0=round( lon_lo )+1;
  i1=round( lon_hi );
  j0=round( (lat_lo+90) )+1;
  j1=round( (lat_hi+90) );
  
  if i1>i0
    ii=i0:i1;
  else
    ii=[i0:360 1:i1];
  end
  ii=mod(ii+360-1,360)+1;
  jj=j0:j1;
  
  nx=prod(size(ii));
  ny=prod(size(jj));
  
  if strcmp(suff,'theta');
    ts=permute(theta(months,:,jj,ii),[4 3 2 1]);
  elseif strcmp(suff,'salt');
    ts=permute(sal(months,:,jj,ii),[4 3 2 1]);
  else
    error(['unknown variable ' suff])
  end
    
  kk = find(isnan(ts));
  ts(kk) = 0;

  x=X(ii);
  ih = find(x>lon_hi); x(ih)=x(ih)-360;
  il = find(x<lon_lo); x(il)=x(il)+360;
  y=Y(jj);
  z=-Z;

  return
1	mlosch	1.1	function [ts,x,y,z]=extract_levitus_monthly(suff,months,lon_lo,lon_hi,lat_lo,lat_hi)
2			%
3			% Extract monthly Levitus data from netcdf files
4			%
5			% e.g.
6			% [T,x,y,z]=extract_levitus_monthly('temp',1,-180,180,-90,90);
7			% [S,x,y,z]=extract_levitus_monthly('salt',1:12,0,360,-90,90);
8
9			% change a name convention:
10			if strcmp(suff,'temp'); suff = 'theta'; end
11
12
13			ncquiet;
14			load DATAPATH datapath
15			fname = fullfile(datapath,['levitus_month_' suff '.cdf']);
16			theResult=ncload(fname);
17			if isempty(theResult);
18			error([fname ' not found'])
19			end
20
21			if size(X,1) > 1; X = X'; end
22			if size(Y,1) > 1; Y = Y'; end
23			if size(Z,1) > 1; Z = Z'; end
24			if size(T,1) > 1; T = T'; end
25
26			i0=round( lon_lo )+1;
27			i1=round( lon_hi );
28			j0=round( (lat_lo+90) )+1;
29			j1=round( (lat_hi+90) );
30
31			if i1>i0
32			ii=i0:i1;
33			else
34			ii=[i0:360 1:i1];
35			end
36			ii=mod(ii+360-1,360)+1;
37			jj=j0:j1;
38
39			nx=prod(size(ii));
40			ny=prod(size(jj));
41
42			if strcmp(suff,'theta');
43			ts=permute(theta(months,:,jj,ii),[4 3 2 1]);
44			elseif strcmp(suff,'salt');
45			ts=permute(sal(months,:,jj,ii),[4 3 2 1]);
46			else
47			error(['unknown variable ' suff])
48			end
49
50			kk = find(isnan(ts));
51			ts(kk) = 0;
52
53			x=X(ii);
54			ih = find(x>lon_hi); x(ih)=x(ih)-360;
55			il = find(x<lon_lo); x(il)=x(il)+360;
56			y=Y(jj);
57			z=-Z;
58
59			return