/[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.1 by heimbach, Mon May 14 22:08:41 2001 UTC revision 1.6 by edhill, Thu Oct 9 04:19:19 2003 UTC
# Line 1  Line 1 
1  c $Header$  c $Header$
2    
3  #include "EXF_CPPOPTIONS.h"  #include "EXF_OPTIONS.h"
4    
5    
6        subroutine exf_Summary(        subroutine exf_Summary(
# 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     changed: Dimitris Menemenlis menemenlis@jpl.nasa.gov 20-Dec-2002
26    c              - modifications for using pkg/exf with pkg/seaice
27    c
28  c     ==================================================================  c     ==================================================================
29  c     SUBROUTINE exf_Summary  c     SUBROUTINE exf_Summary
30  c     ==================================================================  c     ==================================================================
# Line 92  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)')
100       &'// '       &'// ALLOW_ATM_TEMP:                     defined'
       call print_message( msgbuf, standardmessageunit,  
      &                    SQUEEZE_RIGHT , mythid)  
       write(msgbuf,'(a)')  
      &'// 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)')
105       &'// '       &'// ALLOW_ATM_TEMP:                 NOT defined'
       call print_message( msgbuf, standardmessageunit,  
      &                    SQUEEZE_RIGHT , mythid)  
       write(msgbuf,'(a)')  
      &'// 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'
113        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
114       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
115  #else  #else
116        write(msgbuf,'(a)')        write(msgbuf,'(a)')
117       &'// ALLOW_ATM_WIND:                   NOT defined'       &'// ALLOW_ATM_WIND:                 NOT defined'
       call print_message( msgbuf, standardmessageunit,  
      &                    SQUEEZE_RIGHT , mythid)  
 #endif  
 #else  
       write(msgbuf,'(a)')  
      &'// ALLOW_BULKFORMULAE:             NOT defined'  
118        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
119       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
120  #endif  #endif
121    
122  #ifdef ALLOW_KPP  #ifdef ALLOW_DOWNWARD_RADIATION
123        write(msgbuf,'(a)')        write(msgbuf,'(a)')
124       &'// ALLOW_KPP:                          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       &'// ALLOW_KPP:                      NOT defined'       &'// ALLOW_DOWNWARD_RADIATION:       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_BULKFORMULAE  #ifdef ALLOW_BULKFORMULAE
   
 #ifdef ALLOW_ATM_TEMP  
 c     Atmospheric temperature.  
       il = ilnblnk(atempfile)  
       call cal_TimeInterval( atempperiod, 'secs', timeint, mythid )  
   
       write(msgbuf,'(a)')  
      &' '  
       call print_message( msgbuf, standardmessageunit,  
      &                    SQUEEZE_RIGHT , mythid)  
       write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')  
      &'   Atmospheric temperature starts at           ',  
      &    (atempstartdate(i), i=1,2), dayofweek(atempstartdate(4)),  
      &    '.'  
       call print_message( msgbuf, standardmessageunit,  
      &                    SQUEEZE_RIGHT , mythid)  
       write(msgbuf,'(a,i9.8,i7.6)')  
      &'   Atmospheric temperature period is           ',  
      &    (timeint(i), i=1,2)  
       call print_message( msgbuf, standardmessageunit,  
      &                    SQUEEZE_RIGHT , mythid)  
135        write(msgbuf,'(a)')        write(msgbuf,'(a)')
136       &'   Atmospheric temperature is read from file:'       &'// ALLOW_BULKFORMULAE:                 defined'
       call print_message( msgbuf, standardmessageunit,  
      &                    SQUEEZE_RIGHT , mythid)  
       write(msgbuf,'(a,a,a)')  
      &'   >>  ',atempfile(1:il),'  <<'  
       call print_message( msgbuf, standardmessageunit,  
      &                    SQUEEZE_RIGHT , mythid)  
   
 c     Atmospheric specific humidity.  
       il = ilnblnk(aqhfile)  
       call cal_TimeInterval( aqhperiod, 'secs', timeint, mythid )  
   
       write(msgbuf,'(a)')  
      &' '  
       call print_message( msgbuf, standardmessageunit,  
      &                    SQUEEZE_RIGHT , mythid)  
       write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')  
      &'   Atmospheric specific humidity starts at     ',  
      &    (aqhstartdate(i), i=1,2), dayofweek(aqhstartdate(4)),  
      &    '.'  
       call print_message( msgbuf, standardmessageunit,  
      &                    SQUEEZE_RIGHT , mythid)  
       write(msgbuf,'(a,i9.8,i7.6)')  
      &'   Atmospheric specific humidity period is     ',  
      &    (timeint(i), i=1,2)  
137        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
138       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
139    #else
140        write(msgbuf,'(a)')        write(msgbuf,'(a)')
141       &'   Atmospheric specific humidity is read from file:'       &'// ALLOW_BULKFORMULAE:             NOT defined'
       call print_message( msgbuf, standardmessageunit,  
      &                    SQUEEZE_RIGHT , mythid)  
       write(msgbuf,'(a,a,a)')  
      &'   >>  ',aqhfile(1:il),'  <<'  
142        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
143       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
144    #endif
145    
146  c     Short wave radiative flux.  c--   For each data set used the summary prints the calendar data
147        il = ilnblnk(swfluxfile)  c     and the corresponding file from which the data will be read.
       call cal_TimeInterval( swfluxperiod, '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       &'   Short wave rad. flux forcing starts at      ',       &'   Zonal wind stress forcing starts at         ',
159       &    (swfluxstartdate(i), i=1,2), dayofweek(swfluxstartdate(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       &'   Short wave rad. flux forcing 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       &'   Short wave rad. flux forcing 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       &'   >>  ',swfluxfile(1:il),'  <<'       &'   >>  ',ustressfile(1:il),'  <<'
174        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
175       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
176    
177  c     Long wave radiative flux.  c--   Meridional wind stress.
178        il = ilnblnk(lwfluxfile)        il = ilnblnk(vstressfile)
179        call cal_TimeInterval( lwfluxperiod, 'secs', timeint, mythid )        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       &'   Long wave rad. flux forcing starts at       ',       &'   Meridional wind stress forcing starts at    ',
186       &    (lwfluxstartdate(i), i=1,2), dayofweek(lwfluxstartdate(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       &'   Long wave rad. flux forcing period is       ',       &'   Meridional wind stress forcing period is    ',
      &    (timeint(i), i=1,2)  
       call print_message( msgbuf, standardmessageunit,  
      &                    SQUEEZE_RIGHT , mythid)  
       write(msgbuf,'(a)')  
      &'   Long wave rad. flux forcing is read from file:'  
       call print_message( msgbuf, standardmessageunit,  
      &                    SQUEEZE_RIGHT , mythid)  
       write(msgbuf,'(a,a,a)')  
      &'   >>  ',lwfluxfile(1:il),'  <<'  
       call print_message( msgbuf, standardmessageunit,  
      &                    SQUEEZE_RIGHT , mythid)  
   
 c     Precipitation.  
       il = ilnblnk(precipfile)  
       call cal_TimeInterval( precipperiod, 'secs', timeint, mythid )  
   
       write(msgbuf,'(a)')  
      &' '  
       call print_message( msgbuf, standardmessageunit,  
      &                    SQUEEZE_RIGHT , mythid)  
       write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')  
      &'   Precipitation data set starts at            ',  
      &    (precipstartdate(i), i=1,2), dayofweek(precipstartdate(4)),  
      &'.'  
       call print_message( msgbuf, standardmessageunit,  
      &                    SQUEEZE_RIGHT , mythid)  
       write(msgbuf,'(a,i9.8,i7.6)')  
      &'   Precipitation data 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       &'   Precipitation data 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       &'   >>  ',precipfile(1:il),'  <<'       &'   >>  ',vstressfile(1:il),'  <<'
201        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
202       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
203    #endif
204    
205  #else  #ifndef ALLOW_ATM_TEMP
206  c     Heat flux.  c--   Heat flux.
207        il = ilnblnk(hfluxfile)        il = ilnblnk(hfluxfile)
208        call cal_TimeInterval( hfluxperiod, 'secs', timeint, mythid )        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,
# Line 320  c     Heat flux. Line 229  c     Heat flux.
229        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
230       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
231    
232  c     Salt flux.  c--   Salt flux.
233        il = ilnblnk(sfluxfile)        il = ilnblnk(sfluxfile)
234        call cal_TimeInterval( sfluxperiod, '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,
# Line 346  c     Salt flux. Line 254  c     Salt flux.
254       &'   >>  ',sfluxfile(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  #ifdef ALLOW_KPP  #if defined(ALLOW_ATM_TEMP) || defined(SHORTWAVE_HEATING)
260  c     Short wave radiative flux.  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)')
# Line 376  c     Short wave radiative flux. Line 284  c     Short wave radiative flux.
284        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
285       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
286  #endif  #endif
 #endif  
287    
288  #ifdef ALLOW_ATM_WIND  #ifdef ALLOW_ATM_WIND
289  c     Zonal wind.  c--   Zonal wind.
290        il = ilnblnk(uwindfile)        il = ilnblnk(uwindfile)
291        call cal_TimeInterval( uwindperiod, 'secs', timeint, mythid )        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,
# Line 407  c     Zonal wind. Line 313  c     Zonal wind.
313        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
314       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
315    
316  c     Meridional wind.  c--   Meridional wind.
317        il = ilnblnk(vwindfile)        il = ilnblnk(vwindfile)
318        call cal_TimeInterval( vwindperiod, '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,
# Line 434  c     Meridional wind. Line 339  c     Meridional wind.
339       &'   >>  ',vwindfile(1:il),'  <<'       &'   >>  ',vwindfile(1:il),'  <<'
340        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
341       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
342  #else  #endif
 c     Zonal wind stress.  
       il = ilnblnk(ustressfile)  
       call cal_TimeInterval( ustressperiod, 'secs', timeint, mythid )  
343    
344    #ifdef ALLOW_ATM_TEMP
345    c--   Atmospheric temperature.
346          il = ilnblnk(atempfile)
347          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       &'   Zonal wind stress forcing starts at         ',       &'   Atmospheric temperature starts at           ',
354       &    (ustressstartdate(i), i=1,2), dayofweek(ustressstartdate(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       &'   Zonal wind stress 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       &'   Zonal wind stress 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       &'   >>  ',ustressfile(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     Meridional wind stress.  c--   Atmospheric specific humidity.
373        il = ilnblnk(vstressfile)        il = ilnblnk(aqhfile)
374        call cal_TimeInterval( vstressperiod, '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       &'   Meridional wind stress forcing starts at    ',       &'   Atmospheric specific humidity starts at     ',
381       &    (vstressstartdate(i), i=1,2), dayofweek(vstressstartdate(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       &'   Meridional wind stress 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       &'   Meridional wind stress 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       &'   >>  ',vstressfile(1:il),'  <<'       &'   >>  ',aqhfile(1:il),'  <<'
396        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
397       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
 #endif  
   
 #else  
 c     Heat flux.  
       il = ilnblnk(hfluxfile)  
       call cal_TimeInterval( hfluxperiod, 'secs', timeint, mythid )  
398    
399    c--   Net longwave.
400          il = ilnblnk(lwfluxfile)
401          call cal_TimeInterval( lwfluxperiod, '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       &'   Heat flux forcing starts at                 ',       &'   Net longwave flux forcing starts at       ',
408       &    (hfluxstartdate(i), i=1,2), dayofweek(hfluxstartdate(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       &'   Heat 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       &'   Heat 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       &'   >>  ',hfluxfile(1:il),'  <<'       &'   >>  ',lwfluxfile(1:il),'  <<'
423        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
424       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
425    
426  c     Salt flux.  c--   Precipitation.
427        il = ilnblnk(sfluxfile)        il = ilnblnk(precipfile)
428        call cal_TimeInterval( sfluxperiod, 'secs', timeint, mythid )        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       &'   Salt flux forcing starts at                 ',       &'   Precipitation data set starts at            ',
435       &    (sfluxstartdate(i), i=1,2), dayofweek(sfluxstartdate(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       &'   Salt flux 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       &'   Salt flux 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       &'   >>  ',sfluxfile(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     Zonal wind stress.  #ifdef EXF_READ_EVAP
455        il = ilnblnk(ustressfile)  c--   Evaporation.
456        call cal_TimeInterval( ustressperiod, 'secs', timeint, mythid )        write(msgbuf,'(a)')
457         &'// EXF_READ_EVAP:                      defined'
458          call print_message( msgbuf, standardmessageunit,
459         &                    SQUEEZE_RIGHT , mythid)
460          il = ilnblnk(evapfile)
461          call cal_TimeInterval( evapperiod, '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,
485         &                    SQUEEZE_RIGHT , mythid)
486    #else
487          write(msgbuf,'(a)')
488         &'// EXF_READ_EVAP:                  NOT defined'
489        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
490       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
491    #endif
492    
493  c     Meridional wind stress.  #ifdef ALLOW_RUNOFF
494        il = ilnblnk(vstressfile)  c--   Runoff.
495        call cal_TimeInterval( vstressperiod, 'secs', timeint, mythid )        write(msgbuf,'(a)')
496         &'// ALLOW_RUNOFF:                       defined'
497          call print_message( msgbuf, standardmessageunit,
498         &                    SQUEEZE_RIGHT , mythid)
499          il = ilnblnk(runofffile)
500          call cal_TimeInterval( runoffperiod, '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       &'   Meridional wind stress forcing starts at    ',       &'   Runnoff starts at     ',
508       &    (vstressstartdate(i), i=1,2), dayofweek(vstressstartdate(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       &'   Meridional wind stress 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       &'   Meridional wind stress 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       &'   >>  ',vstressfile(1:il),'  <<'       &'   >>  ',runofffile(1:il),'  <<'
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)')
527         &'// ALLOW_RUNOFF:                   NOT defined'
528          call print_message( msgbuf, standardmessageunit,
529         &                    SQUEEZE_RIGHT , mythid)
530    #endif /* ALLOW_RUNOFF */
531    
532  #ifdef ALLOW_KPP  #ifdef DOWNWARD_RADIATION
533  c     Short wave radiative flux.  c--   Downward shortwave.
534        il = ilnblnk(swfluxfile)        il = ilnblnk(swdownfile)
535        call cal_TimeInterval( swfluxperiod, 'secs', timeint, mythid )        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       &'   Short wave rad. flux forcing starts at      ',       &'   Downward shortwave flux forcing starts at      ',
542       &    (swfluxstartdate(i), i=1,2), dayofweek(swfluxstartdate(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       &'   Short wave rad. flux 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       &'   Short wave rad. flux 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       &'   >>  ',swfluxfile(1:il),'  <<'       &'   >>  ',swdownfile(1:il),'  <<'
557        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
558       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
 #endif  
559    
560    c--   Downward longwave.
561          il = ilnblnk(lwdownfile)
562          call cal_TimeInterval( lwdownperiod, 'secs', timeint, mythid )
563          write(msgbuf,'(a)')
564         &' '
565          call print_message( msgbuf, standardmessageunit,
566         &                    SQUEEZE_RIGHT , mythid)
567          write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
568         &'   Downward longwave flux forcing starts at       ',
569         &    (lwdownstartdate(i), i=1,2), dayofweek(lwdownstartdate(4)),
570         &    '.'
571          call print_message( msgbuf, standardmessageunit,
572         &                    SQUEEZE_RIGHT , mythid)
573          write(msgbuf,'(a,i9.8,i7.6)')
574         &'   Downward longwave flux forcing period is       ',
575         &    (timeint(i), i=1,2)
576          call print_message( msgbuf, standardmessageunit,
577         &                    SQUEEZE_RIGHT , mythid)
578          write(msgbuf,'(a)')
579         &'   Downward longwave flux forcing is read from file:'
580          call print_message( msgbuf, standardmessageunit,
581         &                    SQUEEZE_RIGHT , mythid)
582          write(msgbuf,'(a,a,a)')
583         &'   >>  ',lwdownfile(1:il),'  <<'
584          call print_message( msgbuf, standardmessageunit,
585         &                    SQUEEZE_RIGHT , mythid)
586  #endif  #endif
587    
588  #ifdef ALLOW_KPP  #ifdef ATMOSPHERIC_LOADING
589    c--   Atmospheric pressure.
590          il = ilnblnk(apressurefile)
591          call cal_TimeInterval( apressureperiod, 'secs', timeint, mythid )
592        write(msgbuf,'(a)')        write(msgbuf,'(a)')
593       &'// ALLOW_KPP:                          defined'       &' '
594          call print_message( msgbuf, standardmessageunit,
595         &                    SQUEEZE_RIGHT , mythid)
596          write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')
597         &'   Atmospheric pressure forcing starts at      ',
598         &    (apressurestartdate(i), i=1,2),
599         &    dayofweek(apressurestartdate(4)), '.'
600          call print_message( msgbuf, standardmessageunit,
601         &                    SQUEEZE_RIGHT , mythid)
602          write(msgbuf,'(a,i9.8,i7.6)')
603         &'   Atmospheric pressure forcing period is      ',
604         &    (timeint(i), i=1,2)
605        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
606       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
 #else  
607        write(msgbuf,'(a)')        write(msgbuf,'(a)')
608       &'// ALLOW_KPP:                      NOT defined'       &'   Atmospheric pressureforcing is read from file:'
609          call print_message( msgbuf, standardmessageunit,
610         &                    SQUEEZE_RIGHT , mythid)
611          write(msgbuf,'(a,a,a)')
612         &'   >>  ',apressurefile(1:il),'  <<'
613        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
614       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
615  #endif  #endif
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22