/[MITgcm]/MITgcm/pkg/exf/exf_clim_summary.F
ViewVC logotype

Contents of /MITgcm/pkg/exf/exf_clim_summary.F

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

Revision 1.7 - (show annotations) (download)
Thu Oct 9 04:19:19 2003 UTC (20 years, 7 months ago) by edhill
Branch: MAIN
CVS Tags: checkpoint51k_post, checkpoint51i_post, checkpoint51l_pre
Branch point for: tg2-branch
Changes since 1.6: +2 -2 lines 
 o first check-in for the "branch-genmake2" merge
 o verification suite as run on shelley (gcc 3.2.2):

Wed Oct  8 23:42:29 EDT 2003
                T           S           U           V
G D M    c        m  s        m  s        m  s        m  s
E p a R  g  m  m  e  .  m  m  e  .  m  m  e  .  m  m  e  .
N n k u  2  i  a  a  d  i  a  a  d  i  a  a  d  i  a  a  d
2 d e n  d  n  x  n  .  n  x  n  .  n  x  n  .  n  x  n  .

OPTFILE=NONE

Y Y Y Y 13 16 16 16  0 16 16 16 16 16 16 16 16 13 12  0  0 pass  adjustment.128x64x1
Y Y Y Y 16 16 16 16  0 16 16 16 16 16 16  0  0 16 16  0  0 pass  adjustment.cs-32x32x1
Y Y Y Y 16 16 16 16  0 16 16 16 16 16 16 22  0 16 16 22  0 pass  adjust_nlfs.cs-32x32x1
Y Y Y Y -- 13 13 16 16 13 13 13 13 16 16 16 16 16 16 16 16 N/O   advect_cs
Y Y Y Y -- 22 16 16 16 16 16 16 13 16 16 16 16 16 16 16 16 N/O   advect_xy
Y Y Y Y -- 13 16 13 16 16 16 16 16 16 16 22 16 16 16 16 16 N/O   advect_xz
Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 pass  aim.5l_cs
Y Y Y Y 14 16 16 16 16 16 16 16 16 13 16 16 16 16 16 13 16 pass  aim.5l_Equatorial_Channel
Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 13 16 16 13 13 16 pass  aim.5l_LatLon
Y Y Y Y 13 16 16 16 16 16 16 16 16 16 13 12 13 13 16 13 16 pass  exp0
Y Y Y Y 14 16 16 16 16 16 16 16 22 16 16 16 13 16 16 22 16 pass  exp1
Y Y Y Y 13 13 16 13 16 16 16 16 16 13 13 16 16 13 13 13 13 pass  exp2
Y Y Y Y 16 16 16 16 16 16 16 16 22 16 16 16 16 16 16 16 16 pass  exp4
Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 22 16 16 16 22 16 pass  exp5
Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 pass  front_relax
Y Y Y Y 14 16 16 13 13 16 16 13 13 16 13 13 16 12 13 13 16 pass  global_ocean.90x40x15
Y Y Y Y 10 16 16 13 13 16 13 16 16 13 13 13 13 16 16 13 16 FAIL  global_ocean.cs32x15
Y Y Y Y  6 11 12 13 13 12 13 16 13  9  9  9  9 10  9  9 11 FAIL  global_ocean_pressure
Y Y Y Y 14 16 16 13 16 16 16 13 13 13 13 13 16 12 16 13 16 pass  global_with_exf
Y Y Y Y 14 16 16 16 16 16 16 16 16 11 13 22 13 16 16  9 16 pass  hs94.128x64x5
Y Y Y Y 13 16 16 16 16 16 16 16 16 11 16 16 16 13 16 22 13 pass  hs94.1x64x5
Y Y Y Y 14 16 16 16 16 16 16 16 16 13 16 13 13 16 16 22 13 pass  hs94.cs-32x32x5
Y Y Y Y 10 10 16 13 13 16 16 16 22 16 13 13 13 13 13 22 13 FAIL  ideal_2D_oce
Y Y Y Y  8 16 16 16 16 16 16 16 16 13 13  8 16 16 16 16 16 FAIL  internal_wave
Y Y Y Y 14 16 16 16 16 16 16 16 16 13 13 22 13 13 13 22 16 pass  inverted_barometer
Y Y Y Y 12 16 16 16 16 16 16 16 16 16 13 12 13 13 13 13 13 FAIL  lab_sea
Y Y Y Y 11 16 16 16 16 16 16 16 13 13 13 12 13 16 13 12 13 FAIL  natl_box
Y Y Y Y 16 16 16 16 16 16 16 16 22 16 16 16 16 16 16 16 16 pass  plume_on_slope
Y Y Y Y 13 16 16 16 16 13 16 16 16 16 16 16 16 13 16 16 16 pass  solid-body.cs-32x32x1
1 c $Header: /u/u3/gcmpack/MITgcm/pkg/exf/exf_clim_summary.F,v 1.6.2.1 2003/10/02 18:30:07 adcroft Exp $
2
3 #include "EXF_OPTIONS.h"
4
5 subroutine exf_clim_summary(
6 I mythid
7 & )
8
9 c ==================================================================
10 c SUBROUTINE exf_clim_summary
11 c ==================================================================
12 c
13 c o List all the settings of the external forcing.
14 c
15 c started: Ralf.Giering@FastOpt.de 25-Mai-2000
16 c mods for pkg/seaice: menemenlis@jpl.nasa.gov 20-Dec-2002
17 c
18 c ==================================================================
19 c SUBROUTINE exf_clim_summary
20 c ==================================================================
21
22 implicit none
23
24 c == global variables ==
25
26 #include "EEPARAMS.h"
27 #include "SIZE.h"
28 #include "cal.h"
29 #include "exf.h"
30 #include "exf_param.h"
31 #include "exf_clim_param.h"
32
33 c == routine arguments ==
34
35 c mythid - thread number for this instance of the routine.
36
37 integer mythid
38
39 c == local variables ==
40
41 #if defined(ALLOW_CLIMTEMP_RELAXATION) || \
42 defined(ALLOW_CLIMSALT_RELAXATION) || \
43 defined(ALLOW_CLIMSST_RELAXATION ) || \
44 defined(ALLOW_CLIMSSS_RELAXATION )
45 integer i
46 integer il
47 integer timeint(4)
48 #endif
49
50 character*(max_len_mbuf) msgbuf
51
52 c == external ==
53
54 integer ilnblnk
55 external ilnblnk
56
57 c == end of interface ==
58
59 write(msgbuf,'(a)')
60 &' '
61 call print_message( msgbuf, standardmessageunit,
62 & SQUEEZE_RIGHT , mythid)
63 write(msgbuf,'(a)')
64 &'// ======================================================='
65 call print_message( msgbuf, standardmessageunit,
66 & SQUEEZE_RIGHT , mythid)
67 write(msgbuf,'(a)')
68 &'// External forcing climatology configuration >>> START <<<'
69 call print_message( msgbuf, standardmessageunit,
70 & SQUEEZE_RIGHT , mythid)
71 write(msgbuf,'(a)')
72 &'// ======================================================='
73 call print_message( msgbuf, standardmessageunit,
74 & SQUEEZE_RIGHT , mythid)
75 write(msgbuf,'(a)')
76 &' '
77 call print_message( msgbuf, standardmessageunit,
78 & SQUEEZE_RIGHT , mythid)
79 write(msgbuf,'(a,a)')
80 &'External forcing version: ',externalforcingversion
81 call print_message( msgbuf, standardmessageunit,
82 & SQUEEZE_RIGHT , mythid)
83 write(msgbuf,'(a,a)')
84 &'Uses Calendar version: ',usescalendarversion
85 call print_message( msgbuf, standardmessageunit,
86 & SQUEEZE_RIGHT , mythid)
87 write(msgbuf,'(a)')
88 &' '
89 call print_message( msgbuf, standardmessageunit,
90 & SQUEEZE_RIGHT , mythid)
91
92 c For each data set used the summary prints the calendar data
93 c and the corresponding file from which the data will be read.
94
95 #ifdef ALLOW_CLIMTEMP_RELAXATION
96 write(msgbuf,'(a)')
97 &'// ALLOW_CLIMTEMP_RELAXATION: defined'
98 call print_message( msgbuf, standardmessageunit,
99 & SQUEEZE_RIGHT , mythid)
100 #else
101 write(msgbuf,'(a)')
102 &'// ALLOW_CLIMTEMP_RELAXATION: NOT defined'
103 call print_message( msgbuf, standardmessageunit,
104 & SQUEEZE_RIGHT , mythid)
105 #endif
106
107 #ifdef ALLOW_CLIMSALT_RELAXATION
108 write(msgbuf,'(a)')
109 &'// ALLOW_CLIMSALT_RELAXATION: defined'
110 call print_message( msgbuf, standardmessageunit,
111 & SQUEEZE_RIGHT , mythid)
112 #else
113 write(msgbuf,'(a)')
114 &'// ALLOW_CLIMSALT_RELAXATION: NOT defined'
115 call print_message( msgbuf, standardmessageunit,
116 & SQUEEZE_RIGHT , mythid)
117 #endif
118
119 #ifdef ALLOW_CLIMSST_RELAXATION
120 write(msgbuf,'(a)')
121 &'// ALLOW_CLIMSST_RELAXATION: defined'
122 call print_message( msgbuf, standardmessageunit,
123 & SQUEEZE_RIGHT , mythid)
124 #else
125 write(msgbuf,'(a)')
126 &'// ALLOW_CLIMSST_RELAXATION: NOT defined'
127 call print_message( msgbuf, standardmessageunit,
128 & SQUEEZE_RIGHT , mythid)
129 #endif
130
131 #ifdef ALLOW_CLIMSSS_RELAXATION
132 write(msgbuf,'(a)')
133 &'// ALLOW_CLIMSSS_RELAXATION: defined'
134 call print_message( msgbuf, standardmessageunit,
135 & SQUEEZE_RIGHT , mythid)
136 #else
137 write(msgbuf,'(a)')
138 &'// ALLOW_CLIMSSS_RELAXATION: NOT defined'
139 call print_message( msgbuf, standardmessageunit,
140 & SQUEEZE_RIGHT , mythid)
141 #endif
142
143 c The climatological data sets are assumed to contain monthly
144 c data. This can be changed in a later version to an arbitrary
145 c number of intervals during a given year.
146
147 #ifdef ALLOW_CLIMTEMP_RELAXATION
148 c Relaxation to pot. temperature climatology.
149 il = ilnblnk(climtempfile)
150 call cal_TimeInterval( climtempperiod, 'secs', timeint, mythid )
151
152 write(msgbuf,'(a)')
153 &' '
154 call print_message( msgbuf, standardmessageunit,
155 & SQUEEZE_RIGHT , mythid)
156
157 write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
158 &' Climatological pot. temperature starts at ',
159 & (climtempstartdate(i), i=1,2),
160 & dayofweek(climtempstartdate(4)),'.'
161 call print_message( msgbuf, standardmessageunit,
162 & SQUEEZE_RIGHT , mythid)
163 write(msgbuf,'(a,i9.8,i7.6)')
164 &' Climatological pot. temperature period is ',
165 & (timeint(i), i=1,2)
166 call print_message( msgbuf, standardmessageunit,
167 & SQUEEZE_RIGHT , mythid)
168
169 write(msgbuf,'(a)')
170 &' Climatological pot. temperature is read from file:'
171 call print_message( msgbuf, standardmessageunit,
172 & SQUEEZE_RIGHT , mythid)
173 write(msgbuf,'(a,a,a)')
174 &' >> ',climtempfile(1:il),' <<'
175 call print_message( msgbuf, standardmessageunit,
176 & SQUEEZE_RIGHT , mythid)
177 #endif
178
179 #ifdef ALLOW_CLIMSALT_RELAXATION
180 c Relaxation to salinity climatology.
181 il = ilnblnk(climsaltfile)
182 call cal_TimeInterval( climsaltperiod, 'secs', timeint, mythid )
183
184 write(msgbuf,'(a)')
185 &' '
186 call print_message( msgbuf, standardmessageunit,
187 & SQUEEZE_RIGHT , mythid)
188 write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
189 &' Climatological salinity starts at ',
190 & (climsaltstartdate(i), i=1,2),
191 & dayofweek(climsaltstartdate(4)),'.'
192 call print_message( msgbuf, standardmessageunit,
193 & SQUEEZE_RIGHT , mythid)
194 write(msgbuf,'(a,i9.8,i7.6)')
195 &' Climatological salinity period is ',
196 & (timeint(i), i=1,2)
197 call print_message( msgbuf, standardmessageunit,
198 & SQUEEZE_RIGHT , mythid)
199 write(msgbuf,'(a)')
200 &' Climatological salinity is read from file:'
201 call print_message( msgbuf, standardmessageunit,
202 & SQUEEZE_RIGHT , mythid)
203 write(msgbuf,'(a,a,a)')
204 &' >> ',climsaltfile(1:il),' <<'
205 call print_message( msgbuf, standardmessageunit,
206 & SQUEEZE_RIGHT , mythid)
207 #endif
208
209 #ifdef ALLOW_CLIMSST_RELAXATION
210 c Relaxation to SST climatology.
211 il = ilnblnk(climsstfile)
212 call cal_TimeInterval( climsstperiod, 'secs', timeint, mythid )
213
214 write(msgbuf,'(a)')
215 &' '
216 call print_message( msgbuf, standardmessageunit,
217 & SQUEEZE_RIGHT , mythid)
218 write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
219 &' Climatological SST starts at ',
220 & (climsststartdate(i), i=1,2),
221 & dayofweek(climsststartdate(4)),'.'
222 call print_message( msgbuf, standardmessageunit,
223 & SQUEEZE_RIGHT , mythid)
224 write(msgbuf,'(a,i9.8,i7.6)')
225 &' Climatological SST period is ',
226 & (timeint(i), i=1,2)
227 call print_message( msgbuf, standardmessageunit,
228 & SQUEEZE_RIGHT , mythid)
229 write(msgbuf,'(a)')
230 &' Climatological SST is read from file:'
231 call print_message( msgbuf, standardmessageunit,
232 & SQUEEZE_RIGHT , mythid)
233 write(msgbuf,'(a,a,a)')
234 &' >> ',climsstfile(1:il),' <<'
235 call print_message( msgbuf, standardmessageunit,
236 & SQUEEZE_RIGHT , mythid)
237 #endif
238
239 #ifdef ALLOW_CLIMSSS_RELAXATION
240 c Relaxation to SSS climatology.
241 il = ilnblnk(climsssfile)
242 call cal_TimeInterval( climsssperiod, 'secs', timeint, mythid )
243
244 write(msgbuf,'(a)')
245 &' '
246 call print_message( msgbuf, standardmessageunit,
247 & SQUEEZE_RIGHT , mythid)
248 write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
249 &' Climatological SSS starts at ',
250 & (climsssstartdate(i), i=1,2),
251 & dayofweek(climsssstartdate(4)),'.'
252 call print_message( msgbuf, standardmessageunit,
253 & SQUEEZE_RIGHT , mythid)
254 write(msgbuf,'(a,i9.8,i7.6)')
255 &' Climatological SSS period is ',
256 & (timeint(i), i=1,2)
257 call print_message( msgbuf, standardmessageunit,
258 & SQUEEZE_RIGHT , mythid)
259 write(msgbuf,'(a)')
260 &' Climatological SSS is read from file:'
261 call print_message( msgbuf, standardmessageunit,
262 & SQUEEZE_RIGHT , mythid)
263 write(msgbuf,'(a,a,a)')
264 &' >> ',climsssfile(1:il),' <<'
265 call print_message( msgbuf, standardmessageunit,
266 & SQUEEZE_RIGHT , mythid)
267 #endif
268
269 write(msgbuf,'(a)')
270 &' '
271 call print_message( msgbuf, standardmessageunit,
272 & SQUEEZE_RIGHT , mythid)
273 write(msgbuf,'(a)')
274 &'// ======================================================='
275 call print_message( msgbuf, standardmessageunit,
276 & SQUEEZE_RIGHT , mythid)
277 write(msgbuf,'(a)')
278 &'// External forcing climatology configuration >>> END <<<'
279 call print_message( msgbuf, standardmessageunit,
280 & SQUEEZE_RIGHT , mythid)
281 write(msgbuf,'(a)')
282 &'// ======================================================='
283 call print_message( msgbuf, standardmessageunit,
284 & SQUEEZE_RIGHT , mythid)
285 write(msgbuf,'(a)')
286 &' '
287 call print_message( msgbuf, standardmessageunit,
288 & SQUEEZE_RIGHT , mythid)
289
290 end
  ViewVC Help
Powered by ViewVC 1.1.22