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

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

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

Revision 1.1 - (hide annotations) (download)
Mon May 14 22:08:41 2001 UTC (23 years, 1 month 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 heimbach 1.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