/[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.8 by jmc, Mon Apr 16 23:27:21 2007 UTC revision 1.14 by jmc, Sat May 12 23:17:18 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(        SUBROUTINE EXF_SUMMARY( myThid )
      I                        mythid  
      &                      )  
8    
9  c     ==================================================================  c     ==================================================================
10  c     SUBROUTINE exf_Summary  c     SUBROUTINE exf_summary
11  c     ==================================================================  c     ==================================================================
12  c  c
13  c     o List all the settings of the external forcing.  c     o List all the settings of the external forcing.
# Line 27  c     changed: Dimitris Menemenlis menem Line 25  c     changed: Dimitris Menemenlis menem
25  c              - modifications for using pkg/exf with pkg/seaice  c              - modifications for using pkg/exf with pkg/seaice
26  c  c
27  c     ==================================================================  c     ==================================================================
28  c     SUBROUTINE exf_Summary  c     SUBROUTINE exf_summary
29  c     ==================================================================  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.h"  #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 75  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,a)')  
      &'External forcing version: ',externalforcingversion  
       call print_message( msgbuf, standardmessageunit,  
      &                    SQUEEZE_RIGHT , mythid)  
       write(msgbuf,'(a,a)')  
      &'Uses Calendar version:    ',usescalendarversion  
       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( cen2kel, INDEX_NONE, 'cen2kel =',
109         &  ' /* conversion of deg. Centigrade to Kelvin [K] */')
110          CALL WRITE_0D_R8( gravity_mks, INDEX_NONE, 'gravity_mks=',
111         &  ' /* gravitational acceleration [m/s^2] */')
112          CALL WRITE_0D_R8( atmrho, INDEX_NONE, 'atmrho =',
113         & '  /* mean atmospheric density [kg/m^3] */')
114          CALL WRITE_0D_R8( atmcp, INDEX_NONE, 'atmcp =',
115         & '  /* mean atmospheric specific heat [J/kg/K] */')
116          CALL WRITE_0D_R8( flamb, INDEX_NONE, 'flamb =',
117         & '  /* latent heat of evaporation [J/kg] */')
118          CALL WRITE_0D_R8( flami, INDEX_NONE, 'flami =',
119         & '  /* latent heat of pure-ice melting [J/kg] */')
120          CALL WRITE_0D_R8( cvapor_fac, INDEX_NONE, 'cvapor_fac =',
121         &  ' /* const. for Saturation calculation [?] */')
122          CALL WRITE_0D_R8( cvapor_exp, INDEX_NONE, 'cvapor_exp =',
123         &  ' /* const. for Saturation calculation [?] */')
124          CALL WRITE_0D_R8( cvapor_fac_ice, INDEX_NONE, 'cvapor_fac_ice=',
125         &  ' /* const. for Saturation calculation [?] */')
126          CALL WRITE_0D_R8( cvapor_exp_ice, INDEX_NONE, 'cvapor_exp_ice=',
127         &  ' /* const. for Saturation calculation [?] */')
128          CALL WRITE_0D_R8( humid_fac, INDEX_NONE, 'humid_fac =',
129         &  ' /* humidity coef. in virtual temp. [(kg/kg)^-1] */')
130          CALL WRITE_0D_R8( gamma_blk, INDEX_NONE, 'gamma_blk =',
131         &  ' /* adiabatic lapse rate [?] */')
132          CALL WRITE_0D_R8( saltsat, INDEX_NONE, 'saltsat =',
133         &  ' /* reduction of Qsat over salty water [-] */')
134          CALL WRITE_0D_R8( cdrag_1, INDEX_NONE, 'cDrag_1 =',
135         &  ' /* coef used in drag calculation [?] */')
136          CALL WRITE_0D_R8( cdrag_2, INDEX_NONE, 'cDrag_2 =',
137         &  ' /* coef used in drag calculation [?] */')
138          CALL WRITE_0D_R8( cdrag_3, INDEX_NONE, 'cDrag_3 =',
139         &  ' /* coef used in drag calculation [?] */')
140          CALL WRITE_0D_R8( cstanton_1, INDEX_NONE, 'cStanton_1 =',
141         &  ' /* coef used in Stanton number calculation [?] */')
142          CALL WRITE_0D_R8( cstanton_2, INDEX_NONE, 'cStanton_2 =',
143         &  ' /* coef used in Stanton number calculation [?] */')
144          CALL WRITE_0D_R8( cdalton, INDEX_NONE, 'cDalton =',
145         &  ' /* coef used in Dalton number calculation [?] */')
146          CALL WRITE_0D_R8( exf_scal_BulkCdn, INDEX_NONE,
147         &  'exf_scal_BulkCdn=',
148         &  ' /* Drag coefficient scaling factor [-] */')
149          CALL WRITE_0D_R8( zolmin, INDEX_NONE, 'zolmin =',
150         &  ' /* minimum stability parameter [?] */')
151          CALL WRITE_0D_R8( psim_fac, INDEX_NONE, 'psim_fac =',
152         &  ' /* coef used in turbulent fluxes calculation [-] */')
153          CALL WRITE_0D_R8( zref, INDEX_NONE, 'zref =',
154         & '  /* reference height [ m ] */')
155          CALL WRITE_0D_R8( hu, INDEX_NONE, 'hu =',
156         & '  /* height of mean wind [ m ] */')
157          CALL WRITE_0D_R8( ht, INDEX_NONE, 'ht =',
158         & '  /* height of mean temperature [ m ] */')
159          CALL WRITE_0D_R8( hq, INDEX_NONE, 'hq =',
160         & '  /* height of mean spec.humidity [ m ] */')
161          CALL WRITE_0D_R8( umin, INDEX_NONE, 'uMin =',
162         &  ' /* minimum wind speed [m/s] */')
163          CALL WRITE_0D_L( useStabilityFct_overIce, INDEX_NONE,
164         &  'useStabilityFct_overIce=',
165         &  ' /* transfert Coeffs over sea-ice depend on stability */')
166          CALL WRITE_0D_R8( exf_iceCd, INDEX_NONE, 'exf_iceCd =',
167         &  ' /* drag coefficient over sea-ice (fixed) [-] */')
168          CALL WRITE_0D_R8( exf_iceCe, INDEX_NONE, 'exf_iceCe =',
169         &  ' /* transfert coeff. over sea-ice, for Evap (fixed) [-] */')
170          CALL WRITE_0D_R8( exf_iceCh, INDEX_NONE, 'exf_iceCh =',
171         &  ' /* transfert coeff. over sea-ice, Sens.Heat.(fixed)[-] */')
172          CALL WRITE_0D_R8( exf_albedo, INDEX_NONE, 'exf_albedo =',
173         & '  /* Sea-water albedo [-] */')
174          CALL WRITE_0D_R8( ocean_emissivity, INDEX_NONE,
175         &  'ocean_emissivity =',
176         &  ' /* longwave ocean-surface emissivity [-] */')
177          CALL WRITE_0D_R8( ice_emissivity, INDEX_NONE,'ice_emissivity =',
178         &  ' /* longwave seaice emissivity [-] */')
179          CALL WRITE_0D_R8(snow_emissivity, INDEX_NONE,'snow_emissivity =',
180         &  ' /* longwave snow  emissivity [-] */')
181          WRITE(msgBuf,'(A)') ' '
182          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
183         &                    SQUEEZE_RIGHT , myThid )
184    
185    C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
186    C--   Print settings of some CPP flags.
187          WRITE(msgBuf,'(A)') ' EXF main CPP flags:'
188          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
189         &                    SQUEEZE_RIGHT , myThid )
190          WRITE(msgBuf,'(A)') ' '
191          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
192         &                    SQUEEZE_RIGHT , myThid )
193    
194  #ifdef ALLOW_ATM_TEMP  #ifdef ALLOW_ATM_TEMP
195        write(msgbuf,'(a)')        write(msgbuf,'(a)')
# Line 142  c--   Print settings of some CPP flags. Line 239  c--   Print settings of some CPP flags.
239       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
240  #endif  #endif
241    
242  c--   For each data set used the summary prints the calendar data  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
243  c     and the corresponding file from which the data will be read.  
244    C--   For each data set used the summary prints the calendar data
245    C     and the corresponding file from which the data will be read.
246    
247  #ifndef ALLOW_ATM_WIND  #ifndef ALLOW_ATM_WIND
248  c--   Zonal wind stress.  C--   Zonal wind stress.
249        il = ilnblnk(ustressfile)        il = ilnblnk(ustressfile)
250        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
251        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
252       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
253        write(msgbuf,'(a,f12.0)')        write(msgbuf,'(a,f12.0)')
# Line 171  c--   Zonal wind stress. Line 269  c--   Zonal wind stress.
269        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
270       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
271    
272  c--   Meridional wind stress.  C--   Meridional wind stress.
273        il = ilnblnk(vstressfile)        il = ilnblnk(vstressfile)
274        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
275        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
276       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
277        write(msgbuf,'(a,f12.0)')        write(msgbuf,'(a,f12.0)')
# Line 198  c--   Meridional wind stress. Line 295  c--   Meridional wind stress.
295  #endif  #endif
296    
297  #ifndef ALLOW_ATM_TEMP  #ifndef ALLOW_ATM_TEMP
298  c--   Heat flux.  C--   Heat flux.
299        il = ilnblnk(hfluxfile)        il = ilnblnk(hfluxfile)
300        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
301        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
302       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
303        write(msgbuf,'(a,f12.0)')        write(msgbuf,'(a,f12.0)')
# Line 223  c--   Heat flux. Line 319  c--   Heat flux.
319        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
320       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
321    
322  c--   Salt flux.  C--   Salt flux.
323        il = ilnblnk(sfluxfile)        il = ilnblnk(sfluxfile)
324        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
325        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
326       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
327        write(msgbuf,'(a,f12.0)')        write(msgbuf,'(a,f12.0)')
# Line 250  c--   Salt flux. Line 345  c--   Salt flux.
345  #endif  #endif
346    
347  #if defined(ALLOW_ATM_TEMP) || defined(SHORTWAVE_HEATING)  #if defined(ALLOW_ATM_TEMP) || defined(SHORTWAVE_HEATING)
348  c--   Net shortwave.  C--   Net shortwave.
349        il = ilnblnk(swfluxfile)        il = ilnblnk(swfluxfile)
350        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
351        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
352       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
353        write(msgbuf,'(a,f12.0)')        write(msgbuf,'(a,f12.0)')
# Line 277  c--   Net shortwave. Line 371  c--   Net shortwave.
371  #endif  #endif
372    
373  #ifdef ALLOW_ATM_WIND  #ifdef ALLOW_ATM_WIND
374  c--   Zonal wind.  C--   Zonal wind.
375        il = ilnblnk(uwindfile)        il = ilnblnk(uwindfile)
376        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
377        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
378       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
379        write(msgbuf,'(a,f12.0)')        write(msgbuf,'(a,f12.0)')
# Line 302  c--   Zonal wind. Line 395  c--   Zonal wind.
395        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
396       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
397    
398  c--   Meridional wind.  C--   Meridional wind.
399        il = ilnblnk(vwindfile)        il = ilnblnk(vwindfile)
400        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
401        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
402       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
403        write(msgbuf,'(a,f12.0)')        write(msgbuf,'(a,f12.0)')
# Line 329  c--   Meridional wind. Line 421  c--   Meridional wind.
421  #endif  #endif
422    
423  #ifdef ALLOW_ATM_TEMP  #ifdef ALLOW_ATM_TEMP
424  c--   Atmospheric temperature.  C--   Atmospheric temperature.
425        il = ilnblnk(atempfile)        il = ilnblnk(atempfile)
426        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
427        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
428       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
429        write(msgbuf,'(a,f12.0)')        write(msgbuf,'(a,f12.0)')
# Line 354  c--   Atmospheric temperature. Line 445  c--   Atmospheric temperature.
445        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
446       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
447    
448  c--   Atmospheric specific humidity.  C--   Atmospheric specific humidity.
449        il = ilnblnk(aqhfile)        il = ilnblnk(aqhfile)
450        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
451        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
452       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
453        write(msgbuf,'(a,f12.0)')        write(msgbuf,'(a,f12.0)')
# Line 379  c--   Atmospheric specific humidity. Line 469  c--   Atmospheric specific humidity.
469        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
470       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
471    
472  c--   Net longwave.  C--   Net longwave.
473        il = ilnblnk(lwfluxfile)        il = ilnblnk(lwfluxfile)
474        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
475        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
476       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
477        write(msgbuf,'(a,f12.0)')        write(msgbuf,'(a,f12.0)')
# Line 404  c--   Net longwave. Line 493  c--   Net longwave.
493        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
494       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
495    
496  c--   Precipitation.  C--   Precipitation.
497        il = ilnblnk(precipfile)        il = ilnblnk(precipfile)
498        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
499        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
500       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
501        write(msgbuf,'(a,f12.0)')        write(msgbuf,'(a,f12.0)')
# Line 430  c--   Precipitation. Line 518  c--   Precipitation.
518       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
519  #endif  #endif
520    
521    C--   Evaporation.
522          write(msgbuf,'(a)') ' '
523          call print_message( msgbuf, standardmessageunit,
524         &                    SQUEEZE_RIGHT , mythid)
525  #ifdef EXF_READ_EVAP  #ifdef EXF_READ_EVAP
 c--   Evaporation.  
526        write(msgbuf,'(a)')        write(msgbuf,'(a)')
527       &'// EXF_READ_EVAP:                      defined'       &'// EXF_READ_EVAP:                      defined'
528        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
529       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
530        il = ilnblnk(evapfile)        il = ilnblnk(evapfile)
       write(msgbuf,'(a)')  
      &' '  
       call print_message( msgbuf, standardmessageunit,  
      &                    SQUEEZE_RIGHT , mythid)  
531        write(msgbuf,'(a,f12.0)')        write(msgbuf,'(a,f12.0)')
532       &'   Evaporation starts at     ',       &'   Evaporation starts at     ',
533       &    evapstartdate       &    evapstartdate
# Line 466  c--   Evaporation. Line 553  c--   Evaporation.
553       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
554  #endif  #endif
555    
556    C--   Runoff.
557          write(msgbuf,'(a)') ' '
558          call print_message( msgbuf, standardmessageunit,
559         &                    SQUEEZE_RIGHT , mythid)
560  #ifdef ALLOW_RUNOFF  #ifdef ALLOW_RUNOFF
 c--   Runoff.  
561        write(msgbuf,'(a)')        write(msgbuf,'(a)')
562       &'// ALLOW_RUNOFF:                       defined'       &'// ALLOW_RUNOFF:                       defined'
563        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
564       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
565        il = ilnblnk(runofffile)        il = ilnblnk(runofffile)
       write(msgbuf,'(a)')  
      &' '  
       call print_message( msgbuf, standardmessageunit,  
      &                    SQUEEZE_RIGHT , mythid)  
566        write(msgbuf,'(a,f12.0)')        write(msgbuf,'(a,f12.0)')
567       &'   Runnoff starts at     ',       &'   Runnoff starts at     ',
568       &    runoffstartdate       &    runoffstartdate
# Line 503  c--   Runoff. Line 589  c--   Runoff.
589  #endif /* ALLOW_RUNOFF */  #endif /* ALLOW_RUNOFF */
590    
591  #ifdef DOWNWARD_RADIATION  #ifdef DOWNWARD_RADIATION
592  c--   Downward shortwave.  C--   Downward shortwave.
593        il = ilnblnk(swdownfile)        il = ilnblnk(swdownfile)
594        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
595        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
596       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
597        write(msgbuf,'(a,f12.0)')        write(msgbuf,'(a,f12.0)')
# Line 528  c--   Downward shortwave. Line 613  c--   Downward shortwave.
613        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
614       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
615    
616  c--   Downward longwave.  C--   Downward longwave.
617        il = ilnblnk(lwdownfile)        il = ilnblnk(lwdownfile)
618        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
619        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
620       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
621        write(msgbuf,'(a,f12.0)')        write(msgbuf,'(a,f12.0)')
# Line 555  c--   Downward longwave. Line 639  c--   Downward longwave.
639  #endif  #endif
640    
641  #ifdef ATMOSPHERIC_LOADING  #ifdef ATMOSPHERIC_LOADING
642  c--   Atmospheric pressure.  C--   Atmospheric pressure.
643        il = ilnblnk(apressurefile)        il = ilnblnk(apressurefile)
644        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
645        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
646       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
647        write(msgbuf,'(a,f12.0)')        write(msgbuf,'(a,f12.0)')
# Line 581  c--   Atmospheric pressure. Line 664  c--   Atmospheric pressure.
664       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
665  #endif  #endif
666    
667        write(msgbuf,'(a)')        write(msgbuf,'(a)') ' '
      &' '  
668        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
669       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
670        write(msgbuf,'(a)')        write(msgbuf,'(a)')
# Line 597  c--   Atmospheric pressure. Line 679  c--   Atmospheric pressure.
679       &'// ======================================================='       &'// ======================================================='
680        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
681       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
682          write(msgbuf,'(a)') ' '
683          call print_message( msgbuf, standardmessageunit,
684         &                    SQUEEZE_RIGHT , mythid)
685    
686    C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
687    
688          call print_message( msgbuf, standardmessageunit,
689         &                    SQUEEZE_RIGHT , mythid)
690          write(msgbuf,'(a)')
691         &'// ======================================================='
692          call print_message( msgbuf, standardmessageunit,
693         &                    SQUEEZE_RIGHT , mythid)
694          write(msgbuf,'(a)')
695         &'// External forcing climatology configuration >>> START <<<'
696          call print_message( msgbuf, standardmessageunit,
697         &                    SQUEEZE_RIGHT , mythid)
698          write(msgbuf,'(a)')
699         &'// ======================================================='
700          call print_message( msgbuf, standardmessageunit,
701         &                    SQUEEZE_RIGHT , mythid)
702          write(msgbuf,'(a)') ' '
703          call print_message( msgbuf, standardmessageunit,
704         &                    SQUEEZE_RIGHT , mythid)
705    
706    C     For each data set used the summary prints the calendar data
707    C     and the corresponding file from which the data will be read.
708    
709    #ifdef ALLOW_CLIMSST_RELAXATION
710          write(msgbuf,'(a)')
711         &'// ALLOW_CLIMSST_RELAXATION:           defined'
712          call print_message( msgbuf, standardmessageunit,
713         &                    SQUEEZE_RIGHT , mythid)
714    #else
715          write(msgbuf,'(a)')
716         &'// ALLOW_CLIMSST_RELAXATION:       NOT defined'
717          call print_message( msgbuf, standardmessageunit,
718         &                    SQUEEZE_RIGHT , mythid)
719    #endif
720    
721    #ifdef ALLOW_CLIMSSS_RELAXATION
722          write(msgbuf,'(a)')
723         &'// ALLOW_CLIMSSS_RELAXATION:           defined'
724          call print_message( msgbuf, standardmessageunit,
725         &                    SQUEEZE_RIGHT , mythid)
726    #else
727          write(msgbuf,'(a)')
728         &'// ALLOW_CLIMSSS_RELAXATION:       NOT defined'
729          call print_message( msgbuf, standardmessageunit,
730         &                    SQUEEZE_RIGHT , mythid)
731    #endif
732    
733    C     The climatological data sets are assumed to contain monthly
734    C     data. This can be changed in a later version to an arbitrary
735    C     number of intervals during a given year.
736    
737    #ifdef ALLOW_CLIMSST_RELAXATION
738    C     Relaxation to SST climatology.
739          il = ilnblnk(climsstfile)
740          write(msgbuf,'(a)')
741         &' '
742          call print_message( msgbuf, standardmessageunit,
743         &                    SQUEEZE_RIGHT , mythid)
744          write(msgbuf,'(a,f12.0)')
745         &'   Climatological SST starts at         ',
746         &    climsststartdate
747          call print_message( msgbuf, standardmessageunit,
748         &                    SQUEEZE_RIGHT , mythid)
749          write(msgbuf,'(a,f12.0)')
750         &'   Climatological SST period is         ',
751         &    climsstperiod
752          call print_message( msgbuf, standardmessageunit,
753         &                    SQUEEZE_RIGHT , mythid)
754          write(msgbuf,'(a)')
755         &'   Climatological SST is read from file:'
756          call print_message( msgbuf, standardmessageunit,
757         &                    SQUEEZE_RIGHT , mythid)
758          write(msgbuf,'(a,a,a)')
759         &'   >>  ',climsstfile(1:il),'  <<'
760          call print_message( msgbuf, standardmessageunit,
761         &                    SQUEEZE_RIGHT , mythid)
762    #endif
763    
764    #ifdef ALLOW_CLIMSSS_RELAXATION
765    C     Relaxation to SSS climatology.
766          il = ilnblnk(climsssfile)
767        write(msgbuf,'(a)')        write(msgbuf,'(a)')
768       &' '       &' '
769        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
770       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
771          write(msgbuf,'(a,f12.0)')
772         &'   Climatological SSS starts at         ',
773         &    climsssstartdate
774          call print_message( msgbuf, standardmessageunit,
775         &                    SQUEEZE_RIGHT , mythid)
776          write(msgbuf,'(a,f12.0)')
777         &'   Climatological SSS period is         ',
778         &    climsssperiod
779          call print_message( msgbuf, standardmessageunit,
780         &                    SQUEEZE_RIGHT , mythid)
781          write(msgbuf,'(a)')
782         &'   Climatological SSS is read from file:'
783          call print_message( msgbuf, standardmessageunit,
784         &                    SQUEEZE_RIGHT , mythid)
785          write(msgbuf,'(a,a,a)')
786         &'   >>  ',climsssfile(1:il),'  <<'
787          call print_message( msgbuf, standardmessageunit,
788         &                    SQUEEZE_RIGHT , mythid)
789    #endif
790    
791          write(msgbuf,'(a)') ' '
792          call print_message( msgbuf, standardmessageunit,
793         &                    SQUEEZE_RIGHT , mythid)
794          write(msgbuf,'(a)')
795         &'// ======================================================='
796          call print_message( msgbuf, standardmessageunit,
797         &                    SQUEEZE_RIGHT , mythid)
798          write(msgbuf,'(a)')
799         &'// External forcing climatology configuration  >>> END <<<'
800          call print_message( msgbuf, standardmessageunit,
801         &                    SQUEEZE_RIGHT , mythid)
802          write(msgbuf,'(a)')
803         &'// ======================================================='
804          call print_message( msgbuf, standardmessageunit,
805         &                    SQUEEZE_RIGHT , mythid)
806          write(msgbuf,'(a)') ' '
807          call print_message( msgbuf, standardmessageunit,
808         &                    SQUEEZE_RIGHT , mythid)
809    
810          _END_MASTER( myThid )
811    
812        end        RETURN
813          END
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.14
  ViewVC Help
Powered by ViewVC 1.1.22