/[MITgcm]/MITgcm_contrib/ecco_darwin/v4_3deg/input/data.diagnostics
ViewVC logotype

Diff of /MITgcm_contrib/ecco_darwin/v4_3deg/input/data.diagnostics

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

revision 1.1 by dimitri, Tue Jan 14 00:34:16 2020 UTC revision 1.4 by dcarroll, Thu Feb 27 17:51:18 2020 UTC
# Line 1  Line 1 
1  # Diagnostic Package Choices  # Diagnostic Package Choices
2  #--------------------  #-----------------
3  #  dumpAtLast (logical): always write output at the end of simulation (default=F)  # for each output-stream:
4  #  diag_mnc   (logical): write to NetCDF files (default=useMNC)  #  filename(n) : prefix of the output file name (only 8.c long) for outp.stream n
5  #--for each output-stream:  #  frequency(n):< 0 : write snap-shot output every multiple of |frequency| (iter)
6  #  fileName(n) : prefix of the output file name (max 80c long) for outp.stream n  #               > 0 : write time-average output every multiple of frequency (iter)
 #  frequency(n):< 0 : write snap-shot output every |frequency| seconds  
 #               > 0 : write time-average output every frequency seconds  
 #  timePhase(n)     : write at time = timePhase + multiple of |frequency|  
 #    averagingFreq  : frequency (in s) for periodic averaging interval  
 #    averagingPhase : phase     (in s) for periodic averaging interval  
 #    repeatCycle    : number of averaging intervals in 1 cycle  
7  #  levels(:,n) : list of levels to write to file (Notes: declared as REAL)  #  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  #                 when this entry is missing, select all common levels of this list
9  #  fields(:,n) : list of selected diagnostics fields (8.c) in outp.stream n  #  fields(:,n) : list of diagnostics fields (8.c) (see "available_diagnostics" file
10  #                (see "available_diagnostics.log" file for the full list of diags)  #                 for the list of all available diag. in this particular config)
11  #  missing_value(n) : missing value for real-type fields in output file "n"  #--------------------------------------------------------------------
12  #  fileFlags(n)     : specific code (8c string) for output file "n"  #
13  #--------------------   &diagnostics_list
14   &DIAGNOSTICS_LIST  #
15    fields(1:10,1) = 'ETAN    ','ETANSQ  ','DETADT2 ','PHIBOT  ','PHIBOTSQ',     dumpatlast = .TRUE.,
16                     'DICTFLX ','DICOFLX ','DICCFLX ','DICPCO2 ','DICPHAV ',  #---
17     levels(1,1) = 1.,    frequency(1) = 2635200.0,
18     fileName(1) = 'surfDiag',     fields(1:25,1) = 'ETAN    ','SIarea  ','SIheff ','SIhsnow ',
19    frequency(1) = 432000.,  #stuff that is not quite state variables (and may not be quite
20    frequency(1) = 172800.,  #synchroneous) but are added here to reduce number of files
21    fields(1:18,2) = 'UVEL    ','VVEL    ','WVEL    ','PHIHYD  ',                   'DETADT2 ','PHIBOT  ','sIceLoad',
22                     'UVELMASS','VVELMASS',                   'MXLDEPTH','oceSPDep',
23                     'UVELSQ  ','VVELSQ  ','UV_VEL_C','THETASQ ',                   'SIatmQnt','SIatmFW ','oceQnet ','oceFWflx',
24                     'THETA   ','UTHMASS ','VTHMASS ','WVELTH  ',                   'oceTAUX ','oceTAUY ',
25                     'SALT    ','USLTMASS','VSLTMASS','WVELSLT ',                   'ADVxHEFF','ADVyHEFF','DFxEHEFF','DFyEHEFF',
26  #  fileName(2) = 'dynDiag',                   'ADVxSNOW','ADVySNOW','DFxESNOW','DFyESNOW',
27    frequency(2) = 0.,                   'SIuice  ','SIvice  ',
28    fields(1:5,3)  = 'DICBIOA ','DICCARB ',     filename(1) = 'diags/state_2d_set1',
29                     'TRAC05  ','ADVxTr05','DFrITr05',  #---
30  #  fileName(3) = 'dicDiag',    frequency(2) = 2635200.0,
31    frequency(3) = 172800.,     fields(1:3,2) = 'THETA   ','SALT    ',
32    fields(1:8,4)  = 'ADVx_TH ','ADVy_TH ','ADVr_TH ',  #stuff that is not quite state variables (and may not be quite
33                     'DFxE_TH ','DFyE_TH ','DFrE_TH ',  #synchronous) but are added here to reduce file number
34                     'DFrI_TH ',                   'DRHODR  ',
35                     'ADVx_SLT',  #                'KPPviscA','KPPdiffT','KPPdiffTS',
36  #  fileName(4) = 'flxDiag',  #                'GGL90TKE','GGL90Lmx','GGL90Prl',
37    frequency(4) = 0.,  #                '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:8,7) = 'oceFWflx','SIatmFW ','TFLUX   ','SItflux ',
71                       'SFLUX   ','oceQsw  ','oceSPflx','DICTFLX ',
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(18) = -2635200.0,
151       fields(1:2,18) = 'SALT    ','TRAC32  ',
152       timePhase(18)= 0.,
153       filename(18) = 'diags/state_3d_snap_set1',
154    #---
155      frequency(19) = 2635200.0,
156       fields(1:8,19) = 'TRAC32  ','ADVxTr32','ADVyTr32','ADVrTr32',
157            'DFxETr32','DFyETr32','DFrETr32','DFrITr32',
158       filename(19) = 'diags/dic_physics',
159    #---
160      frequency(20) = 2635200.0,
161       fields(1:5,20) = 'BCDIC   ','BCDICPIC','BREMIND ',
162            'BREMINP ','DISSCPIC',
163      filename(20) = 'diags/dic_biology',
164    #---
165      filename(21) = 'diags/CO2_FLUX_BUDGET_DEPTH',
166       frequency(21) = 2635200.0,
167       fields(1,21) = 'TRAC32  ','THETA   ','SALT    ',
168                      'TRAC36  ','TRAC01  ','TRAC04  ',
169    #---
170      filename(22) = 'diags/CO2_FLUX_BUDGET_SURF',
171       frequency(22) = 2635200.0,
172       levels(1,22) = 1.,
173       fields(1,22) = 'DICCFLX ','DICTFLX ',
174    #                  'DICPCO2 ','EXFapco2',
175    
176   &   &
177    
178  #--------------------   /
179    #
180    #
181  # Parameter for Diagnostics of per level statistics:  # Parameter for Diagnostics of per level statistics:
182  #--------------------  #-----------------
183  #  diagSt_mnc (logical): write stat-diags to NetCDF files (default=diag_mnc)  # for each output-stream:
184  #  diagSt_regMaskFile : file containing the region-mask to read-in  #  stat_fname(n) : prefix of the output file name (only 8.c long) for outp.stream n
 #  nSetRegMskFile   : number of region-mask sets within the region-mask file  
 #  set_regMask(i)   : region-mask set-index that identifies the region "i"  
 #  val_regMask(i)   : region "i" identifier value in the region mask  
 #--for each output-stream:  
 #  stat_fName(n) : prefix of the output file name (max 80c long) for outp.stream n  
185  #  stat_freq(n):< 0 : write snap-shot output every |stat_freq| seconds  #  stat_freq(n):< 0 : write snap-shot output every |stat_freq| seconds
186  #               > 0 : write time-average output every stat_freq seconds  #               > 0 : write time-average output every stat_freq seconds
187  #  stat_phase(n)    : write at time = stat_phase + multiple of |stat_freq|  #  stat_phase(n)    : write at time = stat_phase + multiple of |stat_freq|
188  #  stat_region(:,n) : list of "regions" (default: 1 region only=global)  #  stat_region(:,n) : list of "regions" (default: 1 region only=global)
189  #  stat_fields(:,n) : list of selected diagnostics fields (8.c) in outp.stream n  #  stat_fields(:,n) : list of diagnostics fields (8.c) (see "available_diagnostics.log"
190  #                (see "available_diagnostics.log" file for the full list of diags)  #                 file for the list of all available diag. in this particular config)
191  #--------------------  #-----------------
192   &DIAG_STATIS_PARMS   &DIAG_STATIS_PARMS
193  #- an example just to check the agreement with MONITOR output:  # diagSt_regMaskFile='basin_masks_eccollc_90x50.bin',
194  #stat_fields(1:5,1)  = 'ETAN    ','UVEL    ','VVEL    ','WVEL    ', 'THETA   ',  # nSetRegMskFile=1,
195  # stat_fName(1) = 'dynStDiag',  # set_regMask(1)= 1,  1,  1,  1,  1,  1,  1,  1,  1,
196  #  stat_freq(1) = -864000.,  #                 1,  1,  1,  1,  1,  1,  1,  1
197  # stat_phase(1) = 0.,  # val_regMask(1)= 1., 2., 3., 4., 5., 6., 7., 8., 9.,
198   &  #                10.,11.,12.,13.,14.,15.,16.,17.
199    ##---
200    # stat_fields(1,1)= 'ETAN    ','ETANSQ  ','DETADT2 ',
201    #                   'UVEL    ','VVEL    ','WVEL    ',
202    #                   'THETA   ','SALT    ',
203    #    stat_fname(1)= 'dynStDiag',
204    #     stat_freq(1)= 3153600.,
205    # stat_region(1,1)=  1, 2, 3, 4, 5, 6, 7, 8, 9,
206    #                   10,11,12,13,14,15,16,17
207    ##---
208    # stat_fields(1,2)= 'oceTAUX ','oceTAUY ',
209    #                   'surForcT','surForcS','TFLUX   ','SFLUX   ',
210    #                   'oceQnet ','oceSflux','oceFWflx',
211    #    stat_fname(2)= 'surfStDiag',
212    #     stat_freq(2)= 3153600.,
213    # stat_region(1,2)=  1, 2, 3, 4, 5, 6, 7, 8, 9,
214    #                   10,11,12,13,14,15,16,17
215     /

Legend:
Removed from v.1.1  
changed lines
  Added in v.1.4

  ViewVC Help
Powered by ViewVC 1.1.22