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	%userStep is set in basic_diags_ecco in 'user' blocs
2	%this routine is by choice not a function, in
3	%order to have access to entire workspace)
4
5	if ~strcmp(setDiags,'user')
6	error(['unknown setDiags ' setDiags]);
7	end;
8
9	%%%%%%%%%%%%%%%%%%TEMPLATE%%%%%%%%%%%%%%%%%%%
10
11	if strcmp(setDiags,'user')&userStep==1;%diags to be computed
12	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	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	return;
17	end;
18
19	if strcmp(setDiags,'user')&userStep==2;%input files and variables
20	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	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
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	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	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	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	return;
45	end;
46