/[MITgcm]/MITgcm_contrib/gael/comm/course-idma2016/matlab/idma_resample.m
ViewVC logotype

Diff of /MITgcm_contrib/gael/comm/course-idma2016/matlab/idma_resample.m

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

revision 1.1 by gforget, Fri Jan 29 15:29:10 2016 UTC revision 1.2 by gforget, Mon Feb 1 17:56:51 2016 UTC
# Line 1  Line 1 
1    function [profIn]=idma_resample(profIn);
2    %
3    % IDMA_RESAMPLE subsamples ECCO temperature fields to observed locations and times
4    %
5    % Example:
6    %    [p]=idma_float_plot('4900828');
7    %    grid_load;
8    %    p=idma_resample(p);
9    
10  p = genpath('gcmfaces/'); addpath(p);%TBE  p = genpath('gcmfaces/'); addpath(p);%TBE
11  p = genpath('MITprof/'); addpath(p);%TBE  p = genpath('MITprof/'); addpath(p);%TBE
12  p = genpath('matlab/'); addpath(p);%TBE  p = genpath('matlab/'); addpath(p);%TBE
13    
14  grid_load; gcmfaces_global; MITprof_global;  fil='release1/nctiles_climatology/THETA';
   
 profIn=idma_float_plot('4900828');  
 [tmp1,tt]=sort(profIn.prof_date);  
   
 fil=[pwd '/release1/nctiles_climatology/THETA'];  
15  fld=read_nctiles([fil '/THETA'],'THETA');  fld=read_nctiles([fil '/THETA'],'THETA');
16  THETA=mean(fld,4);  THETA=mean(fld,4);
17  THETAfil=fil;  THETAfil=fil;
# Line 25  fldIn.name='prof_THETAclim'; Line 28  fldIn.name='prof_THETAclim';
28  fldIn.fld=[];  fldIn.fld=[];
29  profIn=MITprof_resample(profIn,fldIn);  profIn=MITprof_resample(profIn,fldIn);
30    
31    if isdir('release1/nctiles/THETA/');
32    fldIn.tim='monser';
33    fldIn.fil='release1/nctiles/THETA';
34    fldIn.name='prof_THETA';
35    fldIn.fld=[];
36    profIn=MITprof_resample(profIn,fldIn);
37    end;
38    
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22