/[MITgcm]/MITgcm_contrib/gael/matlab_class/gcmfaces_init.m
ViewVC logotype

Diff of /MITgcm_contrib/gael/matlab_class/gcmfaces_init.m

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

revision 1.7 by gforget, Wed Dec 22 01:14:30 2010 UTC revision 1.10 by gforget, Tue Jun 21 21:33:40 2011 UTC
# Line 6  global gcmfaces_verbose; %to print notes Line 6  global gcmfaces_verbose; %to print notes
6  if isempty(gcmfaces_verbose); gcmfaces_verbose=1; end;  if isempty(gcmfaces_verbose); gcmfaces_verbose=1; end;
7    
8    
9  global mydir; mydir=[pwd '/'];  global gcmfaces_dir; gcmfaces_dir=[pwd '/'];
10    
11  fid=fopen([mydir 'gcmfaces_path.m'],'wt');  fid=fopen([gcmfaces_dir 'gcmfaces_path.m'],'wt');
12  fprintf(fid,['global mydir; mydir=''' mydir ''';\n']);  fprintf(fid,['global gcmfaces_dir; gcmfaces_dir=''' gcmfaces_dir ''';\n']);
13  fprintf(fid,['addpath ' mydir ';\n']);  fprintf(fid,['addpath ' gcmfaces_dir ';\n']);
14  fprintf(fid,['addpath ' mydir 'gcmfaces_IO/;\n']);  fprintf(fid,['addpath ' gcmfaces_dir 'gcmfaces_IO/;\n']);
15  fprintf(fid,['addpath ' mydir 'gcmfaces_convert/;\n']);  fprintf(fid,['addpath ' gcmfaces_dir 'gcmfaces_convert/;\n']);
16  fprintf(fid,['addpath ' mydir 'gcmfaces_exch/;\n']);  fprintf(fid,['addpath ' gcmfaces_dir 'gcmfaces_exch/;\n']);
17  fprintf(fid,['addpath ' mydir 'gcmfaces_maps/;\n']);  fprintf(fid,['addpath ' gcmfaces_dir 'gcmfaces_maps/;\n']);
18  fprintf(fid,['addpath ' mydir 'gcmfaces_misc/;\n']);  fprintf(fid,['addpath ' gcmfaces_dir 'gcmfaces_misc/;\n']);
19  fprintf(fid,['addpath ' mydir 'gcmfaces_calc/;\n']);  fprintf(fid,['addpath ' gcmfaces_dir 'gcmfaces_calc/;\n']);
20  fprintf(fid,['addpath ' mydir 'gcmfaces_smooth/;\n']);  fprintf(fid,['addpath ' gcmfaces_dir 'gcmfaces_smooth/;\n']);
21  %fprintf(fid,['addpath ' mydir 'gcmfaces_specs/;\n']);  %fprintf(fid,['addpath ' gcmfaces_dir 'gcmfaces_specs/;\n']);
22  %fprintf(fid,['addpath ' mydir 'gcmfaces_legacy/;\n']);  %fprintf(fid,['addpath ' gcmfaces_dir 'gcmfaces_legacy/;\n']);
23  %fprintf(fid,['addpath ' mydir 'gcmfaces_devel/;\n']);  fprintf(fid,['addpath ' gcmfaces_dir 'gcmfaces_devel/;\n']);
24  fprintf(fid,['addpath ' mydir 'sample_analysis/;\n']);  fprintf(fid,['addpath ' gcmfaces_dir 'sample_analysis/;\n']);
25  fprintf(fid,['addpath ' mydir 'sample_processing/;\n']);  fprintf(fid,['addpath ' gcmfaces_dir 'sample_processing/;\n']);
26  fprintf(fid,['addpath ' mydir 'ecco_v4/;\n']);  fprintf(fid,['addpath ' gcmfaces_dir 'ecco_v4/;\n']);
27  fclose(fid);  fclose(fid);
28    
29  if gcmfaces_verbose;  if gcmfaces_verbose;
# Line 39  test0=dir('sample_input'); if isempty(te Line 39  test0=dir('sample_input'); if isempty(te
39    
40  fprintf('\n\n basic gcmfaces test: started... \n');  fprintf('\n\n basic gcmfaces test: started... \n');
41    gcmfaces_path;    gcmfaces_path;
42    global mygrid; mygrid=[]; grid_load([mydir '/sample_input/GRIDv4/'],5);    global mygrid; mygrid=[]; grid_load([gcmfaces_dir '/sample_input/GRIDv4/'],5,'compact');
43    nameFld='DDetan'; tt=[53:78]*336; cc=[0 0.10];    nameFld='DDetan'; tt=[53:78]*336; cc=[0 0.10];
44    fld=rdmds2gcmfaces([mydir '/sample_input/SAMPLEv4/' nameFld],tt,5);    fld=rdmds2gcmfaces([gcmfaces_dir '/sample_input/SAMPLEv4/' nameFld],tt);
45    fld=std(fld,[],3); msk=mygrid.hFacC(:,:,1); fld(find(msk==0))=NaN;    fld=std(fld,[],3); msk=mygrid.hFacC(:,:,1); fld(find(msk==0))=NaN;
46  fprintf(' basic gcmfaces test: completed. \n\n');  fprintf(' basic gcmfaces test: completed. \n\n');
47    
# Line 52  fprintf('\n\n\n***********message from g Line 52  fprintf('\n\n\n***********message from g
52  fprintf(' starting 1st example routine: plot_one_field ... \n');  fprintf(' starting 1st example routine: plot_one_field ... \n');
53  end;  end;
54    
55  plot_one_field;  plot_one_field('v4',0);
56    
57  if gcmfaces_verbose;  if gcmfaces_verbose;
58  fprintf('\n\n\n***********message from gcmfaces_init.m************\n ');  fprintf('\n\n\n***********message from gcmfaces_init.m************\n ');
59  fprintf(' starting 2nd example routine: plot_std_field ... \n');  fprintf(' starting 2nd example routine: plot_std_field ... \n');
60  end;  end;
61    
62  plot_std_field;  plot_std_field('v4');
63    
64  end;  end;
65    
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.10
  ViewVC Help
Powered by ViewVC 1.1.22