profilesMatlabProcessing/profiles_misc/model_addinMITprof.m

function model_addinMITprof(ncfile,model,varargin)
% function model_addinMITprof(ncfile,model)
%   add hydrographic profiles from a model interpolated on profile
%   positions available in a MITprof netcdf file.
%
%   fileIn: path (absolute or relative) to MITprof netcdf file
%   model is a string used to select the model to be loaded
%       'OCCA' : ECCOv4 grid + OCCA atlas
%       'SOSE59' : SOSE59 grid + atlas

% process file name
[pathstr, name, ext] = fileparts(ncfile);
if isempty(pathstr) | strcmp(pathstr,'.'), pathstr=pwd; end
if isempty(ext) | ~strcmp(ext,'.nc'), ext='.nc'; end
ncfile=[pathstr '/' name ext];

% load profiles
M=MITprof_load(ncfile);

% interpolate model
M=model_interp(M,model);

% add fields
fillval=double(-9999);
T=getfield(M,['prof_T_' model])'; T(isnan(T))=fillval;
MITprof_addVar(ncfile,['prof_T_' model],'double',{'iDEPTH','iPROF'},T);
S=getfield(M,['prof_S_' model])'; S(isnan(S))=fillval;
MITprof_addVar(ncfile,['prof_S_' model],'double',{'iDEPTH','iPROF'},S);

% add attributes
nc=ncopen(ncfile,'write');
ncaddAtt(nc,['prof_T_' model],'long_name',['pot. temp. model ' model]);
ncaddAtt(nc,['prof_T_' model],'units','degree C');
ncaddAtt(nc,['prof_T_' model],'missing_value',fillval);
ncaddAtt(nc,['prof_T_' model],'_FillValue',fillval);
ncaddAtt(nc,['prof_S_' model],'long_name',['salinity model ' model]);
ncaddAtt(nc,['prof_S_' model],'units','psu');
ncaddAtt(nc,['prof_S_' model],'missing_value',fillval);
ncaddAtt(nc,['prof_S_' model],'_FillValue',fillval);
ncclose(nc);
1	roquet	1.1	function model_addinMITprof(ncfile,model,varargin)
2			% function model_addinMITprof(ncfile,model)
3			% add hydrographic profiles from a model interpolated on profile
4			% positions available in a MITprof netcdf file.
5			%
6			% fileIn: path (absolute or relative) to MITprof netcdf file
7			% model is a string used to select the model to be loaded
8			% 'OCCA' : ECCOv4 grid + OCCA atlas
9			% 'SOSE59' : SOSE59 grid + atlas
10
11			% process file name
12			[pathstr, name, ext] = fileparts(ncfile);
13			if isempty(pathstr) \| strcmp(pathstr,'.'), pathstr=pwd; end
14			if isempty(ext) \| ~strcmp(ext,'.nc'), ext='.nc'; end
15			ncfile=[pathstr '/' name ext];
16
17			% load profiles
18			M=MITprof_load(ncfile);
19
20			% interpolate model
21			M=model_interp(M,model);
22
23			% add fields
24			fillval=double(-9999);
25			T=getfield(M,['prof_T_' model])'; T(isnan(T))=fillval;
26			MITprof_addVar(ncfile,['prof_T_' model],'double',{'iDEPTH','iPROF'},T);
27			S=getfield(M,['prof_S_' model])'; S(isnan(S))=fillval;
28			MITprof_addVar(ncfile,['prof_S_' model],'double',{'iDEPTH','iPROF'},S);
29
30			% add attributes
31			nc=ncopen(ncfile,'write');
32			ncaddAtt(nc,['prof_T_' model],'long_name',['pot. temp. model ' model]);
33			ncaddAtt(nc,['prof_T_' model],'units','degree C');
34			ncaddAtt(nc,['prof_T_' model],'missing_value',fillval);
35			ncaddAtt(nc,['prof_T_' model],'_FillValue',fillval);
36			ncaddAtt(nc,['prof_S_' model],'long_name',['salinity model ' model]);
37			ncaddAtt(nc,['prof_S_' model],'units','psu');
38			ncaddAtt(nc,['prof_S_' model],'missing_value',fillval);
39			ncaddAtt(nc,['prof_S_' model],'_FillValue',fillval);
40			ncclose(nc);