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