/[MITgcm]/MITgcm_contrib/ecco_utils/ecco_v4_release4_devel/input/data.diagnostics.gcmfaces
ViewVC logotype

Annotation of /MITgcm_contrib/ecco_utils/ecco_v4_release4_devel/input/data.diagnostics.gcmfaces

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph


Revision 1.2 - (hide annotations) (download)
Thu Nov 7 22:16:04 2019 UTC (5 years, 8 months ago) by ou.wang
Branch: MAIN
CVS Tags: HEAD
Changes since 1.1: +4 -4 lines
Update code and namelist related to SSH and OBP

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 ou.wang 1.2 fields(1:3,23) = 'SSH', 'SSHIBC', 'SSHNOIBC',
160 ou.wang 1.1 filename(23) = 'diags/state_2d_ssh_mon_mean',
161     #---
162     frequency(24) = 86400.0,
163 ou.wang 1.2 fields(1:3,24) = 'SSH', 'SSHIBC', 'SSHNOIBC',
164 ou.wang 1.1 timePhase(24) = 42300.,
165     filename(24) = 'diags/state_2d_ssh_day_mean',
166     #---
167     frequency(25) = 2635200.0,
168 ou.wang 1.2 fields(1:3,25) = 'OBP', 'OBPGMAP',
169 ou.wang 1.1 filename(25) = 'diags/state_2d_obp_mon_mean',
170     #---
171     frequency(26) = 86400.0,
172 ou.wang 1.2 fields(1:3,26) = 'OBP', 'OBPGMAP',
173 ou.wang 1.1 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     /

  ViewVC Help
Powered by ViewVC 1.1.22