/[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.12 - (hide annotations) (download)
Tue May 8 00:21:31 2007 UTC (17 years ago) by jmc
Branch: MAIN
Changes since 1.11: +173 -94 lines 
print parameters from 1rst namelist.
1 jmc 1.12 C $Header: /u/gcmpack/MITgcm/pkg/exf/exf_summary.F,v 1.11 2007/04/18 19:55:34 heimbach Exp $
2 jmc 1.8 C $Name: $
3 heimbach 1.1
4 edhill 1.6 #include "EXF_OPTIONS.h"
5 heimbach 1.1
6    
7 jmc 1.12 SUBROUTINE EXF_SUMMARY( myThid )
8 heimbach 1.1
9     c ==================================================================
10 heimbach 1.11 c SUBROUTINE exf_summary
11 heimbach 1.1 c ==================================================================
12     c
13     c o List all the settings of the external forcing.
14     c
15 dimitri 1.4 c started: Christian Eckert eckert@mit.edu 11-Jan-1999
16 heimbach 1.1 c
17 dimitri 1.4 c changed: Christian Eckert eckert@mit.edu 12-Feb-2000
18     c - changed routine names (package prefix: exf_)
19 heimbach 1.1 c
20 dimitri 1.4 c changed: Patrick Heimbach heimbach@mit.edu 04-May-2000
21 heimbach 1.1 c - changed the handling of precip and sflux with respect
22     c to CPP options ALLOW_BULKFORMULAE and ALLOW_ATM_TEMP
23     c
24 dimitri 1.4 c changed: Dimitris Menemenlis menemenlis@jpl.nasa.gov 20-Dec-2002
25     c - modifications for using pkg/exf with pkg/seaice
26 dimitri 1.3 c
27 heimbach 1.1 c ==================================================================
28 heimbach 1.11 c SUBROUTINE exf_summary
29 heimbach 1.1 c ==================================================================
30    
31     implicit none
32    
33 jmc 1.12 C == global variables ==
34 heimbach 1.1
35     #include "EEPARAMS.h"
36     #include "SIZE.h"
37     #include "cal.h"
38 jmc 1.12 #include "EXF_CONSTANTS.h"
39 jmc 1.8 #include "EXF_PARAM.h"
40 heimbach 1.1
41 jmc 1.12 C == routine arguments ==
42 heimbach 1.1
43 jmc 1.12 C myThid - thread number for this instance of the routine.
44 heimbach 1.1
45 jmc 1.12 integer myThid
46 heimbach 1.1
47 jmc 1.12 C == local variables ==
48 heimbach 1.1
49     integer il
50    
51     character*(max_len_mbuf) msgbuf
52    
53 jmc 1.12 C == external ==
54 heimbach 1.1
55     integer ilnblnk
56     external ilnblnk
57    
58 jmc 1.12 C == end of interface ==
59 heimbach 1.1
60 jmc 1.12 _BEGIN_MASTER( myThid )
61    
62     write(msgbuf,'(a)') ' '
63 heimbach 1.1 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 jmc 1.12 write(msgbuf,'(a)') ' '
78 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
79     & SQUEEZE_RIGHT , mythid)
80    
81 jmc 1.12 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 heimbach 1.1
183     #ifdef ALLOW_ATM_TEMP
184     write(msgbuf,'(a)')
185 dimitri 1.3 &'// ALLOW_ATM_TEMP: defined'
186 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
187     & SQUEEZE_RIGHT , mythid)
188     #else
189     write(msgbuf,'(a)')
190 dimitri 1.3 &'// ALLOW_ATM_TEMP: NOT defined'
191 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
192     & SQUEEZE_RIGHT , mythid)
193     #endif
194 dimitri 1.4
195 heimbach 1.1 #ifdef ALLOW_ATM_WIND
196     write(msgbuf,'(a)')
197 dimitri 1.3 &'// 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 dimitri 1.4
207     #ifdef ALLOW_DOWNWARD_RADIATION
208 dimitri 1.3 write(msgbuf,'(a)')
209 dimitri 1.4 &'// ALLOW_DOWNWARD_RADIATION: defined'
210 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
211     & SQUEEZE_RIGHT , mythid)
212     #else
213     write(msgbuf,'(a)')
214 dimitri 1.5 &'// ALLOW_DOWNWARD_RADIATION: NOT defined'
215 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
216     & SQUEEZE_RIGHT , mythid)
217     #endif
218    
219 dimitri 1.4 #ifdef ALLOW_BULKFORMULAE
220 heimbach 1.1 write(msgbuf,'(a)')
221 dimitri 1.4 &'// ALLOW_BULKFORMULAE: defined'
222 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
223     & SQUEEZE_RIGHT , mythid)
224     #else
225     write(msgbuf,'(a)')
226 dimitri 1.4 &'// ALLOW_BULKFORMULAE: NOT defined'
227 heimbach 1.2 call print_message( msgbuf, standardmessageunit,
228     & SQUEEZE_RIGHT , mythid)
229     #endif
230    
231 jmc 1.12 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 dimitri 1.3
236 dimitri 1.4 #ifndef ALLOW_ATM_WIND
237 jmc 1.12 C-- Zonal wind stress.
238 dimitri 1.4 il = ilnblnk(ustressfile)
239 jmc 1.12 write(msgbuf,'(a)') ' '
240 dimitri 1.3 call print_message( msgbuf, standardmessageunit,
241     & SQUEEZE_RIGHT , mythid)
242 dimitri 1.7 write(msgbuf,'(a,f12.0)')
243 dimitri 1.4 &' Zonal wind stress forcing starts at ',
244 dimitri 1.7 & ustressstartdate
245 dimitri 1.3 call print_message( msgbuf, standardmessageunit,
246     & SQUEEZE_RIGHT , mythid)
247 dimitri 1.7 write(msgbuf,'(a,f12.0)')
248 dimitri 1.4 &' Zonal wind stress forcing period is ',
249 dimitri 1.7 & ustressperiod
250 dimitri 1.3 call print_message( msgbuf, standardmessageunit,
251     & SQUEEZE_RIGHT , mythid)
252     write(msgbuf,'(a)')
253 dimitri 1.4 &' Zonal wind stress forcing is read from file:'
254 dimitri 1.3 call print_message( msgbuf, standardmessageunit,
255     & SQUEEZE_RIGHT , mythid)
256     write(msgbuf,'(a,a,a)')
257 dimitri 1.4 &' >> ',ustressfile(1:il),' <<'
258 dimitri 1.3 call print_message( msgbuf, standardmessageunit,
259     & SQUEEZE_RIGHT , mythid)
260 heimbach 1.2
261 jmc 1.12 C-- Meridional wind stress.
262 dimitri 1.4 il = ilnblnk(vstressfile)
263 jmc 1.12 write(msgbuf,'(a)') ' '
264 heimbach 1.2 call print_message( msgbuf, standardmessageunit,
265     & SQUEEZE_RIGHT , mythid)
266 dimitri 1.7 write(msgbuf,'(a,f12.0)')
267 dimitri 1.4 &' Meridional wind stress forcing starts at ',
268 dimitri 1.7 & vstressstartdate
269 heimbach 1.2 call print_message( msgbuf, standardmessageunit,
270     & SQUEEZE_RIGHT , mythid)
271 dimitri 1.7 write(msgbuf,'(a,f12.0)')
272 dimitri 1.4 &' Meridional wind stress forcing period is ',
273 dimitri 1.7 & vstressperiod
274 heimbach 1.2 call print_message( msgbuf, standardmessageunit,
275     & SQUEEZE_RIGHT , mythid)
276     write(msgbuf,'(a)')
277 dimitri 1.4 &' Meridional wind stress forcing is read from file:'
278 heimbach 1.2 call print_message( msgbuf, standardmessageunit,
279     & SQUEEZE_RIGHT , mythid)
280     write(msgbuf,'(a,a,a)')
281 dimitri 1.4 &' >> ',vstressfile(1:il),' <<'
282 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
283     & SQUEEZE_RIGHT , mythid)
284 dimitri 1.4 #endif
285 heimbach 1.1
286 dimitri 1.4 #ifndef ALLOW_ATM_TEMP
287 jmc 1.12 C-- Heat flux.
288 dimitri 1.4 il = ilnblnk(hfluxfile)
289 jmc 1.12 write(msgbuf,'(a)') ' '
290 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
291     & SQUEEZE_RIGHT , mythid)
292 dimitri 1.7 write(msgbuf,'(a,f12.0)')
293 dimitri 1.4 &' Heat flux forcing starts at ',
294 dimitri 1.7 & hfluxstartdate
295 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
296     & SQUEEZE_RIGHT , mythid)
297 dimitri 1.7 write(msgbuf,'(a,f12.0)')
298 dimitri 1.4 &' Heat flux forcing period is ',
299 dimitri 1.7 & hfluxperiod
300 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
301     & SQUEEZE_RIGHT , mythid)
302     write(msgbuf,'(a)')
303 dimitri 1.4 &' Heat flux forcing is read from file: '
304 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
305     & SQUEEZE_RIGHT , mythid)
306     write(msgbuf,'(a,a,a)')
307 dimitri 1.4 &' >> ',hfluxfile(1:il),' <<'
308 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
309     & SQUEEZE_RIGHT , mythid)
310    
311 jmc 1.12 C-- Salt flux.
312 dimitri 1.4 il = ilnblnk(sfluxfile)
313 jmc 1.12 write(msgbuf,'(a)') ' '
314 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
315     & SQUEEZE_RIGHT , mythid)
316 dimitri 1.7 write(msgbuf,'(a,f12.0)')
317 dimitri 1.4 &' Salt flux forcing starts at ',
318 dimitri 1.7 & sfluxstartdate
319 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
320     & SQUEEZE_RIGHT , mythid)
321 dimitri 1.7 write(msgbuf,'(a,f12.0)')
322 dimitri 1.4 &' Salt flux forcing period is ',
323 dimitri 1.7 & sfluxperiod
324 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
325     & SQUEEZE_RIGHT , mythid)
326     write(msgbuf,'(a)')
327 dimitri 1.4 &' Salt flux forcing is read from file: '
328 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
329     & SQUEEZE_RIGHT , mythid)
330     write(msgbuf,'(a,a,a)')
331 dimitri 1.4 &' >> ',sfluxfile(1:il),' <<'
332 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
333     & SQUEEZE_RIGHT , mythid)
334 dimitri 1.4 #endif
335 heimbach 1.1
336 dimitri 1.4 #if defined(ALLOW_ATM_TEMP) || defined(SHORTWAVE_HEATING)
337 jmc 1.12 C-- Net shortwave.
338 heimbach 1.1 il = ilnblnk(swfluxfile)
339 jmc 1.12 write(msgbuf,'(a)') ' '
340 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
341     & SQUEEZE_RIGHT , mythid)
342 dimitri 1.7 write(msgbuf,'(a,f12.0)')
343 dimitri 1.4 &' Net shortwave flux forcing starts at ',
344 dimitri 1.7 & swfluxstartdate
345 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
346     & SQUEEZE_RIGHT , mythid)
347 dimitri 1.7 write(msgbuf,'(a,f12.0)')
348 dimitri 1.4 &' Net shortwave flux forcing period is ',
349 dimitri 1.7 & swfluxperiod
350 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
351     & SQUEEZE_RIGHT , mythid)
352     write(msgbuf,'(a)')
353 dimitri 1.4 &' Net shortwave flux forcing is read from file:'
354 heimbach 1.1 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 dimitri 1.4 #endif
361 heimbach 1.1
362 dimitri 1.4 #ifdef ALLOW_ATM_WIND
363 jmc 1.12 C-- Zonal wind.
364 dimitri 1.4 il = ilnblnk(uwindfile)
365 jmc 1.12 write(msgbuf,'(a)') ' '
366 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
367     & SQUEEZE_RIGHT , mythid)
368 dimitri 1.7 write(msgbuf,'(a,f12.0)')
369 dimitri 1.4 &' Zonal wind forcing starts at ',
370 dimitri 1.7 & uwindstartdate
371 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
372     & SQUEEZE_RIGHT , mythid)
373 dimitri 1.7 write(msgbuf,'(a,f12.0)')
374 dimitri 1.4 &' Zonal wind forcing period is ',
375 dimitri 1.7 & uwindperiod
376 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
377     & SQUEEZE_RIGHT , mythid)
378     write(msgbuf,'(a)')
379 dimitri 1.4 &' Zonal wind forcing is read from file:'
380 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
381     & SQUEEZE_RIGHT , mythid)
382     write(msgbuf,'(a,a,a)')
383 dimitri 1.4 &' >> ',uwindfile(1:il),' <<'
384 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
385     & SQUEEZE_RIGHT , mythid)
386    
387 jmc 1.12 C-- Meridional wind.
388 dimitri 1.4 il = ilnblnk(vwindfile)
389 jmc 1.12 write(msgbuf,'(a)') ' '
390 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
391     & SQUEEZE_RIGHT , mythid)
392 dimitri 1.7 write(msgbuf,'(a,f12.0)')
393 dimitri 1.4 &' Meridional wind forcing starts at ',
394 dimitri 1.7 & vwindstartdate
395 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
396     & SQUEEZE_RIGHT , mythid)
397 dimitri 1.7 write(msgbuf,'(a,f12.0)')
398 dimitri 1.4 &' Meridional wind forcing period is ',
399 dimitri 1.7 & vwindperiod
400 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
401     & SQUEEZE_RIGHT , mythid)
402     write(msgbuf,'(a)')
403 dimitri 1.4 &' Meridional wind forcing is read from file:'
404 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
405     & SQUEEZE_RIGHT , mythid)
406     write(msgbuf,'(a,a,a)')
407 dimitri 1.4 &' >> ',vwindfile(1:il),' <<'
408 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
409     & SQUEEZE_RIGHT , mythid)
410 dimitri 1.4 #endif
411 heimbach 1.1
412 dimitri 1.4 #ifdef ALLOW_ATM_TEMP
413 jmc 1.12 C-- Atmospheric temperature.
414 dimitri 1.4 il = ilnblnk(atempfile)
415 jmc 1.12 write(msgbuf,'(a)') ' '
416 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
417     & SQUEEZE_RIGHT , mythid)
418 dimitri 1.7 write(msgbuf,'(a,f12.0)')
419 dimitri 1.4 &' Atmospheric temperature starts at ',
420 dimitri 1.7 & atempstartdate
421 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
422     & SQUEEZE_RIGHT , mythid)
423 dimitri 1.7 write(msgbuf,'(a,f12.0)')
424 dimitri 1.4 &' Atmospheric temperature period is ',
425 dimitri 1.7 & atempperiod
426 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
427     & SQUEEZE_RIGHT , mythid)
428     write(msgbuf,'(a)')
429 dimitri 1.4 &' Atmospheric temperature is read from file:'
430 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
431     & SQUEEZE_RIGHT , mythid)
432     write(msgbuf,'(a,a,a)')
433 dimitri 1.4 &' >> ',atempfile(1:il),' <<'
434 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
435     & SQUEEZE_RIGHT , mythid)
436    
437 jmc 1.12 C-- Atmospheric specific humidity.
438 dimitri 1.4 il = ilnblnk(aqhfile)
439 jmc 1.12 write(msgbuf,'(a)') ' '
440 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
441     & SQUEEZE_RIGHT , mythid)
442 dimitri 1.7 write(msgbuf,'(a,f12.0)')
443 dimitri 1.4 &' Atmospheric specific humidity starts at ',
444 dimitri 1.7 & aqhstartdate
445 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
446     & SQUEEZE_RIGHT , mythid)
447 dimitri 1.7 write(msgbuf,'(a,f12.0)')
448 dimitri 1.4 &' Atmospheric specific humidity period is ',
449 dimitri 1.7 & aqhperiod
450 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
451     & SQUEEZE_RIGHT , mythid)
452     write(msgbuf,'(a)')
453 dimitri 1.4 &' Atmospheric specific humidity is read from file:'
454 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
455     & SQUEEZE_RIGHT , mythid)
456     write(msgbuf,'(a,a,a)')
457 dimitri 1.4 &' >> ',aqhfile(1:il),' <<'
458 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
459     & SQUEEZE_RIGHT , mythid)
460    
461 jmc 1.12 C-- Net longwave.
462 dimitri 1.4 il = ilnblnk(lwfluxfile)
463 jmc 1.12 write(msgbuf,'(a)') ' '
464 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
465     & SQUEEZE_RIGHT , mythid)
466 dimitri 1.7 write(msgbuf,'(a,f12.0)')
467 dimitri 1.4 &' Net longwave flux forcing starts at ',
468 dimitri 1.7 & lwfluxstartdate
469 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
470     & SQUEEZE_RIGHT , mythid)
471 dimitri 1.7 write(msgbuf,'(a,f12.0)')
472 dimitri 1.4 &' Net longwave flux forcing period is ',
473 dimitri 1.7 & lwfluxperiod
474 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
475     & SQUEEZE_RIGHT , mythid)
476     write(msgbuf,'(a)')
477 dimitri 1.4 &' Net longwave flux forcing is read from file:'
478 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
479     & SQUEEZE_RIGHT , mythid)
480     write(msgbuf,'(a,a,a)')
481 dimitri 1.4 &' >> ',lwfluxfile(1:il),' <<'
482 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
483     & SQUEEZE_RIGHT , mythid)
484    
485 jmc 1.12 C-- Precipitation.
486 dimitri 1.4 il = ilnblnk(precipfile)
487 jmc 1.12 write(msgbuf,'(a)') ' '
488 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
489     & SQUEEZE_RIGHT , mythid)
490 dimitri 1.7 write(msgbuf,'(a,f12.0)')
491 dimitri 1.4 &' Precipitation data set starts at ',
492 dimitri 1.7 & precipstartdate
493 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
494     & SQUEEZE_RIGHT , mythid)
495 dimitri 1.7 write(msgbuf,'(a,f12.0)')
496 dimitri 1.4 &' Precipitation data period is ',
497 dimitri 1.7 & precipperiod
498 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
499     & SQUEEZE_RIGHT , mythid)
500     write(msgbuf,'(a)')
501 dimitri 1.4 &' Precipitation data is read from file: '
502 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
503     & SQUEEZE_RIGHT , mythid)
504     write(msgbuf,'(a,a,a)')
505 dimitri 1.4 &' >> ',precipfile(1:il),' <<'
506 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
507     & SQUEEZE_RIGHT , mythid)
508 dimitri 1.4 #endif
509 heimbach 1.1
510 jmc 1.12 C-- Evaporation.
511     write(msgbuf,'(a)') ' '
512     call print_message( msgbuf, standardmessageunit,
513     & SQUEEZE_RIGHT , mythid)
514 dimitri 1.4 #ifdef EXF_READ_EVAP
515 heimbach 1.1 write(msgbuf,'(a)')
516 dimitri 1.4 &'// EXF_READ_EVAP: defined'
517 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
518     & SQUEEZE_RIGHT , mythid)
519 dimitri 1.4 il = ilnblnk(evapfile)
520 dimitri 1.7 write(msgbuf,'(a,f12.0)')
521 dimitri 1.4 &' Evaporation starts at ',
522 dimitri 1.7 & evapstartdate
523 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
524     & SQUEEZE_RIGHT , mythid)
525 dimitri 1.7 write(msgbuf,'(a,f12.0)')
526 dimitri 1.4 &' Evaporation period is ',
527 dimitri 1.7 & evapperiod
528 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
529     & SQUEEZE_RIGHT , mythid)
530     write(msgbuf,'(a)')
531 dimitri 1.4 &' Evaporation is read from file:'
532 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
533     & SQUEEZE_RIGHT , mythid)
534     write(msgbuf,'(a,a,a)')
535 dimitri 1.4 &' >> ',evapfile(1:il),' <<'
536 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
537     & SQUEEZE_RIGHT , mythid)
538 dimitri 1.4 #else
539 heimbach 1.1 write(msgbuf,'(a)')
540 dimitri 1.4 &'// EXF_READ_EVAP: NOT defined'
541 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
542     & SQUEEZE_RIGHT , mythid)
543 dimitri 1.4 #endif
544    
545 jmc 1.12 C-- Runoff.
546     write(msgbuf,'(a)') ' '
547     call print_message( msgbuf, standardmessageunit,
548     & SQUEEZE_RIGHT , mythid)
549 dimitri 1.4 #ifdef ALLOW_RUNOFF
550 heimbach 1.1 write(msgbuf,'(a)')
551 dimitri 1.4 &'// ALLOW_RUNOFF: defined'
552 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
553     & SQUEEZE_RIGHT , mythid)
554 dimitri 1.4 il = ilnblnk(runofffile)
555 dimitri 1.7 write(msgbuf,'(a,f12.0)')
556 dimitri 1.4 &' Runnoff starts at ',
557 dimitri 1.7 & runoffstartdate
558 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
559     & SQUEEZE_RIGHT , mythid)
560 dimitri 1.7 write(msgbuf,'(a,f12.0)')
561 dimitri 1.4 &' Runoff period is ',
562 dimitri 1.7 & runoffperiod
563 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
564     & SQUEEZE_RIGHT , mythid)
565     write(msgbuf,'(a)')
566 dimitri 1.4 &' Runoff is read from file:'
567 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
568     & SQUEEZE_RIGHT , mythid)
569     write(msgbuf,'(a,a,a)')
570 dimitri 1.4 &' >> ',runofffile(1:il),' <<'
571 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
572     & SQUEEZE_RIGHT , mythid)
573 dimitri 1.4 #else /* ALLOW_RUNOFF */
574 heimbach 1.1 write(msgbuf,'(a)')
575 dimitri 1.4 &'// ALLOW_RUNOFF: NOT defined'
576 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
577     & SQUEEZE_RIGHT , mythid)
578 dimitri 1.4 #endif /* ALLOW_RUNOFF */
579 heimbach 1.1
580 dimitri 1.4 #ifdef DOWNWARD_RADIATION
581 jmc 1.12 C-- Downward shortwave.
582 dimitri 1.4 il = ilnblnk(swdownfile)
583 jmc 1.12 write(msgbuf,'(a)') ' '
584 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
585     & SQUEEZE_RIGHT , mythid)
586 dimitri 1.7 write(msgbuf,'(a,f12.0)')
587 dimitri 1.4 &' Downward shortwave flux forcing starts at ',
588 dimitri 1.7 & swdownstartdate
589 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
590     & SQUEEZE_RIGHT , mythid)
591 dimitri 1.7 write(msgbuf,'(a,f12.0)')
592 dimitri 1.4 &' Downward shortwave flux forcing period is ',
593 dimitri 1.7 & swdownperiod
594 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
595     & SQUEEZE_RIGHT , mythid)
596     write(msgbuf,'(a)')
597 dimitri 1.4 &' Downward shortwave flux forcing is read from file:'
598 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
599     & SQUEEZE_RIGHT , mythid)
600     write(msgbuf,'(a,a,a)')
601 dimitri 1.4 &' >> ',swdownfile(1:il),' <<'
602 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
603     & SQUEEZE_RIGHT , mythid)
604    
605 jmc 1.12 C-- Downward longwave.
606 dimitri 1.4 il = ilnblnk(lwdownfile)
607 jmc 1.12 write(msgbuf,'(a)') ' '
608 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
609     & SQUEEZE_RIGHT , mythid)
610 dimitri 1.7 write(msgbuf,'(a,f12.0)')
611 dimitri 1.4 &' Downward longwave flux forcing starts at ',
612 dimitri 1.7 & lwdownstartdate
613 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
614     & SQUEEZE_RIGHT , mythid)
615 dimitri 1.7 write(msgbuf,'(a,f12.0)')
616 dimitri 1.4 &' Downward longwave flux forcing period is ',
617 dimitri 1.7 & lwdownperiod
618 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
619     & SQUEEZE_RIGHT , mythid)
620     write(msgbuf,'(a)')
621 dimitri 1.4 &' Downward longwave flux forcing is read from file:'
622 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
623     & SQUEEZE_RIGHT , mythid)
624     write(msgbuf,'(a,a,a)')
625 dimitri 1.4 &' >> ',lwdownfile(1:il),' <<'
626 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
627     & SQUEEZE_RIGHT , mythid)
628 dimitri 1.4 #endif
629 heimbach 1.1
630 dimitri 1.4 #ifdef ATMOSPHERIC_LOADING
631 jmc 1.12 C-- Atmospheric pressure.
632 dimitri 1.4 il = ilnblnk(apressurefile)
633 jmc 1.12 write(msgbuf,'(a)') ' '
634 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
635     & SQUEEZE_RIGHT , mythid)
636 dimitri 1.7 write(msgbuf,'(a,f12.0)')
637 dimitri 1.4 &' Atmospheric pressure forcing starts at ',
638 dimitri 1.7 & apressurestartdate
639 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
640     & SQUEEZE_RIGHT , mythid)
641 dimitri 1.7 write(msgbuf,'(a,f12.0)')
642 dimitri 1.4 &' Atmospheric pressure forcing period is ',
643 dimitri 1.7 & apressureperiod
644 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
645     & SQUEEZE_RIGHT , mythid)
646     write(msgbuf,'(a)')
647 dimitri 1.4 &' Atmospheric pressureforcing is read from file:'
648 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
649     & SQUEEZE_RIGHT , mythid)
650     write(msgbuf,'(a,a,a)')
651 dimitri 1.4 &' >> ',apressurefile(1:il),' <<'
652 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
653     & SQUEEZE_RIGHT , mythid)
654     #endif
655    
656 jmc 1.12 write(msgbuf,'(a)') ' '
657 heimbach 1.1 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 jmc 1.12 write(msgbuf,'(a)') ' '
672 heimbach 1.1 call print_message( msgbuf, standardmessageunit,
673     & SQUEEZE_RIGHT , mythid)
674    
675 jmc 1.12 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
676    
677 heimbach 1.11 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 jmc 1.12 write(msgbuf,'(a)') ' '
692 heimbach 1.11 call print_message( msgbuf, standardmessageunit,
693     & SQUEEZE_RIGHT , mythid)
694    
695 jmc 1.12 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 heimbach 1.11
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 jmc 1.12 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 heimbach 1.11
726     #ifdef ALLOW_CLIMSST_RELAXATION
727 jmc 1.12 C Relaxation to SST climatology.
728 heimbach 1.11 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 jmc 1.12 C Relaxation to SSS climatology.
755 heimbach 1.11 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 jmc 1.12 write(msgbuf,'(a)') ' '
781 heimbach 1.11 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 jmc 1.12 write(msgbuf,'(a)') ' '
796 heimbach 1.11 call print_message( msgbuf, standardmessageunit,
797     & SQUEEZE_RIGHT , mythid)
798    
799 jmc 1.12 _END_MASTER( myThid )
800    
801     RETURN
802     END
  ViewVC Help
Powered by ViewVC 1.1.22