/[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.1 by gforget, Fri Apr 17 20:39:21 2015 UTC revision 1.3 by gforget, Fri Oct 9 01:30:12 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 48  profOut=NaN*ones(profIn.np,profIn.nr); Line 48  profOut=NaN*ones(profIn.np,profIn.nr);
48    
49  %2) loop over record pairs  %2) loop over record pairs
50  for tt=1:length(rec_fld)-1;  for tt=1:length(rec_fld)-1;
51      fld0=read_bin(fldIn.fil,rec_fld(tt));      fld0=mygrid.mskC.*read_bin(fldIn.fil,rec_fld(tt));
52      fld1=read_bin(fldIn.fil,rec_fld(tt+1));      fld1=mygrid.mskC.*read_bin(fldIn.fil,rec_fld(tt+1));
53      ndim=length(size(fld0{1}));      ndim=length(size(fld0{1}));
54      fld=cat(ndim+1,fld0,fld1);      fld=cat(ndim+1,fld0,fld1);
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            if strcmp(fldIn.tim,'std');
61                arr2=arr2.*randn(size(arr2));
62            end;
63      end;              end;        
64      %          %    
65      k0=floor(tim_prof(ii)); k1=k0+1;      k0=floor(tim_prof(ii)); k1=k0+1;
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22