1 |
|
2 |
# Diagnostic Package Choices |
3 |
#----------------- |
4 |
# for each output-stream: |
5 |
# filename(n) : prefix of the output file name (only 8.c long) for outp.stream n |
6 |
# frequency(n):< 0 : write snap-shot output every multiple of |frequency| (iter) |
7 |
# > 0 : write time-average output every multiple of frequency (iter) |
8 |
# levels(:,n) : list of levels to write to file (Notes: declared as REAL) |
9 |
# when this entry is missing, select all common levels of this list |
10 |
# fields(:,n) : list of diagnostics fields (8.c) (see "available_diagnostics" file |
11 |
# for the list of all available diag. in this particular config) |
12 |
#-------------------------------------------------------------------- |
13 |
# |
14 |
&diagnostics_list |
15 |
# |
16 |
dumpatlast = .TRUE., |
17 |
#--- |
18 |
frequency(1) = 2635200.0, |
19 |
fields(1:24,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 |
'SIuice ','SIvice ', |
29 |
filename(1) = 'diags/state_2d_set1', |
30 |
#--- |
31 |
frequency(2) = 2635200.0, |
32 |
fields(1:17,2) = 'THETA ','SALT ', 'TRAC01 ', |
33 |
'TRAC02 ', 'TRAC03 ', 'TRAC04 ', 'TRAC05 ', |
34 |
'TRAC06 ', 'TRAC07 ', 'TRAC08 ', 'TRAC09 ', |
35 |
'TRAC10 ', 'TRAC11 ', 'TRAC12 ', 'TRAC13 ', |
36 |
'TRAC14 ', 'TRAC15 ', |
37 |
#stuff that is not quite state variables (and may not be quite |
38 |
#synchronous) but are added here to reduce file number |
39 |
# 'DRHODR ', |
40 |
# 'KPPviscA','KPPdiffT','KPPdiffTS', |
41 |
# 'GGL90TKE','GGL90Lmx','GGL90Prl', |
42 |
# 'GGL90ArU','GGL90ArV','GGL90Kr ', |
43 |
filename(2) = 'diags/state_3d_set1', |
44 |
#--- |
45 |
frequency(3) = 2635200.0, |
46 |
fields(1:3,3) = 'UVELMASS','VVELMASS','WVELMASS', |
47 |
# 'GM_PsiX ','GM_PsiY ', |
48 |
#uvel is approximately uvelmass divided by initial hfac, so we dont output the former |
49 |
# 'UVEL ','VVEL ','WVEL ', |
50 |
#full 3D temperature fluxes : 'DFxE_TH ','DFyE_TH ','DFrE_TH ','DFrI_TH ','ADVx_TH ','ADVy_TH ','ADVr_TH ', |
51 |
#but for present computations I only need the vertically integrated horizontal components (see trsp_3d_set2) |
52 |
filename(3) = 'diags/trsp_3d_set1', |
53 |
#--- |
54 |
# frequency(4) = 2635200.0, |
55 |
# fields(1:8,4) = 'DFxE_TH ','DFyE_TH ','ADVx_TH ','ADVy_TH ', |
56 |
# 'DFxE_SLT','DFyE_SLT','ADVx_SLT','ADVy_SLT', |
57 |
# filename(4) = 'diags/trsp_3d_set2', |
58 |
#vertically integrate fields we only use to compute vertically integr. |
59 |
#meridional transports (also omit vertical transports, both to save space) |
60 |
#--- |
61 |
# frequency(5) = -2635200.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) = -2635200.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(4) = 2635200.0, |
75 |
fields(1:6,4) = 'oceFWflx','SIatmFW ','TFLUX ','SItflux ', |
76 |
'SFLUX ','oceQsw ', |
77 |
filename(4) = 'diags/budg2d_zflux_set1', |
78 |
fileFlags(4) = 'D ', |
79 |
#--- |
80 |
frequency(5) = 2635200.0, |
81 |
fields(1:21,5) ='EXFhs ', 'EXFhl ', 'EXFlwnet', 'EXFswnet', 'EXFlwdn ', |
82 |
'EXFswdn ','EXFqnet ','EXFtaux ','EXFtauy ','EXFuwind', |
83 |
'EXFvwind','EXFwspee','EXFatemp','EXFaqh ','EXFevap ', |
84 |
'EXFpreci','EXFsnow ','EXFempmr','EXFpress','EXFroff ', |
85 |
'EXFroft ', |
86 |
filename(5) = 'diags/EXF_forcing_set1', |
87 |
fileFlags(5) = 'D ', |
88 |
#--- |
89 |
frequency(6) = 2635200.0, |
90 |
fields(1:2,6) = 'RHOAnoma','PHIHYD ', |
91 |
filename(6) = 'diags/state_3d_set2', |
92 |
/ |
93 |
# |
94 |
# |
95 |
# Parameter for Diagnostics of per level statistics: |
96 |
#----------------- |
97 |
# for each output-stream: |
98 |
# stat_fname(n) : prefix of the output file name (only 8.c long) for outp.stream n |
99 |
# stat_freq(n):< 0 : write snap-shot output every |stat_freq| seconds |
100 |
# > 0 : write time-average output every stat_freq seconds |
101 |
# stat_phase(n) : write at time = stat_phase + multiple of |stat_freq| |
102 |
# stat_region(:,n) : list of "regions" (default: 1 region only=global) |
103 |
# stat_fields(:,n) : list of diagnostics fields (8.c) (see "available_diagnostics.log" |
104 |
# file for the list of all available diag. in this particular config) |
105 |
#----------------- |
106 |
&DIAG_STATIS_PARMS |
107 |
# diagSt_regMaskFile='basin_masks_eccollc_90x50.bin', |
108 |
# nSetRegMskFile=1, |
109 |
# set_regMask(1)= 1, 1, 1, 1, 1, 1, 1, 1, 1, |
110 |
# 1, 1, 1, 1, 1, 1, 1, 1 |
111 |
# val_regMask(1)= 1., 2., 3., 4., 5., 6., 7., 8., 9., |
112 |
# 10.,11.,12.,13.,14.,15.,16.,17. |
113 |
##--- |
114 |
# stat_fields(1,1)= 'ETAN ','ETANSQ ','DETADT2 ', |
115 |
# 'UVEL ','VVEL ','WVEL ', |
116 |
# 'THETA ','SALT ', |
117 |
# stat_fname(1)= 'dynStDiag', |
118 |
# stat_freq(1)= 3153600., |
119 |
# stat_region(1,1)= 1, 2, 3, 4, 5, 6, 7, 8, 9, |
120 |
# 10,11,12,13,14,15,16,17 |
121 |
##--- |
122 |
# stat_fields(1,2)= 'oceTAUX ','oceTAUY ', |
123 |
# 'surForcT','surForcS','TFLUX ','SFLUX ', |
124 |
# 'oceQnet ','oceSflux','oceFWflx', |
125 |
# stat_fname(2)= 'surfStDiag', |
126 |
# stat_freq(2)= 3153600., |
127 |
# stat_region(1,2)= 1, 2, 3, 4, 5, 6, 7, 8, 9, |
128 |
# 10,11,12,13,14,15,16,17 |
129 |
/ |