matlab_class/gcmfaces_diags/diags_set_user.m

%userStep is set in basic_diags_ecco in 'user' blocs
%this routine is by choice not a function, in
%order to have access to entire workspace)

if ~strcmp(setDiags,'user')
    error(['unknown setDiags ' setDiags]);
end;

%%%%%%%%%%%%%%%%%%TEMPLATE%%%%%%%%%%%%%%%%%%%

if strcmp(setDiags,'user')&userStep==1;%diags to be computed
    fprintf('this is a place-holder for an extra set of diags : to be defined by you if you please\n')
    fprintf('how to proceed : 1) define listDiags HERE 2) define listFlds etc. below \n');
    fprintf('how to proceed : 3) add the computation code below 4) add the display code below\n');
    fprintf('how to proceed : 5) once done, remove this print the ''return'' call below\n');
    return;
end;

if strcmp(setDiags,'user')&userStep==2;%input files and variables
    fprintf('this is a place-holder for an extra set of diags : to be defined by you if you please\n')
    fprintf('how to proceed : 1) define listDiags above 2) define listFlds etc. HERE \n');
    fprintf('how to proceed : 3) add the computation code below 4) add the display code below\n');
    fprintf('how to proceed : 5) once done, remove this print the ''return'' call below\n');
    
    listSubdirs={[dirMat 'BUDG/' ],[dirModel 'diags/BUDG/' ],[dirModel 'diags/OTHER/' ],...
        [dirModel 'diags/STATE/' ],[dirModel 'diags/TRSP/'],[dirModel 'diags/' ]};
    listFiles={listFiles{:},'diags_2d_set1','diags_ice_set1','diags_3d_set1','diags_3d_set2','diags_3d_set3'};%for backward compatibility
    return;
end;

if strcmp(setDiags,'user')&userStep==3;%computational part;
    fprintf('this is a place-holder for an extra set of diags : to be defined by you if you please\n')
    fprintf('how to proceed : 1) define listDiags above 2) define listFlds etc. above \n');
    fprintf('how to proceed : 3) add the computation code HERE 4) add the display code below\n');
    fprintf('how to proceed : 5) once done, remove this print the ''return'' call below\n');
    return;
end;

if strcmp(setDiags,'user')&userStep==-1;%computational part;
    fprintf('this is a place-holder for an extra set of diags : to be defined by you if you please\n')
    fprintf('how to proceed : 1) define listDiags above 2) define listFlds etc. above \n');
    fprintf('how to proceed : 3) add the computation code above 4) add the display code below\n');
    fprintf('how to proceed : 5) once done, remove this print the ''return'' call below\n');
    return;
end;

1	gforget	1.2	%userStep is set in basic_diags_ecco in 'user' blocs
2	gforget	1.1	%this routine is by choice not a function, in
3			%order to have access to entire workspace)
4
5	gforget	1.2	if ~strcmp(setDiags,'user')
6	gforget	1.1	error(['unknown setDiags ' setDiags]);
7			end;
8
9			%%%%%%%%%%%%%%%%%%TEMPLATE%%%%%%%%%%%%%%%%%%%
10
11	gforget	1.2	if strcmp(setDiags,'user')&userStep==1;%diags to be computed
12	gforget	1.1	fprintf('this is a place-holder for an extra set of diags : to be defined by you if you please\n')
13			fprintf('how to proceed : 1) define listDiags HERE 2) define listFlds etc. below \n');
14	gforget	1.2	fprintf('how to proceed : 3) add the computation code below 4) add the display code below\n');
15			fprintf('how to proceed : 5) once done, remove this print the ''return'' call below\n');
16	gforget	1.1	return;
17			end;
18
19	gforget	1.2	if strcmp(setDiags,'user')&userStep==2;%input files and variables
20	gforget	1.1	fprintf('this is a place-holder for an extra set of diags : to be defined by you if you please\n')
21			fprintf('how to proceed : 1) define listDiags above 2) define listFlds etc. HERE \n');
22	gforget	1.2	fprintf('how to proceed : 3) add the computation code below 4) add the display code below\n');
23			fprintf('how to proceed : 5) once done, remove this print the ''return'' call below\n');
24	gforget	1.1
25			listSubdirs={[dirMat 'BUDG/' ],[dirModel 'diags/BUDG/' ],[dirModel 'diags/OTHER/' ],...
26			[dirModel 'diags/STATE/' ],[dirModel 'diags/TRSP/'],[dirModel 'diags/' ]};
27			listFiles={listFiles{:},'diags_2d_set1','diags_ice_set1','diags_3d_set1','diags_3d_set2','diags_3d_set3'};%for backward compatibility
28			return;
29			end;
30
31	gforget	1.2	if strcmp(setDiags,'user')&userStep==3;%computational part;
32			fprintf('this is a place-holder for an extra set of diags : to be defined by you if you please\n')
33			fprintf('how to proceed : 1) define listDiags above 2) define listFlds etc. above \n');
34			fprintf('how to proceed : 3) add the computation code HERE 4) add the display code below\n');
35			fprintf('how to proceed : 5) once done, remove this print the ''return'' call below\n');
36			return;
37			end;
38
39			if strcmp(setDiags,'user')&userStep==-1;%computational part;
40	gforget	1.1	fprintf('this is a place-holder for an extra set of diags : to be defined by you if you please\n')
41			fprintf('how to proceed : 1) define listDiags above 2) define listFlds etc. above \n');
42	gforget	1.2	fprintf('how to proceed : 3) add the computation code above 4) add the display code below\n');
43			fprintf('how to proceed : 5) once done, remove this print the ''return'' call below\n');
44	gforget	1.1	return;
45			end;
46