/[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.5 - (hide annotations) (download)
Tue Feb 25 06:35:46 2003 UTC (21 years, 2 months ago) by dimitri
Branch: MAIN
CVS Tags: checkpoint50c_post, checkpoint50c_pre, checkpoint48i_post, checkpoint51, checkpoint50, checkpoint50d_post, checkpoint50b_pre, checkpoint51f_post, checkpoint51d_post, checkpoint51j_post, checkpoint48h_post, checkpoint51b_pre, checkpoint51h_pre, checkpoint50f_post, checkpoint50a_post, checkpoint50f_pre, branchpoint-genmake2, checkpoint51b_post, checkpoint51c_post, checkpoint50g_post, checkpoint50h_post, checkpoint50e_pre, checkpoint50i_post, checkpoint51i_pre, checkpoint50e_post, checkpoint50d_pre, checkpoint51e_post, checkpoint49, checkpoint51f_pre, checkpoint51g_post, checkpoint50b_post, checkpoint51a_post
Branch point for: branch-genmake2
Changes since 1.4: +2 -2 lines 
o Added net flux to downward flux conversion to pkg/exf/exf_getffields.F
o Added SEAICE_initialHEFF to pkg/seaice
1 dimitri 1.5 c $Header: /u/gcmpack/MITgcm/pkg/exf/exf_summary.F,v 1.1.6.4 2003/02/19 06:00:56 dimitri Exp $
2 heimbach 1.1
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 dimitri 1.4 c started: Christian Eckert eckert@mit.edu 11-Jan-1999
17 heimbach 1.1 c
18 dimitri 1.4 c changed: Christian Eckert eckert@mit.edu 12-Feb-2000
19     c - changed routine names (package prefix: exf_)
20 heimbach 1.1 c
21 dimitri 1.4 c changed: Patrick Heimbach heimbach@mit.edu 04-May-2000
22 heimbach 1.1 c - changed the handling of precip and sflux with respect
23     c to CPP options ALLOW_BULKFORMULAE and ALLOW_ATM_TEMP
24     c
25 dimitri 1.4 c changed: Dimitris Menemenlis menemenlis@jpl.nasa.gov 20-Dec-2002
26     c - modifications for using pkg/exf with pkg/seaice
27 dimitri 1.3 c
28 heimbach 1.1 c ==================================================================
29     c SUBROUTINE exf_Summary
30     c ==================================================================
31    
32     implicit none
33    
34     c == global variables ==
35    
36     #include "EEPARAMS.h"
37     #include "SIZE.h"
38     #include "cal.h"
39     #include "exf.h"
40     #include "exf_param.h"
41    
42     c == routine arguments ==
43    
44     c mythid - thread number for this instance of the routine.
45    
46     integer mythid
47    
48     c == local variables ==
49    
50     integer i
51     integer il
52     integer timeint(4)
53    
54     character*(max_len_mbuf) msgbuf
55    
56     c == external ==
57    
58     integer ilnblnk
59     external ilnblnk
60    
61     c == end of interface ==
62    
63     write(msgbuf,'(a)')
64     &' '
65     call print_message( msgbuf, standardmessageunit,
66     & SQUEEZE_RIGHT , mythid)
67     write(msgbuf,'(a)')
68     &'// ======================================================='
69     call print_message( msgbuf, standardmessageunit,
70     & SQUEEZE_RIGHT , mythid)
71     write(msgbuf,'(a)')
72     &'// External forcing configuration >>> START <<<'
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)')
80     &' '
81     call print_message( msgbuf, standardmessageunit,
82     & SQUEEZE_RIGHT , mythid)
83     write(msgbuf,'(a,a)')
84     &'External forcing version: ',externalforcingversion
85     call print_message( msgbuf, standardmessageunit,
86     & SQUEEZE_RIGHT , mythid)
87     write(msgbuf,'(a,a)')
88     &'Uses Calendar version: ',usescalendarversion
89     call print_message( msgbuf, standardmessageunit,
90     & SQUEEZE_RIGHT , mythid)
91     write(msgbuf,'(a)')
92     &' '
93     call print_message( msgbuf, standardmessageunit,
94     & SQUEEZE_RIGHT , mythid)
95    
96 dimitri 1.4 c-- Print settings of some CPP flags.
97 heimbach 1.1
98     #ifdef ALLOW_ATM_TEMP
99     write(msgbuf,'(a)')
100 dimitri 1.3 &'// ALLOW_ATM_TEMP: defined'
101 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
102     & SQUEEZE_RIGHT , mythid)
103     #else
104     write(msgbuf,'(a)')
105 dimitri 1.3 &'// ALLOW_ATM_TEMP: NOT defined'
106 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
107     & SQUEEZE_RIGHT , mythid)
108     #endif
109 dimitri 1.4
110 heimbach 1.1 #ifdef ALLOW_ATM_WIND
111     write(msgbuf,'(a)')
112 dimitri 1.3 &'// ALLOW_ATM_WIND: defined'
113     call print_message( msgbuf, standardmessageunit,
114     & SQUEEZE_RIGHT , mythid)
115     #else
116     write(msgbuf,'(a)')
117     &'// ALLOW_ATM_WIND: NOT defined'
118     call print_message( msgbuf, standardmessageunit,
119     & SQUEEZE_RIGHT , mythid)
120     #endif
121 dimitri 1.4
122     #ifdef ALLOW_DOWNWARD_RADIATION
123 dimitri 1.3 write(msgbuf,'(a)')
124 dimitri 1.4 &'// ALLOW_DOWNWARD_RADIATION: defined'
125 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
126     & SQUEEZE_RIGHT , mythid)
127     #else
128     write(msgbuf,'(a)')
129 dimitri 1.5 &'// ALLOW_DOWNWARD_RADIATION: NOT defined'
130 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
131     & SQUEEZE_RIGHT , mythid)
132     #endif
133    
134 dimitri 1.4 #ifdef ALLOW_BULKFORMULAE
135 heimbach 1.1 write(msgbuf,'(a)')
136 dimitri 1.4 &'// ALLOW_BULKFORMULAE: defined'
137 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
138     & SQUEEZE_RIGHT , mythid)
139     #else
140     write(msgbuf,'(a)')
141 dimitri 1.4 &'// ALLOW_BULKFORMULAE: NOT defined'
142 heimbach 1.2 call print_message( msgbuf, standardmessageunit,
143     & SQUEEZE_RIGHT , mythid)
144     #endif
145    
146 dimitri 1.4 c-- For each data set used the summary prints the calendar data
147     c and the corresponding file from which the data will be read.
148 dimitri 1.3
149 dimitri 1.4 #ifndef ALLOW_ATM_WIND
150     c-- Zonal wind stress.
151     il = ilnblnk(ustressfile)
152     call cal_TimeInterval( ustressperiod, 'secs', timeint, mythid )
153 dimitri 1.3 write(msgbuf,'(a)')
154     &' '
155     call print_message( msgbuf, standardmessageunit,
156     & SQUEEZE_RIGHT , mythid)
157     write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
158 dimitri 1.4 &' Zonal wind stress forcing starts at ',
159     & (ustressstartdate(i), i=1,2), dayofweek(ustressstartdate(4)),
160 dimitri 1.3 & '.'
161     call print_message( msgbuf, standardmessageunit,
162     & SQUEEZE_RIGHT , mythid)
163     write(msgbuf,'(a,i9.8,i7.6)')
164 dimitri 1.4 &' Zonal wind stress forcing period is ',
165 dimitri 1.3 & (timeint(i), i=1,2)
166     call print_message( msgbuf, standardmessageunit,
167     & SQUEEZE_RIGHT , mythid)
168     write(msgbuf,'(a)')
169 dimitri 1.4 &' Zonal wind stress forcing is read from file:'
170 dimitri 1.3 call print_message( msgbuf, standardmessageunit,
171     & SQUEEZE_RIGHT , mythid)
172     write(msgbuf,'(a,a,a)')
173 dimitri 1.4 &' >> ',ustressfile(1:il),' <<'
174 dimitri 1.3 call print_message( msgbuf, standardmessageunit,
175     & SQUEEZE_RIGHT , mythid)
176 heimbach 1.2
177 dimitri 1.4 c-- Meridional wind stress.
178     il = ilnblnk(vstressfile)
179     call cal_TimeInterval( vstressperiod, 'secs', timeint, mythid )
180 heimbach 1.2 write(msgbuf,'(a)')
181     &' '
182     call print_message( msgbuf, standardmessageunit,
183     & SQUEEZE_RIGHT , mythid)
184     write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
185 dimitri 1.4 &' Meridional wind stress forcing starts at ',
186     & (vstressstartdate(i), i=1,2), dayofweek(vstressstartdate(4)),
187 heimbach 1.2 & '.'
188     call print_message( msgbuf, standardmessageunit,
189     & SQUEEZE_RIGHT , mythid)
190     write(msgbuf,'(a,i9.8,i7.6)')
191 dimitri 1.4 &' Meridional wind stress forcing period is ',
192 heimbach 1.2 & (timeint(i), i=1,2)
193     call print_message( msgbuf, standardmessageunit,
194     & SQUEEZE_RIGHT , mythid)
195     write(msgbuf,'(a)')
196 dimitri 1.4 &' Meridional wind stress forcing is read from file:'
197 heimbach 1.2 call print_message( msgbuf, standardmessageunit,
198     & SQUEEZE_RIGHT , mythid)
199     write(msgbuf,'(a,a,a)')
200 dimitri 1.4 &' >> ',vstressfile(1:il),' <<'
201 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
202     & SQUEEZE_RIGHT , mythid)
203 dimitri 1.4 #endif
204 heimbach 1.1
205 dimitri 1.4 #ifndef ALLOW_ATM_TEMP
206     c-- Heat flux.
207     il = ilnblnk(hfluxfile)
208     call cal_TimeInterval( hfluxperiod, 'secs', timeint, mythid )
209 heimbach 1.1 write(msgbuf,'(a)')
210     &' '
211     call print_message( msgbuf, standardmessageunit,
212     & SQUEEZE_RIGHT , mythid)
213     write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
214 dimitri 1.4 &' Heat flux forcing starts at ',
215     & (hfluxstartdate(i), i=1,2), dayofweek(hfluxstartdate(4)),'.'
216 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
217     & SQUEEZE_RIGHT , mythid)
218     write(msgbuf,'(a,i9.8,i7.6)')
219 dimitri 1.4 &' Heat flux forcing period is ',
220 heimbach 1.1 & (timeint(i), i=1,2)
221     call print_message( msgbuf, standardmessageunit,
222     & SQUEEZE_RIGHT , mythid)
223     write(msgbuf,'(a)')
224 dimitri 1.4 &' Heat flux forcing is read from file: '
225 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
226     & SQUEEZE_RIGHT , mythid)
227     write(msgbuf,'(a,a,a)')
228 dimitri 1.4 &' >> ',hfluxfile(1:il),' <<'
229 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
230     & SQUEEZE_RIGHT , mythid)
231    
232 dimitri 1.4 c-- Salt flux.
233     il = ilnblnk(sfluxfile)
234     call cal_TimeInterval( sfluxperiod, 'secs', timeint, mythid )
235 heimbach 1.1 write(msgbuf,'(a)')
236     &' '
237     call print_message( msgbuf, standardmessageunit,
238     & SQUEEZE_RIGHT , mythid)
239     write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
240 dimitri 1.4 &' Salt flux forcing starts at ',
241     & (sfluxstartdate(i), i=1,2), dayofweek(sfluxstartdate(4)),'.'
242 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
243     & SQUEEZE_RIGHT , mythid)
244     write(msgbuf,'(a,i9.8,i7.6)')
245 dimitri 1.4 &' Salt flux forcing period is ',
246 heimbach 1.1 & (timeint(i), i=1,2)
247     call print_message( msgbuf, standardmessageunit,
248     & SQUEEZE_RIGHT , mythid)
249     write(msgbuf,'(a)')
250 dimitri 1.4 &' Salt flux forcing is read from file: '
251 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
252     & SQUEEZE_RIGHT , mythid)
253     write(msgbuf,'(a,a,a)')
254 dimitri 1.4 &' >> ',sfluxfile(1:il),' <<'
255 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
256     & SQUEEZE_RIGHT , mythid)
257 dimitri 1.4 #endif
258 heimbach 1.1
259 dimitri 1.4 #if defined(ALLOW_ATM_TEMP) || defined(SHORTWAVE_HEATING)
260     c-- Net shortwave.
261 heimbach 1.1 il = ilnblnk(swfluxfile)
262     call cal_TimeInterval( swfluxperiod, 'secs', timeint, mythid )
263     write(msgbuf,'(a)')
264     &' '
265     call print_message( msgbuf, standardmessageunit,
266     & SQUEEZE_RIGHT , mythid)
267     write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
268 dimitri 1.4 &' Net shortwave flux forcing starts at ',
269 heimbach 1.1 & (swfluxstartdate(i), i=1,2), dayofweek(swfluxstartdate(4)),
270     & '.'
271     call print_message( msgbuf, standardmessageunit,
272     & SQUEEZE_RIGHT , mythid)
273     write(msgbuf,'(a,i9.8,i7.6)')
274 dimitri 1.4 &' Net shortwave flux forcing period is ',
275 heimbach 1.1 & (timeint(i), i=1,2)
276     call print_message( msgbuf, standardmessageunit,
277     & SQUEEZE_RIGHT , mythid)
278     write(msgbuf,'(a)')
279 dimitri 1.4 &' Net shortwave flux forcing is read from file:'
280 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
281     & SQUEEZE_RIGHT , mythid)
282     write(msgbuf,'(a,a,a)')
283     &' >> ',swfluxfile(1:il),' <<'
284     call print_message( msgbuf, standardmessageunit,
285     & SQUEEZE_RIGHT , mythid)
286 dimitri 1.4 #endif
287 heimbach 1.1
288 dimitri 1.4 #ifdef ALLOW_ATM_WIND
289     c-- Zonal wind.
290     il = ilnblnk(uwindfile)
291     call cal_TimeInterval( uwindperiod, 'secs', timeint, mythid )
292 heimbach 1.1 write(msgbuf,'(a)')
293     &' '
294     call print_message( msgbuf, standardmessageunit,
295     & SQUEEZE_RIGHT , mythid)
296     write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
297 dimitri 1.4 &' Zonal wind forcing starts at ',
298     & (uwindstartdate(i), i=1,2), dayofweek(uwindstartdate(4)),
299 heimbach 1.1 & '.'
300     call print_message( msgbuf, standardmessageunit,
301     & SQUEEZE_RIGHT , mythid)
302     write(msgbuf,'(a,i9.8,i7.6)')
303 dimitri 1.4 &' Zonal wind forcing period is ',
304 heimbach 1.1 & (timeint(i), i=1,2)
305     call print_message( msgbuf, standardmessageunit,
306     & SQUEEZE_RIGHT , mythid)
307     write(msgbuf,'(a)')
308 dimitri 1.4 &' Zonal wind forcing is read from file:'
309 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
310     & SQUEEZE_RIGHT , mythid)
311     write(msgbuf,'(a,a,a)')
312 dimitri 1.4 &' >> ',uwindfile(1:il),' <<'
313 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
314     & SQUEEZE_RIGHT , mythid)
315    
316 dimitri 1.4 c-- Meridional wind.
317     il = ilnblnk(vwindfile)
318     call cal_TimeInterval( vwindperiod, 'secs', timeint, mythid )
319 heimbach 1.1 write(msgbuf,'(a)')
320     &' '
321     call print_message( msgbuf, standardmessageunit,
322     & SQUEEZE_RIGHT , mythid)
323     write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
324 dimitri 1.4 &' Meridional wind forcing starts at ',
325     & (vwindstartdate(i), i=1,2), dayofweek(vwindstartdate(4)),
326     & '.'
327 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
328     & SQUEEZE_RIGHT , mythid)
329     write(msgbuf,'(a,i9.8,i7.6)')
330 dimitri 1.4 &' Meridional wind forcing period is ',
331 heimbach 1.1 & (timeint(i), i=1,2)
332     call print_message( msgbuf, standardmessageunit,
333     & SQUEEZE_RIGHT , mythid)
334     write(msgbuf,'(a)')
335 dimitri 1.4 &' Meridional wind forcing is read from file:'
336 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
337     & SQUEEZE_RIGHT , mythid)
338     write(msgbuf,'(a,a,a)')
339 dimitri 1.4 &' >> ',vwindfile(1:il),' <<'
340 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
341     & SQUEEZE_RIGHT , mythid)
342 dimitri 1.4 #endif
343 heimbach 1.1
344 dimitri 1.4 #ifdef ALLOW_ATM_TEMP
345     c-- Atmospheric temperature.
346     il = ilnblnk(atempfile)
347     call cal_TimeInterval( atempperiod, 'secs', timeint, mythid )
348 heimbach 1.1 write(msgbuf,'(a)')
349     &' '
350     call print_message( msgbuf, standardmessageunit,
351     & SQUEEZE_RIGHT , mythid)
352     write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
353 dimitri 1.4 &' Atmospheric temperature starts at ',
354     & (atempstartdate(i), i=1,2), dayofweek(atempstartdate(4)),
355     & '.'
356 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
357     & SQUEEZE_RIGHT , mythid)
358     write(msgbuf,'(a,i9.8,i7.6)')
359 dimitri 1.4 &' Atmospheric temperature period is ',
360 heimbach 1.1 & (timeint(i), i=1,2)
361     call print_message( msgbuf, standardmessageunit,
362     & SQUEEZE_RIGHT , mythid)
363     write(msgbuf,'(a)')
364 dimitri 1.4 &' Atmospheric temperature is read from file:'
365 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
366     & SQUEEZE_RIGHT , mythid)
367     write(msgbuf,'(a,a,a)')
368 dimitri 1.4 &' >> ',atempfile(1:il),' <<'
369 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
370     & SQUEEZE_RIGHT , mythid)
371    
372 dimitri 1.4 c-- Atmospheric specific humidity.
373     il = ilnblnk(aqhfile)
374     call cal_TimeInterval( aqhperiod, 'secs', timeint, mythid )
375 heimbach 1.1 write(msgbuf,'(a)')
376     &' '
377     call print_message( msgbuf, standardmessageunit,
378     & SQUEEZE_RIGHT , mythid)
379     write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
380 dimitri 1.4 &' Atmospheric specific humidity starts at ',
381     & (aqhstartdate(i), i=1,2), dayofweek(aqhstartdate(4)),
382     & '.'
383 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
384     & SQUEEZE_RIGHT , mythid)
385     write(msgbuf,'(a,i9.8,i7.6)')
386 dimitri 1.4 &' Atmospheric specific humidity period is ',
387 heimbach 1.1 & (timeint(i), i=1,2)
388     call print_message( msgbuf, standardmessageunit,
389     & SQUEEZE_RIGHT , mythid)
390     write(msgbuf,'(a)')
391 dimitri 1.4 &' Atmospheric specific humidity is read from file:'
392 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
393     & SQUEEZE_RIGHT , mythid)
394     write(msgbuf,'(a,a,a)')
395 dimitri 1.4 &' >> ',aqhfile(1:il),' <<'
396 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
397     & SQUEEZE_RIGHT , mythid)
398    
399 dimitri 1.4 c-- Net longwave.
400     il = ilnblnk(lwfluxfile)
401     call cal_TimeInterval( lwfluxperiod, 'secs', timeint, mythid )
402 heimbach 1.1 write(msgbuf,'(a)')
403     &' '
404     call print_message( msgbuf, standardmessageunit,
405     & SQUEEZE_RIGHT , mythid)
406     write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
407 dimitri 1.4 &' Net longwave flux forcing starts at ',
408     & (lwfluxstartdate(i), i=1,2), dayofweek(lwfluxstartdate(4)),
409 heimbach 1.1 & '.'
410     call print_message( msgbuf, standardmessageunit,
411     & SQUEEZE_RIGHT , mythid)
412     write(msgbuf,'(a,i9.8,i7.6)')
413 dimitri 1.4 &' Net longwave flux forcing period is ',
414 heimbach 1.1 & (timeint(i), i=1,2)
415     call print_message( msgbuf, standardmessageunit,
416     & SQUEEZE_RIGHT , mythid)
417     write(msgbuf,'(a)')
418 dimitri 1.4 &' Net longwave flux forcing is read from file:'
419 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
420     & SQUEEZE_RIGHT , mythid)
421     write(msgbuf,'(a,a,a)')
422 dimitri 1.4 &' >> ',lwfluxfile(1:il),' <<'
423 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
424     & SQUEEZE_RIGHT , mythid)
425    
426 dimitri 1.4 c-- Precipitation.
427     il = ilnblnk(precipfile)
428     call cal_TimeInterval( precipperiod, 'secs', timeint, mythid )
429 heimbach 1.1 write(msgbuf,'(a)')
430     &' '
431     call print_message( msgbuf, standardmessageunit,
432     & SQUEEZE_RIGHT , mythid)
433     write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
434 dimitri 1.4 &' Precipitation data set starts at ',
435     & (precipstartdate(i), i=1,2), dayofweek(precipstartdate(4)),
436     &'.'
437 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
438     & SQUEEZE_RIGHT , mythid)
439     write(msgbuf,'(a,i9.8,i7.6)')
440 dimitri 1.4 &' Precipitation data period is ',
441 heimbach 1.1 & (timeint(i), i=1,2)
442     call print_message( msgbuf, standardmessageunit,
443     & SQUEEZE_RIGHT , mythid)
444     write(msgbuf,'(a)')
445 dimitri 1.4 &' Precipitation data is read from file: '
446 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
447     & SQUEEZE_RIGHT , mythid)
448     write(msgbuf,'(a,a,a)')
449 dimitri 1.4 &' >> ',precipfile(1:il),' <<'
450 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
451     & SQUEEZE_RIGHT , mythid)
452 dimitri 1.4 #endif
453 heimbach 1.1
454 dimitri 1.4 #ifdef EXF_READ_EVAP
455     c-- Evaporation.
456 heimbach 1.1 write(msgbuf,'(a)')
457 dimitri 1.4 &'// EXF_READ_EVAP: defined'
458 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
459     & SQUEEZE_RIGHT , mythid)
460 dimitri 1.4 il = ilnblnk(evapfile)
461     call cal_TimeInterval( evapperiod, 'secs', timeint, mythid )
462 heimbach 1.1
463     write(msgbuf,'(a)')
464     &' '
465     call print_message( msgbuf, standardmessageunit,
466     & SQUEEZE_RIGHT , mythid)
467     write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
468 dimitri 1.4 &' Evaporation starts at ',
469     & (evapstartdate(i), i=1,2), dayofweek(evapstartdate(4)),
470 heimbach 1.1 & '.'
471     call print_message( msgbuf, standardmessageunit,
472     & SQUEEZE_RIGHT , mythid)
473     write(msgbuf,'(a,i9.8,i7.6)')
474 dimitri 1.4 &' Evaporation period is ',
475 heimbach 1.1 & (timeint(i), i=1,2)
476     call print_message( msgbuf, standardmessageunit,
477     & SQUEEZE_RIGHT , mythid)
478     write(msgbuf,'(a)')
479 dimitri 1.4 &' Evaporation is read from file:'
480 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
481     & SQUEEZE_RIGHT , mythid)
482     write(msgbuf,'(a,a,a)')
483 dimitri 1.4 &' >> ',evapfile(1:il),' <<'
484 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
485     & SQUEEZE_RIGHT , mythid)
486 dimitri 1.4 #else
487 heimbach 1.1 write(msgbuf,'(a)')
488 dimitri 1.4 &'// EXF_READ_EVAP: NOT defined'
489 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
490     & SQUEEZE_RIGHT , mythid)
491 dimitri 1.4 #endif
492    
493     #ifdef ALLOW_RUNOFF
494     c-- Runoff.
495 heimbach 1.1 write(msgbuf,'(a)')
496 dimitri 1.4 &'// ALLOW_RUNOFF: defined'
497 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
498     & SQUEEZE_RIGHT , mythid)
499 dimitri 1.4 il = ilnblnk(runofffile)
500     call cal_TimeInterval( runoffperiod, 'secs', timeint, mythid )
501 heimbach 1.1
502     write(msgbuf,'(a)')
503     &' '
504     call print_message( msgbuf, standardmessageunit,
505     & SQUEEZE_RIGHT , mythid)
506     write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
507 dimitri 1.4 &' Runnoff starts at ',
508     & (runoffstartdate(i), i=1,2), dayofweek(runoffstartdate(4)),
509     & '.'
510 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
511     & SQUEEZE_RIGHT , mythid)
512     write(msgbuf,'(a,i9.8,i7.6)')
513 dimitri 1.4 &' Runoff period is ',
514 heimbach 1.1 & (timeint(i), i=1,2)
515     call print_message( msgbuf, standardmessageunit,
516     & SQUEEZE_RIGHT , mythid)
517     write(msgbuf,'(a)')
518 dimitri 1.4 &' Runoff is read from file:'
519 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
520     & SQUEEZE_RIGHT , mythid)
521     write(msgbuf,'(a,a,a)')
522 dimitri 1.4 &' >> ',runofffile(1:il),' <<'
523 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
524     & SQUEEZE_RIGHT , mythid)
525 dimitri 1.4 #else /* ALLOW_RUNOFF */
526 heimbach 1.1 write(msgbuf,'(a)')
527 dimitri 1.4 &'// ALLOW_RUNOFF: NOT defined'
528 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
529     & SQUEEZE_RIGHT , mythid)
530 dimitri 1.4 #endif /* ALLOW_RUNOFF */
531 heimbach 1.1
532 dimitri 1.4 #ifdef DOWNWARD_RADIATION
533     c-- Downward shortwave.
534     il = ilnblnk(swdownfile)
535     call cal_TimeInterval( swdownperiod, 'secs', timeint, mythid )
536 heimbach 1.1 write(msgbuf,'(a)')
537     &' '
538     call print_message( msgbuf, standardmessageunit,
539     & SQUEEZE_RIGHT , mythid)
540     write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
541 dimitri 1.4 &' Downward shortwave flux forcing starts at ',
542     & (swdownstartdate(i), i=1,2), dayofweek(swdownstartdate(4)),
543 heimbach 1.1 & '.'
544     call print_message( msgbuf, standardmessageunit,
545     & SQUEEZE_RIGHT , mythid)
546     write(msgbuf,'(a,i9.8,i7.6)')
547 dimitri 1.4 &' Downward shortwave flux forcing period is ',
548 heimbach 1.1 & (timeint(i), i=1,2)
549     call print_message( msgbuf, standardmessageunit,
550     & SQUEEZE_RIGHT , mythid)
551     write(msgbuf,'(a)')
552 dimitri 1.4 &' Downward shortwave flux forcing is read from file:'
553 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
554     & SQUEEZE_RIGHT , mythid)
555     write(msgbuf,'(a,a,a)')
556 dimitri 1.4 &' >> ',swdownfile(1:il),' <<'
557 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
558     & SQUEEZE_RIGHT , mythid)
559    
560 dimitri 1.4 c-- Downward longwave.
561     il = ilnblnk(lwdownfile)
562     call cal_TimeInterval( lwdownperiod, 'secs', timeint, mythid )
563 heimbach 1.1 write(msgbuf,'(a)')
564     &' '
565     call print_message( msgbuf, standardmessageunit,
566     & SQUEEZE_RIGHT , mythid)
567     write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
568 dimitri 1.4 &' Downward longwave flux forcing starts at ',
569     & (lwdownstartdate(i), i=1,2), dayofweek(lwdownstartdate(4)),
570 heimbach 1.1 & '.'
571     call print_message( msgbuf, standardmessageunit,
572     & SQUEEZE_RIGHT , mythid)
573     write(msgbuf,'(a,i9.8,i7.6)')
574 dimitri 1.4 &' Downward longwave flux forcing period is ',
575 heimbach 1.1 & (timeint(i), i=1,2)
576     call print_message( msgbuf, standardmessageunit,
577     & SQUEEZE_RIGHT , mythid)
578     write(msgbuf,'(a)')
579 dimitri 1.4 &' Downward longwave flux forcing is read from file:'
580 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
581     & SQUEEZE_RIGHT , mythid)
582     write(msgbuf,'(a,a,a)')
583 dimitri 1.4 &' >> ',lwdownfile(1:il),' <<'
584 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
585     & SQUEEZE_RIGHT , mythid)
586 dimitri 1.4 #endif
587 heimbach 1.1
588 dimitri 1.4 #ifdef ATMOSPHERIC_LOADING
589     c-- Atmospheric pressure.
590     il = ilnblnk(apressurefile)
591     call cal_TimeInterval( apressureperiod, 'secs', timeint, mythid )
592 heimbach 1.1 write(msgbuf,'(a)')
593     &' '
594     call print_message( msgbuf, standardmessageunit,
595     & SQUEEZE_RIGHT , mythid)
596     write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
597 dimitri 1.4 &' Atmospheric pressure forcing starts at ',
598     & (apressurestartdate(i), i=1,2),
599     & dayofweek(apressurestartdate(4)), '.'
600 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
601     & SQUEEZE_RIGHT , mythid)
602     write(msgbuf,'(a,i9.8,i7.6)')
603 dimitri 1.4 &' Atmospheric pressure forcing period is ',
604 heimbach 1.1 & (timeint(i), i=1,2)
605     call print_message( msgbuf, standardmessageunit,
606     & SQUEEZE_RIGHT , mythid)
607     write(msgbuf,'(a)')
608 dimitri 1.4 &' Atmospheric pressureforcing is read from file:'
609 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
610     & SQUEEZE_RIGHT , mythid)
611     write(msgbuf,'(a,a,a)')
612 dimitri 1.4 &' >> ',apressurefile(1:il),' <<'
613 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
614     & SQUEEZE_RIGHT , mythid)
615     #endif
616    
617     write(msgbuf,'(a)')
618     &' '
619     call print_message( msgbuf, standardmessageunit,
620     & SQUEEZE_RIGHT , mythid)
621     write(msgbuf,'(a)')
622     &'// ======================================================='
623     call print_message( msgbuf, standardmessageunit,
624     & SQUEEZE_RIGHT , mythid)
625     write(msgbuf,'(a)')
626     &'// External forcing configuration >>> END <<<'
627     call print_message( msgbuf, standardmessageunit,
628     & SQUEEZE_RIGHT , mythid)
629     write(msgbuf,'(a)')
630     &'// ======================================================='
631     call print_message( msgbuf, standardmessageunit,
632     & SQUEEZE_RIGHT , mythid)
633     write(msgbuf,'(a)')
634     &' '
635     call print_message( msgbuf, standardmessageunit,
636     & SQUEEZE_RIGHT , mythid)
637    
638     end
  ViewVC Help
Powered by ViewVC 1.1.22