/[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.19 by jmc, Tue Apr 28 23:27:24 2009 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_RL( exf_monFreq, INDEX_NONE, 'exf_monFreq =',
101         &  ' /* EXF monitor frequency [ s ] */')
102          CALL WRITE_0D_RL( repeatPeriod, INDEX_NONE, 'repeatPeriod =',
103         &  ' /* period for cycling forcing dataset [ s ] */')
104          CALL WRITE_0D_RL( climtempfreeze, INDEX_NONE,'climTempFreeze=',
105         &  ' /* Minimum climatological temperature [deg.C] */')
106          CALL WRITE_0D_RL( windStressMax, INDEX_NONE,'windStressMax =',
107         &  ' /* Maximum absolute windstress [ Pa ] */')
108          CALL WRITE_0D_L( stressIsOnCgrid,INDEX_NONE,'stressIsOnCgrid =',
109         &  ' /* set u,v_stress on Arakawa C-grid */')
110          CALL WRITE_0D_RL( cen2kel, INDEX_NONE, 'cen2kel =',
111         &  ' /* conversion of deg. Centigrade to Kelvin [K] */')
112          CALL WRITE_0D_RL( gravity_mks, INDEX_NONE, 'gravity_mks=',
113         &  ' /* gravitational acceleration [m/s^2] */')
114          CALL WRITE_0D_RL( atmrho, INDEX_NONE, 'atmrho =',
115         & '  /* mean atmospheric density [kg/m^3] */')
116          CALL WRITE_0D_RL( atmcp, INDEX_NONE, 'atmcp =',
117         & '  /* mean atmospheric specific heat [J/kg/K] */')
118          CALL WRITE_0D_RL( flamb, INDEX_NONE, 'flamb =',
119         & '  /* latent heat of evaporation [J/kg] */')
120          CALL WRITE_0D_RL( flami, INDEX_NONE, 'flami =',
121         & '  /* latent heat of pure-ice melting [J/kg] */')
122          CALL WRITE_0D_RL( cvapor_fac, INDEX_NONE, 'cvapor_fac =',
123         &  ' /* const. for Saturation calculation [?] */')
124          CALL WRITE_0D_RL( cvapor_exp, INDEX_NONE, 'cvapor_exp =',
125         &  ' /* const. for Saturation calculation [?] */')
126          CALL WRITE_0D_RL( cvapor_fac_ice, INDEX_NONE, 'cvapor_fac_ice=',
127         &  ' /* const. for Saturation calculation [?] */')
128          CALL WRITE_0D_RL( cvapor_exp_ice, INDEX_NONE, 'cvapor_exp_ice=',
129         &  ' /* const. for Saturation calculation [?] */')
130          CALL WRITE_0D_RL( humid_fac, INDEX_NONE, 'humid_fac =',
131         &  ' /* humidity coef. in virtual temp. [(kg/kg)^-1] */')
132          CALL WRITE_0D_RL( gamma_blk, INDEX_NONE, 'gamma_blk =',
133         &  ' /* adiabatic lapse rate [?] */')
134          CALL WRITE_0D_RL( saltsat, INDEX_NONE, 'saltsat =',
135         &  ' /* reduction of Qsat over salty water [-] */')
136          CALL WRITE_0D_L( noNegativeEvap,INDEX_NONE,'noNegativeEvap =',
137         &  ' /* prevent negative Evaporation */')
138          CALL WRITE_0D_RL( sstExtrapol, INDEX_NONE, 'sstExtrapol =',
139         &  ' /* extrapolation coeff from lev. 1 & 2 to surf [-] */')
140          CALL WRITE_0D_RL( cdrag_1, INDEX_NONE, 'cDrag_1 =',
141         &  ' /* coef used in drag calculation [?] */')
142          CALL WRITE_0D_RL( cdrag_2, INDEX_NONE, 'cDrag_2 =',
143         &  ' /* coef used in drag calculation [?] */')
144          CALL WRITE_0D_RL( cdrag_3, INDEX_NONE, 'cDrag_3 =',
145         &  ' /* coef used in drag calculation [?] */')
146          CALL WRITE_0D_RL( cstanton_1, INDEX_NONE, 'cStanton_1 =',
147         &  ' /* coef used in Stanton number calculation [?] */')
148          CALL WRITE_0D_RL( cstanton_2, INDEX_NONE, 'cStanton_2 =',
149         &  ' /* coef used in Stanton number calculation [?] */')
150          CALL WRITE_0D_RL( cdalton, INDEX_NONE, 'cDalton =',
151         &  ' /* coef used in Dalton number calculation [?] */')
152          CALL WRITE_0D_RL( exf_scal_BulkCdn, INDEX_NONE,
153         &  'exf_scal_BulkCdn=',
154         &  ' /* Drag coefficient scaling factor [-] */')
155          CALL WRITE_0D_RL( zolmin, INDEX_NONE, 'zolmin =',
156         &  ' /* minimum stability parameter [?] */')
157          CALL WRITE_0D_RL( psim_fac, INDEX_NONE, 'psim_fac =',
158         &  ' /* coef used in turbulent fluxes calculation [-] */')
159          CALL WRITE_0D_RL( zref, INDEX_NONE, 'zref =',
160         & '  /* reference height [ m ] */')
161          CALL WRITE_0D_RL( hu, INDEX_NONE, 'hu =',
162         & '  /* height of mean wind [ m ] */')
163          CALL WRITE_0D_RL( ht, INDEX_NONE, 'ht =',
164         & '  /* height of mean temperature [ m ] */')
165          CALL WRITE_0D_RL( hq, INDEX_NONE, 'hq =',
166         & '  /* height of mean spec.humidity [ m ] */')
167          CALL WRITE_0D_RL( umin, INDEX_NONE, 'uMin =',
168         &  ' /* minimum wind speed [m/s] */')
169          CALL WRITE_0D_L( useStabilityFct_overIce, INDEX_NONE,
170         &  'useStabilityFct_overIce=',
171         &  ' /* transfert Coeffs over sea-ice depend on stability */')
172          CALL WRITE_0D_RL( exf_iceCd, INDEX_NONE, 'exf_iceCd =',
173         &  ' /* drag coefficient over sea-ice (fixed) [-] */')
174          CALL WRITE_0D_RL( exf_iceCe, INDEX_NONE, 'exf_iceCe =',
175         &  ' /* transfert coeff. over sea-ice, for Evap (fixed) [-] */')
176          CALL WRITE_0D_RL( exf_iceCh, INDEX_NONE, 'exf_iceCh =',
177         &  ' /* transfert coeff. over sea-ice, Sens.Heat.(fixed)[-] */')
178          CALL WRITE_0D_RL( exf_albedo, INDEX_NONE, 'exf_albedo =',
179         & '  /* Sea-water albedo [-] */')
180          CALL WRITE_0D_RL( ocean_emissivity, INDEX_NONE,
181         &  'ocean_emissivity =',
182         &  ' /* longwave ocean-surface emissivity [-] */')
183          CALL WRITE_0D_RL( ice_emissivity, INDEX_NONE,'ice_emissivity =',
184         &  ' /* longwave seaice emissivity [-] */')
185          CALL WRITE_0D_RL(snow_emissivity, INDEX_NONE,'snow_emissivity =',
186         &  ' /* longwave snow  emissivity [-] */')
187          WRITE(msgBuf,'(A)') ' '
188          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
189         &                    SQUEEZE_RIGHT , myThid )
190    
191    C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
192    C--   Print settings of some CPP flags.
193          WRITE(msgBuf,'(A)') ' EXF main CPP flags:'
194          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
195         &                    SQUEEZE_RIGHT , myThid )
196          WRITE(msgBuf,'(A)') ' '
197          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
198         &                    SQUEEZE_RIGHT , myThid )
199    
200  #ifdef ALLOW_ATM_TEMP  #ifdef ALLOW_ATM_TEMP
201        write(msgbuf,'(a)')        write(msgbuf,'(a)')
# Line 131  c--   Print settings of some CPP flags. Line 245  c--   Print settings of some CPP flags.
245       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
246  #endif  #endif
247    
248  c--   For each data set used the summary prints the calendar data  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
249  c     and the corresponding file from which the data will be read.  
250    C--   For each data set used the summary prints the calendar data
251    C     and the corresponding file from which the data will be read.
252    
253  #ifndef ALLOW_ATM_WIND  #ifndef ALLOW_ATM_WIND
254  c--   Zonal wind stress.  C--   Zonal wind stress.
255        il = ilnblnk(ustressfile)        il = ilnblnk(ustressfile)
256        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
257        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
258       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
259        write(msgbuf,'(a,f12.0)')        write(msgbuf,'(a,f12.0)')
# Line 160  c--   Zonal wind stress. Line 275  c--   Zonal wind stress.
275        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
276       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
277    
278  c--   Meridional wind stress.  C--   Meridional wind stress.
279        il = ilnblnk(vstressfile)        il = ilnblnk(vstressfile)
280        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
281        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
282       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
283        write(msgbuf,'(a,f12.0)')        write(msgbuf,'(a,f12.0)')
# Line 187  c--   Meridional wind stress. Line 301  c--   Meridional wind stress.
301  #endif  #endif
302    
303  #ifndef ALLOW_ATM_TEMP  #ifndef ALLOW_ATM_TEMP
304  c--   Heat flux.  C--   Heat flux.
305        il = ilnblnk(hfluxfile)        il = ilnblnk(hfluxfile)
306        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
307        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
308       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
309        write(msgbuf,'(a,f12.0)')        write(msgbuf,'(a,f12.0)')
# Line 212  c--   Heat flux. Line 325  c--   Heat flux.
325        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
326       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
327    
328  c--   Salt flux.  C--   Salt flux.
329        il = ilnblnk(sfluxfile)        il = ilnblnk(sfluxfile)
330        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
331        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
332       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
333        write(msgbuf,'(a,f12.0)')        write(msgbuf,'(a,f12.0)')
# Line 239  c--   Salt flux. Line 351  c--   Salt flux.
351  #endif  #endif
352    
353  #if defined(ALLOW_ATM_TEMP) || defined(SHORTWAVE_HEATING)  #if defined(ALLOW_ATM_TEMP) || defined(SHORTWAVE_HEATING)
354  c--   Net shortwave.  C--   Net shortwave.
355        il = ilnblnk(swfluxfile)        il = ilnblnk(swfluxfile)
356        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
357        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
358       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
359        write(msgbuf,'(a,f12.0)')        write(msgbuf,'(a,f12.0)')
# Line 266  c--   Net shortwave. Line 377  c--   Net shortwave.
377  #endif  #endif
378    
379  #ifdef ALLOW_ATM_WIND  #ifdef ALLOW_ATM_WIND
380  c--   Zonal wind.  C--   Zonal wind.
381        il = ilnblnk(uwindfile)        il = ilnblnk(uwindfile)
382        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
383        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
384       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
385        write(msgbuf,'(a,f12.0)')        write(msgbuf,'(a,f12.0)')
# Line 291  c--   Zonal wind. Line 401  c--   Zonal wind.
401        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
402       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
403    
404  c--   Meridional wind.  C--   Meridional wind.
405        il = ilnblnk(vwindfile)        il = ilnblnk(vwindfile)
406        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
407        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
408       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
409        write(msgbuf,'(a,f12.0)')        write(msgbuf,'(a,f12.0)')
# Line 318  c--   Meridional wind. Line 427  c--   Meridional wind.
427  #endif  #endif
428    
429  #ifdef ALLOW_ATM_TEMP  #ifdef ALLOW_ATM_TEMP
430  c--   Atmospheric temperature.  C--   Atmospheric temperature.
431        il = ilnblnk(atempfile)        il = ilnblnk(atempfile)
432        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
433        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
434       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
435        write(msgbuf,'(a,f12.0)')        write(msgbuf,'(a,f12.0)')
# Line 343  c--   Atmospheric temperature. Line 451  c--   Atmospheric temperature.
451        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
452       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
453    
454  c--   Atmospheric specific humidity.  C--   Atmospheric specific humidity.
455        il = ilnblnk(aqhfile)        il = ilnblnk(aqhfile)
456        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
457        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
458       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
459        write(msgbuf,'(a,f12.0)')        write(msgbuf,'(a,f12.0)')
# Line 368  c--   Atmospheric specific humidity. Line 475  c--   Atmospheric specific humidity.
475        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
476       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
477    
478  c--   Net longwave.  C--   Net longwave.
479        il = ilnblnk(lwfluxfile)        il = ilnblnk(lwfluxfile)
480        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
481        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
482       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
483        write(msgbuf,'(a,f12.0)')        write(msgbuf,'(a,f12.0)')
# Line 393  c--   Net longwave. Line 499  c--   Net longwave.
499        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
500       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
501    
502  c--   Precipitation.  C--   Precipitation.
503        il = ilnblnk(precipfile)        il = ilnblnk(precipfile)
504        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
505        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
506       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
507        write(msgbuf,'(a,f12.0)')        write(msgbuf,'(a,f12.0)')
# Line 419  c--   Precipitation. Line 524  c--   Precipitation.
524       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
525  #endif  #endif
526    
527    C--   Evaporation.
528          write(msgbuf,'(a)') ' '
529          call print_message( msgbuf, standardmessageunit,
530         &                    SQUEEZE_RIGHT , mythid)
531  #ifdef EXF_READ_EVAP  #ifdef EXF_READ_EVAP
 c--   Evaporation.  
532        write(msgbuf,'(a)')        write(msgbuf,'(a)')
533       &'// EXF_READ_EVAP:                      defined'       &'// EXF_READ_EVAP:                      defined'
534        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
535       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
536        il = ilnblnk(evapfile)        il = ilnblnk(evapfile)
       write(msgbuf,'(a)')  
      &' '  
       call print_message( msgbuf, standardmessageunit,  
      &                    SQUEEZE_RIGHT , mythid)  
537        write(msgbuf,'(a,f12.0)')        write(msgbuf,'(a,f12.0)')
538       &'   Evaporation starts at     ',       &'   Evaporation starts at     ',
539       &    evapstartdate       &    evapstartdate
# Line 455  c--   Evaporation. Line 559  c--   Evaporation.
559       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
560  #endif  #endif
561    
562    C--   Runoff.
563          write(msgbuf,'(a)') ' '
564          call print_message( msgbuf, standardmessageunit,
565         &                    SQUEEZE_RIGHT , mythid)
566  #ifdef ALLOW_RUNOFF  #ifdef ALLOW_RUNOFF
 c--   Runoff.  
567        write(msgbuf,'(a)')        write(msgbuf,'(a)')
568       &'// ALLOW_RUNOFF:                       defined'       &'// ALLOW_RUNOFF:                       defined'
569        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
570       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
571        il = ilnblnk(runofffile)        il = ilnblnk(runofffile)
       write(msgbuf,'(a)')  
      &' '  
       call print_message( msgbuf, standardmessageunit,  
      &                    SQUEEZE_RIGHT , mythid)  
572        write(msgbuf,'(a,f12.0)')        write(msgbuf,'(a,f12.0)')
573       &'   Runnoff starts at     ',       &'   Runoff starts at     ',
574       &    runoffstartdate       &    runoffstartdate
575        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
576       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
# Line 492  c--   Runoff. Line 595  c--   Runoff.
595  #endif /* ALLOW_RUNOFF */  #endif /* ALLOW_RUNOFF */
596    
597  #ifdef DOWNWARD_RADIATION  #ifdef DOWNWARD_RADIATION
598  c--   Downward shortwave.  C--   Downward shortwave.
599        il = ilnblnk(swdownfile)        il = ilnblnk(swdownfile)
600        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
601        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
602       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
603        write(msgbuf,'(a,f12.0)')        write(msgbuf,'(a,f12.0)')
# Line 517  c--   Downward shortwave. Line 619  c--   Downward shortwave.
619        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
620       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
621    
622  c--   Downward longwave.  C--   Downward longwave.
623        il = ilnblnk(lwdownfile)        il = ilnblnk(lwdownfile)
624        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
625        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
626       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
627        write(msgbuf,'(a,f12.0)')        write(msgbuf,'(a,f12.0)')
# Line 544  c--   Downward longwave. Line 645  c--   Downward longwave.
645  #endif  #endif
646    
647  #ifdef ATMOSPHERIC_LOADING  #ifdef ATMOSPHERIC_LOADING
648  c--   Atmospheric pressure.  C--   Atmospheric pressure.
649        il = ilnblnk(apressurefile)        il = ilnblnk(apressurefile)
650        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
651        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
652       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
653        write(msgbuf,'(a,f12.0)')        write(msgbuf,'(a,f12.0)')
# Line 570  c--   Atmospheric pressure. Line 670  c--   Atmospheric pressure.
670       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
671  #endif  #endif
672    
673    #ifdef ALLOW_ICE_AREAMASK
674    C--   fractional ice-covered area MASK.
675          il = ilnblnk(areamaskfile)
676          write(msgbuf,'(a)') ' '
677          call print_message( msgbuf, standardmessageunit,
678         &                    SQUEEZE_RIGHT , mythid)
679          write(msgbuf,'(a,f12.0)')
680         &'   fractional ice-covered area MASK starts at      ',
681         &    areamaskstartdate
682          call print_message( msgbuf, standardmessageunit,
683         &                    SQUEEZE_RIGHT , mythid)
684          write(msgbuf,'(a,f12.0)')
685         &'   fractional ice-covered area MASK period is      ',
686         &    areamaskperiod
687          call print_message( msgbuf, standardmessageunit,
688         &                    SQUEEZE_RIGHT , mythid)
689        write(msgbuf,'(a)')        write(msgbuf,'(a)')
690       &' '       &'   fractional ice-covered area MASK is read from file:'
691          call print_message( msgbuf, standardmessageunit,
692         &                    SQUEEZE_RIGHT , mythid)
693          write(msgbuf,'(a,a,a)')
694         &'   >>  ',areamaskfile(1:il),'  <<'
695          call print_message( msgbuf, standardmessageunit,
696         &                    SQUEEZE_RIGHT , mythid)
697    #endif
698    
699          write(msgbuf,'(a)') ' '
700        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
701       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
702        write(msgbuf,'(a)')        write(msgbuf,'(a)')
# Line 586  c--   Atmospheric pressure. Line 711  c--   Atmospheric pressure.
711       &'// ======================================================='       &'// ======================================================='
712        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
713       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
714        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
715        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
716       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
717    
718        write(msgbuf,'(a)')  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
719       &' '  
720        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
721       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
722        write(msgbuf,'(a)')        write(msgbuf,'(a)')
# Line 607  c--   Atmospheric pressure. Line 731  c--   Atmospheric pressure.
731       &'// ======================================================='       &'// ======================================================='
732        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
733       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
734        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
       call print_message( msgbuf, standardmessageunit,  
      &                    SQUEEZE_RIGHT , mythid)  
735        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
736       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
737    
738  c     For each data set used the summary prints the calendar data  C     For each data set used the summary prints the calendar data
739  c     and the corresponding file from which the data will be read.  C     and the corresponding file from which the data will be read.
740    
741  #ifdef ALLOW_CLIMSST_RELAXATION  #ifdef ALLOW_CLIMSST_RELAXATION
742        write(msgbuf,'(a)')        write(msgbuf,'(a)')
# Line 641  c     and the corresponding file from wh Line 762  c     and the corresponding file from wh
762       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
763  #endif  #endif
764    
765  c     The climatological data sets are assumed to contain monthly  C     The climatological data sets are assumed to contain monthly
766  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
767  c     number of intervals during a given year.  C     number of intervals during a given year.
768    
769  #ifdef ALLOW_CLIMSST_RELAXATION  #ifdef ALLOW_CLIMSST_RELAXATION
770  c     Relaxation to SST climatology.  C     Relaxation to SST climatology.
771        il = ilnblnk(climsstfile)        il = ilnblnk(climsstfile)
772        write(msgbuf,'(a)')        write(msgbuf,'(a)')
773       &' '       &' '
# Line 673  c     Relaxation to SST climatology. Line 794  c     Relaxation to SST climatology.
794  #endif  #endif
795    
796  #ifdef ALLOW_CLIMSSS_RELAXATION  #ifdef ALLOW_CLIMSSS_RELAXATION
797  c     Relaxation to SSS climatology.  C     Relaxation to SSS climatology.
798        il = ilnblnk(climsssfile)        il = ilnblnk(climsssfile)
799        write(msgbuf,'(a)')        write(msgbuf,'(a)')
800       &' '       &' '
# Line 699  c     Relaxation to SSS climatology. Line 820  c     Relaxation to SSS climatology.
820       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
821  #endif  #endif
822    
823        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
824        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
825       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
826        write(msgbuf,'(a)')        write(msgbuf,'(a)')
# Line 715  c     Relaxation to SSS climatology. Line 835  c     Relaxation to SSS climatology.
835       &'// ======================================================='       &'// ======================================================='
836        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
837       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
838        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
839        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
840       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
841    
842        end        _END_MASTER( myThid )
843    
844          RETURN
845          END
Legend:
Removed from v.1.11  
changed lines
  Added in v.1.19
  ViewVC Help
Powered by ViewVC 1.1.22