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

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

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


Revision 1.6 - (show annotations) (download)
Thu Nov 4 15:22:40 2010 UTC (14 years, 8 months ago) by gforget
Branch: MAIN
Changes since 1.5: +30 -6 lines
- introduce basic comp. test (load a field and compute std).
- introduce gcmfaces_skipplottest to only do the basic comp. test.
- add gcmfaces_verbose brackets.

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

  ViewVC Help
Powered by ViewVC 1.1.22