/[MITgcm]/MITgcm/pkg/exf/exf_summary.F
ViewVC logotype

Diff of /MITgcm/pkg/exf/exf_summary.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.6 by edhill, Thu Oct 9 04:19:19 2003 UTC revision 1.14 by jmc, Sat May 12 23:17:18 2007 UTC
# Line 1  Line 1 
1  c $Header$  C $Header$
2    C $Name$
3    
4  #include "EXF_OPTIONS.h"  #include "EXF_OPTIONS.h"
5    
6    
7        subroutine exf_Summary(        SUBROUTINE EXF_SUMMARY( myThid )
      I                        mythid  
      &                      )  
8    
9  c     ==================================================================  c     ==================================================================
10  c     SUBROUTINE exf_Summary  c     SUBROUTINE exf_summary
11  c     ==================================================================  c     ==================================================================
12  c  c
13  c     o List all the settings of the external forcing.  c     o List all the settings of the external forcing.
# Line 26  c     changed: Dimitris Menemenlis menem Line 25  c     changed: Dimitris Menemenlis menem
25  c              - modifications for using pkg/exf with pkg/seaice  c              - modifications for using pkg/exf with pkg/seaice
26  c  c
27  c     ==================================================================  c     ==================================================================
28  c     SUBROUTINE exf_Summary  c     SUBROUTINE exf_summary
29  c     ==================================================================  c     ==================================================================
30    
31        implicit none        implicit none
32    
33  c     == global variables ==  C     == global variables ==
34    
35  #include "EEPARAMS.h"  #include "EEPARAMS.h"
36  #include "SIZE.h"  #include "SIZE.h"
37  #include "cal.h"  #include "cal.h"
38  #include "exf.h"  #include "EXF_CONSTANTS.h"
39  #include "exf_param.h"  #include "EXF_PARAM.h"
40    
41  c     == routine arguments ==  C     == routine arguments ==
42    
43  c     mythid  - thread number for this instance of the routine.  C     myThid  - thread number for this instance of the routine.
44    
45        integer mythid        integer myThid
46    
47  c     == local variables ==  C     == local variables ==
48    
       integer i  
49        integer il        integer il
       integer timeint(4)  
50    
51        character*(max_len_mbuf) msgbuf        character*(max_len_mbuf) msgbuf
52    
53  c     == external ==  C     == external ==
54    
55        integer  ilnblnk        integer  ilnblnk
56        external ilnblnk        external ilnblnk
57    
58  c     == end of interface ==  C     == end of interface ==
59    
60        write(msgbuf,'(a)')        _BEGIN_MASTER( myThid )
61       &' '  
62          write(msgbuf,'(a)') ' '
63        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
64       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
65        write(msgbuf,'(a)')        write(msgbuf,'(a)')
# Line 76  c     == end of interface == Line 74  c     == end of interface ==
74       &'// ======================================================='       &'// ======================================================='
75        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
76       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
77        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
       call print_message( msgbuf, standardmessageunit,  
      &                    SQUEEZE_RIGHT , mythid)  
       write(msgbuf,'(a,a)')  
      &'External forcing version: ',externalforcingversion  
       call print_message( msgbuf, standardmessageunit,  
      &                    SQUEEZE_RIGHT , mythid)  
       write(msgbuf,'(a,a)')  
      &'Uses Calendar version:    ',usescalendarversion  
       call print_message( msgbuf, standardmessageunit,  
      &                    SQUEEZE_RIGHT , mythid)  
       write(msgbuf,'(a)')  
      &' '  
78        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
79       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
80    
81  c--   Print settings of some CPP flags.  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( 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          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          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          CALL WRITE_0D_R8( exf_scal_BulkCdn, INDEX_NONE,
147         &  'exf_scal_BulkCdn=',
148         &  ' /* Drag coefficient scaling factor [-] */')
149          CALL WRITE_0D_R8( zolmin, INDEX_NONE, 'zolmin =',
150         &  ' /* minimum stability parameter [?] */')
151          CALL WRITE_0D_R8( psim_fac, INDEX_NONE, 'psim_fac =',
152         &  ' /* coef used in turbulent fluxes calculation [-] */')
153          CALL WRITE_0D_R8( zref, INDEX_NONE, 'zref =',
154         & '  /* reference height [ m ] */')
155          CALL WRITE_0D_R8( hu, INDEX_NONE, 'hu =',
156         & '  /* height of mean wind [ m ] */')
157          CALL WRITE_0D_R8( ht, INDEX_NONE, 'ht =',
158         & '  /* height of mean temperature [ m ] */')
159          CALL WRITE_0D_R8( hq, INDEX_NONE, 'hq =',
160         & '  /* height of mean spec.humidity [ m ] */')
161          CALL WRITE_0D_R8( umin, INDEX_NONE, 'uMin =',
162         &  ' /* minimum wind speed [m/s] */')
163          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          CALL WRITE_0D_R8( exf_iceCh, INDEX_NONE, 'exf_iceCh =',
171         &  ' /* 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          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    
194  #ifdef ALLOW_ATM_TEMP  #ifdef ALLOW_ATM_TEMP
195        write(msgbuf,'(a)')        write(msgbuf,'(a)')
# Line 143  c--   Print settings of some CPP flags. Line 239  c--   Print settings of some CPP flags.
239       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
240  #endif  #endif
241    
242  c--   For each data set used the summary prints the calendar data  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
243  c     and the corresponding file from which the data will be read.  
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    
247  #ifndef ALLOW_ATM_WIND  #ifndef ALLOW_ATM_WIND
248  c--   Zonal wind stress.  C--   Zonal wind stress.
249        il = ilnblnk(ustressfile)        il = ilnblnk(ustressfile)
250        call cal_TimeInterval( ustressperiod, 'secs', timeint, mythid )        write(msgbuf,'(a)') ' '
       write(msgbuf,'(a)')  
      &' '  
251        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
252       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
253        write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')        write(msgbuf,'(a,f12.0)')
254       &'   Zonal wind stress forcing starts at         ',       &'   Zonal wind stress forcing starts at         ',
255       &    (ustressstartdate(i), i=1,2), dayofweek(ustressstartdate(4)),       &    ustressstartdate
      &    '.'  
256        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
257       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
258        write(msgbuf,'(a,i9.8,i7.6)')        write(msgbuf,'(a,f12.0)')
259       &'   Zonal wind stress forcing period is         ',       &'   Zonal wind stress forcing period is         ',
260       &    (timeint(i), i=1,2)       &    ustressperiod
261        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
262       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
263        write(msgbuf,'(a)')        write(msgbuf,'(a)')
# Line 174  c--   Zonal wind stress. Line 269  c--   Zonal wind stress.
269        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
270       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
271    
272  c--   Meridional wind stress.  C--   Meridional wind stress.
273        il = ilnblnk(vstressfile)        il = ilnblnk(vstressfile)
274        call cal_TimeInterval( vstressperiod, 'secs', timeint, mythid )        write(msgbuf,'(a)') ' '
       write(msgbuf,'(a)')  
      &' '  
275        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
276       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
277        write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')        write(msgbuf,'(a,f12.0)')
278       &'   Meridional wind stress forcing starts at    ',       &'   Meridional wind stress forcing starts at    ',
279       &    (vstressstartdate(i), i=1,2), dayofweek(vstressstartdate(4)),       &    vstressstartdate
      &    '.'  
280        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
281       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
282        write(msgbuf,'(a,i9.8,i7.6)')        write(msgbuf,'(a,f12.0)')
283       &'   Meridional wind stress forcing period is    ',       &'   Meridional wind stress forcing period is    ',
284       &    (timeint(i), i=1,2)       &    vstressperiod
285        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
286       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
287        write(msgbuf,'(a)')        write(msgbuf,'(a)')
# Line 203  c--   Meridional wind stress. Line 295  c--   Meridional wind stress.
295  #endif  #endif
296    
297  #ifndef ALLOW_ATM_TEMP  #ifndef ALLOW_ATM_TEMP
298  c--   Heat flux.  C--   Heat flux.
299        il = ilnblnk(hfluxfile)        il = ilnblnk(hfluxfile)
300        call cal_TimeInterval( hfluxperiod, 'secs', timeint, mythid )        write(msgbuf,'(a)') ' '
       write(msgbuf,'(a)')  
      &' '  
301        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
302       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
303        write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')        write(msgbuf,'(a,f12.0)')
304       &'   Heat flux forcing starts at                ',       &'   Heat flux forcing starts at                ',
305       &    (hfluxstartdate(i), i=1,2), dayofweek(hfluxstartdate(4)),'.'       &    hfluxstartdate
306        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
307       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
308        write(msgbuf,'(a,i9.8,i7.6)')        write(msgbuf,'(a,f12.0)')
309       &'   Heat flux forcing period is                 ',       &'   Heat flux forcing period is                 ',
310       &    (timeint(i), i=1,2)       &    hfluxperiod
311        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
312       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
313        write(msgbuf,'(a)')        write(msgbuf,'(a)')
# Line 229  c--   Heat flux. Line 319  c--   Heat flux.
319        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
320       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
321    
322  c--   Salt flux.  C--   Salt flux.
323        il = ilnblnk(sfluxfile)        il = ilnblnk(sfluxfile)
324        call cal_TimeInterval( sfluxperiod, 'secs', timeint, mythid )        write(msgbuf,'(a)') ' '
       write(msgbuf,'(a)')  
      &' '  
325        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
326       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
327        write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')        write(msgbuf,'(a,f12.0)')
328       &'   Salt flux forcing starts at                 ',       &'   Salt flux forcing starts at                 ',
329       &    (sfluxstartdate(i), i=1,2), dayofweek(sfluxstartdate(4)),'.'       &    sfluxstartdate
330        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
331       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
332        write(msgbuf,'(a,i9.8,i7.6)')        write(msgbuf,'(a,f12.0)')
333       &'   Salt flux forcing period is                 ',       &'   Salt flux forcing period is                 ',
334       &    (timeint(i), i=1,2)       &    sfluxperiod
335        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
336       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
337        write(msgbuf,'(a)')        write(msgbuf,'(a)')
# Line 257  c--   Salt flux. Line 345  c--   Salt flux.
345  #endif  #endif
346    
347  #if defined(ALLOW_ATM_TEMP) || defined(SHORTWAVE_HEATING)  #if defined(ALLOW_ATM_TEMP) || defined(SHORTWAVE_HEATING)
348  c--   Net shortwave.  C--   Net shortwave.
349        il = ilnblnk(swfluxfile)        il = ilnblnk(swfluxfile)
350        call cal_TimeInterval( swfluxperiod, 'secs', timeint, mythid )        write(msgbuf,'(a)') ' '
       write(msgbuf,'(a)')  
      &' '  
351        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
352       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
353        write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')        write(msgbuf,'(a,f12.0)')
354       &'   Net shortwave flux forcing starts at      ',       &'   Net shortwave flux forcing starts at      ',
355       &    (swfluxstartdate(i), i=1,2), dayofweek(swfluxstartdate(4)),       &    swfluxstartdate
      &    '.'  
356        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
357       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
358        write(msgbuf,'(a,i9.8,i7.6)')        write(msgbuf,'(a,f12.0)')
359       &'   Net shortwave flux forcing period is      ',       &'   Net shortwave flux forcing period is      ',
360       &    (timeint(i), i=1,2)       &    swfluxperiod
361        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
362       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
363        write(msgbuf,'(a)')        write(msgbuf,'(a)')
# Line 286  c--   Net shortwave. Line 371  c--   Net shortwave.
371  #endif  #endif
372    
373  #ifdef ALLOW_ATM_WIND  #ifdef ALLOW_ATM_WIND
374  c--   Zonal wind.  C--   Zonal wind.
375        il = ilnblnk(uwindfile)        il = ilnblnk(uwindfile)
376        call cal_TimeInterval( uwindperiod, 'secs', timeint, mythid )        write(msgbuf,'(a)') ' '
       write(msgbuf,'(a)')  
      &' '  
377        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
378       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
379        write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')        write(msgbuf,'(a,f12.0)')
380       &'   Zonal wind forcing starts at                ',       &'   Zonal wind forcing starts at                ',
381       &    (uwindstartdate(i), i=1,2), dayofweek(uwindstartdate(4)),       &    uwindstartdate
      &    '.'  
382        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
383       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
384        write(msgbuf,'(a,i9.8,i7.6)')        write(msgbuf,'(a,f12.0)')
385       &'   Zonal wind forcing period is                ',       &'   Zonal wind forcing period is                ',
386       &    (timeint(i), i=1,2)       &    uwindperiod
387        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
388       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
389        write(msgbuf,'(a)')        write(msgbuf,'(a)')
# Line 313  c--   Zonal wind. Line 395  c--   Zonal wind.
395        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
396       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
397    
398  c--   Meridional wind.  C--   Meridional wind.
399        il = ilnblnk(vwindfile)        il = ilnblnk(vwindfile)
400        call cal_TimeInterval( vwindperiod, 'secs', timeint, mythid )        write(msgbuf,'(a)') ' '
       write(msgbuf,'(a)')  
      &' '  
401        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
402       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
403        write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')        write(msgbuf,'(a,f12.0)')
404       &'   Meridional wind forcing starts at           ',       &'   Meridional wind forcing starts at           ',
405       &    (vwindstartdate(i), i=1,2), dayofweek(vwindstartdate(4)),       &    vwindstartdate
      &    '.'  
406        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
407       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
408        write(msgbuf,'(a,i9.8,i7.6)')        write(msgbuf,'(a,f12.0)')
409       &'   Meridional wind forcing period is           ',       &'   Meridional wind forcing period is           ',
410       &    (timeint(i), i=1,2)       &    vwindperiod
411        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
412       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
413        write(msgbuf,'(a)')        write(msgbuf,'(a)')
# Line 342  c--   Meridional wind. Line 421  c--   Meridional wind.
421  #endif  #endif
422    
423  #ifdef ALLOW_ATM_TEMP  #ifdef ALLOW_ATM_TEMP
424  c--   Atmospheric temperature.  C--   Atmospheric temperature.
425        il = ilnblnk(atempfile)        il = ilnblnk(atempfile)
426        call cal_TimeInterval( atempperiod, 'secs', timeint, mythid )        write(msgbuf,'(a)') ' '
       write(msgbuf,'(a)')  
      &' '  
427        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
428       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
429        write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')        write(msgbuf,'(a,f12.0)')
430       &'   Atmospheric temperature starts at           ',       &'   Atmospheric temperature starts at           ',
431       &    (atempstartdate(i), i=1,2), dayofweek(atempstartdate(4)),       &    atempstartdate
      &    '.'  
432        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
433       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
434        write(msgbuf,'(a,i9.8,i7.6)')        write(msgbuf,'(a,f12.0)')
435       &'   Atmospheric temperature period is           ',       &'   Atmospheric temperature period is           ',
436       &    (timeint(i), i=1,2)       &    atempperiod
437        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
438       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
439        write(msgbuf,'(a)')        write(msgbuf,'(a)')
# Line 369  c--   Atmospheric temperature. Line 445  c--   Atmospheric temperature.
445        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
446       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
447    
448  c--   Atmospheric specific humidity.  C--   Atmospheric specific humidity.
449        il = ilnblnk(aqhfile)        il = ilnblnk(aqhfile)
450        call cal_TimeInterval( aqhperiod, 'secs', timeint, mythid )        write(msgbuf,'(a)') ' '
       write(msgbuf,'(a)')  
      &' '  
451        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
452       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
453        write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')        write(msgbuf,'(a,f12.0)')
454       &'   Atmospheric specific humidity starts at     ',       &'   Atmospheric specific humidity starts at     ',
455       &    (aqhstartdate(i), i=1,2), dayofweek(aqhstartdate(4)),       &    aqhstartdate
      &    '.'  
456        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
457       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
458        write(msgbuf,'(a,i9.8,i7.6)')        write(msgbuf,'(a,f12.0)')
459       &'   Atmospheric specific humidity period is     ',       &'   Atmospheric specific humidity period is     ',
460       &    (timeint(i), i=1,2)       &    aqhperiod
461        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
462       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
463        write(msgbuf,'(a)')        write(msgbuf,'(a)')
# Line 396  c--   Atmospheric specific humidity. Line 469  c--   Atmospheric specific humidity.
469        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
470       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
471    
472  c--   Net longwave.  C--   Net longwave.
473        il = ilnblnk(lwfluxfile)        il = ilnblnk(lwfluxfile)
474        call cal_TimeInterval( lwfluxperiod, 'secs', timeint, mythid )        write(msgbuf,'(a)') ' '
       write(msgbuf,'(a)')  
      &' '  
475        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
476       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
477        write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')        write(msgbuf,'(a,f12.0)')
478       &'   Net longwave flux forcing starts at       ',       &'   Net longwave flux forcing starts at       ',
479       &    (lwfluxstartdate(i), i=1,2), dayofweek(lwfluxstartdate(4)),       &    lwfluxstartdate
      &    '.'  
480        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
481       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
482        write(msgbuf,'(a,i9.8,i7.6)')        write(msgbuf,'(a,f12.0)')
483       &'   Net longwave flux forcing period is       ',       &'   Net longwave flux forcing period is       ',
484       &    (timeint(i), i=1,2)       &    lwfluxperiod
485        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
486       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
487        write(msgbuf,'(a)')        write(msgbuf,'(a)')
# Line 423  c--   Net longwave. Line 493  c--   Net longwave.
493        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
494       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
495    
496  c--   Precipitation.  C--   Precipitation.
497        il = ilnblnk(precipfile)        il = ilnblnk(precipfile)
498        call cal_TimeInterval( precipperiod, 'secs', timeint, mythid )        write(msgbuf,'(a)') ' '
       write(msgbuf,'(a)')  
      &' '  
499        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
500       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
501        write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')        write(msgbuf,'(a,f12.0)')
502       &'   Precipitation data set starts at            ',       &'   Precipitation data set starts at            ',
503       &    (precipstartdate(i), i=1,2), dayofweek(precipstartdate(4)),       &    precipstartdate
      &'.'  
504        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
505       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
506        write(msgbuf,'(a,i9.8,i7.6)')        write(msgbuf,'(a,f12.0)')
507       &'   Precipitation data period is                ',       &'   Precipitation data period is                ',
508       &    (timeint(i), i=1,2)       &    precipperiod
509        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
510       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
511        write(msgbuf,'(a)')        write(msgbuf,'(a)')
# Line 451  c--   Precipitation. Line 518  c--   Precipitation.
518       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
519  #endif  #endif
520    
521    C--   Evaporation.
522          write(msgbuf,'(a)') ' '
523          call print_message( msgbuf, standardmessageunit,
524         &                    SQUEEZE_RIGHT , mythid)
525  #ifdef EXF_READ_EVAP  #ifdef EXF_READ_EVAP
 c--   Evaporation.  
526        write(msgbuf,'(a)')        write(msgbuf,'(a)')
527       &'// EXF_READ_EVAP:                      defined'       &'// EXF_READ_EVAP:                      defined'
528        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
529       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
530        il = ilnblnk(evapfile)        il = ilnblnk(evapfile)
531        call cal_TimeInterval( evapperiod, 'secs', timeint, mythid )        write(msgbuf,'(a,f12.0)')
   
       write(msgbuf,'(a)')  
      &' '  
       call print_message( msgbuf, standardmessageunit,  
      &                    SQUEEZE_RIGHT , mythid)  
       write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')  
532       &'   Evaporation starts at     ',       &'   Evaporation starts at     ',
533       &    (evapstartdate(i), i=1,2), dayofweek(evapstartdate(4)),       &    evapstartdate
      &    '.'  
534        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
535       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
536        write(msgbuf,'(a,i9.8,i7.6)')        write(msgbuf,'(a,f12.0)')
537       &'   Evaporation period is     ',       &'   Evaporation period is     ',
538       &    (timeint(i), i=1,2)       &    evapperiod
539        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
540       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
541        write(msgbuf,'(a)')        write(msgbuf,'(a)')
# Line 490  c--   Evaporation. Line 553  c--   Evaporation.
553       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
554  #endif  #endif
555    
556    C--   Runoff.
557          write(msgbuf,'(a)') ' '
558          call print_message( msgbuf, standardmessageunit,
559         &                    SQUEEZE_RIGHT , mythid)
560  #ifdef ALLOW_RUNOFF  #ifdef ALLOW_RUNOFF
 c--   Runoff.  
561        write(msgbuf,'(a)')        write(msgbuf,'(a)')
562       &'// ALLOW_RUNOFF:                       defined'       &'// ALLOW_RUNOFF:                       defined'
563        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
564       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
565        il = ilnblnk(runofffile)        il = ilnblnk(runofffile)
566        call cal_TimeInterval( runoffperiod, 'secs', timeint, mythid )        write(msgbuf,'(a,f12.0)')
   
       write(msgbuf,'(a)')  
      &' '  
       call print_message( msgbuf, standardmessageunit,  
      &                    SQUEEZE_RIGHT , mythid)  
       write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')  
567       &'   Runnoff starts at     ',       &'   Runnoff starts at     ',
568       &    (runoffstartdate(i), i=1,2), dayofweek(runoffstartdate(4)),       &    runoffstartdate
      &    '.'  
569        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
570       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
571        write(msgbuf,'(a,i9.8,i7.6)')        write(msgbuf,'(a,f12.0)')
572       &'   Runoff period is     ',       &'   Runoff period is     ',
573       &    (timeint(i), i=1,2)       &    runoffperiod
574        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
575       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
576        write(msgbuf,'(a)')        write(msgbuf,'(a)')
# Line 530  c--   Runoff. Line 589  c--   Runoff.
589  #endif /* ALLOW_RUNOFF */  #endif /* ALLOW_RUNOFF */
590    
591  #ifdef DOWNWARD_RADIATION  #ifdef DOWNWARD_RADIATION
592  c--   Downward shortwave.  C--   Downward shortwave.
593        il = ilnblnk(swdownfile)        il = ilnblnk(swdownfile)
594        call cal_TimeInterval( swdownperiod, 'secs', timeint, mythid )        write(msgbuf,'(a)') ' '
       write(msgbuf,'(a)')  
      &' '  
595        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
596       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
597        write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')        write(msgbuf,'(a,f12.0)')
598       &'   Downward shortwave flux forcing starts at      ',       &'   Downward shortwave flux forcing starts at      ',
599       &    (swdownstartdate(i), i=1,2), dayofweek(swdownstartdate(4)),       &    swdownstartdate
      &    '.'  
600        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
601       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
602        write(msgbuf,'(a,i9.8,i7.6)')        write(msgbuf,'(a,f12.0)')
603       &'   Downward shortwave flux forcing period is      ',       &'   Downward shortwave flux forcing period is      ',
604       &    (timeint(i), i=1,2)       &    swdownperiod
605        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
606       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
607        write(msgbuf,'(a)')        write(msgbuf,'(a)')
# Line 557  c--   Downward shortwave. Line 613  c--   Downward shortwave.
613        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
614       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
615    
616  c--   Downward longwave.  C--   Downward longwave.
617        il = ilnblnk(lwdownfile)        il = ilnblnk(lwdownfile)
618        call cal_TimeInterval( lwdownperiod, 'secs', timeint, mythid )        write(msgbuf,'(a)') ' '
       write(msgbuf,'(a)')  
      &' '  
619        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
620       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
621        write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')        write(msgbuf,'(a,f12.0)')
622       &'   Downward longwave flux forcing starts at       ',       &'   Downward longwave flux forcing starts at       ',
623       &    (lwdownstartdate(i), i=1,2), dayofweek(lwdownstartdate(4)),       &    lwdownstartdate
      &    '.'  
624        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
625       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
626        write(msgbuf,'(a,i9.8,i7.6)')        write(msgbuf,'(a,f12.0)')
627       &'   Downward longwave flux forcing period is       ',       &'   Downward longwave flux forcing period is       ',
628       &    (timeint(i), i=1,2)       &    lwdownperiod
629        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
630       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
631        write(msgbuf,'(a)')        write(msgbuf,'(a)')
# Line 586  c--   Downward longwave. Line 639  c--   Downward longwave.
639  #endif  #endif
640    
641  #ifdef ATMOSPHERIC_LOADING  #ifdef ATMOSPHERIC_LOADING
642  c--   Atmospheric pressure.  C--   Atmospheric pressure.
643        il = ilnblnk(apressurefile)        il = ilnblnk(apressurefile)
644        call cal_TimeInterval( apressureperiod, 'secs', timeint, mythid )        write(msgbuf,'(a)') ' '
       write(msgbuf,'(a)')  
      &' '  
645        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
646       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
647        write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')        write(msgbuf,'(a,f12.0)')
648       &'   Atmospheric pressure forcing starts at      ',       &'   Atmospheric pressure forcing starts at      ',
649       &    (apressurestartdate(i), i=1,2),       &    apressurestartdate
      &    dayofweek(apressurestartdate(4)), '.'  
650        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
651       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
652        write(msgbuf,'(a,i9.8,i7.6)')        write(msgbuf,'(a,f12.0)')
653       &'   Atmospheric pressure forcing period is      ',       &'   Atmospheric pressure forcing period is      ',
654       &    (timeint(i), i=1,2)       &    apressureperiod
655        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
656       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
657        write(msgbuf,'(a)')        write(msgbuf,'(a)')
# Line 614  c--   Atmospheric pressure. Line 664  c--   Atmospheric pressure.
664       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
665  #endif  #endif
666    
667        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
668        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
669       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
670        write(msgbuf,'(a)')        write(msgbuf,'(a)')
# Line 630  c--   Atmospheric pressure. Line 679  c--   Atmospheric pressure.
679       &'// ======================================================='       &'// ======================================================='
680        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
681       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
682          write(msgbuf,'(a)') ' '
683          call print_message( msgbuf, standardmessageunit,
684         &                    SQUEEZE_RIGHT , mythid)
685    
686    C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
687    
688          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          write(msgbuf,'(a)') ' '
703          call print_message( msgbuf, standardmessageunit,
704         &                    SQUEEZE_RIGHT , mythid)
705    
706    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    
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    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    
737    #ifdef ALLOW_CLIMSST_RELAXATION
738    C     Relaxation to SST climatology.
739          il = ilnblnk(climsstfile)
740        write(msgbuf,'(a)')        write(msgbuf,'(a)')
741       &' '       &' '
742        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
743       &                    SQUEEZE_RIGHT , mythid)       &                    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    C     Relaxation to SSS climatology.
766          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          write(msgbuf,'(a)') ' '
792          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          write(msgbuf,'(a)') ' '
807          call print_message( msgbuf, standardmessageunit,
808         &                    SQUEEZE_RIGHT , mythid)
809    
810          _END_MASTER( myThid )
811    
812        end        RETURN
813          END
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.14
  ViewVC Help
Powered by ViewVC 1.1.22