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