/[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.3 by dimitri, Sat Dec 28 10:11:11 2002 UTC revision 1.4 by dimitri, Tue Feb 18 05:33:55 2003 UTC
# Line 13  c     ================================== Line 13  c     ==================================
13  c  c
14  c     o List all the settings of the external forcing.  c     o List all the settings of the external forcing.
15  c  c
16  c     started: Christian Eckert eckert@mit.edu  11-Jan-1999  c     started: Christian Eckert eckert@mit.edu 11-Jan-1999
17  c  c
18  c     changed: Christian Eckert eckert@mit.edu  12-Feb-2000  c     changed: Christian Eckert eckert@mit.edu 12-Feb-2000
19  c  c              - changed routine names (package prefix: exf_)
 c              - Changed Routine names (package prefix: exf_)  
 c  
 c              Patrick Heimbach, heimbach@mit.edu  04-May-2000  
20  c  c
21    c     changed: Patrick Heimbach heimbach@mit.edu 04-May-2000
22  c              - changed the handling of precip and sflux with respect  c              - changed the handling of precip and sflux with respect
23  c                to CPP options ALLOW_BULKFORMULAE and ALLOW_ATM_TEMP  c                to CPP options ALLOW_BULKFORMULAE and ALLOW_ATM_TEMP
24  c  c
25  c              Dimitris Menemenlis menemenlis@jpl.nasa.gov  20-Dec-2002  c     changed: Dimitris Menemenlis menemenlis@jpl.nasa.gov 20-Dec-2002
26  c  c              - modifications for using pkg/exf with pkg/seaice
 c              - Added EXF_READ_EVAP and EXF_NO_BULK_COMPUTATIONS.  
27  c  c
28  c     ==================================================================  c     ==================================================================
29  c     SUBROUTINE exf_Summary  c     SUBROUTINE exf_Summary
# Line 96  c     == end of interface == Line 93  c     == end of interface ==
93        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
94       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
95    
96  c     For each data set used the summary prints the calendar data  c--   Print settings of some CPP flags.
 c     and the corresponding file from which the data will be read.  
97    
 #ifdef ALLOW_BULKFORMULAE  
       write(msgbuf,'(a)')  
      &'// ALLOW_BULKFORMULAE:                 defined'  
       call print_message( msgbuf, standardmessageunit,  
      &                    SQUEEZE_RIGHT , mythid)  
98  #ifdef ALLOW_ATM_TEMP  #ifdef ALLOW_ATM_TEMP
99        write(msgbuf,'(a)')        write(msgbuf,'(a)')
      &'// '  
       call print_message( msgbuf, standardmessageunit,  
      &                    SQUEEZE_RIGHT , mythid)  
       write(msgbuf,'(a)')  
100       &'// ALLOW_ATM_TEMP:                     defined'       &'// ALLOW_ATM_TEMP:                     defined'
101        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
102       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
103  #else  #else
104        write(msgbuf,'(a)')        write(msgbuf,'(a)')
      &'// '  
       call print_message( msgbuf, standardmessageunit,  
      &                    SQUEEZE_RIGHT , mythid)  
       write(msgbuf,'(a)')  
105       &'// ALLOW_ATM_TEMP:                 NOT defined'       &'// ALLOW_ATM_TEMP:                 NOT defined'
106        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
107       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
108  #endif  #endif
109    
110  #ifdef ALLOW_ATM_WIND  #ifdef ALLOW_ATM_WIND
111        write(msgbuf,'(a)')        write(msgbuf,'(a)')
112       &'// ALLOW_ATM_WIND:                     defined'       &'// ALLOW_ATM_WIND:                     defined'
# Line 134  c     and the corresponding file from wh Line 118  c     and the corresponding file from wh
118        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
119       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
120  #endif  #endif
121  #ifdef EXF_NO_BULK_COMPUTATIONS  
122    #ifdef ALLOW_DOWNWARD_RADIATION
123        write(msgbuf,'(a)')        write(msgbuf,'(a)')
124       &'// EXF_NO_BULK_COMPUTATIONS:           defined'       &'// ALLOW_DOWNWARD_RADIATION:           defined'
125        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
126       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
127  #else  #else
128        write(msgbuf,'(a)')        write(msgbuf,'(a)')
129       &'// EXF_NO_BULK_COMPUTATIONS:       NOT defined'       &'// DOWNWARD_RADIATION:             NOT defined'
       call print_message( msgbuf, standardmessageunit,  
      &                    SQUEEZE_RIGHT , mythid)  
 #endif  
 #else  
       write(msgbuf,'(a)')  
      &'// ALLOW_BULKFORMULAE:             NOT defined'  
130        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
131       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
132  #endif  #endif
133    
134  #ifdef ALLOW_KPP  #ifdef ALLOW_BULKFORMULAE
135        write(msgbuf,'(a)')        write(msgbuf,'(a)')
136       &'// ALLOW_KPP:                          defined'       &'// ALLOW_BULKFORMULAE:                 defined'
137        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
138       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
139  #else  #else
140        write(msgbuf,'(a)')        write(msgbuf,'(a)')
141       &'// ALLOW_KPP:                      NOT defined'       &'// ALLOW_BULKFORMULAE:             NOT defined'
142        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
143       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
144  #endif  #endif
145    
146  #ifdef EXF_READ_EVAP  c--   For each data set used the summary prints the calendar data
147        write(msgbuf,'(a)')  c     and the corresponding file from which the data will be read.
      &'// EXF_READ_EVAP:                      defined'  
       call print_message( msgbuf, standardmessageunit,  
      &                    SQUEEZE_RIGHT , mythid)  
       il = ilnblnk(evapfile)  
       call cal_TimeInterval( evapperiod, 'secs', timeint, mythid )  
148    
149    #ifndef ALLOW_ATM_WIND
150    c--   Zonal wind stress.
151          il = ilnblnk(ustressfile)
152          call cal_TimeInterval( ustressperiod, 'secs', timeint, mythid )
153        write(msgbuf,'(a)')        write(msgbuf,'(a)')
154       &' '       &' '
155        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
156       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
157        write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')        write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
158       &'   Evaporation starts at     ',       &'   Zonal wind stress forcing starts at         ',
159       &    (evapstartdate(i), i=1,2), dayofweek(evapstartdate(4)),       &    (ustressstartdate(i), i=1,2), dayofweek(ustressstartdate(4)),
160       &    '.'       &    '.'
161        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
162       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
163        write(msgbuf,'(a,i9.8,i7.6)')        write(msgbuf,'(a,i9.8,i7.6)')
164       &'   Evaporation period is     ',       &'   Zonal wind stress forcing period is         ',
165       &    (timeint(i), i=1,2)       &    (timeint(i), i=1,2)
166        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
167       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
168        write(msgbuf,'(a)')        write(msgbuf,'(a)')
169       &'   Evaporation is read from file:'       &'   Zonal wind stress forcing is read from file:'
170        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
171       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
172        write(msgbuf,'(a,a,a)')        write(msgbuf,'(a,a,a)')
173       &'   >>  ',evapfile(1:il),'  <<'       &'   >>  ',ustressfile(1:il),'  <<'
       call print_message( msgbuf, standardmessageunit,  
      &                    SQUEEZE_RIGHT , mythid)  
 #else EXF_READ_EVAP  
       write(msgbuf,'(a)')  
      &'// EXF_READ_EVAP:                  NOT defined'  
       call print_message( msgbuf, standardmessageunit,  
      &                    SQUEEZE_RIGHT , mythid)  
 #endif EXF_READ_EVAP  
   
 #ifdef ALLOW_RUNOFF  
       write(msgbuf,'(a)')  
      &'// ALLOW_RUNOFF:                       defined'  
174        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
175       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
       il = ilnblnk(runofffile)  
       call cal_TimeInterval( runoffperiod, 'secs', timeint, mythid )  
176    
177    c--   Meridional wind stress.
178          il = ilnblnk(vstressfile)
179          call cal_TimeInterval( vstressperiod, 'secs', timeint, mythid )
180        write(msgbuf,'(a)')        write(msgbuf,'(a)')
181       &' '       &' '
182        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
183       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
184        write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')        write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
185       &'   Runnoff starts at     ',       &'   Meridional wind stress forcing starts at    ',
186       &    (runoffstartdate(i), i=1,2), dayofweek(runoffstartdate(4)),       &    (vstressstartdate(i), i=1,2), dayofweek(vstressstartdate(4)),
187       &    '.'       &    '.'
188        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
189       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
190        write(msgbuf,'(a,i9.8,i7.6)')        write(msgbuf,'(a,i9.8,i7.6)')
191       &'   Runoff period is     ',       &'   Meridional wind stress forcing period is    ',
192       &    (timeint(i), i=1,2)       &    (timeint(i), i=1,2)
193        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
194       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
195        write(msgbuf,'(a)')        write(msgbuf,'(a)')
196       &'   Runoff is read from file:'       &'   Meridional wind stress forcing is read from file:'
197        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
198       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
199        write(msgbuf,'(a,a,a)')        write(msgbuf,'(a,a,a)')
200       &'   >>  ',runofffile(1:il),'  <<'       &'   >>  ',vstressfile(1:il),'  <<'
       call print_message( msgbuf, standardmessageunit,  
      &                    SQUEEZE_RIGHT , mythid)  
 #else ALLOW_RUNOFF  
       write(msgbuf,'(a)')  
      &'// ALLOW_RUNOFF:                   NOT defined'  
201        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
202       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
203  #endif ALLOW_RUNOFF  #endif
   
 #ifdef ALLOW_BULKFORMULAE  
   
 #ifdef ALLOW_ATM_TEMP  
 c     Atmospheric temperature.  
       il = ilnblnk(atempfile)  
       call cal_TimeInterval( atempperiod, 'secs', timeint, mythid )  
204    
205    #ifndef ALLOW_ATM_TEMP
206    c--   Heat flux.
207          il = ilnblnk(hfluxfile)
208          call cal_TimeInterval( hfluxperiod, 'secs', timeint, mythid )
209        write(msgbuf,'(a)')        write(msgbuf,'(a)')
210       &' '       &' '
211        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
212       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
213        write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')        write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
214       &'   Atmospheric temperature starts at           ',       &'   Heat flux forcing starts at                ',
215       &    (atempstartdate(i), i=1,2), dayofweek(atempstartdate(4)),       &    (hfluxstartdate(i), i=1,2), dayofweek(hfluxstartdate(4)),'.'
      &    '.'  
216        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
217       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
218        write(msgbuf,'(a,i9.8,i7.6)')        write(msgbuf,'(a,i9.8,i7.6)')
219       &'   Atmospheric temperature period is           ',       &'   Heat flux forcing period is                 ',
220       &    (timeint(i), i=1,2)       &    (timeint(i), i=1,2)
221        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
222       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
223        write(msgbuf,'(a)')        write(msgbuf,'(a)')
224       &'   Atmospheric temperature is read from file:'       &'   Heat flux forcing is read from file:        '
225        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
226       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
227        write(msgbuf,'(a,a,a)')        write(msgbuf,'(a,a,a)')
228       &'   >>  ',atempfile(1:il),'  <<'       &'   >>  ',hfluxfile(1:il),'  <<'
229        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
230       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
231    
232  c     Atmospheric specific humidity.  c--   Salt flux.
233        il = ilnblnk(aqhfile)        il = ilnblnk(sfluxfile)
234        call cal_TimeInterval( aqhperiod, 'secs', timeint, mythid )        call cal_TimeInterval( sfluxperiod, 'secs', timeint, mythid )
   
235        write(msgbuf,'(a)')        write(msgbuf,'(a)')
236       &' '       &' '
237        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
238       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
239        write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')        write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
240       &'   Atmospheric specific humidity starts at     ',       &'   Salt flux forcing starts at                 ',
241       &    (aqhstartdate(i), i=1,2), dayofweek(aqhstartdate(4)),       &    (sfluxstartdate(i), i=1,2), dayofweek(sfluxstartdate(4)),'.'
      &    '.'  
242        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
243       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
244        write(msgbuf,'(a,i9.8,i7.6)')        write(msgbuf,'(a,i9.8,i7.6)')
245       &'   Atmospheric specific humidity period is     ',       &'   Salt flux forcing period is                 ',
246       &    (timeint(i), i=1,2)       &    (timeint(i), i=1,2)
247        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
248       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
249        write(msgbuf,'(a)')        write(msgbuf,'(a)')
250       &'   Atmospheric specific humidity is read from file:'       &'   Salt flux forcing is read from file:        '
251        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
252       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
253        write(msgbuf,'(a,a,a)')        write(msgbuf,'(a,a,a)')
254       &'   >>  ',aqhfile(1:il),'  <<'       &'   >>  ',sfluxfile(1:il),'  <<'
255        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
256       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
257    #endif
258    
259  c     Short wave radiative flux.  #if defined(ALLOW_ATM_TEMP) || defined(SHORTWAVE_HEATING)
260    c--   Net shortwave.
261        il = ilnblnk(swfluxfile)        il = ilnblnk(swfluxfile)
262        call cal_TimeInterval( swfluxperiod, 'secs', timeint, mythid )        call cal_TimeInterval( swfluxperiod, 'secs', timeint, mythid )
   
263        write(msgbuf,'(a)')        write(msgbuf,'(a)')
264       &' '       &' '
265        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
266       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
267        write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')        write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
268       &'   Short wave rad. flux forcing starts at      ',       &'   Net shortwave flux forcing starts at      ',
269       &    (swfluxstartdate(i), i=1,2), dayofweek(swfluxstartdate(4)),       &    (swfluxstartdate(i), i=1,2), dayofweek(swfluxstartdate(4)),
270       &    '.'       &    '.'
271        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
272       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
273        write(msgbuf,'(a,i9.8,i7.6)')        write(msgbuf,'(a,i9.8,i7.6)')
274       &'   Short wave rad. flux forcing period is      ',       &'   Net shortwave flux forcing period is      ',
275       &    (timeint(i), i=1,2)       &    (timeint(i), i=1,2)
276        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
277       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
278        write(msgbuf,'(a)')        write(msgbuf,'(a)')
279       &'   Short wave rad. flux forcing is read from file:'       &'   Net shortwave flux forcing is read from file:'
280        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
281       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
282        write(msgbuf,'(a,a,a)')        write(msgbuf,'(a,a,a)')
283       &'   >>  ',swfluxfile(1:il),'  <<'       &'   >>  ',swfluxfile(1:il),'  <<'
284        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
285       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
286    #endif
287    
288  c     Long wave radiative flux.  #ifdef ALLOW_ATM_WIND
289        il = ilnblnk(lwfluxfile)  c--   Zonal wind.
290        call cal_TimeInterval( lwfluxperiod, 'secs', timeint, mythid )        il = ilnblnk(uwindfile)
291          call cal_TimeInterval( uwindperiod, 'secs', timeint, mythid )
292        write(msgbuf,'(a)')        write(msgbuf,'(a)')
293       &' '       &' '
294        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
295       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
296        write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')        write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
297       &'   Long wave rad. flux forcing starts at       ',       &'   Zonal wind forcing starts at                ',
298       &    (lwfluxstartdate(i), i=1,2), dayofweek(lwfluxstartdate(4)),       &    (uwindstartdate(i), i=1,2), dayofweek(uwindstartdate(4)),
299       &    '.'       &    '.'
300        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
301       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
302        write(msgbuf,'(a,i9.8,i7.6)')        write(msgbuf,'(a,i9.8,i7.6)')
303       &'   Long wave rad. flux forcing period is       ',       &'   Zonal wind forcing period is                ',
304       &    (timeint(i), i=1,2)       &    (timeint(i), i=1,2)
305        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
306       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
307        write(msgbuf,'(a)')        write(msgbuf,'(a)')
308       &'   Long wave rad. flux forcing is read from file:'       &'   Zonal wind forcing is read from file:'
309        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
310       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
311        write(msgbuf,'(a,a,a)')        write(msgbuf,'(a,a,a)')
312       &'   >>  ',lwfluxfile(1:il),'  <<'       &'   >>  ',uwindfile(1:il),'  <<'
313        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
314       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
315    
316  c     Precipitation.  c--   Meridional wind.
317        il = ilnblnk(precipfile)        il = ilnblnk(vwindfile)
318        call cal_TimeInterval( precipperiod, 'secs', timeint, mythid )        call cal_TimeInterval( vwindperiod, 'secs', timeint, mythid )
   
319        write(msgbuf,'(a)')        write(msgbuf,'(a)')
320       &' '       &' '
321        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
322       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
323        write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')        write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
324       &'   Precipitation data set starts at            ',       &'   Meridional wind forcing starts at           ',
325       &    (precipstartdate(i), i=1,2), dayofweek(precipstartdate(4)),       &    (vwindstartdate(i), i=1,2), dayofweek(vwindstartdate(4)),
326       &'.'       &    '.'
327        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
328       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
329        write(msgbuf,'(a,i9.8,i7.6)')        write(msgbuf,'(a,i9.8,i7.6)')
330       &'   Precipitation data period is                ',       &'   Meridional wind forcing period is           ',
331       &    (timeint(i), i=1,2)       &    (timeint(i), i=1,2)
332        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
333       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
334        write(msgbuf,'(a)')        write(msgbuf,'(a)')
335       &'   Precipitation data is read from file:       '       &'   Meridional wind forcing is read from file:'
336        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
337       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
338        write(msgbuf,'(a,a,a)')        write(msgbuf,'(a,a,a)')
339       &'   >>  ',precipfile(1:il),'  <<'       &'   >>  ',vwindfile(1:il),'  <<'
340        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
341       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
342    #endif
343    
344  #else  #ifdef ALLOW_ATM_TEMP
345  c     Heat flux.  c--   Atmospheric temperature.
346        il = ilnblnk(hfluxfile)        il = ilnblnk(atempfile)
347        call cal_TimeInterval( hfluxperiod, 'secs', timeint, mythid )        call cal_TimeInterval( atempperiod, 'secs', timeint, mythid )
   
348        write(msgbuf,'(a)')        write(msgbuf,'(a)')
349       &' '       &' '
350        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
351       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
352        write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')        write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
353       &'   Heat flux forcing starts at                ',       &'   Atmospheric temperature starts at           ',
354       &    (hfluxstartdate(i), i=1,2), dayofweek(hfluxstartdate(4)),'.'       &    (atempstartdate(i), i=1,2), dayofweek(atempstartdate(4)),
355         &    '.'
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,i9.8,i7.6)')
359       &'   Heat flux forcing period is                 ',       &'   Atmospheric temperature period is           ',
360       &    (timeint(i), i=1,2)       &    (timeint(i), i=1,2)
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)')
364       &'   Heat flux forcing is read from file:        '       &'   Atmospheric temperature is read from file:'
365        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
366       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
367        write(msgbuf,'(a,a,a)')        write(msgbuf,'(a,a,a)')
368       &'   >>  ',hfluxfile(1:il),'  <<'       &'   >>  ',atempfile(1:il),'  <<'
369        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
370       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
371    
372  c     Salt flux.  c--   Atmospheric specific humidity.
373        il = ilnblnk(sfluxfile)        il = ilnblnk(aqhfile)
374        call cal_TimeInterval( sfluxperiod, 'secs', timeint, mythid )        call cal_TimeInterval( aqhperiod, 'secs', timeint, mythid )
   
375        write(msgbuf,'(a)')        write(msgbuf,'(a)')
376       &' '       &' '
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,i9.8,i7.6,1x,a,a)')
380       &'   Salt flux forcing starts at                 ',       &'   Atmospheric specific humidity starts at     ',
381       &    (sfluxstartdate(i), i=1,2), dayofweek(sfluxstartdate(4)),'.'       &    (aqhstartdate(i), i=1,2), dayofweek(aqhstartdate(4)),
382         &    '.'
383        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
384       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
385        write(msgbuf,'(a,i9.8,i7.6)')        write(msgbuf,'(a,i9.8,i7.6)')
386       &'   Salt flux forcing period is                 ',       &'   Atmospheric specific humidity period is     ',
387       &    (timeint(i), i=1,2)       &    (timeint(i), i=1,2)
388        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
389       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
390        write(msgbuf,'(a)')        write(msgbuf,'(a)')
391       &'   Salt flux forcing is read from file:        '       &'   Atmospheric specific humidity is read from file:'
392        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
393       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
394        write(msgbuf,'(a,a,a)')        write(msgbuf,'(a,a,a)')
395       &'   >>  ',sfluxfile(1:il),'  <<'       &'   >>  ',aqhfile(1:il),'  <<'
396        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
397       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
398    
399  #ifdef ALLOW_KPP  c--   Net longwave.
400  c     Short wave radiative flux.        il = ilnblnk(lwfluxfile)
401        il = ilnblnk(swfluxfile)        call cal_TimeInterval( lwfluxperiod, 'secs', timeint, mythid )
       call cal_TimeInterval( swfluxperiod, 'secs', timeint, mythid )  
   
402        write(msgbuf,'(a)')        write(msgbuf,'(a)')
403       &' '       &' '
404        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
405       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
406        write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')        write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
407       &'   Short wave rad. flux forcing starts at      ',       &'   Net longwave flux forcing starts at       ',
408       &    (swfluxstartdate(i), i=1,2), dayofweek(swfluxstartdate(4)),       &    (lwfluxstartdate(i), i=1,2), dayofweek(lwfluxstartdate(4)),
409       &    '.'       &    '.'
410        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
411       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
412        write(msgbuf,'(a,i9.8,i7.6)')        write(msgbuf,'(a,i9.8,i7.6)')
413       &'   Short wave rad. flux forcing period is      ',       &'   Net longwave flux forcing period is       ',
414       &    (timeint(i), i=1,2)       &    (timeint(i), i=1,2)
415        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
416       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
417        write(msgbuf,'(a)')        write(msgbuf,'(a)')
418       &'   Short wave rad. flux forcing is read from file:'       &'   Net longwave flux forcing is read from file:'
419        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
420       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
421        write(msgbuf,'(a,a,a)')        write(msgbuf,'(a,a,a)')
422       &'   >>  ',swfluxfile(1:il),'  <<'       &'   >>  ',lwfluxfile(1:il),'  <<'
423        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
424       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
 #endif  
 #endif  
   
 #ifdef ALLOW_ATM_WIND  
 c     Zonal wind.  
       il = ilnblnk(uwindfile)  
       call cal_TimeInterval( uwindperiod, 'secs', timeint, mythid )  
425    
426    c--   Precipitation.
427          il = ilnblnk(precipfile)
428          call cal_TimeInterval( precipperiod, 'secs', timeint, mythid )
429        write(msgbuf,'(a)')        write(msgbuf,'(a)')
430       &' '       &' '
431        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
432       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
433        write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')        write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
434       &'   Zonal wind forcing starts at                ',       &'   Precipitation data set starts at            ',
435       &    (uwindstartdate(i), i=1,2), dayofweek(uwindstartdate(4)),       &    (precipstartdate(i), i=1,2), dayofweek(precipstartdate(4)),
436       &    '.'       &'.'
437        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
438       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
439        write(msgbuf,'(a,i9.8,i7.6)')        write(msgbuf,'(a,i9.8,i7.6)')
440       &'   Zonal wind forcing period is                ',       &'   Precipitation data period is                ',
441       &    (timeint(i), i=1,2)       &    (timeint(i), i=1,2)
442        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
443       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
444        write(msgbuf,'(a)')        write(msgbuf,'(a)')
445       &'   Zonal wind forcing is read from file:'       &'   Precipitation data is read from file:       '
446        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
447       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
448        write(msgbuf,'(a,a,a)')        write(msgbuf,'(a,a,a)')
449       &'   >>  ',uwindfile(1:il),'  <<'       &'   >>  ',precipfile(1:il),'  <<'
450        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
451       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
452    #endif
453    
454  c     Meridional wind.  #ifdef EXF_READ_EVAP
455        il = ilnblnk(vwindfile)  c--   Evaporation.
       call cal_TimeInterval( vwindperiod, 'secs', timeint, mythid )  
   
       write(msgbuf,'(a)')  
      &' '  
       call print_message( msgbuf, standardmessageunit,  
      &                    SQUEEZE_RIGHT , mythid)  
       write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')  
      &'   Meridional wind forcing starts at           ',  
      &    (vwindstartdate(i), i=1,2), dayofweek(vwindstartdate(4)),  
      &    '.'  
       call print_message( msgbuf, standardmessageunit,  
      &                    SQUEEZE_RIGHT , mythid)  
       write(msgbuf,'(a,i9.8,i7.6)')  
      &'   Meridional wind forcing period is           ',  
      &    (timeint(i), i=1,2)  
       call print_message( msgbuf, standardmessageunit,  
      &                    SQUEEZE_RIGHT , mythid)  
456        write(msgbuf,'(a)')        write(msgbuf,'(a)')
457       &'   Meridional wind forcing is read from file:'       &'// EXF_READ_EVAP:                      defined'
       call print_message( msgbuf, standardmessageunit,  
      &                    SQUEEZE_RIGHT , mythid)  
       write(msgbuf,'(a,a,a)')  
      &'   >>  ',vwindfile(1:il),'  <<'  
458        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
459       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
460  #else        il = ilnblnk(evapfile)
461  c     Zonal wind stress.        call cal_TimeInterval( evapperiod, 'secs', timeint, mythid )
       il = ilnblnk(ustressfile)  
       call cal_TimeInterval( ustressperiod, 'secs', timeint, mythid )  
462    
463        write(msgbuf,'(a)')        write(msgbuf,'(a)')
464       &' '       &' '
465        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
466       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
467        write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')        write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
468       &'   Zonal wind stress forcing starts at         ',       &'   Evaporation starts at     ',
469       &    (ustressstartdate(i), i=1,2), dayofweek(ustressstartdate(4)),       &    (evapstartdate(i), i=1,2), dayofweek(evapstartdate(4)),
470       &    '.'       &    '.'
471        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
472       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
473        write(msgbuf,'(a,i9.8,i7.6)')        write(msgbuf,'(a,i9.8,i7.6)')
474       &'   Zonal wind stress forcing period is         ',       &'   Evaporation period is     ',
475       &    (timeint(i), i=1,2)       &    (timeint(i), i=1,2)
476        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
477       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
478        write(msgbuf,'(a)')        write(msgbuf,'(a)')
479       &'   Zonal wind stress forcing is read from file:'       &'   Evaporation is read from file:'
480        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
481       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
482        write(msgbuf,'(a,a,a)')        write(msgbuf,'(a,a,a)')
483       &'   >>  ',ustressfile(1:il),'  <<'       &'   >>  ',evapfile(1:il),'  <<'
484        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
485       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
486    #else
 c     Meridional wind stress.  
       il = ilnblnk(vstressfile)  
       call cal_TimeInterval( vstressperiod, 'secs', timeint, mythid )  
   
487        write(msgbuf,'(a)')        write(msgbuf,'(a)')
488       &' '       &'// EXF_READ_EVAP:                  NOT defined'
       call print_message( msgbuf, standardmessageunit,  
      &                    SQUEEZE_RIGHT , mythid)  
       write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')  
      &'   Meridional wind stress forcing starts at    ',  
      &    (vstressstartdate(i), i=1,2), dayofweek(vstressstartdate(4)),  
      &    '.'  
       call print_message( msgbuf, standardmessageunit,  
      &                    SQUEEZE_RIGHT , mythid)  
       write(msgbuf,'(a,i9.8,i7.6)')  
      &'   Meridional wind stress forcing period is    ',  
      &    (timeint(i), i=1,2)  
489        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
490       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
491    #endif
492    
493    #ifdef ALLOW_RUNOFF
494    c--   Runoff.
495        write(msgbuf,'(a)')        write(msgbuf,'(a)')
496       &'   Meridional wind stress forcing is read from file:'       &'// ALLOW_RUNOFF:                       defined'
       call print_message( msgbuf, standardmessageunit,  
      &                    SQUEEZE_RIGHT , mythid)  
       write(msgbuf,'(a,a,a)')  
      &'   >>  ',vstressfile(1:il),'  <<'  
497        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
498       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
499  #endif        il = ilnblnk(runofffile)
500          call cal_TimeInterval( runoffperiod, 'secs', timeint, mythid )
 #else  
 c     Heat flux.  
       il = ilnblnk(hfluxfile)  
       call cal_TimeInterval( hfluxperiod, 'secs', timeint, mythid )  
501    
502        write(msgbuf,'(a)')        write(msgbuf,'(a)')
503       &' '       &' '
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,1x,a,a)')        write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
507       &'   Heat flux forcing starts at                 ',       &'   Runnoff starts at     ',
508       &    (hfluxstartdate(i), i=1,2), dayofweek(hfluxstartdate(4)),'.'       &    (runoffstartdate(i), i=1,2), dayofweek(runoffstartdate(4)),
509         &    '.'
510        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
511       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
512        write(msgbuf,'(a,i9.8,i7.6)')        write(msgbuf,'(a,i9.8,i7.6)')
513       &'   Heat flux forcing period is                 ',       &'   Runoff period is     ',
514       &    (timeint(i), i=1,2)       &    (timeint(i), i=1,2)
515        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
516       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
517        write(msgbuf,'(a)')        write(msgbuf,'(a)')
518       &'   Heat flux forcing is read from file:        '       &'   Runoff is read from file:'
519        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
520       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
521        write(msgbuf,'(a,a,a)')        write(msgbuf,'(a,a,a)')
522       &'   >>  ',hfluxfile(1:il),'  <<'       &'   >>  ',runofffile(1:il),'  <<'
       call print_message( msgbuf, standardmessageunit,  
      &                    SQUEEZE_RIGHT , mythid)  
   
 c     Salt flux.  
       il = ilnblnk(sfluxfile)  
       call cal_TimeInterval( sfluxperiod, 'secs', timeint, mythid )  
   
       write(msgbuf,'(a)')  
      &' '  
       call print_message( msgbuf, standardmessageunit,  
      &                    SQUEEZE_RIGHT , mythid)  
       write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')  
      &'   Salt flux forcing starts at                 ',  
      &    (sfluxstartdate(i), i=1,2), dayofweek(sfluxstartdate(4)),'.'  
       call print_message( msgbuf, standardmessageunit,  
      &                    SQUEEZE_RIGHT , mythid)  
       write(msgbuf,'(a,i9.8,i7.6)')  
      &'   Salt flux forcing period is                 ',  
      &    (timeint(i), i=1,2)  
523        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
524       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
525    #else /* ALLOW_RUNOFF */
526        write(msgbuf,'(a)')        write(msgbuf,'(a)')
527       &'   Salt flux forcing is read from file:        '       &'// ALLOW_RUNOFF:                   NOT defined'
       call print_message( msgbuf, standardmessageunit,  
      &                    SQUEEZE_RIGHT , mythid)  
       write(msgbuf,'(a,a,a)')  
      &'   >>  ',sfluxfile(1:il),'  <<'  
528        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
529       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
530    #endif /* ALLOW_RUNOFF */
531    
532  c     Zonal wind stress.  #ifdef DOWNWARD_RADIATION
533        il = ilnblnk(ustressfile)  c--   Downward shortwave.
534        call cal_TimeInterval( ustressperiod, 'secs', timeint, mythid )        il = ilnblnk(swdownfile)
535          call cal_TimeInterval( swdownperiod, 'secs', timeint, mythid )
536        write(msgbuf,'(a)')        write(msgbuf,'(a)')
537       &' '       &' '
538        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
539       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
540        write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')        write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
541       &'   Zonal wind stress forcing starts at         ',       &'   Downward shortwave flux forcing starts at      ',
542       &    (ustressstartdate(i), i=1,2), dayofweek(ustressstartdate(4)),       &    (swdownstartdate(i), i=1,2), dayofweek(swdownstartdate(4)),
543       &    '.'       &    '.'
544        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
545       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
546        write(msgbuf,'(a,i9.8,i7.6)')        write(msgbuf,'(a,i9.8,i7.6)')
547       &'   Zonal wind stress forcing period is         ',       &'   Downward shortwave flux forcing period is      ',
548       &    (timeint(i), i=1,2)       &    (timeint(i), i=1,2)
549        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
550       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
551        write(msgbuf,'(a)')        write(msgbuf,'(a)')
552       &'   Zonal wind stress forcing is read from file:'       &'   Downward shortwave flux forcing is read from file:'
553        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
554       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
555        write(msgbuf,'(a,a,a)')        write(msgbuf,'(a,a,a)')
556       &'   >>  ',ustressfile(1:il),'  <<'       &'   >>  ',swdownfile(1:il),'  <<'
557        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
558       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
559    
560  c     Meridional wind stress.  c--   Downward longwave.
561        il = ilnblnk(vstressfile)        il = ilnblnk(lwdownfile)
562        call cal_TimeInterval( vstressperiod, 'secs', timeint, mythid )        call cal_TimeInterval( lwdownperiod, 'secs', timeint, mythid )
   
563        write(msgbuf,'(a)')        write(msgbuf,'(a)')
564       &' '       &' '
565        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
566       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
567        write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')        write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
568       &'   Meridional wind stress forcing starts at    ',       &'   Downward longwave flux forcing starts at       ',
569       &    (vstressstartdate(i), i=1,2), dayofweek(vstressstartdate(4)),       &    (lwdownstartdate(i), i=1,2), dayofweek(lwdownstartdate(4)),
570       &    '.'       &    '.'
571        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
572       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
573        write(msgbuf,'(a,i9.8,i7.6)')        write(msgbuf,'(a,i9.8,i7.6)')
574       &'   Meridional wind stress forcing period is    ',       &'   Downward longwave flux forcing period is       ',
575       &    (timeint(i), i=1,2)       &    (timeint(i), i=1,2)
576        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
577       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
578        write(msgbuf,'(a)')        write(msgbuf,'(a)')
579       &'   Meridional wind stress forcing is read from file:'       &'   Downward longwave flux forcing is read from file:'
580        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
581       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
582        write(msgbuf,'(a,a,a)')        write(msgbuf,'(a,a,a)')
583       &'   >>  ',vstressfile(1:il),'  <<'       &'   >>  ',lwdownfile(1:il),'  <<'
584        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
585       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
586    #endif
587    
588  #ifdef ALLOW_KPP  #ifdef ATMOSPHERIC_LOADING
589  c     Short wave radiative flux.  c--   Atmospheric pressure.
590        il = ilnblnk(swfluxfile)        il = ilnblnk(apressurefile)
591        call cal_TimeInterval( swfluxperiod, 'secs', timeint, mythid )        call cal_TimeInterval( apressureperiod, 'secs', timeint, mythid )
   
592        write(msgbuf,'(a)')        write(msgbuf,'(a)')
593       &' '       &' '
594        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
595       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
596        write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')        write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
597       &'   Short wave rad. flux forcing starts at      ',       &'   Atmospheric pressure forcing starts at      ',
598       &    (swfluxstartdate(i), i=1,2), dayofweek(swfluxstartdate(4)),       &    (apressurestartdate(i), i=1,2),
599       &    '.'       &    dayofweek(apressurestartdate(4)), '.'
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,i9.8,i7.6)')
603       &'   Short wave rad. flux forcing period is      ',       &'   Atmospheric pressure forcing period is      ',
604       &    (timeint(i), i=1,2)       &    (timeint(i), i=1,2)
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)')
608       &'   Short wave rad. flux forcing is read from file:'       &'   Atmospheric pressureforcing is read from file:'
609        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
610       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
611        write(msgbuf,'(a,a,a)')        write(msgbuf,'(a,a,a)')
612       &'   >>  ',swfluxfile(1:il),'  <<'       &'   >>  ',apressurefile(1:il),'  <<'
613        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
614       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
615  #endif  #endif
616    
 #endif  
   
617        write(msgbuf,'(a)')        write(msgbuf,'(a)')
618       &' '       &' '
619        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22