/[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.2 by gforget, Sat Apr 18 12:54:20 2015 UTC revision 1.4 by gforget, Fri Oct 9 23:43:34 2015 UTC
# Line 35  if strcmp(fldIn.tim,'monclim'); Line 35  if strcmp(fldIn.tim,'monclim');
35      tim_prof=(profIn.prof_date-tmp2);      tim_prof=(profIn.prof_date-tmp2);
36      tim_prof(tim_prof>365)=365;      tim_prof(tim_prof>365)=365;
37      tim_prof=tim_prof/365*12;%neglecting differences in months length      tim_prof=tim_prof/365*12;%neglecting differences in months length
38  elseif strcmp(fldIn.tim,'const');  elseif strcmp(fldIn.tim,'const')|strcmp(fldIn.tim,'std');
39      tim_fld=[1 2]; rec_fld=[1 2];      tim_fld=[1 2]; rec_fld=[1 1];
40      tim_prof=1.5*ones(profIn.np,1);      tim_prof=1.5*ones(profIn.np,1);
41  else;  else;
42      error('this case remains to be implemented');      error('this case remains to be implemented');
# Line 55  for tt=1:length(rec_fld)-1; Line 55  for tt=1:length(rec_fld)-1;
55      %      %
56      ii=find(tim_prof>=tim_fld(tt)&tim_prof<tim_fld(tt+1));      ii=find(tim_prof>=tim_fld(tt)&tim_prof<tim_fld(tt+1));
57      if ~isempty(ii);      if ~isempty(ii);
58          arr=gcmfaces_interp(fld,lon(ii),lat(ii),'polygons');          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);
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.2  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22