/[MITgcm]/MITgcm_contrib/gael/profilesMatlabProcessing/profiles_misc/MITprof_resample.m
ViewVC logotype

Diff of /MITgcm_contrib/gael/profilesMatlabProcessing/profiles_misc/MITprof_resample.m

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.3 by gforget, Fri Oct 9 01:30:12 2015 UTC revision 1.4 by gforget, Fri Oct 9 23:43:34 2015 UTC
# Line 57  for tt=1:length(rec_fld)-1; Line 57  for tt=1:length(rec_fld)-1;
57      if ~isempty(ii);      if ~isempty(ii);
58          arr=gcmfaces_interp(fld,lon(ii),lat(ii),method);          arr=gcmfaces_interp(fld,lon(ii),lat(ii),method);
59          arr2=gcmfaces_interp_1d(2,depIn,arr,depOut);          arr2=gcmfaces_interp_1d(2,depIn,arr,depOut);
         if strcmp(fldIn.tim,'std');  
             arr2=arr2.*randn(size(arr2));  
         end;  
60      end;              end;        
61      %          %    
62      k0=floor(tim_prof(ii)); k1=k0+1;      k0=floor(tim_prof(ii)); k1=k0+1;
63      a0=tim_prof(ii)-k0; a0=a0*ones(1,profIn.nr);      a0=tim_prof(ii)-k0; a0=a0*ones(1,profIn.nr);
64      profOut(ii,:)=(1-a0).*arr2(:,:,1)+a0.*arr2(:,:,2);      profOut(ii,:)=(1-a0).*arr2(:,:,1)+a0.*arr2(:,:,2);
65        if strcmp(fldIn.tim,'std');
66            profOut(ii,:)=profOut(ii,:).*randn(size(profOut(ii,:)));
67        end;
68  end;  end;
69    
70  %3) deal with file output  %3) deal with file output
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22