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