/[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.11 by heimbach, Wed Apr 18 19:55:34 2007 UTC revision 1.12 by jmc, Tue May 8 00:21:31 2007 UTC
# Line 4  C $Name$ Line 4  C $Name$
4  #include "EXF_OPTIONS.h"  #include "EXF_OPTIONS.h"
5    
6    
7        subroutine exf_summary( mythid )        SUBROUTINE EXF_SUMMARY( myThid )
8    
9  c     ==================================================================  c     ==================================================================
10  c     SUBROUTINE exf_summary  c     SUBROUTINE exf_summary
# Line 30  c     ================================== Line 30  c     ==================================
30    
31        implicit none        implicit none
32    
33  c     == global variables ==  C     == global variables ==
34    
35  #include "EEPARAMS.h"  #include "EEPARAMS.h"
36  #include "SIZE.h"  #include "SIZE.h"
37  #include "cal.h"  #include "cal.h"
38    #include "EXF_CONSTANTS.h"
39  #include "EXF_PARAM.h"  #include "EXF_PARAM.h"
40    
41  c     == routine arguments ==  C     == routine arguments ==
42    
43  c     mythid  - thread number for this instance of the routine.  C     myThid  - thread number for this instance of the routine.
44    
45        integer mythid        integer myThid
46    
47  c     == local variables ==  C     == local variables ==
48    
49        integer il        integer il
50    
51        character*(max_len_mbuf) msgbuf        character*(max_len_mbuf) msgbuf
52    
53  c     == external ==  C     == external ==
54    
55        integer  ilnblnk        integer  ilnblnk
56        external ilnblnk        external ilnblnk
57    
58  c     == end of interface ==  C     == end of interface ==
59    
60        write(msgbuf,'(a)')        _BEGIN_MASTER( myThid )
61       &' '  
62          write(msgbuf,'(a)') ' '
63        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
64       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
65        write(msgbuf,'(a)')        write(msgbuf,'(a)')
# Line 72  c     == end of interface == Line 74  c     == end of interface ==
74       &'// ======================================================='       &'// ======================================================='
75        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
76       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
77        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
       call print_message( msgbuf, standardmessageunit,  
      &                    SQUEEZE_RIGHT , mythid)  
       write(msgbuf,'(a)')  
      &' '  
78        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
79       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
80    
81  c--   Print settings of some CPP flags.  C--   Print general parameters:
82          WRITE(msgBuf,'(A)') ' EXF general parameters:'
83          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
84         &                    SQUEEZE_RIGHT , myThid )
85          WRITE(msgBuf,'(A)') ' '
86          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
87         &                    SQUEEZE_RIGHT , myThid )
88          CALL WRITE_0D_C( exf_yftype, 0, INDEX_NONE, 'exf_yftype = ',
89         &  ' /* ? */')
90          CALL WRITE_0D_I( exf_iprec, INDEX_NONE, 'exf_iprec =',
91         &  ' /* exf file precision */')
92          CALL WRITE_0D_L( useExfYearlyFields, INDEX_NONE,
93         &  'useExfYearlyFields =',
94         &  ' /* add extension _YEAR to input file names */')
95          CALL WRITE_0D_L( twoDigitYear, INDEX_NONE, 'twoDigitYear =',
96         &  ' /* use 2-digit year extension */')
97          CALL WRITE_0D_L( useExfCheckRange, INDEX_NONE,
98         &  'useExfCheckRange =',
99         &  ' /* check for fields range */')
100          CALL WRITE_0D_R8( exf_monFreq, INDEX_NONE, 'exf_monFreq =',
101         &  ' /* EXF monitor frequency [ s ] */')
102          CALL WRITE_0D_R8( repeatPeriod, INDEX_NONE, 'repeatPeriod =',
103         &  ' /* period for cycling forcing dataset [ s ] */')
104          CALL WRITE_0D_R8( climtempfreeze, INDEX_NONE,'climTempFreeze=',
105         &  ' /* Minimum climatological temperature [deg.C] */')
106          CALL WRITE_0D_R8( windStressMax, INDEX_NONE,'windStressMax =',
107         &  ' /* Maximum absolute windstress [ Pa ] */')
108          CALL WRITE_0D_R8( exf_albedo, INDEX_NONE, 'exf_albedo =',
109         & '  /* Sea-water albedo [-] */')
110          CALL WRITE_0D_R8( ocean_emissivity, INDEX_NONE,
111         &  'ocean_emissivity =',
112         &  ' /* longwave ocean-surface emissivity [-] */')
113          CALL WRITE_0D_R8( ice_emissivity, INDEX_NONE,'ice_emissivity =',
114         &  ' /* longwave seaice emissivity [-] */')
115          CALL WRITE_0D_R8(snow_emissivity, INDEX_NONE,'snow_emissivity =',
116         &  ' /* longwave snow  emissivity [-] */')
117          CALL WRITE_0D_R8( atmrho, INDEX_NONE, 'atmrho =',
118         & '  /* mean atmospheric density [kg/m^3] */')
119          CALL WRITE_0D_R8( atmcp, INDEX_NONE, 'atmcp =',
120         & '  /* mean atmospheric specific heat [J/kg/K] */')
121          CALL WRITE_0D_R8( flamb, INDEX_NONE, 'flamb =',
122         & '  /* latent heat of evaporation [J/kg] */')
123          CALL WRITE_0D_R8( flami, INDEX_NONE, 'flami =',
124         & '  /* latent heat of pure-ice melting [J/kg] */')
125          CALL WRITE_0D_R8( cvapor_fac, INDEX_NONE, 'cvapor_fac =',
126         &  ' /* const. for Saturation calculation [?] */')
127          CALL WRITE_0D_R8( cvapor_exp, INDEX_NONE, 'cvapor_exp =',
128         &  ' /* const. for Saturation calculation [?] */')
129          CALL WRITE_0D_R8( cvapor_fac_ice, INDEX_NONE, 'cvapor_fac_ice=',
130         &  ' /* const. for Saturation calculation [?] */')
131          CALL WRITE_0D_R8( cvapor_exp_ice, INDEX_NONE, 'cvapor_exp_ice=',
132         &  ' /* const. for Saturation calculation [?] */')
133          CALL WRITE_0D_R8( humid_fac, INDEX_NONE, 'humid_fac =',
134         &  ' /* humidity coef. in virtual temp. [(kg/kg)^-1] */')
135          CALL WRITE_0D_R8( gamma_blk, INDEX_NONE, 'gamma_blk =',
136         &  ' /* adiabatic lapse rate [?] */')
137          CALL WRITE_0D_R8( saltsat, INDEX_NONE, 'saltsat =',
138         &  ' /* reduction of Qsat over salty water [-] */')
139          CALL WRITE_0D_R8( cen2kel, INDEX_NONE, 'cen2kel =',
140         &  ' /* conversion of deg. Centigrade to Kelvin [K] */')
141          CALL WRITE_0D_R8( exf_scal_BulkCdn, INDEX_NONE,
142         &  'exf_scal_BulkCdn=',
143         &  ' /* Drag coefficient scaling factor [-] */')
144          CALL WRITE_0D_R8( psim_fac, INDEX_NONE, 'psim_fac =',
145         &  ' /* coef used in turbulent fluxes calculation [-] */')
146          CALL WRITE_0D_R8( zolmin, INDEX_NONE, 'zolmin =',
147         &  ' /* minimum stability parameter [?] */')
148          CALL WRITE_0D_R8( hu, INDEX_NONE, 'zref =',
149         & '  /* reference height [ m ] */')
150          CALL WRITE_0D_R8( hu, INDEX_NONE, 'hu =',
151         & '  /* height of mean wind [ m ] */')
152          CALL WRITE_0D_R8( hu, INDEX_NONE, 'ht =',
153         & '  /* height of mean temperature [ m ] */')
154          CALL WRITE_0D_R8( hu, INDEX_NONE, 'hq =',
155         & '  /* height of mean spec.humidity [ m ] */')
156          CALL WRITE_0D_R8( umin, INDEX_NONE, 'uMin =',
157         &  ' /* minimum wind speed [m/s] */')
158          CALL WRITE_0D_R8( cdrag_1, INDEX_NONE, 'cDrag_1 =',
159         &  ' /* coef used in drag calculation [?] */')
160          CALL WRITE_0D_R8( cdrag_2, INDEX_NONE, 'cDrag_2 =',
161         &  ' /* coef used in drag calculation [?] */')
162          CALL WRITE_0D_R8( cdrag_3, INDEX_NONE, 'cDrag_3 =',
163         &  ' /* coef used in drag calculation [?] */')
164          CALL WRITE_0D_R8( cstanton_1, INDEX_NONE, 'cStanton_1 =',
165         &  ' /* coef used in Stanton number calculation [?] */')
166          CALL WRITE_0D_R8( cstanton_2, INDEX_NONE, 'cStanton_2 =',
167         &  ' /* coef used in Stanton number calculation [?] */')
168          CALL WRITE_0D_R8( cdalton, INDEX_NONE, 'cDalton =',
169         &  ' /* coef used in Dalton number calculation [?] */')
170          WRITE(msgBuf,'(A)') ' '
171          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
172         &                    SQUEEZE_RIGHT , myThid )
173    
174    C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
175    C--   Print settings of some CPP flags.
176          WRITE(msgBuf,'(A)') ' EXF main CPP flags:'
177          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
178         &                    SQUEEZE_RIGHT , myThid )
179          WRITE(msgBuf,'(A)') ' '
180          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
181         &                    SQUEEZE_RIGHT , myThid )
182    
183  #ifdef ALLOW_ATM_TEMP  #ifdef ALLOW_ATM_TEMP
184        write(msgbuf,'(a)')        write(msgbuf,'(a)')
# Line 131  c--   Print settings of some CPP flags. Line 228  c--   Print settings of some CPP flags.
228       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
229  #endif  #endif
230    
231  c--   For each data set used the summary prints the calendar data  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
232  c     and the corresponding file from which the data will be read.  
233    C--   For each data set used the summary prints the calendar data
234    C     and the corresponding file from which the data will be read.
235    
236  #ifndef ALLOW_ATM_WIND  #ifndef ALLOW_ATM_WIND
237  c--   Zonal wind stress.  C--   Zonal wind stress.
238        il = ilnblnk(ustressfile)        il = ilnblnk(ustressfile)
239        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
240        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
241       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
242        write(msgbuf,'(a,f12.0)')        write(msgbuf,'(a,f12.0)')
# Line 160  c--   Zonal wind stress. Line 258  c--   Zonal wind stress.
258        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
259       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
260    
261  c--   Meridional wind stress.  C--   Meridional wind stress.
262        il = ilnblnk(vstressfile)        il = ilnblnk(vstressfile)
263        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
264        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
265       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
266        write(msgbuf,'(a,f12.0)')        write(msgbuf,'(a,f12.0)')
# Line 187  c--   Meridional wind stress. Line 284  c--   Meridional wind stress.
284  #endif  #endif
285    
286  #ifndef ALLOW_ATM_TEMP  #ifndef ALLOW_ATM_TEMP
287  c--   Heat flux.  C--   Heat flux.
288        il = ilnblnk(hfluxfile)        il = ilnblnk(hfluxfile)
289        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
290        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
291       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
292        write(msgbuf,'(a,f12.0)')        write(msgbuf,'(a,f12.0)')
# Line 212  c--   Heat flux. Line 308  c--   Heat flux.
308        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
309       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
310    
311  c--   Salt flux.  C--   Salt flux.
312        il = ilnblnk(sfluxfile)        il = ilnblnk(sfluxfile)
313        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
314        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
315       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
316        write(msgbuf,'(a,f12.0)')        write(msgbuf,'(a,f12.0)')
# Line 239  c--   Salt flux. Line 334  c--   Salt flux.
334  #endif  #endif
335    
336  #if defined(ALLOW_ATM_TEMP) || defined(SHORTWAVE_HEATING)  #if defined(ALLOW_ATM_TEMP) || defined(SHORTWAVE_HEATING)
337  c--   Net shortwave.  C--   Net shortwave.
338        il = ilnblnk(swfluxfile)        il = ilnblnk(swfluxfile)
339        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
340        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
341       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
342        write(msgbuf,'(a,f12.0)')        write(msgbuf,'(a,f12.0)')
# Line 266  c--   Net shortwave. Line 360  c--   Net shortwave.
360  #endif  #endif
361    
362  #ifdef ALLOW_ATM_WIND  #ifdef ALLOW_ATM_WIND
363  c--   Zonal wind.  C--   Zonal wind.
364        il = ilnblnk(uwindfile)        il = ilnblnk(uwindfile)
365        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
366        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
367       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
368        write(msgbuf,'(a,f12.0)')        write(msgbuf,'(a,f12.0)')
# Line 291  c--   Zonal wind. Line 384  c--   Zonal wind.
384        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
385       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
386    
387  c--   Meridional wind.  C--   Meridional wind.
388        il = ilnblnk(vwindfile)        il = ilnblnk(vwindfile)
389        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
390        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
391       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
392        write(msgbuf,'(a,f12.0)')        write(msgbuf,'(a,f12.0)')
# Line 318  c--   Meridional wind. Line 410  c--   Meridional wind.
410  #endif  #endif
411    
412  #ifdef ALLOW_ATM_TEMP  #ifdef ALLOW_ATM_TEMP
413  c--   Atmospheric temperature.  C--   Atmospheric temperature.
414        il = ilnblnk(atempfile)        il = ilnblnk(atempfile)
415        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
416        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
417       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
418        write(msgbuf,'(a,f12.0)')        write(msgbuf,'(a,f12.0)')
# Line 343  c--   Atmospheric temperature. Line 434  c--   Atmospheric temperature.
434        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
435       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
436    
437  c--   Atmospheric specific humidity.  C--   Atmospheric specific humidity.
438        il = ilnblnk(aqhfile)        il = ilnblnk(aqhfile)
439        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
440        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
441       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
442        write(msgbuf,'(a,f12.0)')        write(msgbuf,'(a,f12.0)')
# Line 368  c--   Atmospheric specific humidity. Line 458  c--   Atmospheric specific humidity.
458        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
459       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
460    
461  c--   Net longwave.  C--   Net longwave.
462        il = ilnblnk(lwfluxfile)        il = ilnblnk(lwfluxfile)
463        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
464        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
465       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
466        write(msgbuf,'(a,f12.0)')        write(msgbuf,'(a,f12.0)')
# Line 393  c--   Net longwave. Line 482  c--   Net longwave.
482        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
483       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
484    
485  c--   Precipitation.  C--   Precipitation.
486        il = ilnblnk(precipfile)        il = ilnblnk(precipfile)
487        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
488        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
489       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
490        write(msgbuf,'(a,f12.0)')        write(msgbuf,'(a,f12.0)')
# Line 419  c--   Precipitation. Line 507  c--   Precipitation.
507       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
508  #endif  #endif
509    
510    C--   Evaporation.
511          write(msgbuf,'(a)') ' '
512          call print_message( msgbuf, standardmessageunit,
513         &                    SQUEEZE_RIGHT , mythid)
514  #ifdef EXF_READ_EVAP  #ifdef EXF_READ_EVAP
 c--   Evaporation.  
515        write(msgbuf,'(a)')        write(msgbuf,'(a)')
516       &'// EXF_READ_EVAP:                      defined'       &'// EXF_READ_EVAP:                      defined'
517        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
518       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
519        il = ilnblnk(evapfile)        il = ilnblnk(evapfile)
       write(msgbuf,'(a)')  
      &' '  
       call print_message( msgbuf, standardmessageunit,  
      &                    SQUEEZE_RIGHT , mythid)  
520        write(msgbuf,'(a,f12.0)')        write(msgbuf,'(a,f12.0)')
521       &'   Evaporation starts at     ',       &'   Evaporation starts at     ',
522       &    evapstartdate       &    evapstartdate
# Line 455  c--   Evaporation. Line 542  c--   Evaporation.
542       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
543  #endif  #endif
544    
545    C--   Runoff.
546          write(msgbuf,'(a)') ' '
547          call print_message( msgbuf, standardmessageunit,
548         &                    SQUEEZE_RIGHT , mythid)
549  #ifdef ALLOW_RUNOFF  #ifdef ALLOW_RUNOFF
 c--   Runoff.  
550        write(msgbuf,'(a)')        write(msgbuf,'(a)')
551       &'// ALLOW_RUNOFF:                       defined'       &'// ALLOW_RUNOFF:                       defined'
552        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
553       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
554        il = ilnblnk(runofffile)        il = ilnblnk(runofffile)
       write(msgbuf,'(a)')  
      &' '  
       call print_message( msgbuf, standardmessageunit,  
      &                    SQUEEZE_RIGHT , mythid)  
555        write(msgbuf,'(a,f12.0)')        write(msgbuf,'(a,f12.0)')
556       &'   Runnoff starts at     ',       &'   Runnoff starts at     ',
557       &    runoffstartdate       &    runoffstartdate
# Line 492  c--   Runoff. Line 578  c--   Runoff.
578  #endif /* ALLOW_RUNOFF */  #endif /* ALLOW_RUNOFF */
579    
580  #ifdef DOWNWARD_RADIATION  #ifdef DOWNWARD_RADIATION
581  c--   Downward shortwave.  C--   Downward shortwave.
582        il = ilnblnk(swdownfile)        il = ilnblnk(swdownfile)
583        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
584        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
585       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
586        write(msgbuf,'(a,f12.0)')        write(msgbuf,'(a,f12.0)')
# Line 517  c--   Downward shortwave. Line 602  c--   Downward shortwave.
602        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
603       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
604    
605  c--   Downward longwave.  C--   Downward longwave.
606        il = ilnblnk(lwdownfile)        il = ilnblnk(lwdownfile)
607        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
608        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
609       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
610        write(msgbuf,'(a,f12.0)')        write(msgbuf,'(a,f12.0)')
# Line 544  c--   Downward longwave. Line 628  c--   Downward longwave.
628  #endif  #endif
629    
630  #ifdef ATMOSPHERIC_LOADING  #ifdef ATMOSPHERIC_LOADING
631  c--   Atmospheric pressure.  C--   Atmospheric pressure.
632        il = ilnblnk(apressurefile)        il = ilnblnk(apressurefile)
633        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
634        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
635       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
636        write(msgbuf,'(a,f12.0)')        write(msgbuf,'(a,f12.0)')
# Line 570  c--   Atmospheric pressure. Line 653  c--   Atmospheric pressure.
653       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
654  #endif  #endif
655    
656        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
657        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
658       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
659        write(msgbuf,'(a)')        write(msgbuf,'(a)')
# Line 586  c--   Atmospheric pressure. Line 668  c--   Atmospheric pressure.
668       &'// ======================================================='       &'// ======================================================='
669        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
670       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
671        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
672        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
673       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
674    
675        write(msgbuf,'(a)')  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
676       &' '  
677        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
678       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
679        write(msgbuf,'(a)')        write(msgbuf,'(a)')
# Line 607  c--   Atmospheric pressure. Line 688  c--   Atmospheric pressure.
688       &'// ======================================================='       &'// ======================================================='
689        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
690       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
691        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
       call print_message( msgbuf, standardmessageunit,  
      &                    SQUEEZE_RIGHT , mythid)  
692        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
693       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
694    
695  c     For each data set used the summary prints the calendar data  C     For each data set used the summary prints the calendar data
696  c     and the corresponding file from which the data will be read.  C     and the corresponding file from which the data will be read.
697    
698  #ifdef ALLOW_CLIMSST_RELAXATION  #ifdef ALLOW_CLIMSST_RELAXATION
699        write(msgbuf,'(a)')        write(msgbuf,'(a)')
# Line 641  c     and the corresponding file from wh Line 719  c     and the corresponding file from wh
719       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
720  #endif  #endif
721    
722  c     The climatological data sets are assumed to contain monthly  C     The climatological data sets are assumed to contain monthly
723  c     data. This can be changed in a later version to an arbitrary  C     data. This can be changed in a later version to an arbitrary
724  c     number of intervals during a given year.  C     number of intervals during a given year.
725    
726  #ifdef ALLOW_CLIMSST_RELAXATION  #ifdef ALLOW_CLIMSST_RELAXATION
727  c     Relaxation to SST climatology.  C     Relaxation to SST climatology.
728        il = ilnblnk(climsstfile)        il = ilnblnk(climsstfile)
729        write(msgbuf,'(a)')        write(msgbuf,'(a)')
730       &' '       &' '
# Line 673  c     Relaxation to SST climatology. Line 751  c     Relaxation to SST climatology.
751  #endif  #endif
752    
753  #ifdef ALLOW_CLIMSSS_RELAXATION  #ifdef ALLOW_CLIMSSS_RELAXATION
754  c     Relaxation to SSS climatology.  C     Relaxation to SSS climatology.
755        il = ilnblnk(climsssfile)        il = ilnblnk(climsssfile)
756        write(msgbuf,'(a)')        write(msgbuf,'(a)')
757       &' '       &' '
# Line 699  c     Relaxation to SSS climatology. Line 777  c     Relaxation to SSS climatology.
777       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
778  #endif  #endif
779    
780        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
781        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
782       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
783        write(msgbuf,'(a)')        write(msgbuf,'(a)')
# Line 715  c     Relaxation to SSS climatology. Line 792  c     Relaxation to SSS climatology.
792       &'// ======================================================='       &'// ======================================================='
793        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
794       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
795        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
796        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
797       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
798    
799        end        _END_MASTER( myThid )
800    
801          RETURN
802          END
Legend:
Removed from v.1.11  
changed lines
  Added in v.1.12
  ViewVC Help
Powered by ViewVC 1.1.22