gael/matlab/myverticaldeform.m

function [varout]=myverticaldeform(varin)
%varin : profondeur en m (positif)
%varout : z correspondant (negatif), avec 
%       deformation selon le niveau

%var de sortie :
varout=zeros(size(varin));
%valeurs limite :
%val_lim=500;fprintf('deform v:500m\n');
%val_lim=1000;
val_lim=250; fprintf('deform v:250m\n');
%val_lim=-0.0001;
val_lim2=5500;

%valeurs inferieures a 500m :
tmp1=find(varin<=val_lim);
tmp2=-100*varin(tmp1)/val_lim;
varout(tmp1)=tmp2;

%valeurs superieurs a 500m :
tmp1=find(varin>val_lim);
tmp2=-100-200*(varin(tmp1)-val_lim)./(val_lim2-val_lim);
varout(tmp1)=tmp2;


1	function [varout]=myverticaldeform(varin)
2	%varin : profondeur en m (positif)
3	%varout : z correspondant (negatif), avec
4	% deformation selon le niveau
5
6	%var de sortie :
7	varout=zeros(size(varin));
8	%valeurs limite :
9	%val_lim=500;fprintf('deform v:500m\n');
10	%val_lim=1000;
11	val_lim=250; fprintf('deform v:250m\n');
12	%val_lim=-0.0001;
13	val_lim2=5500;
14
15	%valeurs inferieures a 500m :
16	tmp1=find(varin<=val_lim);
17	tmp2=-100*varin(tmp1)/val_lim;
18	varout(tmp1)=tmp2;
19
20	%valeurs superieurs a 500m :
21	tmp1=find(varin>val_lim);
22	tmp2=-100-200*(varin(tmp1)-val_lim)./(val_lim2-val_lim);
23	varout(tmp1)=tmp2;
24
25
26