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

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

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

revision 1.6 by edhill, Thu Oct 9 04:19:19 2003 UTC revision 1.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_OPTIONS.h"  #include "EXF_OPTIONS.h"
5    
# Line 36  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 47  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 149  c     and the corresponding file from wh Line 148  c     and the corresponding file from wh
148  #ifndef ALLOW_ATM_WIND  #ifndef ALLOW_ATM_WIND
149  c--   Zonal wind stress.  c--   Zonal wind stress.
150        il = ilnblnk(ustressfile)        il = ilnblnk(ustressfile)
       call cal_TimeInterval( ustressperiod, 'secs', timeint, mythid )  
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       &'   Zonal wind stress forcing starts at         ',       &'   Zonal wind stress forcing starts at         ',
157       &    (ustressstartdate(i), i=1,2), dayofweek(ustressstartdate(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       &'   Zonal wind stress 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)')
# Line 176  c--   Zonal wind stress. Line 173  c--   Zonal wind stress.
173    
174  c--   Meridional wind stress.  c--   Meridional wind stress.
175        il = ilnblnk(vstressfile)        il = ilnblnk(vstressfile)
       call cal_TimeInterval( vstressperiod, '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       &'   Meridional wind stress forcing starts at    ',       &'   Meridional wind stress forcing starts at    ',
182       &    (vstressstartdate(i), i=1,2), dayofweek(vstressstartdate(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       &'   Meridional wind stress forcing 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)')
# Line 205  c--   Meridional wind stress. Line 200  c--   Meridional wind stress.
200  #ifndef ALLOW_ATM_TEMP  #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 231  c--   Heat flux. Line 225  c--   Heat flux.
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 259  c--   Salt flux. Line 252  c--   Salt flux.
252  #if defined(ALLOW_ATM_TEMP) || defined(SHORTWAVE_HEATING)  #if defined(ALLOW_ATM_TEMP) || defined(SHORTWAVE_HEATING)
253  c--   Net shortwave.  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       &'   Net shortwave 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       &'   Net shortwave 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)')
# Line 288  c--   Net shortwave. Line 279  c--   Net shortwave.
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 315  c--   Zonal wind. Line 304  c--   Zonal wind.
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 344  c--   Meridional wind. Line 331  c--   Meridional wind.
331  #ifdef ALLOW_ATM_TEMP  #ifdef ALLOW_ATM_TEMP
332  c--   Atmospheric temperature.  c--   Atmospheric temperature.
333        il = ilnblnk(atempfile)        il = ilnblnk(atempfile)
       call cal_TimeInterval( atempperiod, 'secs', timeint, mythid )  
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       &'   Atmospheric temperature starts at           ',       &'   Atmospheric temperature starts at           ',
340       &    (atempstartdate(i), i=1,2), dayofweek(atempstartdate(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       &'   Atmospheric temperature 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)')
# Line 371  c--   Atmospheric temperature. Line 356  c--   Atmospheric temperature.
356    
357  c--   Atmospheric specific humidity.  c--   Atmospheric specific humidity.
358        il = ilnblnk(aqhfile)        il = ilnblnk(aqhfile)
       call cal_TimeInterval( aqhperiod, '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       &'   Atmospheric specific humidity starts at     ',       &'   Atmospheric specific humidity starts at     ',
365       &    (aqhstartdate(i), i=1,2), dayofweek(aqhstartdate(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       &'   Atmospheric specific humidity 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)')
# Line 398  c--   Atmospheric specific humidity. Line 381  c--   Atmospheric specific humidity.
381    
382  c--   Net longwave.  c--   Net longwave.
383        il = ilnblnk(lwfluxfile)        il = ilnblnk(lwfluxfile)
       call cal_TimeInterval( lwfluxperiod, 'secs', timeint, mythid )  
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       &'   Net longwave flux forcing starts at       ',       &'   Net longwave flux forcing starts at       ',
390       &    (lwfluxstartdate(i), i=1,2), dayofweek(lwfluxstartdate(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       &'   Net longwave 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)')
# Line 425  c--   Net longwave. Line 406  c--   Net longwave.
406    
407  c--   Precipitation.  c--   Precipitation.
408        il = ilnblnk(precipfile)        il = ilnblnk(precipfile)
       call cal_TimeInterval( precipperiod, '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       &'   Precipitation data set starts at            ',       &'   Precipitation data set starts at            ',
415       &    (precipstartdate(i), i=1,2), dayofweek(precipstartdate(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       &'   Precipitation data 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)')
# Line 458  c--   Evaporation. Line 437  c--   Evaporation.
437        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
438       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
439        il = ilnblnk(evapfile)        il = ilnblnk(evapfile)
       call cal_TimeInterval( evapperiod, 'secs', timeint, mythid )  
   
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       &'   Evaporation starts at     ',       &'   Evaporation starts at     ',
446       &    (evapstartdate(i), i=1,2), dayofweek(evapstartdate(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       &'   Evaporation 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)')
# Line 497  c--   Runoff. Line 473  c--   Runoff.
473        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
474       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
475        il = ilnblnk(runofffile)        il = ilnblnk(runofffile)
       call cal_TimeInterval( runoffperiod, 'secs', timeint, mythid )  
   
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       &'   Runnoff starts at     ',       &'   Runnoff starts at     ',
482       &    (runoffstartdate(i), i=1,2), dayofweek(runoffstartdate(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       &'   Runoff 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)')
# Line 532  c--   Runoff. Line 505  c--   Runoff.
505  #ifdef DOWNWARD_RADIATION  #ifdef DOWNWARD_RADIATION
506  c--   Downward shortwave.  c--   Downward shortwave.
507        il = ilnblnk(swdownfile)        il = ilnblnk(swdownfile)
       call cal_TimeInterval( swdownperiod, '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       &'   Downward shortwave flux forcing starts at      ',       &'   Downward shortwave flux forcing starts at      ',
514       &    (swdownstartdate(i), i=1,2), dayofweek(swdownstartdate(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       &'   Downward shortwave 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)')
# Line 559  c--   Downward shortwave. Line 530  c--   Downward shortwave.
530    
531  c--   Downward longwave.  c--   Downward longwave.
532        il = ilnblnk(lwdownfile)        il = ilnblnk(lwdownfile)
       call cal_TimeInterval( lwdownperiod, 'secs', timeint, mythid )  
533        write(msgbuf,'(a)')        write(msgbuf,'(a)')
534       &' '       &' '
535        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
536       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
537        write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')        write(msgbuf,'(a,f12.0)')
538       &'   Downward longwave flux forcing starts at       ',       &'   Downward longwave flux forcing starts at       ',
539       &    (lwdownstartdate(i), i=1,2), dayofweek(lwdownstartdate(4)),       &    lwdownstartdate
      &    '.'  
540        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
541       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
542        write(msgbuf,'(a,i9.8,i7.6)')        write(msgbuf,'(a,f12.0)')
543       &'   Downward longwave flux forcing period is       ',       &'   Downward longwave flux forcing period is       ',
544       &    (timeint(i), i=1,2)       &    lwdownperiod
545        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
546       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
547        write(msgbuf,'(a)')        write(msgbuf,'(a)')
# Line 588  c--   Downward longwave. Line 557  c--   Downward longwave.
557  #ifdef ATMOSPHERIC_LOADING  #ifdef ATMOSPHERIC_LOADING
558  c--   Atmospheric pressure.  c--   Atmospheric pressure.
559        il = ilnblnk(apressurefile)        il = ilnblnk(apressurefile)
       call cal_TimeInterval( apressureperiod, 'secs', timeint, mythid )  
560        write(msgbuf,'(a)')        write(msgbuf,'(a)')
561       &' '       &' '
562        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
563       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
564        write(msgbuf,'(a,i9.8,i7.6,1x,a,a)')        write(msgbuf,'(a,f12.0)')
565       &'   Atmospheric pressure forcing starts at      ',       &'   Atmospheric pressure forcing starts at      ',
566       &    (apressurestartdate(i), i=1,2),       &    apressurestartdate
      &    dayofweek(apressurestartdate(4)), '.'  
567        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
568       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
569        write(msgbuf,'(a,i9.8,i7.6)')        write(msgbuf,'(a,f12.0)')
570       &'   Atmospheric pressure forcing period is      ',       &'   Atmospheric pressure forcing period is      ',
571       &    (timeint(i), i=1,2)       &    apressureperiod
572        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
573       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
574        write(msgbuf,'(a)')        write(msgbuf,'(a)')
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22