/[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.2 by gforget, Sun May 2 21:28:52 2010 UTC revision 1.7 by gforget, Wed Dec 22 01:14:30 2010 UTC
# Line 1  Line 1 
1    
2    global gcmfaces_skipplottest; %externally set key that will bypass the interactive plotting test
3    if isempty(gcmfaces_skipplottest); gcmfaces_skipplottest=0; end;
4    
5    global gcmfaces_verbose; %to print notes to screen
6    if isempty(gcmfaces_verbose); gcmfaces_verbose=1; end;
7    
8    
9  global mydir; mydir=[pwd '/'];  global mydir; mydir=[pwd '/'];
10    
11  fid=fopen([mydir 'gcmfaces_path.m'],'wt');  fid=fopen([mydir 'gcmfaces_path.m'],'wt');
12  fprintf(fid,['global mydir; mydir=''' mydir ''';\n']);  fprintf(fid,['global mydir; mydir=''' mydir ''';\n']);
13  fprintf(fid,['addpath ' mydir ';\n']);  fprintf(fid,['addpath ' mydir ';\n']);
 fprintf(fid,['addpath ' mydir 'gcmfaces_specs/;\n']);  
14  fprintf(fid,['addpath ' mydir 'gcmfaces_IO/;\n']);  fprintf(fid,['addpath ' mydir 'gcmfaces_IO/;\n']);
15  fprintf(fid,['addpath ' mydir 'gcmfaces_convert/;\n']);  fprintf(fid,['addpath ' mydir 'gcmfaces_convert/;\n']);
16  fprintf(fid,['addpath ' mydir 'gcmfaces_exch/;\n']);  fprintf(fid,['addpath ' mydir 'gcmfaces_exch/;\n']);
# Line 12  fprintf(fid,['addpath ' mydir 'gcmfaces_ Line 18  fprintf(fid,['addpath ' mydir 'gcmfaces_
18  fprintf(fid,['addpath ' mydir 'gcmfaces_misc/;\n']);  fprintf(fid,['addpath ' mydir 'gcmfaces_misc/;\n']);
19  fprintf(fid,['addpath ' mydir 'gcmfaces_calc/;\n']);  fprintf(fid,['addpath ' mydir 'gcmfaces_calc/;\n']);
20  fprintf(fid,['addpath ' mydir 'gcmfaces_smooth/;\n']);  fprintf(fid,['addpath ' mydir 'gcmfaces_smooth/;\n']);
21    %fprintf(fid,['addpath ' mydir 'gcmfaces_specs/;\n']);
22  %fprintf(fid,['addpath ' mydir 'gcmfaces_legacy/;\n']);  %fprintf(fid,['addpath ' mydir 'gcmfaces_legacy/;\n']);
23  %fprintf(fid,['addpath ' mydir 'gcmfaces_devel/;\n']);  %fprintf(fid,['addpath ' mydir 'gcmfaces_devel/;\n']);
24  fprintf(fid,['addpath ' mydir 'sample_analysis/;\n']);  fprintf(fid,['addpath ' mydir 'sample_analysis/;\n']);
25  fprintf(fid,['addpath ' mydir 'sample_processing/;\n']);  fprintf(fid,['addpath ' mydir 'sample_processing/;\n']);
26    fprintf(fid,['addpath ' mydir 'ecco_v4/;\n']);
27  fclose(fid);  fclose(fid);
28    
29    if gcmfaces_verbose;
30  fprintf('\n\n\n***********message from gcmfaces_init.m************ \n');  fprintf('\n\n\n***********message from gcmfaces_init.m************ \n');
31  fprintf(' gcmfaces_path.m was created that, when executed, \n');  fprintf(' gcmfaces_path.m was created that, when executed, \n');
32  fprintf(' adds the gcmfaces directories to your path  \n\n\n');  fprintf(' adds the gcmfaces directories to your path  \n\n\n');
33  fprintf(' ultimately you may want to copy it in your startup.m \n');  fprintf(' ultimately you may want to copy it in your startup.m \n');
34    end;
35    
36  gcmfaces_path;  gcmfaces_path;
37    
38    test0=dir('sample_input'); if isempty(test0); fprintf('no sample input data found\n'); return; end;
39    
40    fprintf('\n\n basic gcmfaces test: started... \n');
41      gcmfaces_path;
42      global mygrid; mygrid=[]; grid_load([mydir '/sample_input/GRIDv4/'],5);
43      nameFld='DDetan'; tt=[53:78]*336; cc=[0 0.10];
44      fld=rdmds2gcmfaces([mydir '/sample_input/SAMPLEv4/' nameFld],tt,5);
45      fld=std(fld,[],3); msk=mygrid.hFacC(:,:,1); fld(find(msk==0))=NaN;
46    fprintf(' basic gcmfaces test: completed. \n\n');
47    
48    if ~gcmfaces_skipplottest;
49    
50    if gcmfaces_verbose;
51  fprintf('\n\n\n***********message from gcmfaces_init.m************\n ');  fprintf('\n\n\n***********message from gcmfaces_init.m************\n ');
52  fprintf(' starting 1st example routine: plot_one_field ... \n');  fprintf(' starting 1st example routine: plot_one_field ... \n');
53    end;
54    
55  plot_one_field;  plot_one_field;
56    
57    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;
61    
62  plot_std_field;  plot_std_field;
63    
64  fprintf('\n\n\n***********message from gcmfaces_init.m************\n ');  end;
 fprintf(' >>> initialization of gcmfaces completed correctly \n');  
 fprintf(' >>> you are all set and may now use the gcmfaces package \n');  
 fprintf(' >>> eventually, to avoid running gcmfaces_init over again \n');  
 fprintf('  you may want to copy the gcmfaces_path.m code to your startup.m \n\n\n');  
65    
66    if gcmfaces_verbose;
67    fprintf('\n\n\n***********message from gcmfaces_init.m************\n ');
68    fprintf(' --- initialization of gcmfaces completed correctly \n');
69    fprintf(' --- you are all set and may now use the gcmfaces package \n');
70    fprintf(' --- eventually, to avoid running gcmfaces_init over again \n');
71    fprintf(' --- you may want to copy the gcmfaces_path.m code to your startup.m \n\n\n');
72    end;
73    
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22