/[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.12 - (show annotations) (download)
Tue May 8 00:21:31 2007 UTC (17 years, 1 month ago) by jmc
Branch: MAIN
Changes since 1.11: +173 -94 lines
print parameters from 1rst namelist.

1 C $Header: /u/gcmpack/MITgcm/pkg/exf/exf_summary.F,v 1.11 2007/04/18 19:55:34 heimbach Exp $
2 C $Name: $
3
4 #include "EXF_OPTIONS.h"
5
6
7 SUBROUTINE EXF_SUMMARY( myThid )
8
9 c ==================================================================
10 c SUBROUTINE exf_summary
11 c ==================================================================
12 c
13 c o List all the settings of the external forcing.
14 c
15 c started: Christian Eckert eckert@mit.edu 11-Jan-1999
16 c
17 c changed: Christian Eckert eckert@mit.edu 12-Feb-2000
18 c - changed routine names (package prefix: exf_)
19 c
20 c changed: Patrick Heimbach heimbach@mit.edu 04-May-2000
21 c - changed the handling of precip and sflux with respect
22 c to CPP options ALLOW_BULKFORMULAE and ALLOW_ATM_TEMP
23 c
24 c changed: Dimitris Menemenlis menemenlis@jpl.nasa.gov 20-Dec-2002
25 c - modifications for using pkg/exf with pkg/seaice
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_CONSTANTS.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 il
50
51 character*(max_len_mbuf) msgbuf
52
53 C == external ==
54
55 integer ilnblnk
56 external ilnblnk
57
58 C == end of interface ==
59
60 _BEGIN_MASTER( myThid )
61
62 write(msgbuf,'(a)') ' '
63 call print_message( msgbuf, standardmessageunit,
64 & SQUEEZE_RIGHT , mythid)
65 write(msgbuf,'(a)')
66 &'// ======================================================='
67 call print_message( msgbuf, standardmessageunit,
68 & SQUEEZE_RIGHT , mythid)
69 write(msgbuf,'(a)')
70 &'// External forcing configuration >>> START <<<'
71 call print_message( msgbuf, standardmessageunit,
72 & SQUEEZE_RIGHT , mythid)
73 write(msgbuf,'(a)')
74 &'// ======================================================='
75 call print_message( msgbuf, standardmessageunit,
76 & SQUEEZE_RIGHT , mythid)
77 write(msgbuf,'(a)') ' '
78 call print_message( msgbuf, standardmessageunit,
79 & SQUEEZE_RIGHT , mythid)
80
81 C-- Print general parameters:
82 WRITE(msgBuf,'(A)') ' EXF general parameters:'
83 CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
84 & SQUEEZE_RIGHT , myThid )
85 WRITE(msgBuf,'(A)') ' '
86 CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
87 & SQUEEZE_RIGHT , myThid )
88 CALL WRITE_0D_C( exf_yftype, 0, INDEX_NONE, 'exf_yftype = ',
89 & ' /* ? */')
90 CALL WRITE_0D_I( exf_iprec, INDEX_NONE, 'exf_iprec =',
91 & ' /* exf file precision */')
92 CALL WRITE_0D_L( useExfYearlyFields, INDEX_NONE,
93 & 'useExfYearlyFields =',
94 & ' /* add extension _YEAR to input file names */')
95 CALL WRITE_0D_L( twoDigitYear, INDEX_NONE, 'twoDigitYear =',
96 & ' /* use 2-digit year extension */')
97 CALL WRITE_0D_L( useExfCheckRange, INDEX_NONE,
98 & 'useExfCheckRange =',
99 & ' /* check for fields range */')
100 CALL WRITE_0D_R8( exf_monFreq, INDEX_NONE, 'exf_monFreq =',
101 & ' /* EXF monitor frequency [ s ] */')
102 CALL WRITE_0D_R8( repeatPeriod, INDEX_NONE, 'repeatPeriod =',
103 & ' /* period for cycling forcing dataset [ s ] */')
104 CALL WRITE_0D_R8( climtempfreeze, INDEX_NONE,'climTempFreeze=',
105 & ' /* Minimum climatological temperature [deg.C] */')
106 CALL WRITE_0D_R8( windStressMax, INDEX_NONE,'windStressMax =',
107 & ' /* Maximum absolute windstress [ Pa ] */')
108 CALL WRITE_0D_R8( exf_albedo, INDEX_NONE, 'exf_albedo =',
109 & ' /* Sea-water albedo [-] */')
110 CALL WRITE_0D_R8( ocean_emissivity, INDEX_NONE,
111 & 'ocean_emissivity =',
112 & ' /* longwave ocean-surface emissivity [-] */')
113 CALL WRITE_0D_R8( ice_emissivity, INDEX_NONE,'ice_emissivity =',
114 & ' /* longwave seaice emissivity [-] */')
115 CALL WRITE_0D_R8(snow_emissivity, INDEX_NONE,'snow_emissivity =',
116 & ' /* longwave snow emissivity [-] */')
117 CALL WRITE_0D_R8( atmrho, INDEX_NONE, 'atmrho =',
118 & ' /* mean atmospheric density [kg/m^3] */')
119 CALL WRITE_0D_R8( atmcp, INDEX_NONE, 'atmcp =',
120 & ' /* mean atmospheric specific heat [J/kg/K] */')
121 CALL WRITE_0D_R8( flamb, INDEX_NONE, 'flamb =',
122 & ' /* latent heat of evaporation [J/kg] */')
123 CALL WRITE_0D_R8( flami, INDEX_NONE, 'flami =',
124 & ' /* latent heat of pure-ice melting [J/kg] */')
125 CALL WRITE_0D_R8( cvapor_fac, INDEX_NONE, 'cvapor_fac =',
126 & ' /* const. for Saturation calculation [?] */')
127 CALL WRITE_0D_R8( cvapor_exp, INDEX_NONE, 'cvapor_exp =',
128 & ' /* const. for Saturation calculation [?] */')
129 CALL WRITE_0D_R8( cvapor_fac_ice, INDEX_NONE, 'cvapor_fac_ice=',
130 & ' /* const. for Saturation calculation [?] */')
131 CALL WRITE_0D_R8( cvapor_exp_ice, INDEX_NONE, 'cvapor_exp_ice=',
132 & ' /* const. for Saturation calculation [?] */')
133 CALL WRITE_0D_R8( humid_fac, INDEX_NONE, 'humid_fac =',
134 & ' /* humidity coef. in virtual temp. [(kg/kg)^-1] */')
135 CALL WRITE_0D_R8( gamma_blk, INDEX_NONE, 'gamma_blk =',
136 & ' /* adiabatic lapse rate [?] */')
137 CALL WRITE_0D_R8( saltsat, INDEX_NONE, 'saltsat =',
138 & ' /* reduction of Qsat over salty water [-] */')
139 CALL WRITE_0D_R8( cen2kel, INDEX_NONE, 'cen2kel =',
140 & ' /* conversion of deg. Centigrade to Kelvin [K] */')
141 CALL WRITE_0D_R8( exf_scal_BulkCdn, INDEX_NONE,
142 & 'exf_scal_BulkCdn=',
143 & ' /* Drag coefficient scaling factor [-] */')
144 CALL WRITE_0D_R8( psim_fac, INDEX_NONE, 'psim_fac =',
145 & ' /* coef used in turbulent fluxes calculation [-] */')
146 CALL WRITE_0D_R8( zolmin, INDEX_NONE, 'zolmin =',
147 & ' /* minimum stability parameter [?] */')
148 CALL WRITE_0D_R8( hu, INDEX_NONE, 'zref =',
149 & ' /* reference height [ m ] */')
150 CALL WRITE_0D_R8( hu, INDEX_NONE, 'hu =',
151 & ' /* height of mean wind [ m ] */')
152 CALL WRITE_0D_R8( hu, INDEX_NONE, 'ht =',
153 & ' /* height of mean temperature [ m ] */')
154 CALL WRITE_0D_R8( hu, INDEX_NONE, 'hq =',
155 & ' /* height of mean spec.humidity [ m ] */')
156 CALL WRITE_0D_R8( umin, INDEX_NONE, 'uMin =',
157 & ' /* minimum wind speed [m/s] */')
158 CALL WRITE_0D_R8( cdrag_1, INDEX_NONE, 'cDrag_1 =',
159 & ' /* coef used in drag calculation [?] */')
160 CALL WRITE_0D_R8( cdrag_2, INDEX_NONE, 'cDrag_2 =',
161 & ' /* coef used in drag calculation [?] */')
162 CALL WRITE_0D_R8( cdrag_3, INDEX_NONE, 'cDrag_3 =',
163 & ' /* coef used in drag calculation [?] */')
164 CALL WRITE_0D_R8( cstanton_1, INDEX_NONE, 'cStanton_1 =',
165 & ' /* coef used in Stanton number calculation [?] */')
166 CALL WRITE_0D_R8( cstanton_2, INDEX_NONE, 'cStanton_2 =',
167 & ' /* coef used in Stanton number calculation [?] */')
168 CALL WRITE_0D_R8( cdalton, INDEX_NONE, 'cDalton =',
169 & ' /* coef used in Dalton number calculation [?] */')
170 WRITE(msgBuf,'(A)') ' '
171 CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
172 & SQUEEZE_RIGHT , myThid )
173
174 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
175 C-- Print settings of some CPP flags.
176 WRITE(msgBuf,'(A)') ' EXF main CPP flags:'
177 CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
178 & SQUEEZE_RIGHT , myThid )
179 WRITE(msgBuf,'(A)') ' '
180 CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
181 & SQUEEZE_RIGHT , myThid )
182
183 #ifdef ALLOW_ATM_TEMP
184 write(msgbuf,'(a)')
185 &'// ALLOW_ATM_TEMP: defined'
186 call print_message( msgbuf, standardmessageunit,
187 & SQUEEZE_RIGHT , mythid)
188 #else
189 write(msgbuf,'(a)')
190 &'// ALLOW_ATM_TEMP: NOT defined'
191 call print_message( msgbuf, standardmessageunit,
192 & SQUEEZE_RIGHT , mythid)
193 #endif
194
195 #ifdef ALLOW_ATM_WIND
196 write(msgbuf,'(a)')
197 &'// ALLOW_ATM_WIND: defined'
198 call print_message( msgbuf, standardmessageunit,
199 & SQUEEZE_RIGHT , mythid)
200 #else
201 write(msgbuf,'(a)')
202 &'// ALLOW_ATM_WIND: NOT defined'
203 call print_message( msgbuf, standardmessageunit,
204 & SQUEEZE_RIGHT , mythid)
205 #endif
206
207 #ifdef ALLOW_DOWNWARD_RADIATION
208 write(msgbuf,'(a)')
209 &'// ALLOW_DOWNWARD_RADIATION: defined'
210 call print_message( msgbuf, standardmessageunit,
211 & SQUEEZE_RIGHT , mythid)
212 #else
213 write(msgbuf,'(a)')
214 &'// ALLOW_DOWNWARD_RADIATION: NOT defined'
215 call print_message( msgbuf, standardmessageunit,
216 & SQUEEZE_RIGHT , mythid)
217 #endif
218
219 #ifdef ALLOW_BULKFORMULAE
220 write(msgbuf,'(a)')
221 &'// ALLOW_BULKFORMULAE: defined'
222 call print_message( msgbuf, standardmessageunit,
223 & SQUEEZE_RIGHT , mythid)
224 #else
225 write(msgbuf,'(a)')
226 &'// ALLOW_BULKFORMULAE: NOT defined'
227 call print_message( msgbuf, standardmessageunit,
228 & SQUEEZE_RIGHT , mythid)
229 #endif
230
231 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
232
233 C-- For each data set used the summary prints the calendar data
234 C and the corresponding file from which the data will be read.
235
236 #ifndef ALLOW_ATM_WIND
237 C-- Zonal wind stress.
238 il = ilnblnk(ustressfile)
239 write(msgbuf,'(a)') ' '
240 call print_message( msgbuf, standardmessageunit,
241 & SQUEEZE_RIGHT , mythid)
242 write(msgbuf,'(a,f12.0)')
243 &' Zonal wind stress forcing starts at ',
244 & ustressstartdate
245 call print_message( msgbuf, standardmessageunit,
246 & SQUEEZE_RIGHT , mythid)
247 write(msgbuf,'(a,f12.0)')
248 &' Zonal wind stress forcing period is ',
249 & ustressperiod
250 call print_message( msgbuf, standardmessageunit,
251 & SQUEEZE_RIGHT , mythid)
252 write(msgbuf,'(a)')
253 &' Zonal wind stress forcing is read from file:'
254 call print_message( msgbuf, standardmessageunit,
255 & SQUEEZE_RIGHT , mythid)
256 write(msgbuf,'(a,a,a)')
257 &' >> ',ustressfile(1:il),' <<'
258 call print_message( msgbuf, standardmessageunit,
259 & SQUEEZE_RIGHT , mythid)
260
261 C-- Meridional wind stress.
262 il = ilnblnk(vstressfile)
263 write(msgbuf,'(a)') ' '
264 call print_message( msgbuf, standardmessageunit,
265 & SQUEEZE_RIGHT , mythid)
266 write(msgbuf,'(a,f12.0)')
267 &' Meridional wind stress forcing starts at ',
268 & vstressstartdate
269 call print_message( msgbuf, standardmessageunit,
270 & SQUEEZE_RIGHT , mythid)
271 write(msgbuf,'(a,f12.0)')
272 &' Meridional wind stress forcing period is ',
273 & vstressperiod
274 call print_message( msgbuf, standardmessageunit,
275 & SQUEEZE_RIGHT , mythid)
276 write(msgbuf,'(a)')
277 &' Meridional wind stress forcing is read from file:'
278 call print_message( msgbuf, standardmessageunit,
279 & SQUEEZE_RIGHT , mythid)
280 write(msgbuf,'(a,a,a)')
281 &' >> ',vstressfile(1:il),' <<'
282 call print_message( msgbuf, standardmessageunit,
283 & SQUEEZE_RIGHT , mythid)
284 #endif
285
286 #ifndef ALLOW_ATM_TEMP
287 C-- Heat flux.
288 il = ilnblnk(hfluxfile)
289 write(msgbuf,'(a)') ' '
290 call print_message( msgbuf, standardmessageunit,
291 & SQUEEZE_RIGHT , mythid)
292 write(msgbuf,'(a,f12.0)')
293 &' Heat flux forcing starts at ',
294 & hfluxstartdate
295 call print_message( msgbuf, standardmessageunit,
296 & SQUEEZE_RIGHT , mythid)
297 write(msgbuf,'(a,f12.0)')
298 &' Heat flux forcing period is ',
299 & hfluxperiod
300 call print_message( msgbuf, standardmessageunit,
301 & SQUEEZE_RIGHT , mythid)
302 write(msgbuf,'(a)')
303 &' Heat flux forcing is read from file: '
304 call print_message( msgbuf, standardmessageunit,
305 & SQUEEZE_RIGHT , mythid)
306 write(msgbuf,'(a,a,a)')
307 &' >> ',hfluxfile(1:il),' <<'
308 call print_message( msgbuf, standardmessageunit,
309 & SQUEEZE_RIGHT , mythid)
310
311 C-- Salt flux.
312 il = ilnblnk(sfluxfile)
313 write(msgbuf,'(a)') ' '
314 call print_message( msgbuf, standardmessageunit,
315 & SQUEEZE_RIGHT , mythid)
316 write(msgbuf,'(a,f12.0)')
317 &' Salt flux forcing starts at ',
318 & sfluxstartdate
319 call print_message( msgbuf, standardmessageunit,
320 & SQUEEZE_RIGHT , mythid)
321 write(msgbuf,'(a,f12.0)')
322 &' Salt flux forcing period is ',
323 & sfluxperiod
324 call print_message( msgbuf, standardmessageunit,
325 & SQUEEZE_RIGHT , mythid)
326 write(msgbuf,'(a)')
327 &' Salt flux forcing is read from file: '
328 call print_message( msgbuf, standardmessageunit,
329 & SQUEEZE_RIGHT , mythid)
330 write(msgbuf,'(a,a,a)')
331 &' >> ',sfluxfile(1:il),' <<'
332 call print_message( msgbuf, standardmessageunit,
333 & SQUEEZE_RIGHT , mythid)
334 #endif
335
336 #if defined(ALLOW_ATM_TEMP) || defined(SHORTWAVE_HEATING)
337 C-- Net shortwave.
338 il = ilnblnk(swfluxfile)
339 write(msgbuf,'(a)') ' '
340 call print_message( msgbuf, standardmessageunit,
341 & SQUEEZE_RIGHT , mythid)
342 write(msgbuf,'(a,f12.0)')
343 &' Net shortwave flux forcing starts at ',
344 & swfluxstartdate
345 call print_message( msgbuf, standardmessageunit,
346 & SQUEEZE_RIGHT , mythid)
347 write(msgbuf,'(a,f12.0)')
348 &' Net shortwave flux forcing period is ',
349 & swfluxperiod
350 call print_message( msgbuf, standardmessageunit,
351 & SQUEEZE_RIGHT , mythid)
352 write(msgbuf,'(a)')
353 &' Net shortwave flux forcing is read from file:'
354 call print_message( msgbuf, standardmessageunit,
355 & SQUEEZE_RIGHT , mythid)
356 write(msgbuf,'(a,a,a)')
357 &' >> ',swfluxfile(1:il),' <<'
358 call print_message( msgbuf, standardmessageunit,
359 & SQUEEZE_RIGHT , mythid)
360 #endif
361
362 #ifdef ALLOW_ATM_WIND
363 C-- Zonal wind.
364 il = ilnblnk(uwindfile)
365 write(msgbuf,'(a)') ' '
366 call print_message( msgbuf, standardmessageunit,
367 & SQUEEZE_RIGHT , mythid)
368 write(msgbuf,'(a,f12.0)')
369 &' Zonal wind forcing starts at ',
370 & uwindstartdate
371 call print_message( msgbuf, standardmessageunit,
372 & SQUEEZE_RIGHT , mythid)
373 write(msgbuf,'(a,f12.0)')
374 &' Zonal wind forcing period is ',
375 & uwindperiod
376 call print_message( msgbuf, standardmessageunit,
377 & SQUEEZE_RIGHT , mythid)
378 write(msgbuf,'(a)')
379 &' Zonal wind forcing is read from file:'
380 call print_message( msgbuf, standardmessageunit,
381 & SQUEEZE_RIGHT , mythid)
382 write(msgbuf,'(a,a,a)')
383 &' >> ',uwindfile(1:il),' <<'
384 call print_message( msgbuf, standardmessageunit,
385 & SQUEEZE_RIGHT , mythid)
386
387 C-- Meridional wind.
388 il = ilnblnk(vwindfile)
389 write(msgbuf,'(a)') ' '
390 call print_message( msgbuf, standardmessageunit,
391 & SQUEEZE_RIGHT , mythid)
392 write(msgbuf,'(a,f12.0)')
393 &' Meridional wind forcing starts at ',
394 & vwindstartdate
395 call print_message( msgbuf, standardmessageunit,
396 & SQUEEZE_RIGHT , mythid)
397 write(msgbuf,'(a,f12.0)')
398 &' Meridional wind forcing period is ',
399 & vwindperiod
400 call print_message( msgbuf, standardmessageunit,
401 & SQUEEZE_RIGHT , mythid)
402 write(msgbuf,'(a)')
403 &' Meridional wind forcing is read from file:'
404 call print_message( msgbuf, standardmessageunit,
405 & SQUEEZE_RIGHT , mythid)
406 write(msgbuf,'(a,a,a)')
407 &' >> ',vwindfile(1:il),' <<'
408 call print_message( msgbuf, standardmessageunit,
409 & SQUEEZE_RIGHT , mythid)
410 #endif
411
412 #ifdef ALLOW_ATM_TEMP
413 C-- Atmospheric temperature.
414 il = ilnblnk(atempfile)
415 write(msgbuf,'(a)') ' '
416 call print_message( msgbuf, standardmessageunit,
417 & SQUEEZE_RIGHT , mythid)
418 write(msgbuf,'(a,f12.0)')
419 &' Atmospheric temperature starts at ',
420 & atempstartdate
421 call print_message( msgbuf, standardmessageunit,
422 & SQUEEZE_RIGHT , mythid)
423 write(msgbuf,'(a,f12.0)')
424 &' Atmospheric temperature period is ',
425 & atempperiod
426 call print_message( msgbuf, standardmessageunit,
427 & SQUEEZE_RIGHT , mythid)
428 write(msgbuf,'(a)')
429 &' Atmospheric temperature is read from file:'
430 call print_message( msgbuf, standardmessageunit,
431 & SQUEEZE_RIGHT , mythid)
432 write(msgbuf,'(a,a,a)')
433 &' >> ',atempfile(1:il),' <<'
434 call print_message( msgbuf, standardmessageunit,
435 & SQUEEZE_RIGHT , mythid)
436
437 C-- Atmospheric specific humidity.
438 il = ilnblnk(aqhfile)
439 write(msgbuf,'(a)') ' '
440 call print_message( msgbuf, standardmessageunit,
441 & SQUEEZE_RIGHT , mythid)
442 write(msgbuf,'(a,f12.0)')
443 &' Atmospheric specific humidity starts at ',
444 & aqhstartdate
445 call print_message( msgbuf, standardmessageunit,
446 & SQUEEZE_RIGHT , mythid)
447 write(msgbuf,'(a,f12.0)')
448 &' Atmospheric specific humidity period is ',
449 & aqhperiod
450 call print_message( msgbuf, standardmessageunit,
451 & SQUEEZE_RIGHT , mythid)
452 write(msgbuf,'(a)')
453 &' Atmospheric specific humidity is read from file:'
454 call print_message( msgbuf, standardmessageunit,
455 & SQUEEZE_RIGHT , mythid)
456 write(msgbuf,'(a,a,a)')
457 &' >> ',aqhfile(1:il),' <<'
458 call print_message( msgbuf, standardmessageunit,
459 & SQUEEZE_RIGHT , mythid)
460
461 C-- Net longwave.
462 il = ilnblnk(lwfluxfile)
463 write(msgbuf,'(a)') ' '
464 call print_message( msgbuf, standardmessageunit,
465 & SQUEEZE_RIGHT , mythid)
466 write(msgbuf,'(a,f12.0)')
467 &' Net longwave flux forcing starts at ',
468 & lwfluxstartdate
469 call print_message( msgbuf, standardmessageunit,
470 & SQUEEZE_RIGHT , mythid)
471 write(msgbuf,'(a,f12.0)')
472 &' Net longwave flux forcing period is ',
473 & lwfluxperiod
474 call print_message( msgbuf, standardmessageunit,
475 & SQUEEZE_RIGHT , mythid)
476 write(msgbuf,'(a)')
477 &' Net longwave flux forcing is read from file:'
478 call print_message( msgbuf, standardmessageunit,
479 & SQUEEZE_RIGHT , mythid)
480 write(msgbuf,'(a,a,a)')
481 &' >> ',lwfluxfile(1:il),' <<'
482 call print_message( msgbuf, standardmessageunit,
483 & SQUEEZE_RIGHT , mythid)
484
485 C-- Precipitation.
486 il = ilnblnk(precipfile)
487 write(msgbuf,'(a)') ' '
488 call print_message( msgbuf, standardmessageunit,
489 & SQUEEZE_RIGHT , mythid)
490 write(msgbuf,'(a,f12.0)')
491 &' Precipitation data set starts at ',
492 & precipstartdate
493 call print_message( msgbuf, standardmessageunit,
494 & SQUEEZE_RIGHT , mythid)
495 write(msgbuf,'(a,f12.0)')
496 &' Precipitation data period is ',
497 & precipperiod
498 call print_message( msgbuf, standardmessageunit,
499 & SQUEEZE_RIGHT , mythid)
500 write(msgbuf,'(a)')
501 &' Precipitation data is read from file: '
502 call print_message( msgbuf, standardmessageunit,
503 & SQUEEZE_RIGHT , mythid)
504 write(msgbuf,'(a,a,a)')
505 &' >> ',precipfile(1:il),' <<'
506 call print_message( msgbuf, standardmessageunit,
507 & SQUEEZE_RIGHT , mythid)
508 #endif
509
510 C-- Evaporation.
511 write(msgbuf,'(a)') ' '
512 call print_message( msgbuf, standardmessageunit,
513 & SQUEEZE_RIGHT , mythid)
514 #ifdef EXF_READ_EVAP
515 write(msgbuf,'(a)')
516 &'// EXF_READ_EVAP: defined'
517 call print_message( msgbuf, standardmessageunit,
518 & SQUEEZE_RIGHT , mythid)
519 il = ilnblnk(evapfile)
520 write(msgbuf,'(a,f12.0)')
521 &' Evaporation starts at ',
522 & evapstartdate
523 call print_message( msgbuf, standardmessageunit,
524 & SQUEEZE_RIGHT , mythid)
525 write(msgbuf,'(a,f12.0)')
526 &' Evaporation period is ',
527 & evapperiod
528 call print_message( msgbuf, standardmessageunit,
529 & SQUEEZE_RIGHT , mythid)
530 write(msgbuf,'(a)')
531 &' Evaporation is read from file:'
532 call print_message( msgbuf, standardmessageunit,
533 & SQUEEZE_RIGHT , mythid)
534 write(msgbuf,'(a,a,a)')
535 &' >> ',evapfile(1:il),' <<'
536 call print_message( msgbuf, standardmessageunit,
537 & SQUEEZE_RIGHT , mythid)
538 #else
539 write(msgbuf,'(a)')
540 &'// EXF_READ_EVAP: NOT defined'
541 call print_message( msgbuf, standardmessageunit,
542 & SQUEEZE_RIGHT , mythid)
543 #endif
544
545 C-- Runoff.
546 write(msgbuf,'(a)') ' '
547 call print_message( msgbuf, standardmessageunit,
548 & SQUEEZE_RIGHT , mythid)
549 #ifdef ALLOW_RUNOFF
550 write(msgbuf,'(a)')
551 &'// ALLOW_RUNOFF: defined'
552 call print_message( msgbuf, standardmessageunit,
553 & SQUEEZE_RIGHT , mythid)
554 il = ilnblnk(runofffile)
555 write(msgbuf,'(a,f12.0)')
556 &' Runnoff starts at ',
557 & runoffstartdate
558 call print_message( msgbuf, standardmessageunit,
559 & SQUEEZE_RIGHT , mythid)
560 write(msgbuf,'(a,f12.0)')
561 &' Runoff period is ',
562 & runoffperiod
563 call print_message( msgbuf, standardmessageunit,
564 & SQUEEZE_RIGHT , mythid)
565 write(msgbuf,'(a)')
566 &' Runoff is read from file:'
567 call print_message( msgbuf, standardmessageunit,
568 & SQUEEZE_RIGHT , mythid)
569 write(msgbuf,'(a,a,a)')
570 &' >> ',runofffile(1:il),' <<'
571 call print_message( msgbuf, standardmessageunit,
572 & SQUEEZE_RIGHT , mythid)
573 #else /* ALLOW_RUNOFF */
574 write(msgbuf,'(a)')
575 &'// ALLOW_RUNOFF: NOT defined'
576 call print_message( msgbuf, standardmessageunit,
577 & SQUEEZE_RIGHT , mythid)
578 #endif /* ALLOW_RUNOFF */
579
580 #ifdef DOWNWARD_RADIATION
581 C-- Downward shortwave.
582 il = ilnblnk(swdownfile)
583 write(msgbuf,'(a)') ' '
584 call print_message( msgbuf, standardmessageunit,
585 & SQUEEZE_RIGHT , mythid)
586 write(msgbuf,'(a,f12.0)')
587 &' Downward shortwave flux forcing starts at ',
588 & swdownstartdate
589 call print_message( msgbuf, standardmessageunit,
590 & SQUEEZE_RIGHT , mythid)
591 write(msgbuf,'(a,f12.0)')
592 &' Downward shortwave flux forcing period is ',
593 & swdownperiod
594 call print_message( msgbuf, standardmessageunit,
595 & SQUEEZE_RIGHT , mythid)
596 write(msgbuf,'(a)')
597 &' Downward shortwave flux forcing is read from file:'
598 call print_message( msgbuf, standardmessageunit,
599 & SQUEEZE_RIGHT , mythid)
600 write(msgbuf,'(a,a,a)')
601 &' >> ',swdownfile(1:il),' <<'
602 call print_message( msgbuf, standardmessageunit,
603 & SQUEEZE_RIGHT , mythid)
604
605 C-- Downward longwave.
606 il = ilnblnk(lwdownfile)
607 write(msgbuf,'(a)') ' '
608 call print_message( msgbuf, standardmessageunit,
609 & SQUEEZE_RIGHT , mythid)
610 write(msgbuf,'(a,f12.0)')
611 &' Downward longwave flux forcing starts at ',
612 & lwdownstartdate
613 call print_message( msgbuf, standardmessageunit,
614 & SQUEEZE_RIGHT , mythid)
615 write(msgbuf,'(a,f12.0)')
616 &' Downward longwave flux forcing period is ',
617 & lwdownperiod
618 call print_message( msgbuf, standardmessageunit,
619 & SQUEEZE_RIGHT , mythid)
620 write(msgbuf,'(a)')
621 &' Downward longwave flux forcing is read from file:'
622 call print_message( msgbuf, standardmessageunit,
623 & SQUEEZE_RIGHT , mythid)
624 write(msgbuf,'(a,a,a)')
625 &' >> ',lwdownfile(1:il),' <<'
626 call print_message( msgbuf, standardmessageunit,
627 & SQUEEZE_RIGHT , mythid)
628 #endif
629
630 #ifdef ATMOSPHERIC_LOADING
631 C-- Atmospheric pressure.
632 il = ilnblnk(apressurefile)
633 write(msgbuf,'(a)') ' '
634 call print_message( msgbuf, standardmessageunit,
635 & SQUEEZE_RIGHT , mythid)
636 write(msgbuf,'(a,f12.0)')
637 &' Atmospheric pressure forcing starts at ',
638 & apressurestartdate
639 call print_message( msgbuf, standardmessageunit,
640 & SQUEEZE_RIGHT , mythid)
641 write(msgbuf,'(a,f12.0)')
642 &' Atmospheric pressure forcing period is ',
643 & apressureperiod
644 call print_message( msgbuf, standardmessageunit,
645 & SQUEEZE_RIGHT , mythid)
646 write(msgbuf,'(a)')
647 &' Atmospheric pressureforcing is read from file:'
648 call print_message( msgbuf, standardmessageunit,
649 & SQUEEZE_RIGHT , mythid)
650 write(msgbuf,'(a,a,a)')
651 &' >> ',apressurefile(1:il),' <<'
652 call print_message( msgbuf, standardmessageunit,
653 & SQUEEZE_RIGHT , mythid)
654 #endif
655
656 write(msgbuf,'(a)') ' '
657 call print_message( msgbuf, standardmessageunit,
658 & SQUEEZE_RIGHT , mythid)
659 write(msgbuf,'(a)')
660 &'// ======================================================='
661 call print_message( msgbuf, standardmessageunit,
662 & SQUEEZE_RIGHT , mythid)
663 write(msgbuf,'(a)')
664 &'// External forcing configuration >>> END <<<'
665 call print_message( msgbuf, standardmessageunit,
666 & SQUEEZE_RIGHT , mythid)
667 write(msgbuf,'(a)')
668 &'// ======================================================='
669 call print_message( msgbuf, standardmessageunit,
670 & SQUEEZE_RIGHT , mythid)
671 write(msgbuf,'(a)') ' '
672 call print_message( msgbuf, standardmessageunit,
673 & SQUEEZE_RIGHT , mythid)
674
675 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
676
677 call print_message( msgbuf, standardmessageunit,
678 & SQUEEZE_RIGHT , mythid)
679 write(msgbuf,'(a)')
680 &'// ======================================================='
681 call print_message( msgbuf, standardmessageunit,
682 & SQUEEZE_RIGHT , mythid)
683 write(msgbuf,'(a)')
684 &'// External forcing climatology configuration >>> START <<<'
685 call print_message( msgbuf, standardmessageunit,
686 & SQUEEZE_RIGHT , mythid)
687 write(msgbuf,'(a)')
688 &'// ======================================================='
689 call print_message( msgbuf, standardmessageunit,
690 & SQUEEZE_RIGHT , mythid)
691 write(msgbuf,'(a)') ' '
692 call print_message( msgbuf, standardmessageunit,
693 & SQUEEZE_RIGHT , mythid)
694
695 C For each data set used the summary prints the calendar data
696 C and the corresponding file from which the data will be read.
697
698 #ifdef ALLOW_CLIMSST_RELAXATION
699 write(msgbuf,'(a)')
700 &'// ALLOW_CLIMSST_RELAXATION: defined'
701 call print_message( msgbuf, standardmessageunit,
702 & SQUEEZE_RIGHT , mythid)
703 #else
704 write(msgbuf,'(a)')
705 &'// ALLOW_CLIMSST_RELAXATION: NOT defined'
706 call print_message( msgbuf, standardmessageunit,
707 & SQUEEZE_RIGHT , mythid)
708 #endif
709
710 #ifdef ALLOW_CLIMSSS_RELAXATION
711 write(msgbuf,'(a)')
712 &'// ALLOW_CLIMSSS_RELAXATION: defined'
713 call print_message( msgbuf, standardmessageunit,
714 & SQUEEZE_RIGHT , mythid)
715 #else
716 write(msgbuf,'(a)')
717 &'// ALLOW_CLIMSSS_RELAXATION: NOT defined'
718 call print_message( msgbuf, standardmessageunit,
719 & SQUEEZE_RIGHT , mythid)
720 #endif
721
722 C The climatological data sets are assumed to contain monthly
723 C data. This can be changed in a later version to an arbitrary
724 C number of intervals during a given year.
725
726 #ifdef ALLOW_CLIMSST_RELAXATION
727 C Relaxation to SST climatology.
728 il = ilnblnk(climsstfile)
729 write(msgbuf,'(a)')
730 &' '
731 call print_message( msgbuf, standardmessageunit,
732 & SQUEEZE_RIGHT , mythid)
733 write(msgbuf,'(a,f12.0)')
734 &' Climatological SST starts at ',
735 & climsststartdate
736 call print_message( msgbuf, standardmessageunit,
737 & SQUEEZE_RIGHT , mythid)
738 write(msgbuf,'(a,f12.0)')
739 &' Climatological SST period is ',
740 & climsstperiod
741 call print_message( msgbuf, standardmessageunit,
742 & SQUEEZE_RIGHT , mythid)
743 write(msgbuf,'(a)')
744 &' Climatological SST is read from file:'
745 call print_message( msgbuf, standardmessageunit,
746 & SQUEEZE_RIGHT , mythid)
747 write(msgbuf,'(a,a,a)')
748 &' >> ',climsstfile(1:il),' <<'
749 call print_message( msgbuf, standardmessageunit,
750 & SQUEEZE_RIGHT , mythid)
751 #endif
752
753 #ifdef ALLOW_CLIMSSS_RELAXATION
754 C Relaxation to SSS climatology.
755 il = ilnblnk(climsssfile)
756 write(msgbuf,'(a)')
757 &' '
758 call print_message( msgbuf, standardmessageunit,
759 & SQUEEZE_RIGHT , mythid)
760 write(msgbuf,'(a,f12.0)')
761 &' Climatological SSS starts at ',
762 & climsssstartdate
763 call print_message( msgbuf, standardmessageunit,
764 & SQUEEZE_RIGHT , mythid)
765 write(msgbuf,'(a,f12.0)')
766 &' Climatological SSS period is ',
767 & climsssperiod
768 call print_message( msgbuf, standardmessageunit,
769 & SQUEEZE_RIGHT , mythid)
770 write(msgbuf,'(a)')
771 &' Climatological SSS is read from file:'
772 call print_message( msgbuf, standardmessageunit,
773 & SQUEEZE_RIGHT , mythid)
774 write(msgbuf,'(a,a,a)')
775 &' >> ',climsssfile(1:il),' <<'
776 call print_message( msgbuf, standardmessageunit,
777 & SQUEEZE_RIGHT , mythid)
778 #endif
779
780 write(msgbuf,'(a)') ' '
781 call print_message( msgbuf, standardmessageunit,
782 & SQUEEZE_RIGHT , mythid)
783 write(msgbuf,'(a)')
784 &'// ======================================================='
785 call print_message( msgbuf, standardmessageunit,
786 & SQUEEZE_RIGHT , mythid)
787 write(msgbuf,'(a)')
788 &'// External forcing climatology configuration >>> END <<<'
789 call print_message( msgbuf, standardmessageunit,
790 & SQUEEZE_RIGHT , mythid)
791 write(msgbuf,'(a)')
792 &'// ======================================================='
793 call print_message( msgbuf, standardmessageunit,
794 & SQUEEZE_RIGHT , mythid)
795 write(msgbuf,'(a)') ' '
796 call print_message( msgbuf, standardmessageunit,
797 & SQUEEZE_RIGHT , mythid)
798
799 _END_MASTER( myThid )
800
801 RETURN
802 END

  ViewVC Help
Powered by ViewVC 1.1.22