Changeset 153

Timestamp:

07/08/15 17:07:32 (9 years ago)

Author:

jbrlod

Message:

modif interp back to handle integer time steps

File:

• altifloat/matlab_toolbox/interp_bck.m (modified) (4 diffs)

altifloat/matlab_toolbox/interp_bck.m

@@ -1 @@
-function [UV_bck,mask_back]=interp_bck(uv,it,dt,meshmask,outfile,outmask)
+function [UV_bck,mask_back]=interp_bck(uv,it,dt,meshmask,outfile,outmask,flag)
  %[]=interp_bck(uv_b1,uv_b2,it1,it2,dt)
  %Time unit : hour
@@ -21 +21 @@
 UV_bck=[];
 mask_back=ones(Nlat*Nlong,1);
-
+yao_time=1;
 
 for idf=1:length(uv)-1
@@ -32 +32 @@
     %     uv1=load(uv_b1);
     %     uv2=load(uv_b2);
-    [uv1,mask1]=aviso2yao(uv_b1,meshmask,[],[],true);
-    [uv2,mask2]=aviso2yao(uv_b2,meshmask,[],[],true);
+    [uv1,mask1]=aviso2yao(uv_b1,meshmask,[],[],flag);
+    [uv2,mask2]=aviso2yao(uv_b2,meshmask,[],[],flag);
     
     mask_back=mask_back.*mask1.*mask2;
@@ -62 +62 @@
        vint=vint';
        %end
-       tt=repmat(t,size(uv1,1),1);
+       
+       tt=repmat(yao_time:1:yao_time+length(t)-1,size(uv1,1),1);
+       yao_time=yao_time+length(t);
        tt=tt(:);
        
        UV_bck=[UV_bck;[tt,repmat(uv1(:,2:3),length(t),1),uint(:),vint(:)]];
-
+fprintf(1,'size(t)=%d (it1=%d,it2=%d)\n',length(t),it1,it2);
 end %for idf