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

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

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

Revision 1.1 - (show annotations) (download)
Mon May 14 22:08:41 2001 UTC (23 years ago) by heimbach
Branch: MAIN
CVS Tags: checkpoint44e_post, checkpoint46l_post, checkpoint46g_pre, checkpoint46f_post, checkpoint44f_post, checkpoint46b_post, checkpoint43a-release1mods, checkpoint40pre3, checkpoint40pre1, checkpoint40pre7, checkpoint40pre6, checkpoint40pre9, checkpoint40pre8, checkpoint46l_pre, chkpt44d_post, release1_p1, release1_p2, release1_p3, release1_p4, release1_p5, release1_p6, checkpoint44e_pre, release1_b1, checkpoint43, release1_chkpt44d_post, checkpoint46d_pre, checkpoint40pre2, release1-branch_tutorials, checkpoint45d_post, checkpoint46j_pre, chkpt44a_post, checkpoint44h_pre, checkpoint40pre4, checkpoint46a_post, checkpoint46j_post, checkpoint46k_post, chkpt44c_pre, checkpoint45a_post, ecco_c44_e19, ecco_c44_e18, ecco_c44_e17, ecco_c44_e16, checkpoint44g_post, checkpoint46e_pre, checkpoint45b_post, checkpoint46b_pre, release1-branch-end, release1_final_v1, checkpoint46c_pre, checkpoint46, checkpoint44b_post, checkpoint46h_pre, checkpoint46m_post, checkpoint46a_pre, checkpoint45c_post, checkpoint44h_post, checkpoint46g_post, checkpoint39, checkpoint40pre5, chkpt44a_pre, checkpoint46i_post, ecco_c44_e20, ecco_c44_e21, checkpoint46c_post, ecco-branch-mod1, ecco-branch-mod2, ecco-branch-mod3, ecco-branch-mod4, ecco-branch-mod5, checkpoint46e_post, release1_beta1, checkpoint44b_pre, checkpoint42, checkpoint40, checkpoint41, checkpoint44, checkpoint45, checkpoint46h_post, chkpt44c_post, checkpoint44f_pre, checkpoint46d_post, release1-branch_branchpoint
Branch point for: release1_final, release1-branch, release1, ecco-branch, release1_coupled 
Added external forcing package.
Not presently supported by mitgcm, i.e. disabled by default.
1 c $Header: /u/gcmpack/development/heimbach/ecco_env/pkg/exf/exf_summary.F,v 1.4 2001/02/02 19:43:48 heimbach Exp $
2
3 #include "EXF_CPPOPTIONS.h"
4
5
6 subroutine exf_Summary(
7 I mythid
8 & )
9
10 c ==================================================================
11 c SUBROUTINE exf_Summary
12 c ==================================================================
13 c
14 c o List all the settings of the external forcing.
15 c
16 c started: Christian Eckert eckert@mit.edu 11-Jan-1999
17 c
18 c changed: Christian Eckert eckert@mit.edu 12-Feb-2000
19 c
20 c - Changed Routine names (package prefix: exf_)
21 c
22 c Patrick Heimbach, heimbach@mit.edu 04-May-2000
23 c
24 c - changed the handling of precip and sflux with respect
25 c to CPP options ALLOW_BULKFORMULAE and ALLOW_ATM_TEMP
26 c
27 c ==================================================================
28 c SUBROUTINE exf_Summary
29 c ==================================================================
30
31 implicit none
32
33 c == global variables ==
34
35 #include "EEPARAMS.h"
36 #include "SIZE.h"
37 #include "cal.h"
38 #include "exf.h"
39 #include "exf_param.h"
40
41 c == routine arguments ==
42
43 c mythid - thread number for this instance of the routine.
44
45 integer mythid
46
47 c == local variables ==
48
49 integer i
50 integer il
51 integer timeint(4)
52
53 character*(max_len_mbuf) msgbuf
54
55 c == external ==
56
57 integer ilnblnk
58 external ilnblnk
59
60 c == end of interface ==
61
62 write(msgbuf,'(a)')
63 &' '
64 call print_message( msgbuf, standardmessageunit,
65 & SQUEEZE_RIGHT , mythid)
66 write(msgbuf,'(a)')
67 &'// ======================================================='
68 call print_message( msgbuf, standardmessageunit,
69 & SQUEEZE_RIGHT , mythid)
70 write(msgbuf,'(a)')
71 &'// External forcing configuration >>> START <<<'
72 call print_message( msgbuf, standardmessageunit,
73 & SQUEEZE_RIGHT , mythid)
74 write(msgbuf,'(a)')
75 &'// ======================================================='
76 call print_message( msgbuf, standardmessageunit,
77 & SQUEEZE_RIGHT , mythid)
78 write(msgbuf,'(a)')
79 &' '
80 call print_message( msgbuf, standardmessageunit,
81 & SQUEEZE_RIGHT , mythid)
82 write(msgbuf,'(a,a)')
83 &'External forcing version: ',externalforcingversion
84 call print_message( msgbuf, standardmessageunit,
85 & SQUEEZE_RIGHT , mythid)
86 write(msgbuf,'(a,a)')
87 &'Uses Calendar version: ',usescalendarversion
88 call print_message( msgbuf, standardmessageunit,
89 & SQUEEZE_RIGHT , mythid)
90 write(msgbuf,'(a)')
91 &' '
92 call print_message( msgbuf, standardmessageunit,
93 & SQUEEZE_RIGHT , mythid)
94
95 c For each data set used the summary prints the calendar data
96 c and the corresponding file from which the data will be read.
97
98 #ifdef ALLOW_BULKFORMULAE
99 write(msgbuf,'(a)')
100 &'// ALLOW_BULKFORMULAE: defined'
101 call print_message( msgbuf, standardmessageunit,
102 & SQUEEZE_RIGHT , mythid)
103 #ifdef ALLOW_ATM_TEMP
104 write(msgbuf,'(a)')
105 &'// '
106 call print_message( msgbuf, standardmessageunit,
107 & SQUEEZE_RIGHT , mythid)
108 write(msgbuf,'(a)')
109 &'// ALLOW_ATM_TEMP: defined'
110 call print_message( msgbuf, standardmessageunit,
111 & SQUEEZE_RIGHT , mythid)
112 #else
113 write(msgbuf,'(a)')
114 &'// '
115 call print_message( msgbuf, standardmessageunit,
116 & SQUEEZE_RIGHT , mythid)
117 write(msgbuf,'(a)')
118 &'// ALLOW_ATM_TEMP: NOT defined'
119 call print_message( msgbuf, standardmessageunit,
120 & SQUEEZE_RIGHT , mythid)
121 #endif
122 #ifdef ALLOW_ATM_WIND
123 write(msgbuf,'(a)')
124 &'// ALLOW_ATM_WIND: defined'
125 call print_message( msgbuf, standardmessageunit,
126 & SQUEEZE_RIGHT , mythid)
127 #else
128 write(msgbuf,'(a)')
129 &'// ALLOW_ATM_WIND: NOT defined'
130 call print_message( msgbuf, standardmessageunit,
131 & SQUEEZE_RIGHT , mythid)
132 #endif
133 #else
134 write(msgbuf,'(a)')
135 &'// ALLOW_BULKFORMULAE: NOT defined'
136 call print_message( msgbuf, standardmessageunit,
137 & SQUEEZE_RIGHT , mythid)
138 #endif
139
140 #ifdef ALLOW_KPP
141 write(msgbuf,'(a)')
142 &'// ALLOW_KPP: defined'
143 call print_message( msgbuf, standardmessageunit,
144 & SQUEEZE_RIGHT , mythid)
145 #else
146 write(msgbuf,'(a)')
147 &'// ALLOW_KPP: NOT defined'
148 call print_message( msgbuf, standardmessageunit,
149 & SQUEEZE_RIGHT , mythid)
150 #endif
151
152 #ifdef ALLOW_BULKFORMULAE
153
154 #ifdef ALLOW_ATM_TEMP
155 c Atmospheric temperature.
156 il = ilnblnk(atempfile)
157 call cal_TimeInterval( atempperiod, 'secs', timeint, mythid )
158
159 write(msgbuf,'(a)')
160 &' '
161 call print_message( msgbuf, standardmessageunit,
162 & SQUEEZE_RIGHT , mythid)
163 write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
164 &' Atmospheric temperature starts at ',
165 & (atempstartdate(i), i=1,2), dayofweek(atempstartdate(4)),
166 & '.'
167 call print_message( msgbuf, standardmessageunit,
168 & SQUEEZE_RIGHT , mythid)
169 write(msgbuf,'(a,i9.8,i7.6)')
170 &' Atmospheric temperature period is ',
171 & (timeint(i), i=1,2)
172 call print_message( msgbuf, standardmessageunit,
173 & SQUEEZE_RIGHT , mythid)
174 write(msgbuf,'(a)')
175 &' Atmospheric temperature is read from file:'
176 call print_message( msgbuf, standardmessageunit,
177 & SQUEEZE_RIGHT , mythid)
178 write(msgbuf,'(a,a,a)')
179 &' >> ',atempfile(1:il),' <<'
180 call print_message( msgbuf, standardmessageunit,
181 & SQUEEZE_RIGHT , mythid)
182
183 c Atmospheric specific humidity.
184 il = ilnblnk(aqhfile)
185 call cal_TimeInterval( aqhperiod, 'secs', timeint, mythid )
186
187 write(msgbuf,'(a)')
188 &' '
189 call print_message( msgbuf, standardmessageunit,
190 & SQUEEZE_RIGHT , mythid)
191 write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
192 &' Atmospheric specific humidity starts at ',
193 & (aqhstartdate(i), i=1,2), dayofweek(aqhstartdate(4)),
194 & '.'
195 call print_message( msgbuf, standardmessageunit,
196 & SQUEEZE_RIGHT , mythid)
197 write(msgbuf,'(a,i9.8,i7.6)')
198 &' Atmospheric specific humidity period is ',
199 & (timeint(i), i=1,2)
200 call print_message( msgbuf, standardmessageunit,
201 & SQUEEZE_RIGHT , mythid)
202 write(msgbuf,'(a)')
203 &' Atmospheric specific humidity is read from file:'
204 call print_message( msgbuf, standardmessageunit,
205 & SQUEEZE_RIGHT , mythid)
206 write(msgbuf,'(a,a,a)')
207 &' >> ',aqhfile(1:il),' <<'
208 call print_message( msgbuf, standardmessageunit,
209 & SQUEEZE_RIGHT , mythid)
210
211 c Short wave radiative flux.
212 il = ilnblnk(swfluxfile)
213 call cal_TimeInterval( swfluxperiod, 'secs', timeint, mythid )
214
215 write(msgbuf,'(a)')
216 &' '
217 call print_message( msgbuf, standardmessageunit,
218 & SQUEEZE_RIGHT , mythid)
219 write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
220 &' Short wave rad. flux forcing starts at ',
221 & (swfluxstartdate(i), i=1,2), dayofweek(swfluxstartdate(4)),
222 & '.'
223 call print_message( msgbuf, standardmessageunit,
224 & SQUEEZE_RIGHT , mythid)
225 write(msgbuf,'(a,i9.8,i7.6)')
226 &' Short wave rad. flux forcing period is ',
227 & (timeint(i), i=1,2)
228 call print_message( msgbuf, standardmessageunit,
229 & SQUEEZE_RIGHT , mythid)
230 write(msgbuf,'(a)')
231 &' Short wave rad. flux forcing is read from file:'
232 call print_message( msgbuf, standardmessageunit,
233 & SQUEEZE_RIGHT , mythid)
234 write(msgbuf,'(a,a,a)')
235 &' >> ',swfluxfile(1:il),' <<'
236 call print_message( msgbuf, standardmessageunit,
237 & SQUEEZE_RIGHT , mythid)
238
239 c Long wave radiative flux.
240 il = ilnblnk(lwfluxfile)
241 call cal_TimeInterval( lwfluxperiod, 'secs', timeint, mythid )
242
243 write(msgbuf,'(a)')
244 &' '
245 call print_message( msgbuf, standardmessageunit,
246 & SQUEEZE_RIGHT , mythid)
247 write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
248 &' Long wave rad. flux forcing starts at ',
249 & (lwfluxstartdate(i), i=1,2), dayofweek(lwfluxstartdate(4)),
250 & '.'
251 call print_message( msgbuf, standardmessageunit,
252 & SQUEEZE_RIGHT , mythid)
253 write(msgbuf,'(a,i9.8,i7.6)')
254 &' Long wave rad. flux forcing period is ',
255 & (timeint(i), i=1,2)
256 call print_message( msgbuf, standardmessageunit,
257 & SQUEEZE_RIGHT , mythid)
258 write(msgbuf,'(a)')
259 &' Long wave rad. flux forcing is read from file:'
260 call print_message( msgbuf, standardmessageunit,
261 & SQUEEZE_RIGHT , mythid)
262 write(msgbuf,'(a,a,a)')
263 &' >> ',lwfluxfile(1:il),' <<'
264 call print_message( msgbuf, standardmessageunit,
265 & SQUEEZE_RIGHT , mythid)
266
267 c Precipitation.
268 il = ilnblnk(precipfile)
269 call cal_TimeInterval( precipperiod, 'secs', timeint, mythid )
270
271 write(msgbuf,'(a)')
272 &' '
273 call print_message( msgbuf, standardmessageunit,
274 & SQUEEZE_RIGHT , mythid)
275 write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
276 &' Precipitation data set starts at ',
277 & (precipstartdate(i), i=1,2), dayofweek(precipstartdate(4)),
278 &'.'
279 call print_message( msgbuf, standardmessageunit,
280 & SQUEEZE_RIGHT , mythid)
281 write(msgbuf,'(a,i9.8,i7.6)')
282 &' Precipitation data period is ',
283 & (timeint(i), i=1,2)
284 call print_message( msgbuf, standardmessageunit,
285 & SQUEEZE_RIGHT , mythid)
286 write(msgbuf,'(a)')
287 &' Precipitation data is read from file: '
288 call print_message( msgbuf, standardmessageunit,
289 & SQUEEZE_RIGHT , mythid)
290 write(msgbuf,'(a,a,a)')
291 &' >> ',precipfile(1:il),' <<'
292 call print_message( msgbuf, standardmessageunit,
293 & SQUEEZE_RIGHT , mythid)
294
295 #else
296 c Heat flux.
297 il = ilnblnk(hfluxfile)
298 call cal_TimeInterval( hfluxperiod, 'secs', timeint, mythid )
299
300 write(msgbuf,'(a)')
301 &' '
302 call print_message( msgbuf, standardmessageunit,
303 & SQUEEZE_RIGHT , mythid)
304 write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
305 &' Heat flux forcing starts at ',
306 & (hfluxstartdate(i), i=1,2), dayofweek(hfluxstartdate(4)),'.'
307 call print_message( msgbuf, standardmessageunit,
308 & SQUEEZE_RIGHT , mythid)
309 write(msgbuf,'(a,i9.8,i7.6)')
310 &' Heat flux forcing period is ',
311 & (timeint(i), i=1,2)
312 call print_message( msgbuf, standardmessageunit,
313 & SQUEEZE_RIGHT , mythid)
314 write(msgbuf,'(a)')
315 &' Heat flux forcing is read from file: '
316 call print_message( msgbuf, standardmessageunit,
317 & SQUEEZE_RIGHT , mythid)
318 write(msgbuf,'(a,a,a)')
319 &' >> ',hfluxfile(1:il),' <<'
320 call print_message( msgbuf, standardmessageunit,
321 & SQUEEZE_RIGHT , mythid)
322
323 c Salt flux.
324 il = ilnblnk(sfluxfile)
325 call cal_TimeInterval( sfluxperiod, 'secs', timeint, mythid )
326
327 write(msgbuf,'(a)')
328 &' '
329 call print_message( msgbuf, standardmessageunit,
330 & SQUEEZE_RIGHT , mythid)
331 write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
332 &' Salt flux forcing starts at ',
333 & (sfluxstartdate(i), i=1,2), dayofweek(sfluxstartdate(4)),'.'
334 call print_message( msgbuf, standardmessageunit,
335 & SQUEEZE_RIGHT , mythid)
336 write(msgbuf,'(a,i9.8,i7.6)')
337 &' Salt flux forcing period is ',
338 & (timeint(i), i=1,2)
339 call print_message( msgbuf, standardmessageunit,
340 & SQUEEZE_RIGHT , mythid)
341 write(msgbuf,'(a)')
342 &' Salt flux forcing is read from file: '
343 call print_message( msgbuf, standardmessageunit,
344 & SQUEEZE_RIGHT , mythid)
345 write(msgbuf,'(a,a,a)')
346 &' >> ',sfluxfile(1:il),' <<'
347 call print_message( msgbuf, standardmessageunit,
348 & SQUEEZE_RIGHT , mythid)
349
350 #ifdef ALLOW_KPP
351 c Short wave radiative flux.
352 il = ilnblnk(swfluxfile)
353 call cal_TimeInterval( swfluxperiod, 'secs', timeint, mythid )
354
355 write(msgbuf,'(a)')
356 &' '
357 call print_message( msgbuf, standardmessageunit,
358 & SQUEEZE_RIGHT , mythid)
359 write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
360 &' Short wave rad. flux forcing starts at ',
361 & (swfluxstartdate(i), i=1,2), dayofweek(swfluxstartdate(4)),
362 & '.'
363 call print_message( msgbuf, standardmessageunit,
364 & SQUEEZE_RIGHT , mythid)
365 write(msgbuf,'(a,i9.8,i7.6)')
366 &' Short wave rad. flux forcing period is ',
367 & (timeint(i), i=1,2)
368 call print_message( msgbuf, standardmessageunit,
369 & SQUEEZE_RIGHT , mythid)
370 write(msgbuf,'(a)')
371 &' Short wave rad. flux forcing is read from file:'
372 call print_message( msgbuf, standardmessageunit,
373 & SQUEEZE_RIGHT , mythid)
374 write(msgbuf,'(a,a,a)')
375 &' >> ',swfluxfile(1:il),' <<'
376 call print_message( msgbuf, standardmessageunit,
377 & SQUEEZE_RIGHT , mythid)
378 #endif
379 #endif
380
381 #ifdef ALLOW_ATM_WIND
382 c Zonal wind.
383 il = ilnblnk(uwindfile)
384 call cal_TimeInterval( uwindperiod, 'secs', timeint, mythid )
385
386 write(msgbuf,'(a)')
387 &' '
388 call print_message( msgbuf, standardmessageunit,
389 & SQUEEZE_RIGHT , mythid)
390 write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
391 &' Zonal wind forcing starts at ',
392 & (uwindstartdate(i), i=1,2), dayofweek(uwindstartdate(4)),
393 & '.'
394 call print_message( msgbuf, standardmessageunit,
395 & SQUEEZE_RIGHT , mythid)
396 write(msgbuf,'(a,i9.8,i7.6)')
397 &' Zonal wind forcing period is ',
398 & (timeint(i), i=1,2)
399 call print_message( msgbuf, standardmessageunit,
400 & SQUEEZE_RIGHT , mythid)
401 write(msgbuf,'(a)')
402 &' Zonal wind forcing is read from file:'
403 call print_message( msgbuf, standardmessageunit,
404 & SQUEEZE_RIGHT , mythid)
405 write(msgbuf,'(a,a,a)')
406 &' >> ',uwindfile(1:il),' <<'
407 call print_message( msgbuf, standardmessageunit,
408 & SQUEEZE_RIGHT , mythid)
409
410 c Meridional wind.
411 il = ilnblnk(vwindfile)
412 call cal_TimeInterval( vwindperiod, 'secs', timeint, mythid )
413
414 write(msgbuf,'(a)')
415 &' '
416 call print_message( msgbuf, standardmessageunit,
417 & SQUEEZE_RIGHT , mythid)
418 write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
419 &' Meridional wind forcing starts at ',
420 & (vwindstartdate(i), i=1,2), dayofweek(vwindstartdate(4)),
421 & '.'
422 call print_message( msgbuf, standardmessageunit,
423 & SQUEEZE_RIGHT , mythid)
424 write(msgbuf,'(a,i9.8,i7.6)')
425 &' Meridional wind forcing period is ',
426 & (timeint(i), i=1,2)
427 call print_message( msgbuf, standardmessageunit,
428 & SQUEEZE_RIGHT , mythid)
429 write(msgbuf,'(a)')
430 &' Meridional wind forcing is read from file:'
431 call print_message( msgbuf, standardmessageunit,
432 & SQUEEZE_RIGHT , mythid)
433 write(msgbuf,'(a,a,a)')
434 &' >> ',vwindfile(1:il),' <<'
435 call print_message( msgbuf, standardmessageunit,
436 & SQUEEZE_RIGHT , mythid)
437 #else
438 c Zonal wind stress.
439 il = ilnblnk(ustressfile)
440 call cal_TimeInterval( ustressperiod, 'secs', timeint, mythid )
441
442 write(msgbuf,'(a)')
443 &' '
444 call print_message( msgbuf, standardmessageunit,
445 & SQUEEZE_RIGHT , mythid)
446 write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
447 &' Zonal wind stress forcing starts at ',
448 & (ustressstartdate(i), i=1,2), dayofweek(ustressstartdate(4)),
449 & '.'
450 call print_message( msgbuf, standardmessageunit,
451 & SQUEEZE_RIGHT , mythid)
452 write(msgbuf,'(a,i9.8,i7.6)')
453 &' Zonal wind stress forcing period is ',
454 & (timeint(i), i=1,2)
455 call print_message( msgbuf, standardmessageunit,
456 & SQUEEZE_RIGHT , mythid)
457 write(msgbuf,'(a)')
458 &' Zonal wind stress forcing is read from file:'
459 call print_message( msgbuf, standardmessageunit,
460 & SQUEEZE_RIGHT , mythid)
461 write(msgbuf,'(a,a,a)')
462 &' >> ',ustressfile(1:il),' <<'
463 call print_message( msgbuf, standardmessageunit,
464 & SQUEEZE_RIGHT , mythid)
465
466 c Meridional wind stress.
467 il = ilnblnk(vstressfile)
468 call cal_TimeInterval( vstressperiod, 'secs', timeint, mythid )
469
470 write(msgbuf,'(a)')
471 &' '
472 call print_message( msgbuf, standardmessageunit,
473 & SQUEEZE_RIGHT , mythid)
474 write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
475 &' Meridional wind stress forcing starts at ',
476 & (vstressstartdate(i), i=1,2), dayofweek(vstressstartdate(4)),
477 & '.'
478 call print_message( msgbuf, standardmessageunit,
479 & SQUEEZE_RIGHT , mythid)
480 write(msgbuf,'(a,i9.8,i7.6)')
481 &' Meridional wind stress forcing period is ',
482 & (timeint(i), i=1,2)
483 call print_message( msgbuf, standardmessageunit,
484 & SQUEEZE_RIGHT , mythid)
485 write(msgbuf,'(a)')
486 &' Meridional wind stress forcing is read from file:'
487 call print_message( msgbuf, standardmessageunit,
488 & SQUEEZE_RIGHT , mythid)
489 write(msgbuf,'(a,a,a)')
490 &' >> ',vstressfile(1:il),' <<'
491 call print_message( msgbuf, standardmessageunit,
492 & SQUEEZE_RIGHT , mythid)
493 #endif
494
495 #else
496 c Heat flux.
497 il = ilnblnk(hfluxfile)
498 call cal_TimeInterval( hfluxperiod, 'secs', timeint, mythid )
499
500 write(msgbuf,'(a)')
501 &' '
502 call print_message( msgbuf, standardmessageunit,
503 & SQUEEZE_RIGHT , mythid)
504 write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
505 &' Heat flux forcing starts at ',
506 & (hfluxstartdate(i), i=1,2), dayofweek(hfluxstartdate(4)),'.'
507 call print_message( msgbuf, standardmessageunit,
508 & SQUEEZE_RIGHT , mythid)
509 write(msgbuf,'(a,i9.8,i7.6)')
510 &' Heat flux forcing period is ',
511 & (timeint(i), i=1,2)
512 call print_message( msgbuf, standardmessageunit,
513 & SQUEEZE_RIGHT , mythid)
514 write(msgbuf,'(a)')
515 &' Heat flux forcing is read from file: '
516 call print_message( msgbuf, standardmessageunit,
517 & SQUEEZE_RIGHT , mythid)
518 write(msgbuf,'(a,a,a)')
519 &' >> ',hfluxfile(1:il),' <<'
520 call print_message( msgbuf, standardmessageunit,
521 & SQUEEZE_RIGHT , mythid)
522
523 c Salt flux.
524 il = ilnblnk(sfluxfile)
525 call cal_TimeInterval( sfluxperiod, 'secs', timeint, mythid )
526
527 write(msgbuf,'(a)')
528 &' '
529 call print_message( msgbuf, standardmessageunit,
530 & SQUEEZE_RIGHT , mythid)
531 write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
532 &' Salt flux forcing starts at ',
533 & (sfluxstartdate(i), i=1,2), dayofweek(sfluxstartdate(4)),'.'
534 call print_message( msgbuf, standardmessageunit,
535 & SQUEEZE_RIGHT , mythid)
536 write(msgbuf,'(a,i9.8,i7.6)')
537 &' Salt flux forcing period is ',
538 & (timeint(i), i=1,2)
539 call print_message( msgbuf, standardmessageunit,
540 & SQUEEZE_RIGHT , mythid)
541 write(msgbuf,'(a)')
542 &' Salt flux forcing is read from file: '
543 call print_message( msgbuf, standardmessageunit,
544 & SQUEEZE_RIGHT , mythid)
545 write(msgbuf,'(a,a,a)')
546 &' >> ',sfluxfile(1:il),' <<'
547 call print_message( msgbuf, standardmessageunit,
548 & SQUEEZE_RIGHT , mythid)
549
550 c Zonal wind stress.
551 il = ilnblnk(ustressfile)
552 call cal_TimeInterval( ustressperiod, 'secs', timeint, mythid )
553
554 write(msgbuf,'(a)')
555 &' '
556 call print_message( msgbuf, standardmessageunit,
557 & SQUEEZE_RIGHT , mythid)
558 write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
559 &' Zonal wind stress forcing starts at ',
560 & (ustressstartdate(i), i=1,2), dayofweek(ustressstartdate(4)),
561 & '.'
562 call print_message( msgbuf, standardmessageunit,
563 & SQUEEZE_RIGHT , mythid)
564 write(msgbuf,'(a,i9.8,i7.6)')
565 &' Zonal wind stress forcing period is ',
566 & (timeint(i), i=1,2)
567 call print_message( msgbuf, standardmessageunit,
568 & SQUEEZE_RIGHT , mythid)
569 write(msgbuf,'(a)')
570 &' Zonal wind stress forcing is read from file:'
571 call print_message( msgbuf, standardmessageunit,
572 & SQUEEZE_RIGHT , mythid)
573 write(msgbuf,'(a,a,a)')
574 &' >> ',ustressfile(1:il),' <<'
575 call print_message( msgbuf, standardmessageunit,
576 & SQUEEZE_RIGHT , mythid)
577
578 c Meridional wind stress.
579 il = ilnblnk(vstressfile)
580 call cal_TimeInterval( vstressperiod, 'secs', timeint, mythid )
581
582 write(msgbuf,'(a)')
583 &' '
584 call print_message( msgbuf, standardmessageunit,
585 & SQUEEZE_RIGHT , mythid)
586 write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
587 &' Meridional wind stress forcing starts at ',
588 & (vstressstartdate(i), i=1,2), dayofweek(vstressstartdate(4)),
589 & '.'
590 call print_message( msgbuf, standardmessageunit,
591 & SQUEEZE_RIGHT , mythid)
592 write(msgbuf,'(a,i9.8,i7.6)')
593 &' Meridional wind stress forcing period is ',
594 & (timeint(i), i=1,2)
595 call print_message( msgbuf, standardmessageunit,
596 & SQUEEZE_RIGHT , mythid)
597 write(msgbuf,'(a)')
598 &' Meridional wind stress forcing is read from file:'
599 call print_message( msgbuf, standardmessageunit,
600 & SQUEEZE_RIGHT , mythid)
601 write(msgbuf,'(a,a,a)')
602 &' >> ',vstressfile(1:il),' <<'
603 call print_message( msgbuf, standardmessageunit,
604 & SQUEEZE_RIGHT , mythid)
605
606 #ifdef ALLOW_KPP
607 c Short wave radiative flux.
608 il = ilnblnk(swfluxfile)
609 call cal_TimeInterval( swfluxperiod, 'secs', timeint, mythid )
610
611 write(msgbuf,'(a)')
612 &' '
613 call print_message( msgbuf, standardmessageunit,
614 & SQUEEZE_RIGHT , mythid)
615 write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
616 &' Short wave rad. flux forcing starts at ',
617 & (swfluxstartdate(i), i=1,2), dayofweek(swfluxstartdate(4)),
618 & '.'
619 call print_message( msgbuf, standardmessageunit,
620 & SQUEEZE_RIGHT , mythid)
621 write(msgbuf,'(a,i9.8,i7.6)')
622 &' Short wave rad. flux forcing period is ',
623 & (timeint(i), i=1,2)
624 call print_message( msgbuf, standardmessageunit,
625 & SQUEEZE_RIGHT , mythid)
626 write(msgbuf,'(a)')
627 &' Short wave rad. flux forcing is read from file:'
628 call print_message( msgbuf, standardmessageunit,
629 & SQUEEZE_RIGHT , mythid)
630 write(msgbuf,'(a,a,a)')
631 &' >> ',swfluxfile(1:il),' <<'
632 call print_message( msgbuf, standardmessageunit,
633 & SQUEEZE_RIGHT , mythid)
634 #endif
635
636 #endif
637
638 #ifdef ALLOW_KPP
639 write(msgbuf,'(a)')
640 &'// ALLOW_KPP: defined'
641 call print_message( msgbuf, standardmessageunit,
642 & SQUEEZE_RIGHT , mythid)
643 #else
644 write(msgbuf,'(a)')
645 &'// ALLOW_KPP: NOT defined'
646 call print_message( msgbuf, standardmessageunit,
647 & SQUEEZE_RIGHT , mythid)
648 #endif
649
650 write(msgbuf,'(a)')
651 &' '
652 call print_message( msgbuf, standardmessageunit,
653 & SQUEEZE_RIGHT , mythid)
654 write(msgbuf,'(a)')
655 &'// ======================================================='
656 call print_message( msgbuf, standardmessageunit,
657 & SQUEEZE_RIGHT , mythid)
658 write(msgbuf,'(a)')
659 &'// External forcing configuration >>> END <<<'
660 call print_message( msgbuf, standardmessageunit,
661 & SQUEEZE_RIGHT , mythid)
662 write(msgbuf,'(a)')
663 &'// ======================================================='
664 call print_message( msgbuf, standardmessageunit,
665 & SQUEEZE_RIGHT , mythid)
666 write(msgbuf,'(a)')
667 &' '
668 call print_message( msgbuf, standardmessageunit,
669 & SQUEEZE_RIGHT , mythid)
670
671 end
  ViewVC Help
Powered by ViewVC 1.1.22