1 |
jmc |
1.1 |
# Diagnostic Package Choices |
2 |
jmc |
1.13 |
#-------------------- |
3 |
|
|
# dumpAtLast (logical): always write output at the end of simulation (default=F) |
4 |
|
|
# diag_mnc (logical): write to NetCDF files (default=useMNC) |
5 |
|
|
#--for each output-stream: |
6 |
|
|
# fileName(n) : prefix of the output file name (max 80c long) for outp.stream n |
7 |
jmc |
1.3 |
# frequency(n):< 0 : write snap-shot output every |frequency| seconds |
8 |
|
|
# > 0 : write time-average output every frequency seconds |
9 |
|
|
# timePhase(n) : write at time = timePhase + multiple of |frequency| |
10 |
jmc |
1.13 |
# averagingFreq : frequency (in s) for periodic averaging interval |
11 |
|
|
# averagingPhase : phase (in s) for periodic averaging interval |
12 |
|
|
# repeatCycle : number of averaging intervals in 1 cycle |
13 |
jmc |
1.3 |
# levels(:,n) : list of levels to write to file (Notes: declared as REAL) |
14 |
jmc |
1.13 |
# when this entry is missing, select all common levels of this list |
15 |
|
|
# fields(:,n) : list of selected diagnostics fields (8.c) in outp.stream n |
16 |
|
|
# (see "available_diagnostics.log" file for the full list of diags) |
17 |
|
|
# missing_value(n) : missing value for real-type fields in output file "n" |
18 |
|
|
# fileFlags(n) : specific code (8c string) for output file "n" |
19 |
|
|
#-------------------- |
20 |
|
|
&DIAGNOSTICS_LIST |
21 |
|
|
fields(1:3,1) = 'ETAN ','ETANSQ ','DETADT2 ', |
22 |
|
|
levels(1,1) = 1., |
23 |
|
|
fileName(1) = 'surfDiag', |
24 |
jmc |
1.5 |
frequency(1) = 4500., |
25 |
jmc |
1.13 |
fields(1:20,2) = 'UVEL ','VVEL ','WVEL ','THETA ','SALT ', |
26 |
|
|
'UVELSQ ','VVELSQ ','WVELSQ ','THETASQ ', |
27 |
|
|
'UV_VEL_Z','PHIHYD ','PHIHYDSQ', |
28 |
|
|
'UVELMASS','VVELMASS', |
29 |
|
|
'UTHMASS ','VTHMASS ','WVELTH ', |
30 |
|
|
'USLTMASS','VSLTMASS','WVELSLT ', |
31 |
|
|
levels(1:5,2) = 1.,2.,3.,4.,5., |
32 |
|
|
fileName(2) = 'dynDiag', |
33 |
jmc |
1.3 |
frequency(2) = 4500., |
34 |
jmc |
1.13 |
fields(1:20,3) = 'UFLUX ','VFLUX ', |
35 |
|
|
'TSR ','OLR ','RADSWG ','RADLWG ', |
36 |
|
|
'HFLUX ','EVAP ','PRECON ','PRECLS ', |
37 |
|
|
'CLDFRC ','CLDPRS ','CLDMAS ', |
38 |
|
|
'DRAG ','WINDS ','TS ','QS ', |
39 |
|
|
'ENPREC ','ALBVISDF','DTSIMPL ', |
40 |
|
|
# 'DWNLWG ', |
41 |
|
|
# 'TSRCLR ','OLRCLR ','SWGCLR ','LWGCLR ', |
42 |
|
|
fileName(3) = 'aimDiag', |
43 |
jmc |
1.3 |
frequency(3) = 4500., |
44 |
jmc |
1.13 |
fields(1:9,4) = 'RELHUM ','DIABT ','DIABQ ', |
45 |
|
|
'RADSW ','RADLW ','DTCONV ','DTLS ', |
46 |
|
|
'DQCONV ','DQLS ', |
47 |
|
|
# 'SWCLR ', 'LWCLR ', |
48 |
|
|
fileName(4) = 'aimDDT', |
49 |
jmc |
1.3 |
frequency(4) = 4500., |
50 |
jmc |
1.13 |
fields(1:14,5) = 'SI_Fract','SI_Thick','SI_SnowH', |
51 |
|
|
'SI_Tsrf ','SI_Tice1','SI_Tice2', |
52 |
|
|
'SI_Qice1','SI_Qice2','SIsnwAge', |
53 |
|
|
'SIsnwPrc','SIalbedo', |
54 |
|
|
'SIflx2oc','SIfrw2oc','SIsaltFx', |
55 |
|
|
# 'SIflxAtm','SIfrwAtm', |
56 |
|
|
# 'SItOcMxL','SIsOcMxL', |
57 |
|
|
fileName(5) = 'thSIceDiag', |
58 |
jmc |
1.3 |
frequency(5) = 4500., |
59 |
jmc |
1.4 |
& |
60 |
|
|
|
61 |
jmc |
1.13 |
#-------------------- |
62 |
jmc |
1.4 |
# Parameter for Diagnostics of per level statistics: |
63 |
jmc |
1.13 |
#-------------------- |
64 |
|
|
# diagSt_mnc (logical): write stat-diags to NetCDF files (default=diag_mnc) |
65 |
|
|
# diagSt_regMaskFile : file containing the region-mask to read-in |
66 |
|
|
# nSetRegMskFile : number of region-mask sets within the region-mask file |
67 |
|
|
# set_regMask(i) : region-mask set-index that identifies the region "i" |
68 |
|
|
# val_regMask(i) : region "i" identifier value in the region mask |
69 |
|
|
#--for each output-stream: |
70 |
|
|
# stat_fName(n) : prefix of the output file name (max 80c long) for outp.stream n |
71 |
jmc |
1.4 |
# stat_freq(n):< 0 : write snap-shot output every |stat_freq| seconds |
72 |
|
|
# > 0 : write time-average output every stat_freq seconds |
73 |
|
|
# stat_phase(n) : write at time = stat_phase + multiple of |stat_freq| |
74 |
|
|
# stat_region(:,n) : list of "regions" (default: 1 region only=global) |
75 |
jmc |
1.13 |
# stat_fields(:,n) : list of selected diagnostics fields (8.c) in outp.stream n |
76 |
|
|
# (see "available_diagnostics.log" file for the full list of diags) |
77 |
|
|
#-------------------- |
78 |
jmc |
1.4 |
&DIAG_STATIS_PARMS |
79 |
jmc |
1.12 |
#- regional mask: 3 lat. band: 1 : y <= -24 ; 2 : -24<y<24 ; 3 : 24 <= y |
80 |
|
|
diagSt_regMaskFile='regMask_lat24.bin', |
81 |
jmc |
1.13 |
nSetRegMskFile = 1, |
82 |
|
|
set_regMask(1:3) = 1, 1, 1, |
83 |
|
|
val_regMask(1:3) = 1., 2., 3., |
84 |
jmc |
1.4 |
#- an example just to check the agreement with MONITOR output: |
85 |
jmc |
1.13 |
stat_fields(1:5,1) = 'ETAN ','UVEL ','VVEL ','WVEL ', 'THETA ', |
86 |
|
|
stat_fName(1) = 'dynStDiag', |
87 |
|
|
stat_freq(1) = -3600., |
88 |
|
|
stat_phase(1) = 0., |
89 |
|
|
stat_fields(1:4,4) = 'LdSnowH ', 'GrdSurfT', 'GrdTemp ', 'GrdWater', |
90 |
|
|
stat_fName(4) = 'landStDiag', |
91 |
|
|
stat_freq(4) = -3600., |
92 |
|
|
stat_phase(4) = 0., |
93 |
|
|
stat_fields(1:12,5) = 'SI_Fract', 'SI_Thick', 'SI_SnowH', |
94 |
jmc |
1.7 |
'SI_Tsrf ', 'SI_Tice1', 'SI_Tice2', |
95 |
jmc |
1.12 |
# 'SI_Qice1', 'SI_Qice2', |
96 |
|
|
'SIsnwPrc', 'SIalbedo', 'SIsnwAge', |
97 |
|
|
'SIflx2oc', 'SIfrw2oc', 'SIsaltFx', |
98 |
jmc |
1.13 |
stat_region(1:3,5) = 1, 3, 0, |
99 |
|
|
stat_fName(5) = 'thSIceStDiag', |
100 |
|
|
stat_freq(5) = -3600., |
101 |
|
|
stat_phase(5) = 0., |
102 |
jmc |
1.4 |
& |