26 |
% where netcdf_domain and netcdf_suff are also in global |
% where netcdf_domain and netcdf_suff are also in global |
27 |
% THE DOT IS ADDED IN SUB-PROG, SO AVOID IT IN DEFINITIONS |
% THE DOT IS ADDED IN SUB-PROG, SO AVOID IT IN DEFINITIONS |
28 |
% |
% |
29 |
% Note that Q is not defined with the ratio by SIGMA |
% Note that Q is not initialy defined with the ratio by RHO (see below). |
30 |
% |
% |
31 |
% A simple potential vorticity (splQ) computing is also available. |
% A simple potential vorticity (splQ) computing is also available. |
32 |
% It is defined as: splQ = f. dSIGMATHETA/dz |
% It is defined as: splQ = f. dSIGMATHETA/dz |
33 |
% |
% |
34 |
% 06/07/2006 |
% It's also possible to add a real last step 4 to compute PV as: |
35 |
|
% Q = -1/RHO * [Ox.dSTdx + Oy.dSTdy + (f+Oz).dSTdz] |
36 |
|
% Note that in this case, program loads the PV output from the |
37 |
|
% routine C_compute_potential_vorticity and simply multiply it |
38 |
|
% by: -1/RHO. |
39 |
|
% RHO may be computed with the routine compute_density.m |
40 |
|
% |
41 |
|
% 06/21/2006 |
42 |
% gmaze@mit.edu |
% gmaze@mit.edu |
43 |
% |
% |
44 |
clear |
clear |
54 |
|
|
55 |
% Files are looked for in subdirectory defined by: ./netcdf-files/<snapshot>/ |
% Files are looked for in subdirectory defined by: ./netcdf-files/<snapshot>/ |
56 |
% So let's define the snapshot: |
% So let's define the snapshot: |
57 |
snapshot = '200103'; |
snapshot = ''; |
58 |
|
|
59 |
|
|
60 |
% File's name: |
% File's name: |
65 |
netcdf_THETA = 'THETA'; |
netcdf_THETA = 'THETA'; |
66 |
netcdf_SALTanom = 'SALTanom'; |
netcdf_SALTanom = 'SALTanom'; |
67 |
netcdf_suff = 'nc'; |
netcdf_suff = 'nc'; |
68 |
netcdf_domain = 'north_atlantic'; |
netcdf_domain = 'north_atlantic'; % Must not be empty ! |
69 |
|
|
|
% FLAGS: |
|
70 |
|
|
71 |
|
% FLAGS: |
72 |
% Turn 0/1 the following flag to determine which PV to compute: |
% Turn 0/1 the following flag to determine which PV to compute: |
73 |
wantsplPV = 0; % (turn 1 for simple PV computing) |
wantsplPV = 0; % (turn 1 for simple PV computing) |
|
|
|
74 |
% Turn 0/1 this flag to get online computing informations: |
% Turn 0/1 this flag to get online computing informations: |
75 |
global toshow |
global toshow |
76 |
toshow = 0; |
toshow = 0; |
77 |
|
|
78 |
|
|
79 |
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% COMPUTING PV: |
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% COMPUTING PV: |
80 |
|
if isempty('snapshot') |
81 |
|
disp(' THIS FILE IS AN EXAMPLE OF USE, COPY IT INTO YOUR') |
82 |
|
disp(' WORKING DIRECTORY WITH A DIFFERENT NAME AND') |
83 |
|
disp(' CUSTOMIZE IT !'); |
84 |
|
end |
85 |
|
|
86 |
% STEP 1: |
% STEP 1: |
87 |
% Output netcdf file is: |
% Output netcdf file is: |
88 |
% ./netcdf-files/<snapshot>/SIGMATHETA.<netcdf_domain>.<netcdf_suff> |
% ./netcdf-files/<snapshot>/SIGMATHETA.<netcdf_domain>.<netcdf_suff> |
104 |
% ./netcdf-files/<snapshot>/PV.<netcdf_domain>.<netcdf_suff> |
% ./netcdf-files/<snapshot>/PV.<netcdf_domain>.<netcdf_suff> |
105 |
C_compute_potential_vorticity(snapshot,wantsplPV) |
C_compute_potential_vorticity(snapshot,wantsplPV) |
106 |
|
|
107 |
|
% STEP 4: |
108 |
|
% Output netcdf file is (replace last one): |
109 |
|
% ./netcdf-files/<snapshot>/PV.<netcdf_domain>.<netcdf_suff> |
110 |
|
global netcdf_RHO netcdf_PV |
111 |
|
netcdf_RHO = 'RHO'; |
112 |
|
netcdf_PV = 'PV'; |
113 |
|
D_compute_potential_vorticity(snapshot) |
114 |
|
|
115 |
|
|
116 |
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% THAT'S IT ! |
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% THAT'S IT ! |