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