/[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.5 - (show 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 c $Header: /u/gcmpack/MITgcm/pkg/exf/exf_summary.F,v 1.1.6.4 2003/02/19 06:00:56 dimitri 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 - changed routine names (package prefix: exf_)
20 c
21 c changed: Patrick Heimbach heimbach@mit.edu 04-May-2000
22 c - changed the handling of precip and sflux with respect
23 c to CPP options ALLOW_BULKFORMULAE and ALLOW_ATM_TEMP
24 c
25 c changed: Dimitris Menemenlis menemenlis@jpl.nasa.gov 20-Dec-2002
26 c - modifications for using pkg/exf with pkg/seaice
27 c
28 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 c-- Print settings of some CPP flags.
97
98 #ifdef ALLOW_ATM_TEMP
99 write(msgbuf,'(a)')
100 &'// ALLOW_ATM_TEMP: defined'
101 call print_message( msgbuf, standardmessageunit,
102 & SQUEEZE_RIGHT , mythid)
103 #else
104 write(msgbuf,'(a)')
105 &'// ALLOW_ATM_TEMP: NOT defined'
106 call print_message( msgbuf, standardmessageunit,
107 & SQUEEZE_RIGHT , mythid)
108 #endif
109
110 #ifdef ALLOW_ATM_WIND
111 write(msgbuf,'(a)')
112 &'// 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
122 #ifdef ALLOW_DOWNWARD_RADIATION
123 write(msgbuf,'(a)')
124 &'// ALLOW_DOWNWARD_RADIATION: defined'
125 call print_message( msgbuf, standardmessageunit,
126 & SQUEEZE_RIGHT , mythid)
127 #else
128 write(msgbuf,'(a)')
129 &'// ALLOW_DOWNWARD_RADIATION: NOT defined'
130 call print_message( msgbuf, standardmessageunit,
131 & SQUEEZE_RIGHT , mythid)
132 #endif
133
134 #ifdef ALLOW_BULKFORMULAE
135 write(msgbuf,'(a)')
136 &'// ALLOW_BULKFORMULAE: defined'
137 call print_message( msgbuf, standardmessageunit,
138 & SQUEEZE_RIGHT , mythid)
139 #else
140 write(msgbuf,'(a)')
141 &'// ALLOW_BULKFORMULAE: NOT defined'
142 call print_message( msgbuf, standardmessageunit,
143 & SQUEEZE_RIGHT , mythid)
144 #endif
145
146 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
149 #ifndef ALLOW_ATM_WIND
150 c-- Zonal wind stress.
151 il = ilnblnk(ustressfile)
152 call cal_TimeInterval( ustressperiod, 'secs', timeint, mythid )
153 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 &' Zonal wind stress forcing starts at ',
159 & (ustressstartdate(i), i=1,2), dayofweek(ustressstartdate(4)),
160 & '.'
161 call print_message( msgbuf, standardmessageunit,
162 & SQUEEZE_RIGHT , mythid)
163 write(msgbuf,'(a,i9.8,i7.6)')
164 &' Zonal wind stress forcing period is ',
165 & (timeint(i), i=1,2)
166 call print_message( msgbuf, standardmessageunit,
167 & SQUEEZE_RIGHT , mythid)
168 write(msgbuf,'(a)')
169 &' Zonal wind stress forcing is read from file:'
170 call print_message( msgbuf, standardmessageunit,
171 & SQUEEZE_RIGHT , mythid)
172 write(msgbuf,'(a,a,a)')
173 &' >> ',ustressfile(1:il),' <<'
174 call print_message( msgbuf, standardmessageunit,
175 & SQUEEZE_RIGHT , mythid)
176
177 c-- Meridional wind stress.
178 il = ilnblnk(vstressfile)
179 call cal_TimeInterval( vstressperiod, 'secs', timeint, mythid )
180 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 &' Meridional wind stress forcing starts at ',
186 & (vstressstartdate(i), i=1,2), dayofweek(vstressstartdate(4)),
187 & '.'
188 call print_message( msgbuf, standardmessageunit,
189 & SQUEEZE_RIGHT , mythid)
190 write(msgbuf,'(a,i9.8,i7.6)')
191 &' Meridional wind stress forcing period is ',
192 & (timeint(i), i=1,2)
193 call print_message( msgbuf, standardmessageunit,
194 & SQUEEZE_RIGHT , mythid)
195 write(msgbuf,'(a)')
196 &' Meridional wind stress forcing is read from file:'
197 call print_message( msgbuf, standardmessageunit,
198 & SQUEEZE_RIGHT , mythid)
199 write(msgbuf,'(a,a,a)')
200 &' >> ',vstressfile(1:il),' <<'
201 call print_message( msgbuf, standardmessageunit,
202 & SQUEEZE_RIGHT , mythid)
203 #endif
204
205 #ifndef ALLOW_ATM_TEMP
206 c-- Heat flux.
207 il = ilnblnk(hfluxfile)
208 call cal_TimeInterval( hfluxperiod, 'secs', timeint, mythid )
209 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 &' Heat flux forcing starts at ',
215 & (hfluxstartdate(i), i=1,2), dayofweek(hfluxstartdate(4)),'.'
216 call print_message( msgbuf, standardmessageunit,
217 & SQUEEZE_RIGHT , mythid)
218 write(msgbuf,'(a,i9.8,i7.6)')
219 &' Heat flux forcing period is ',
220 & (timeint(i), i=1,2)
221 call print_message( msgbuf, standardmessageunit,
222 & SQUEEZE_RIGHT , mythid)
223 write(msgbuf,'(a)')
224 &' Heat flux forcing is read from file: '
225 call print_message( msgbuf, standardmessageunit,
226 & SQUEEZE_RIGHT , mythid)
227 write(msgbuf,'(a,a,a)')
228 &' >> ',hfluxfile(1:il),' <<'
229 call print_message( msgbuf, standardmessageunit,
230 & SQUEEZE_RIGHT , mythid)
231
232 c-- Salt flux.
233 il = ilnblnk(sfluxfile)
234 call cal_TimeInterval( sfluxperiod, 'secs', timeint, mythid )
235 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 &' Salt flux forcing starts at ',
241 & (sfluxstartdate(i), i=1,2), dayofweek(sfluxstartdate(4)),'.'
242 call print_message( msgbuf, standardmessageunit,
243 & SQUEEZE_RIGHT , mythid)
244 write(msgbuf,'(a,i9.8,i7.6)')
245 &' Salt flux forcing period is ',
246 & (timeint(i), i=1,2)
247 call print_message( msgbuf, standardmessageunit,
248 & SQUEEZE_RIGHT , mythid)
249 write(msgbuf,'(a)')
250 &' Salt flux forcing is read from file: '
251 call print_message( msgbuf, standardmessageunit,
252 & SQUEEZE_RIGHT , mythid)
253 write(msgbuf,'(a,a,a)')
254 &' >> ',sfluxfile(1:il),' <<'
255 call print_message( msgbuf, standardmessageunit,
256 & SQUEEZE_RIGHT , mythid)
257 #endif
258
259 #if defined(ALLOW_ATM_TEMP) || defined(SHORTWAVE_HEATING)
260 c-- Net shortwave.
261 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 &' Net shortwave flux forcing starts at ',
269 & (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 &' Net shortwave flux forcing period is ',
275 & (timeint(i), i=1,2)
276 call print_message( msgbuf, standardmessageunit,
277 & SQUEEZE_RIGHT , mythid)
278 write(msgbuf,'(a)')
279 &' Net shortwave flux forcing is read from file:'
280 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 #endif
287
288 #ifdef ALLOW_ATM_WIND
289 c-- Zonal wind.
290 il = ilnblnk(uwindfile)
291 call cal_TimeInterval( uwindperiod, 'secs', timeint, mythid )
292 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 &' Zonal wind forcing starts at ',
298 & (uwindstartdate(i), i=1,2), dayofweek(uwindstartdate(4)),
299 & '.'
300 call print_message( msgbuf, standardmessageunit,
301 & SQUEEZE_RIGHT , mythid)
302 write(msgbuf,'(a,i9.8,i7.6)')
303 &' Zonal wind forcing period is ',
304 & (timeint(i), i=1,2)
305 call print_message( msgbuf, standardmessageunit,
306 & SQUEEZE_RIGHT , mythid)
307 write(msgbuf,'(a)')
308 &' Zonal wind forcing is read from file:'
309 call print_message( msgbuf, standardmessageunit,
310 & SQUEEZE_RIGHT , mythid)
311 write(msgbuf,'(a,a,a)')
312 &' >> ',uwindfile(1:il),' <<'
313 call print_message( msgbuf, standardmessageunit,
314 & SQUEEZE_RIGHT , mythid)
315
316 c-- Meridional wind.
317 il = ilnblnk(vwindfile)
318 call cal_TimeInterval( vwindperiod, 'secs', timeint, mythid )
319 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 &' Meridional wind forcing starts at ',
325 & (vwindstartdate(i), i=1,2), dayofweek(vwindstartdate(4)),
326 & '.'
327 call print_message( msgbuf, standardmessageunit,
328 & SQUEEZE_RIGHT , mythid)
329 write(msgbuf,'(a,i9.8,i7.6)')
330 &' Meridional wind forcing period is ',
331 & (timeint(i), i=1,2)
332 call print_message( msgbuf, standardmessageunit,
333 & SQUEEZE_RIGHT , mythid)
334 write(msgbuf,'(a)')
335 &' Meridional wind forcing is read from file:'
336 call print_message( msgbuf, standardmessageunit,
337 & SQUEEZE_RIGHT , mythid)
338 write(msgbuf,'(a,a,a)')
339 &' >> ',vwindfile(1:il),' <<'
340 call print_message( msgbuf, standardmessageunit,
341 & SQUEEZE_RIGHT , mythid)
342 #endif
343
344 #ifdef ALLOW_ATM_TEMP
345 c-- Atmospheric temperature.
346 il = ilnblnk(atempfile)
347 call cal_TimeInterval( atempperiod, 'secs', timeint, mythid )
348 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 &' Atmospheric temperature starts at ',
354 & (atempstartdate(i), i=1,2), dayofweek(atempstartdate(4)),
355 & '.'
356 call print_message( msgbuf, standardmessageunit,
357 & SQUEEZE_RIGHT , mythid)
358 write(msgbuf,'(a,i9.8,i7.6)')
359 &' Atmospheric temperature period is ',
360 & (timeint(i), i=1,2)
361 call print_message( msgbuf, standardmessageunit,
362 & SQUEEZE_RIGHT , mythid)
363 write(msgbuf,'(a)')
364 &' Atmospheric temperature is read from file:'
365 call print_message( msgbuf, standardmessageunit,
366 & SQUEEZE_RIGHT , mythid)
367 write(msgbuf,'(a,a,a)')
368 &' >> ',atempfile(1:il),' <<'
369 call print_message( msgbuf, standardmessageunit,
370 & SQUEEZE_RIGHT , mythid)
371
372 c-- Atmospheric specific humidity.
373 il = ilnblnk(aqhfile)
374 call cal_TimeInterval( aqhperiod, 'secs', timeint, mythid )
375 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 &' Atmospheric specific humidity starts at ',
381 & (aqhstartdate(i), i=1,2), dayofweek(aqhstartdate(4)),
382 & '.'
383 call print_message( msgbuf, standardmessageunit,
384 & SQUEEZE_RIGHT , mythid)
385 write(msgbuf,'(a,i9.8,i7.6)')
386 &' Atmospheric specific humidity period is ',
387 & (timeint(i), i=1,2)
388 call print_message( msgbuf, standardmessageunit,
389 & SQUEEZE_RIGHT , mythid)
390 write(msgbuf,'(a)')
391 &' Atmospheric specific humidity is read from file:'
392 call print_message( msgbuf, standardmessageunit,
393 & SQUEEZE_RIGHT , mythid)
394 write(msgbuf,'(a,a,a)')
395 &' >> ',aqhfile(1:il),' <<'
396 call print_message( msgbuf, standardmessageunit,
397 & SQUEEZE_RIGHT , mythid)
398
399 c-- Net longwave.
400 il = ilnblnk(lwfluxfile)
401 call cal_TimeInterval( lwfluxperiod, 'secs', timeint, mythid )
402 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 &' Net longwave flux forcing starts at ',
408 & (lwfluxstartdate(i), i=1,2), dayofweek(lwfluxstartdate(4)),
409 & '.'
410 call print_message( msgbuf, standardmessageunit,
411 & SQUEEZE_RIGHT , mythid)
412 write(msgbuf,'(a,i9.8,i7.6)')
413 &' Net longwave flux forcing period is ',
414 & (timeint(i), i=1,2)
415 call print_message( msgbuf, standardmessageunit,
416 & SQUEEZE_RIGHT , mythid)
417 write(msgbuf,'(a)')
418 &' Net longwave flux forcing is read from file:'
419 call print_message( msgbuf, standardmessageunit,
420 & SQUEEZE_RIGHT , mythid)
421 write(msgbuf,'(a,a,a)')
422 &' >> ',lwfluxfile(1:il),' <<'
423 call print_message( msgbuf, standardmessageunit,
424 & SQUEEZE_RIGHT , mythid)
425
426 c-- Precipitation.
427 il = ilnblnk(precipfile)
428 call cal_TimeInterval( precipperiod, 'secs', timeint, mythid )
429 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 &' Precipitation data set starts at ',
435 & (precipstartdate(i), i=1,2), dayofweek(precipstartdate(4)),
436 &'.'
437 call print_message( msgbuf, standardmessageunit,
438 & SQUEEZE_RIGHT , mythid)
439 write(msgbuf,'(a,i9.8,i7.6)')
440 &' Precipitation data period is ',
441 & (timeint(i), i=1,2)
442 call print_message( msgbuf, standardmessageunit,
443 & SQUEEZE_RIGHT , mythid)
444 write(msgbuf,'(a)')
445 &' Precipitation data is read from file: '
446 call print_message( msgbuf, standardmessageunit,
447 & SQUEEZE_RIGHT , mythid)
448 write(msgbuf,'(a,a,a)')
449 &' >> ',precipfile(1:il),' <<'
450 call print_message( msgbuf, standardmessageunit,
451 & SQUEEZE_RIGHT , mythid)
452 #endif
453
454 #ifdef EXF_READ_EVAP
455 c-- Evaporation.
456 write(msgbuf,'(a)')
457 &'// EXF_READ_EVAP: defined'
458 call print_message( msgbuf, standardmessageunit,
459 & SQUEEZE_RIGHT , mythid)
460 il = ilnblnk(evapfile)
461 call cal_TimeInterval( evapperiod, 'secs', timeint, mythid )
462
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 &' Evaporation starts at ',
469 & (evapstartdate(i), i=1,2), dayofweek(evapstartdate(4)),
470 & '.'
471 call print_message( msgbuf, standardmessageunit,
472 & SQUEEZE_RIGHT , mythid)
473 write(msgbuf,'(a,i9.8,i7.6)')
474 &' Evaporation period is ',
475 & (timeint(i), i=1,2)
476 call print_message( msgbuf, standardmessageunit,
477 & SQUEEZE_RIGHT , mythid)
478 write(msgbuf,'(a)')
479 &' Evaporation is read from file:'
480 call print_message( msgbuf, standardmessageunit,
481 & SQUEEZE_RIGHT , mythid)
482 write(msgbuf,'(a,a,a)')
483 &' >> ',evapfile(1:il),' <<'
484 call print_message( msgbuf, standardmessageunit,
485 & SQUEEZE_RIGHT , mythid)
486 #else
487 write(msgbuf,'(a)')
488 &'// EXF_READ_EVAP: NOT defined'
489 call print_message( msgbuf, standardmessageunit,
490 & SQUEEZE_RIGHT , mythid)
491 #endif
492
493 #ifdef ALLOW_RUNOFF
494 c-- Runoff.
495 write(msgbuf,'(a)')
496 &'// ALLOW_RUNOFF: defined'
497 call print_message( msgbuf, standardmessageunit,
498 & SQUEEZE_RIGHT , mythid)
499 il = ilnblnk(runofffile)
500 call cal_TimeInterval( runoffperiod, 'secs', timeint, mythid )
501
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 &' Runnoff starts at ',
508 & (runoffstartdate(i), i=1,2), dayofweek(runoffstartdate(4)),
509 & '.'
510 call print_message( msgbuf, standardmessageunit,
511 & SQUEEZE_RIGHT , mythid)
512 write(msgbuf,'(a,i9.8,i7.6)')
513 &' Runoff period is ',
514 & (timeint(i), i=1,2)
515 call print_message( msgbuf, standardmessageunit,
516 & SQUEEZE_RIGHT , mythid)
517 write(msgbuf,'(a)')
518 &' Runoff is read from file:'
519 call print_message( msgbuf, standardmessageunit,
520 & SQUEEZE_RIGHT , mythid)
521 write(msgbuf,'(a,a,a)')
522 &' >> ',runofffile(1:il),' <<'
523 call print_message( msgbuf, standardmessageunit,
524 & SQUEEZE_RIGHT , mythid)
525 #else /* ALLOW_RUNOFF */
526 write(msgbuf,'(a)')
527 &'// ALLOW_RUNOFF: NOT defined'
528 call print_message( msgbuf, standardmessageunit,
529 & SQUEEZE_RIGHT , mythid)
530 #endif /* ALLOW_RUNOFF */
531
532 #ifdef DOWNWARD_RADIATION
533 c-- Downward shortwave.
534 il = ilnblnk(swdownfile)
535 call cal_TimeInterval( swdownperiod, 'secs', timeint, mythid )
536 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 &' Downward shortwave flux forcing starts at ',
542 & (swdownstartdate(i), i=1,2), dayofweek(swdownstartdate(4)),
543 & '.'
544 call print_message( msgbuf, standardmessageunit,
545 & SQUEEZE_RIGHT , mythid)
546 write(msgbuf,'(a,i9.8,i7.6)')
547 &' Downward shortwave flux forcing period is ',
548 & (timeint(i), i=1,2)
549 call print_message( msgbuf, standardmessageunit,
550 & SQUEEZE_RIGHT , mythid)
551 write(msgbuf,'(a)')
552 &' Downward shortwave flux forcing is read from file:'
553 call print_message( msgbuf, standardmessageunit,
554 & SQUEEZE_RIGHT , mythid)
555 write(msgbuf,'(a,a,a)')
556 &' >> ',swdownfile(1:il),' <<'
557 call print_message( msgbuf, standardmessageunit,
558 & SQUEEZE_RIGHT , mythid)
559
560 c-- Downward longwave.
561 il = ilnblnk(lwdownfile)
562 call cal_TimeInterval( lwdownperiod, 'secs', timeint, mythid )
563 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 &' Downward longwave flux forcing starts at ',
569 & (lwdownstartdate(i), i=1,2), dayofweek(lwdownstartdate(4)),
570 & '.'
571 call print_message( msgbuf, standardmessageunit,
572 & SQUEEZE_RIGHT , mythid)
573 write(msgbuf,'(a,i9.8,i7.6)')
574 &' Downward longwave flux forcing period is ',
575 & (timeint(i), i=1,2)
576 call print_message( msgbuf, standardmessageunit,
577 & SQUEEZE_RIGHT , mythid)
578 write(msgbuf,'(a)')
579 &' Downward longwave flux forcing is read from file:'
580 call print_message( msgbuf, standardmessageunit,
581 & SQUEEZE_RIGHT , mythid)
582 write(msgbuf,'(a,a,a)')
583 &' >> ',lwdownfile(1:il),' <<'
584 call print_message( msgbuf, standardmessageunit,
585 & SQUEEZE_RIGHT , mythid)
586 #endif
587
588 #ifdef ATMOSPHERIC_LOADING
589 c-- Atmospheric pressure.
590 il = ilnblnk(apressurefile)
591 call cal_TimeInterval( apressureperiod, 'secs', timeint, mythid )
592 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 &' Atmospheric pressure forcing starts at ',
598 & (apressurestartdate(i), i=1,2),
599 & dayofweek(apressurestartdate(4)), '.'
600 call print_message( msgbuf, standardmessageunit,
601 & SQUEEZE_RIGHT , mythid)
602 write(msgbuf,'(a,i9.8,i7.6)')
603 &' Atmospheric pressure forcing period is ',
604 & (timeint(i), i=1,2)
605 call print_message( msgbuf, standardmessageunit,
606 & SQUEEZE_RIGHT , mythid)
607 write(msgbuf,'(a)')
608 &' Atmospheric pressureforcing is read from file:'
609 call print_message( msgbuf, standardmessageunit,
610 & SQUEEZE_RIGHT , mythid)
611 write(msgbuf,'(a,a,a)')
612 &' >> ',apressurefile(1:il),' <<'
613 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