1 |
ou.wang |
1.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', 'SSHIB', 'SSHDYN', |
160 |
|
|
filename(23) = 'diags/state_2d_ssh_mon_mean', |
161 |
|
|
#--- |
162 |
|
|
frequency(24) = 86400.0, |
163 |
|
|
fields(1:3,24) = 'SSH', 'SSHIB', 'SSHDYN', |
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', 'OBPNOPAB', |
169 |
|
|
filename(25) = 'diags/state_2d_obp_mon_mean', |
170 |
|
|
#--- |
171 |
|
|
frequency(26) = 86400.0, |
172 |
|
|
fields(1:3,26) = 'OBP', 'OBPNOPAB', |
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 |
|
|
/ |