1 |
# Diagnostic Package Choices |
2 |
#----------------- |
3 |
# for each output-stream: |
4 |
# filename(n) : prefix of the output file name (only 8.c long) for outp.stream n |
5 |
# frequency(n):< 0 : write snap-shot output every multiple of |frequency| (iter) |
6 |
# > 0 : write time-average output every multiple of frequency (iter) |
7 |
# levels(:,n) : list of levels to write to file (Notes: declared as REAL) |
8 |
# when this entry is missing, select all common levels of this list |
9 |
# fields(:,n) : list of diagnostics fields (8.c) (see "available_diagnostics" file |
10 |
# for the list of all available diag. in this particular config) |
11 |
#-------------------------------------------------------------------- |
12 |
# |
13 |
&diagnostics_list |
14 |
# |
15 |
dumpatlast = .TRUE., |
16 |
#--- |
17 |
frequency(1) = 2635200.0, |
18 |
fields(1:25,1) = 'ETAN ','SIarea ','SIheff ','SIhsnow ', |
19 |
#stuff that is not quite state variables (and may not be quite |
20 |
#synchroneous) but are added here to reduce number of files |
21 |
'DETADT2 ','PHIBOT ','sIceLoad', |
22 |
'MXLDEPTH','oceSPDep', |
23 |
'SIatmQnt','SIatmFW ','oceQnet ','oceFWflx', |
24 |
'oceTAUX ','oceTAUY ', |
25 |
'ADVxHEFF','ADVyHEFF','DFxEHEFF','DFyEHEFF', |
26 |
'ADVxSNOW','ADVySNOW','DFxESNOW','DFyESNOW', |
27 |
'SIuice ','SIvice ', |
28 |
filename(1) = 'diags/state_2d_set1', |
29 |
#--- |
30 |
frequency(2) = 2635200.0, |
31 |
fields(1:3,2) = 'THETA ','SALT ', |
32 |
#stuff that is not quite state variables (and may not be quite |
33 |
#synchroneous) but are added here to reduce file number |
34 |
'DRHODR ', |
35 |
# 'KPPviscA','KPPdiffT','KPPdiffTS', |
36 |
# 'GGL90TKE','GGL90Lmx','GGL90Prl', |
37 |
# 'GGL90ArU','GGL90ArV','GGL90Kr ', |
38 |
filename(2) = 'diags/state_3d_set1', |
39 |
#--- |
40 |
frequency(3) = 2635200.0, |
41 |
fields(1:5,3) = 'UVELMASS','VVELMASS','WVELMASS', |
42 |
'GM_PsiX ','GM_PsiY ', |
43 |
#uvel is approximately uvelmass divided by initial hfac, so we dont output the former |
44 |
# 'UVEL ','VVEL ','WVEL ', |
45 |
#full 3D temperature fluxes : 'DFxE_TH ','DFyE_TH ','DFrE_TH ','DFrI_TH ','ADVx_TH ','ADVy_TH ','ADVr_TH ', |
46 |
#but for present computations I only need the vertically integrated horizontal components (see trsp_3d_set2) |
47 |
filename(3) = 'diags/trsp_3d_set1', |
48 |
#--- |
49 |
frequency(4) = 2635200.0, |
50 |
fields(1:8,4) = 'DFxE_TH ','DFyE_TH ','ADVx_TH ','ADVy_TH ', |
51 |
'DFxE_SLT','DFyE_SLT','ADVx_SLT','ADVy_SLT', |
52 |
filename(4) = 'diags/trsp_3d_set2', |
53 |
#vertically integrate fields we only use to compute vertically integr. |
54 |
#meridional transports (also omit vertical transports, both to save space) |
55 |
#--- |
56 |
frequency(5) = -2635200.0, |
57 |
fields(1:6,5) = 'ETAN ','SIheff ','SIhsnow ', |
58 |
'SIarea ','sIceLoad','PHIBOT ', |
59 |
filename(5) = 'diags/budg2d_snap_set1', |
60 |
timePhase(5)= 0., |
61 |
fileFlags(5) = 'D ', |
62 |
#--- |
63 |
frequency(6) = -2635200.0, |
64 |
fields(1:2,6) = 'THETA ','SALT ', |
65 |
filename(6) = 'diags/budg2d_snap_set2', |
66 |
timePhase(6)= 0., |
67 |
fileFlags(6) = 'DI ', |
68 |
#--- |
69 |
frequency(7) = 2635200.0, |
70 |
fields(1:7,7) = 'oceFWflx','SIatmFW ','TFLUX ','SItflux ', |
71 |
'SFLUX ','oceQsw ','oceSPflx', |
72 |
filename(7) = 'diags/budg2d_zflux_set1', |
73 |
fileFlags(7) = 'D ', |
74 |
#--- |
75 |
frequency(8) = 2635200.0, |
76 |
fields(1:13,8) ='UVELMASS','VVELMASS', |
77 |
'ADVx_TH ','ADVy_TH ','DFxE_TH ','DFyE_TH ', |
78 |
'ADVx_SLT','ADVy_SLT','DFxE_SLT','DFyE_SLT', |
79 |
#the following are not transports but tendencies |
80 |
'oceSPtnd','AB_gT ','AB_gS ', |
81 |
filename(8) = 'diags/budg2d_hflux_set2', |
82 |
fileFlags(8) = 'DI ', |
83 |
#--- |
84 |
frequency(9) = 2635200.0, |
85 |
fields(1:8,9) ='ADVxHEFF','ADVyHEFF','DFxEHEFF','DFyEHEFF', |
86 |
'ADVxSNOW','ADVySNOW','DFxESNOW','DFyESNOW', |
87 |
filename(9) = 'diags/budg2d_hflux_set1', |
88 |
fileFlags(9) = 'D ', |
89 |
#--- |
90 |
# this one is important because it activates the vertical advection diags |
91 |
frequency(10) = 2635200.0, |
92 |
fields(1:7,10) = 'ADVr_TH ','DFrE_TH ','DFrI_TH ', |
93 |
'ADVr_SLT','DFrE_SLT','DFrI_SLT', |
94 |
'WVELMASS', |
95 |
filename(10) = 'diags/budg2d_zflux_set3_11', |
96 |
levels(1, 10)= 11., |
97 |
fileFlags(10) = 'D ', |
98 |
#--- |
99 |
frequency(11) = 2635200.0, |
100 |
fields(1:10,11) ='SRELAX ','TRELAX ','WTHMASS ','WSLTMASS', |
101 |
'oceSflux','oceQnet ','SIatmQnt', |
102 |
'SIaaflux','SIsnPrcp','SIacSubl', |
103 |
filename(11) = 'diags/budg2d_zflux_set2', |
104 |
fileFlags(11) = 'D ', |
105 |
#--- |
106 |
frequency(12) = 2635200.0, |
107 |
fields(1:13,12) ='UVELMASS','VVELMASS', |
108 |
'ADVx_TH ','ADVy_TH ','DFxE_TH ','DFyE_TH ', |
109 |
'ADVx_SLT','ADVy_SLT','DFxE_SLT','DFyE_SLT', |
110 |
'oceSPtnd','AB_gT ','AB_gS ', |
111 |
filename(12) = 'diags/budg2d_hflux_set3_11', |
112 |
levels(1:40,12) = 11.,12.,13.,14.,15.,16.,17.,18.,19.,20., |
113 |
21.,22.,23.,24.,25.,26.,27.,28.,29.,30., |
114 |
31.,32.,33.,34.,35.,36.,37.,38.,39.,40., |
115 |
41.,42.,43.,44.,45.,46.,47.,48.,49.,50., |
116 |
fileFlags(12) = 'DI ', |
117 |
#--- |
118 |
frequency(13) = -2635200.0, |
119 |
fields(1:2,13) = 'THETA ','SALT ', |
120 |
filename(13) = 'diags/budg2d_snap_set3_11', |
121 |
timePhase(13)= 0., |
122 |
levels(1:40,13) = 11.,12.,13.,14.,15.,16.,17.,18.,19.,20., |
123 |
21.,22.,23.,24.,25.,26.,27.,28.,29.,30., |
124 |
31.,32.,33.,34.,35.,36.,37.,38.,39.,40., |
125 |
41.,42.,43.,44.,45.,46.,47.,48.,49.,50., |
126 |
fileFlags(13) = 'DI ', |
127 |
#--- |
128 |
frequency(14) = 2635200.0, |
129 |
fields(1:8,14) = 'DFxE_TH ','DFyE_TH ','ADVx_TH ','ADVy_TH ', |
130 |
'DFxE_SLT','DFyE_SLT','ADVx_SLT','ADVy_SLT', |
131 |
filename(14) = 'diags/trsp_2d_set1', |
132 |
#vertically integrate fields we only use to compute vertically integr. |
133 |
#meridional transports (also omit vertical transports, both to save space) |
134 |
fileFlags(14) = ' I ', |
135 |
#--- |
136 |
frequency(15) = 2635200.0, |
137 |
fields(1:3,15) = 'RHOAnoma','PHIHYD ','oceSPtnd', |
138 |
filename(15) = 'diags/state_3d_set2', |
139 |
#--- |
140 |
frequency(16) = 2635200.0, |
141 |
fields(1:7,16) = 'ADVr_TH ','DFrE_TH ','DFrI_TH ', |
142 |
'ADVr_SLT','DFrE_SLT','DFrI_SLT','WVELMASS', |
143 |
filename(16) = 'diags/trsp_3d_set3', |
144 |
#--- |
145 |
frequency(17) = 2635200.0, |
146 |
fields(1:7,17) = 'TFLUX ','SItflux ','SFLUX ', |
147 |
'oceQsw ','oceSPflx','SIsnPrcp','SIacSubl', |
148 |
filename(17) = 'diags/state_2d_set2', |
149 |
#--- |
150 |
frequency(20) = 2635200.0, |
151 |
fields(1:7,20) = 'EXFswdn', 'EXFlwdn', 'EXFatemp', 'EXFaqh','EXFpreci', 'EXFtaux','EXFtauy', |
152 |
filename(20) = 'diags/exf_zflux_set1', |
153 |
#--- |
154 |
frequency(22) = 2635200.0, |
155 |
fields(1:8,22) = 'EXFhs', 'EXFhl', 'EXFlwnet', 'EXFswnet', 'EXFlwdn', 'EXFswdn', 'EXFqnet', 'EXFpress', |
156 |
filename(22) = 'diags/exf_zflux_set2', |
157 |
#--- |
158 |
frequency(23) = 2635200.0, |
159 |
fields(1:3,23) = 'SSH', 'SSHIBC', 'SSHNOIBC', |
160 |
filename(23) = 'diags/state_2d_ssh_mon_mean', |
161 |
#--- |
162 |
frequency(24) = 86400.0, |
163 |
fields(1:3,24) = 'SSH', 'SSHIBC', 'SSHNOIBC', |
164 |
timePhase(24) = 42300., |
165 |
filename(24) = 'diags/state_2d_ssh_day_mean', |
166 |
#--- |
167 |
frequency(25) = 2635200.0, |
168 |
fields(1:3,25) = 'OBP', 'OBPGMAP', |
169 |
filename(25) = 'diags/state_2d_obp_mon_mean', |
170 |
#--- |
171 |
frequency(26) = 86400.0, |
172 |
fields(1:3,26) = 'OBP', 'OBPGMAP', |
173 |
timePhase(26) = 42300., |
174 |
filename(26) = 'diags/state_2d_obp_day_mean', |
175 |
|
176 |
|
177 |
|
178 |
|
179 |
#--- |
180 |
# frequency(1) = 2635200.0, |
181 |
# fields(1:7,1) = 'SIpress', |
182 |
# filename(1) = 'diags/SIpress', |
183 |
#--- |
184 |
# frequency(18) = 2635200.0, |
185 |
# fields(1:1,18) = 'PHIHYDcR' |
186 |
# filename(18) = 'diags/state_3d_set3', |
187 |
#--- |
188 |
# frequency(19) = -2635200.0, |
189 |
# fields(1:4,19) = 'THETA ','SALT ', |
190 |
# 'PHIHYDcR ', 'RHOAnoma', |
191 |
# timePhase(19)= 0., |
192 |
# filename(19) = 'diags/state_3d_snap_set1', |
193 |
|
194 |
|
195 |
#--- |
196 |
/ |
197 |
# |
198 |
# |
199 |
# Parameter for Diagnostics of per level statistics: |
200 |
#----------------- |
201 |
# for each output-stream: |
202 |
# stat_fname(n) : prefix of the output file name (only 8.c long) for outp.stream n |
203 |
# stat_freq(n):< 0 : write snap-shot output every |stat_freq| seconds |
204 |
# > 0 : write time-average output every stat_freq seconds |
205 |
# stat_phase(n) : write at time = stat_phase + multiple of |stat_freq| |
206 |
# stat_region(:,n) : list of "regions" (default: 1 region only=global) |
207 |
# stat_fields(:,n) : list of diagnostics fields (8.c) (see "available_diagnostics.log" |
208 |
# file for the list of all available diag. in this particular config) |
209 |
#----------------- |
210 |
&DIAG_STATIS_PARMS |
211 |
# diagSt_regMaskFile='basin_masks_eccollc_90x50.bin', |
212 |
# nSetRegMskFile=1, |
213 |
# set_regMask(1)= 1, 1, 1, 1, 1, 1, 1, 1, 1, |
214 |
# 1, 1, 1, 1, 1, 1, 1, 1 |
215 |
# val_regMask(1)= 1., 2., 3., 4., 5., 6., 7., 8., 9., |
216 |
# 10.,11.,12.,13.,14.,15.,16.,17. |
217 |
##--- |
218 |
# stat_fields(1,1)= 'ETAN ','ETANSQ ','DETADT2 ', |
219 |
# 'UVEL ','VVEL ','WVEL ', |
220 |
# 'THETA ','SALT ', |
221 |
# stat_fname(1)= 'dynStDiag', |
222 |
# stat_freq(1)= 3153600., |
223 |
# stat_region(1,1)= 1, 2, 3, 4, 5, 6, 7, 8, 9, |
224 |
# 10,11,12,13,14,15,16,17 |
225 |
##--- |
226 |
# stat_fields(1,2)= 'oceTAUX ','oceTAUY ', |
227 |
# 'surForcT','surForcS','TFLUX ','SFLUX ', |
228 |
# 'oceQnet ','oceSflux','oceFWflx', |
229 |
# stat_fname(2)= 'surfStDiag', |
230 |
# stat_freq(2)= 3153600., |
231 |
# stat_region(1,2)= 1, 2, 3, 4, 5, 6, 7, 8, 9, |
232 |
# 10,11,12,13,14,15,16,17 |
233 |
/ |