/[MITgcm]/MITgcm/model/src/config_summary.F
ViewVC logotype

Diff of /MITgcm/model/src/config_summary.F

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

revision 1.41 by jmc, Thu Apr 17 13:16:23 2003 UTC revision 1.69 by baylor, Thu Mar 10 03:45:11 2005 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2  C $Name$  C $Name$
3    
4    #include "PACKAGES_CONFIG.h"
5  #include "CPP_OPTIONS.h"  #include "CPP_OPTIONS.h"
6    
7    C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
8  CBOP  CBOP
9  C     !ROUTINE: CONFIG_SUMMARY  C     !ROUTINE: CONFIG_SUMMARY
10    
11  C     !INTERFACE:  C     !INTERFACE:
12        SUBROUTINE CONFIG_SUMMARY( myThid )        SUBROUTINE CONFIG_SUMMARY( myThid )
13  C     !DESCRIPTION: \bv  
14  C     *=========================================================*  C     !DESCRIPTION:
15  C     | SUBROUTINE CONFIG_SUMMARY                                  C     This routine summarizes the model parameter settings by writing a
16  C     | o Summarize model parameter settings.                      C     tabulated list of the kernel model configuration variables.  It
17  C     *=========================================================*  C     describes all the parameter settings in force and the meaning and
18  C     | This routine writes a tabulated summary of the kernel      C     units of those parameters. Individal packages report a similar
19  C     | model configuration. Information describes all the  C     table for each package using the same format as employed here. If
20  C     | parameter setting in force and the meaning and units of  C     parameters are missing or incorrectly described or dimensioned
21  C     | those parameters. Individal packages report a similar  C     please contact <MITgcm-support@mitgcm.org>
 C     | table for each package using the same format as employed  
 C     | here. If parameters are missing or incorrectly described  
 C     | or dimensioned please contact support@mitgcm.org  
 C     *=========================================================*  
 C     \ev  
22    
23  C     !USES:  C     !USES:
24        IMPLICIT NONE        IMPLICIT NONE
 C     === Global variables ===  
25  #include "SIZE.h"  #include "SIZE.h"
26  #include "EEPARAMS.h"  #include "EEPARAMS.h"
27  #include "PARAMS.h"  #include "PARAMS.h"
28  #include "EOS.h"  #include "EOS.h"
29  #include "GRID.h"  #include "GRID.h"
30  #include "DYNVARS.h"  #include "DYNVARS.h"
31    #ifdef ALLOW_MNC
32    #include "MNC_PARAMS.h"
33    #endif
34    
35  C     !INPUT/OUTPUT PARAMETERS:  C     !INPUT/OUTPUT PARAMETERS:
36  C     == Routine arguments ==  C     myThid ::  Number of this instance of CONFIG_SUMMARY
 C     myThid -  Number of this instance of CONFIG_SUMMARY  
37        INTEGER myThid        INTEGER myThid
38  CEndOfInterface  CEOP
39    
40  C     !LOCAL VARIABLES:  C     !LOCAL VARIABLES:
 C     == Local variables ==  
41  C     msgBuf :: Temp. for building output string.  C     msgBuf :: Temp. for building output string.
42  C     I,J,K  :: Loop counters.  C     I,J,K  :: Loop counters.
43  C     bi,bj  :: Tile loop counters.  C     bi,bj  :: Tile loop counters.
# Line 54  C     zcoord :: Line 52  C     zcoord ::
52        _RL                     rcoord(Nr+1)        _RL                     rcoord(Nr+1)
53        INTEGER coordLine        INTEGER coordLine
54        INTEGER tileLine        INTEGER tileLine
 CEOP  
55    
56    
57        _BARRIER        _BARRIER
# Line 85  CEOP Line 82  CEOP
82        WRITE(msgBuf,'(A,A40)') 'buoyancyRelation = ', buoyancyRelation        WRITE(msgBuf,'(A,A40)') 'buoyancyRelation = ', buoyancyRelation
83        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
84       &                    SQUEEZE_RIGHT , 1)       &                    SQUEEZE_RIGHT , 1)
85          CALL WRITE_0D_L( fluidIsAir,   INDEX_NONE,
86         & 'fluidIsAir  =', '  /* fluid major constituent is Air */')
87          CALL WRITE_0D_L( fluidIsWater, INDEX_NONE,
88         & 'fluidIsWater=', '  /* fuild major constituent is Water */')
89          CALL WRITE_0D_L( usingPCoords, INDEX_NONE,
90         & 'usingPCoords =', '  /* use p (or p*) vertical coordinate */')
91          CALL WRITE_0D_L( usingZCoords, INDEX_NONE,
92         & 'usingZCoords =', '  /* use z (or z*) vertical coordinate */')
93        CALL WRITE_1D_R8( tRef, Nr, INDEX_K,'tRef =',        CALL WRITE_1D_R8( tRef, Nr, INDEX_K,'tRef =',
94       &'   /* Reference temperature profile ( oC or oK ) */')       &'   /* Reference temperature profile ( oC or oK ) */')
95        CALL WRITE_1D_R8( sRef, Nr, INDEX_K,'sRef =',        CALL WRITE_1D_R8( sRef, Nr, INDEX_K,'sRef =',
96       &'   /* Reference salinity profile ( ppt ) */')       &'   /* Reference salinity profile ( ppt ) */')
97        CALL WRITE_0D_R8( viscAh, INDEX_NONE,'viscAh =',        CALL WRITE_0D_R8( viscAh, INDEX_NONE,'viscAh =',
98       &'   /* Lateral eddy viscosity ( m^2/s ) */')       &'   /* Lateral eddy viscosity ( m^2/s ) */')
99        CALL WRITE_0D_R8( viscA4, INDEX_NONE,'viscAh =',        IF ( viscAhD.NE.viscAh )
100         & CALL WRITE_0D_R8( viscAhD, INDEX_NONE,'viscAhD =',
101         & '  /* Lateral eddy viscosity (Divergence)( m^2/s ) */')
102          IF ( viscAhZ.NE.viscAh )
103         & CALL WRITE_0D_R8( viscAhZ, INDEX_NONE,'viscAhZ =',
104         & '  /* Lateral eddy viscosity (Vorticity) ( m^2/s ) */')
105          CALL WRITE_0D_R8( viscAhMax, INDEX_NONE,'viscAhMax =',
106         &'   /* Maximum lateral eddy viscosity ( m^2/s ) */')
107          CALL WRITE_0D_R8( viscAhGrid, INDEX_NONE,'viscAhGrid =',
108         &'   /* Grid dependent lateral eddy viscosity ( non-dim. ) */')
109          CALL WRITE_0D_L( useFullLeith, INDEX_NONE,
110         &'useFullLeith =',
111         &'   /* Use Full Form of Leith Viscosity on/off flag*/')
112          CALL WRITE_0D_R8( viscC2leith, INDEX_NONE,'viscC2leith =',
113         &'   /* Leith harmonic viscosity factor ( on grad(vort), non-dim. ) */')
114          CALL WRITE_0D_R8( viscC2leithD, INDEX_NONE,'viscC2leithD =',
115         &'   /* Leith harmonic viscosity factor ( on grad(div), non-dim. ) */')
116          CALL WRITE_0D_R8( viscC2smag, INDEX_NONE,'viscC2smag =',
117         &'   /* Smagorinsky harmonic viscosity factor (non-dim.) */')
118          CALL WRITE_0D_R8( viscA4, INDEX_NONE,'viscA4 =',
119       &'   /* Lateral biharmonic viscosity ( m^4/s ) */')       &'   /* Lateral biharmonic viscosity ( m^4/s ) */')
120          IF ( viscA4D.NE.viscA4 )
121         & CALL WRITE_0D_R8( viscA4D, INDEX_NONE,'viscA4D =',
122         & '  /* Lateral biharmonic viscosity (Divergence)( m^4/s ) */')
123          IF ( viscA4Z.NE.viscA4 )
124         & CALL WRITE_0D_R8( viscA4Z, INDEX_NONE,'viscA4Z =',
125         & '  /* Lateral biharmonic viscosity (Vorticity) ( m^4/s ) */')
126          CALL WRITE_0D_R8( viscA4Max, INDEX_NONE,'viscA4Max =',
127         &'   /* Maximum biharmonic viscosity ( m^2/s ) */')
128          CALL WRITE_0D_R8( viscA4Grid, INDEX_NONE,'viscA4Grid =',
129         &'   /* Grid dependent biharmonic viscosity ( non-dim. ) */')
130          CALL WRITE_0D_R8( viscC4leith, INDEX_NONE,'viscC4leith =',
131         &'   /* Leith biharmonic viscosity factor ( on grad(vort), non-dom. ) */')
132          CALL WRITE_0D_R8( viscC4leithD, INDEX_NONE,'viscC4leithD =',
133         &'   /* Leith biharmonic viscosity factor ( on grad(div), non-dom. ) */')
134        CALL WRITE_0D_L( no_slip_sides, INDEX_NONE,        CALL WRITE_0D_L( no_slip_sides, INDEX_NONE,
135       & 'no_slip_sides =', '  /* Viscous BCs: No-slip sides */')       & 'no_slip_sides =', '  /* Viscous BCs: No-slip sides */')
       IF ( viscAz .NE. UNSET_RL ) THEN  
        CALL WRITE_0D_R8( viscAz, INDEX_NONE,'viscAz =',  
      & '   /* Vertical eddy viscosity ( m^2/s ) */')  
       ENDIF  
       IF ( viscAp .NE. UNSET_RL ) THEN  
        CALL WRITE_0D_R8( viscAp, INDEX_NONE,'viscAp =',  
      & '   /* Vertical eddy viscosity ( Pa^2/s ) */')  
       ENDIF  
136        CALL WRITE_0D_R8( viscAr,  INDEX_NONE,'viscAr =',        CALL WRITE_0D_R8( viscAr,  INDEX_NONE,'viscAr =',
137       &'   /* Vertical eddy viscosity ( units of r^2/s ) */')       &'   /* Vertical eddy viscosity ( units of r^2/s ) */')
138          CALL WRITE_0D_L( no_slip_bottom, INDEX_NONE,
139         & 'no_slip_bottom =', '  /* Viscous BCs: No-slip bottom */')
140        CALL WRITE_0D_R8( diffKhT, INDEX_NONE,'diffKhT =',        CALL WRITE_0D_R8( diffKhT, INDEX_NONE,'diffKhT =',
141       &'   /* Laplacian diffusion of heat laterally ( m^2/s ) */')       &'   /* Laplacian diffusion of heat laterally ( m^2/s ) */')
142        CALL WRITE_0D_R8( diffK4T, INDEX_NONE,'diffK4T =',        CALL WRITE_0D_R8( diffK4T, INDEX_NONE,'diffK4T =',
143       &'   /* Bihaarmonic diffusion of heat laterally ( m^4/s ) */')       &'   /* Bihaarmonic diffusion of heat laterally ( m^4/s ) */')
       CALL WRITE_0D_R8( diffKzT, INDEX_NONE,'diffKzT =',  
      &'   /* Laplacian diffusion of heat vertically ( m^2/s ) */')  
       CALL WRITE_0D_R8( diffKrT, INDEX_NONE,'diffKrT =',  
      &'   /* Laplacian diffusion of heat vertically ( m^2/s ) */')  
144        CALL WRITE_0D_R8( diffKhS, INDEX_NONE,'diffKhS =',        CALL WRITE_0D_R8( diffKhS, INDEX_NONE,'diffKhS =',
145       &'   /* Laplacian diffusion of salt laterally ( m^2/s ) */')       &'   /* Laplacian diffusion of salt laterally ( m^2/s ) */')
146        CALL WRITE_0D_R8( diffK4S, INDEX_NONE,'diffK4S =',        CALL WRITE_0D_R8( diffK4S, INDEX_NONE,'diffK4S =',
147       &'   /* Bihaarmonic diffusion of salt laterally ( m^4/s ) */')       &'   /* Bihaarmonic diffusion of salt laterally ( m^4/s ) */')
148        CALL WRITE_0D_R8( diffKzS, INDEX_NONE,'diffKzS =',        CALL WRITE_1D_R8( diffKrNrT, Nr, INDEX_K,'diffKrNrT =',
149       &'   /* Laplacian diffusion of salt vertically ( m^2/s ) */')       & ' /* vertical profile of vertical diffusion of Temp ( m^2/s )*/')
150        CALL WRITE_0D_R8( diffKrS, INDEX_NONE,'diffKrS =',        CALL WRITE_1D_R8( diffKrNrS, Nr, INDEX_K,'diffKrNrS =',
151       &'   /* Laplacian diffusion of salt vertically ( m^2/s ) */')       & ' /* vertical profile of vertical diffusion of Salt ( m^2/s )*/')
152          CALL WRITE_0D_R8( diffKrBL79surf, INDEX_NONE,'diffKrBL79surf =',
153         &'   /* Surface diffusion for Bryan and Lewis 1979 ( m^2/s ) */')
154          CALL WRITE_0D_R8( diffKrBL79deep, INDEX_NONE,'diffKrBL79deep =',
155         &'   /* Deep diffusion for Bryan and Lewis 1979 ( m^2/s ) */')
156          CALL WRITE_0D_R8( diffKrBL79scl, INDEX_NONE,'diffKrBL79scl =',
157         &'   /* Depth scale for Bryan and Lewis 1979 ( m ) */')
158          CALL WRITE_0D_R8( diffKrBL79Ho, INDEX_NONE,'diffKrBL79Ho =',
159         &'   /* Turning depth for Bryan and Lewis 1979 ( m ) */')
160        WRITE(msgBuf,'(2A)') ' Equation of State : eosType = ', eosType        WRITE(msgBuf,'(2A)') ' Equation of State : eosType = ', eosType
161        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
162       &                    SQUEEZE_RIGHT , 1)       &                    SQUEEZE_RIGHT , 1)
# Line 138  CEOP Line 174  CEOP
174       &                    SQUEEZE_RIGHT , 1)       &                    SQUEEZE_RIGHT , 1)
175          ENDDO          ENDDO
176        ENDIF        ENDIF
177        IF ( buoyancyRelation .eq. 'ATMOSPHERIC' ) THEN        IF ( fluidIsAir ) THEN
178         CALL WRITE_0D_R8( atm_Rd, INDEX_NONE, 'atm_Rd =',         CALL WRITE_0D_R8( atm_Rd, INDEX_NONE, 'atm_Rd =',
179       & '  /* gas constant for dry air ( J/kg/K ) */')       & '  /* gas constant for dry air ( J/kg/K ) */')
180        CALL WRITE_0D_R8( atm_Cp, INDEX_NONE, 'atm_Cp =',        CALL WRITE_0D_R8( atm_Cp, INDEX_NONE, 'atm_Cp =',
181       & '  /* specific heat (Cp) of dry air ( J/kg/K ) */')       & '  /* specific heat (Cp) of dry air ( J/kg/K ) */')
182        CALL WRITE_0D_R8( atm_kappa, INDEX_NONE, 'atm_kappa =',        CALL WRITE_0D_R8( atm_kappa, INDEX_NONE, 'atm_kappa =',
183       & '  /* kappa (=Rd/Cp ) of dry air */')       & '  /* kappa (=Rd/Cp ) of dry air */')
184           CALL WRITE_0D_R8( atm_Rq, INDEX_NONE, 'atm_Rq =',
185         &  ' /* water vap. specific vol. anomaly relative to dry air */')
186        CALL WRITE_0D_R8( atm_Po, INDEX_NONE, 'atm_Po =',        CALL WRITE_0D_R8( atm_Po, INDEX_NONE, 'atm_Po =',
187       & '  /* standard reference pressure ( Pa ) */')       & '  /* standard reference pressure ( Pa ) */')
188        CALL WRITE_0D_I( integr_GeoPot, INDEX_NONE, 'integr_GeoPot =',        CALL WRITE_0D_I( integr_GeoPot, INDEX_NONE, 'integr_GeoPot =',
# Line 218  CEOP Line 256  CEOP
256        CALL WRITE_0D_R8( salt_EvPrRn, INDEX_NONE,        CALL WRITE_0D_R8( salt_EvPrRn, INDEX_NONE,
257       &'salt_EvPrRn =',       &'salt_EvPrRn =',
258       &' /* Salin. of Evap/Prec/R (UNSET=use local S)(ppt)*/')       &' /* Salin. of Evap/Prec/R (UNSET=use local S)(ppt)*/')
       CALL WRITE_0D_R8( trac_EvPrRn, INDEX_NONE,  
      &'trac_EvPrRn =',  
      &' /* Tracer in Evap/Prec/R (UNSET=use local Tr)*/')  
259        ELSE        ELSE
260        CALL WRITE_0D_R8( convertFW2Salt, INDEX_NONE,        CALL WRITE_0D_R8( convertFW2Salt, INDEX_NONE,
261       &'convertFW2Salt =',       &'convertFW2Salt =',
262       &' /* convert F.W. Flux to Salt Flux (-1=use local S)(ppt)*/')       &' /* convert F.W. Flux to Salt Flux (-1=use local S)(ppt)*/')
263        ENDIF        ENDIF
264    
265        CALL WRITE_0D_L( multiDimAdvection, INDEX_NONE,        CALL WRITE_0D_L( nonHydrostatic, INDEX_NONE,
266       & 'multiDimAdvection =',       & 'nonHydrostatic =', '  /* Non-Hydrostatic on/off flag */')
      &'   /* enable/disable Multi-Dim Advection */')  
       CALL WRITE_0D_L( staggerTimeStep, INDEX_NONE,  
      &                 'staggerTimeStep =',  
      &'   /* Stagger time stepping on/off flag */')  
267        CALL WRITE_0D_L( momStepping,  INDEX_NONE,        CALL WRITE_0D_L( momStepping,  INDEX_NONE,
268       & 'momStepping =', '  /* Momentum equation on/off flag */')       & 'momStepping =', '  /* Momentum equation on/off flag */')
269        CALL WRITE_0D_L( momAdvection, INDEX_NONE,        CALL WRITE_0D_L( momAdvection, INDEX_NONE,
270       & 'momAdvection =', '  /* Momentum advection on/off flag */')       & 'momAdvection =', '  /* Momentum advection on/off flag */')
271        CALL WRITE_0D_L( momViscosity, INDEX_NONE,        CALL WRITE_0D_L( momViscosity, INDEX_NONE,
272       & 'momViscosity =', '  /* Momentum viscosity on/off flag */')       & 'momViscosity =', '  /* Momentum viscosity on/off flag */')
273          CALL WRITE_0D_L( momImplVertAdv, INDEX_NONE, 'momImplVertAdv =',
274         &                '/* Momentum implicit vert. advection on/off*/')
275          CALL WRITE_0D_L( implicitViscosity, INDEX_NONE,
276         & 'implicitViscosity =', ' /* Implicit viscosity on/off flag */')
277        CALL WRITE_0D_L( useCoriolis,  INDEX_NONE,        CALL WRITE_0D_L( useCoriolis,  INDEX_NONE,
278       & 'useCoriolis =', '  /* Coriolis on/off flag */')       & 'useCoriolis =', '  /* Coriolis on/off flag */')
279          CALL WRITE_0D_L( useCDscheme,  INDEX_NONE,
280         & 'useCDscheme =', '  /* CD scheme on/off flag */')
281          CALL WRITE_0D_L( useJamartWetPoints,  INDEX_NONE,
282         & 'useJamartWetPoints=',' /* Coriolis WetPoints method flag */')
283          CALL WRITE_0D_L( useJamartMomAdv,  INDEX_NONE,
284         & 'useJamartMomAdv=',' /* V.I. Non-linear terms Jamart flag */')
285          CALL WRITE_0D_L( SadournyCoriolis,  INDEX_NONE,
286         & 'SadournyCoriolis=',' /* Sadourny Coriolis discr. flag */')
287          CALL WRITE_0D_L( upwindVorticity,  INDEX_NONE,
288         & 'upwindVorticity=',' /* Upwind bias vorticity flag */')
289          CALL WRITE_0D_L( useAbsVorticity,  INDEX_NONE,
290         & 'useAbsVorticity=',' /* Work with f+zeta in Coriolis */')
291          CALL WRITE_0D_L( highOrderVorticity,  INDEX_NONE,
292         & 'highOrderVorticity=',' /* High order interp. of vort. flag */')
293        CALL WRITE_0D_L( momForcing,   INDEX_NONE,        CALL WRITE_0D_L( momForcing,   INDEX_NONE,
294       & 'momForcing =', '  /* Momentum forcing on/off flag */')       & 'momForcing =', '  /* Momentum forcing on/off flag */')
295        CALL WRITE_0D_L( momPressureForcing, INDEX_NONE,        CALL WRITE_0D_L( momPressureForcing, INDEX_NONE,
296       & 'momPressureForcing =',         & 'momPressureForcing =',  
297       & '  /* Momentum pressure term on/off flag */')       & '  /* Momentum pressure term on/off flag */')
298          CALL WRITE_0D_L( staggerTimeStep, INDEX_NONE,
299         &                 'staggerTimeStep =',
300         &'   /* Stagger time stepping on/off flag */')
301          CALL WRITE_0D_L( multiDimAdvection, INDEX_NONE,
302         & 'multiDimAdvection =',
303         &'   /* enable/disable Multi-Dim Advection */')
304          CALL WRITE_0D_L( useMultiDimAdvec, INDEX_NONE,
305         & 'useMultiDimAdvec =',
306         &'   /* Multi-Dim Advection is/is-not used */')
307          CALL WRITE_0D_L( implicitDiffusion, INDEX_NONE,
308         & 'implicitDiffusion =','/* Implicit Diffusion on/off flag */')
309        CALL WRITE_0D_L( tempStepping,  INDEX_NONE,        CALL WRITE_0D_L( tempStepping,  INDEX_NONE,
310       & 'tempStepping =', '  /* Temperature equation on/off flag */')       & 'tempStepping =', '  /* Temperature equation on/off flag */')
311        CALL WRITE_0D_L( tempAdvection,  INDEX_NONE,        CALL WRITE_0D_L( tempAdvection,  INDEX_NONE,
312       & 'tempAdvection=', '  /* Temperature advection on/off flag */')       & 'tempAdvection=', '  /* Temperature advection on/off flag */')
313          CALL WRITE_0D_L( tempImplVertAdv,INDEX_NONE,'tempImplVertAdv =',
314         &                '/* Temp. implicit vert. advection on/off */')
315        CALL WRITE_0D_L( tempForcing,  INDEX_NONE,        CALL WRITE_0D_L( tempForcing,  INDEX_NONE,
316       & 'tempForcing  =', '  /* Temperature forcing on/off flag */')       & 'tempForcing  =', '  /* Temperature forcing on/off flag */')
317        CALL WRITE_0D_L( saltStepping,  INDEX_NONE,        CALL WRITE_0D_L( saltStepping,  INDEX_NONE,
318       & 'saltStepping =', '  /* Salinity equation on/off flag */')       & 'saltStepping =', '  /* Salinity equation on/off flag */')
319        CALL WRITE_0D_L( saltAdvection,  INDEX_NONE,        CALL WRITE_0D_L( saltAdvection,  INDEX_NONE,
320       & 'saltAdvection=', '  /* Salinity advection on/off flag */')       & 'saltAdvection=', '  /* Salinity advection on/off flag */')
321          CALL WRITE_0D_L( saltImplVertAdv,INDEX_NONE,'saltImplVertAdv =',
322         &                '/* Sali. implicit vert. advection on/off */')
323        CALL WRITE_0D_L( saltForcing,  INDEX_NONE,        CALL WRITE_0D_L( saltForcing,  INDEX_NONE,
324       & 'saltForcing  =', '  /* Salinity forcing on/off flag */')       & 'saltForcing  =', '  /* Salinity forcing on/off flag */')
       CALL WRITE_0D_L( nonHydrostatic, INDEX_NONE,  
      & 'nonHydrostatic =', '  /* Non-Hydrostatic on/off flag */')  
       CALL WRITE_0D_L( useCDscheme,  INDEX_NONE,  
      & 'useCDscheme =', '  /* CD scheme on/off flag */')  
       CALL WRITE_0D_L( useJamartWetPoints,  INDEX_NONE,  
      & 'useJamartWetPoints=',' /* Coriolis WetPoints method flag */')  
325        WRITE(msgBuf,'(A)') '//  '        WRITE(msgBuf,'(A)') '//  '
326        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
327       &                    SQUEEZE_RIGHT , 1)       &                    SQUEEZE_RIGHT , 1)
# Line 282  CEOP Line 340  CEOP
340        CALL WRITE_0D_R8( cg2dTargetResidual, INDEX_NONE,        CALL WRITE_0D_R8( cg2dTargetResidual, INDEX_NONE,
341       & 'cg2dTargetResidual =',       & 'cg2dTargetResidual =',
342       &'   /* 2d con. grad target residual  */')       &'   /* 2d con. grad target residual  */')
343          CALL WRITE_0D_R8( cg2dTargetResWunit, INDEX_NONE,
344         & 'cg2dTargetResWunit =',
345         &'   /* CG2d target residual [W units] */')
346          CALL WRITE_0D_I( cg2dPreCondFreq, INDEX_NONE,'cg2dPreCondFreq =',
347         &'   /* Freq. for updating cg2d preconditioner */')
348    
349        WRITE(msgBuf,'(A)') '//  '        WRITE(msgBuf,'(A)') '//  '
350        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
# Line 301  CEOP Line 364  CEOP
364       &'   /* Momentum equation timestep ( s ) */')       &'   /* Momentum equation timestep ( s ) */')
365        CALL WRITE_0D_R8( deltaTfreesurf,INDEX_NONE,'deltaTfreesurf =',        CALL WRITE_0D_R8( deltaTfreesurf,INDEX_NONE,'deltaTfreesurf =',
366       &'   /* FreeSurface equation timestep ( s ) */')       &'   /* FreeSurface equation timestep ( s ) */')
367        CALL WRITE_0D_R8( deltaTtracer, INDEX_NONE,'deltatTtracer =',        CALL WRITE_1D_R8( dTtracerLev, Nr, INDEX_K, 'dTtracerLev =',
368       &'   /* Tracer equation timestep ( s ) */')       &'   /* Tracer equation timestep ( s ) */')
369        CALL WRITE_0D_R8( deltaTClock, INDEX_NONE,'deltatTClock  =',        CALL WRITE_0D_R8( deltaTClock, INDEX_NONE,'deltatTClock  =',
370       &'   /* Model clock timestep ( s ) */')       &'   /* Model clock timestep ( s ) */')
# Line 325  CEOP Line 388  CEOP
388       &'   /* Permanent restart/checkpoint file interval ( s ). */')       &'   /* Permanent restart/checkpoint file interval ( s ). */')
389        CALL WRITE_0D_R8( chkPtFreq, INDEX_NONE,'chkPtFreq =',        CALL WRITE_0D_R8( chkPtFreq, INDEX_NONE,'chkPtFreq =',
390       &'   /* Rolling restart/checkpoint file interval ( s ). */')       &'   /* Rolling restart/checkpoint file interval ( s ). */')
391          CALL WRITE_0D_L(pickup_write_mdsio,INDEX_NONE,
392         &     'pickup_write_mdsio =', '   /* Model IO flag. */')
393          CALL WRITE_0D_L(pickup_read_mdsio,INDEX_NONE,
394         &     'pickup_read_mdsio =', '   /* Model IO flag. */')
395    #ifdef ALLOW_MNC
396          CALL WRITE_0D_L(pickup_write_mnc,INDEX_NONE,
397         &     'pickup_write_mnc =', '   /* Model IO flag. */')
398          CALL WRITE_0D_L(pickup_read_mnc,INDEX_NONE,
399         &     'pickup_read_mnc =', '   /* Model IO flag. */')
400    #endif
401          CALL WRITE_0D_L(pickup_write_immed,INDEX_NONE,
402         &     'pickup_write_immed =','   /* Model IO flag. */')
403        CALL WRITE_0D_R8( dumpFreq, INDEX_NONE,'dumpFreq =',        CALL WRITE_0D_R8( dumpFreq, INDEX_NONE,'dumpFreq =',
404       &'   /* Model state write out interval ( s ). */')       &'   /* Model state write out interval ( s ). */')
405          CALL WRITE_0D_L(snapshot_mdsio,INDEX_NONE,
406         &     'snapshot_mdsio =', '   /* Model IO flag. */')
407    #ifdef ALLOW_MNC
408          CALL WRITE_0D_L(snapshot_mnc,INDEX_NONE,
409         &     'snapshot_mnc =', '   /* Model IO flag. */')
410    #endif
411          CALL WRITE_0D_R8( monitorFreq, INDEX_NONE,'monitorFreq =',
412         &'   /* Monitor output interval ( s ). */')
413          CALL WRITE_0D_L(monitor_stdio,INDEX_NONE,
414         &     'monitor_stdio =', '   /* Model IO flag. */')
415    #ifdef ALLOW_MNC
416          CALL WRITE_0D_L(monitor_mnc,INDEX_NONE,
417         &     'monitor_mnc =', '   /* Model IO flag. */')
418    #endif
419          CALL WRITE_0D_R8( externForcingPeriod, INDEX_NONE,
420         &   'externForcingPeriod =', '   /* forcing period (s) */')
421          CALL WRITE_0D_R8( externForcingCycle, INDEX_NONE,
422         &   'externForcingCycle =', '   /* period of the cyle (s). */')
423          CALL WRITE_0D_R8( tauThetaClimRelax, INDEX_NONE,
424         &   'tauThetaClimRelax =', '   /* relaxation time scale (s) */')
425          CALL WRITE_0D_R8( tauSaltClimRelax, INDEX_NONE,
426         &   'tauSaltClimRelax =',  '   /* relaxation time scale (s) */')
427          CALL WRITE_0D_R8( latBandClimRelax, INDEX_NONE,
428         &   'latBandClimRelax =', '   /* max. Lat. where relaxation */')
429        WRITE(msgBuf,'(A)') '//  '        WRITE(msgBuf,'(A)') '//  '
430        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
431       &                    SQUEEZE_RIGHT , 1)       &                    SQUEEZE_RIGHT , 1)
# Line 344  CEOP Line 442  CEOP
442        CALL WRITE_0D_L( usingSphericalPolarGrid, INDEX_NONE,        CALL WRITE_0D_L( usingSphericalPolarGrid, INDEX_NONE,
443       & 'usingSphericalPolarGrid =',       & 'usingSphericalPolarGrid =',
444       &'   /* Spherical coordinates flag ( True / False ) */')       &'   /* Spherical coordinates flag ( True / False ) */')
445          CALL WRITE_0D_L( usingCylindricalGrid, INDEX_NONE,
446         & 'usingCylindricalGrid =',
447         &'   /* Spherical coordinates flag ( True / False ) */')
448        CALL WRITE_0D_L( groundAtK1, INDEX_NONE, 'groundAtK1 =',        CALL WRITE_0D_L( groundAtK1, INDEX_NONE, 'groundAtK1 =',
449       &'   /* Lower Boundary (ground) at the surface(k=1) ( T / F ) */')       &'   /* Lower Boundary (ground) at the surface(k=1) ( T / F ) */')
450        CALL WRITE_0D_R8( Ro_SeaLevel, INDEX_NONE,'Ro_SeaLevel =',        CALL WRITE_0D_R8( Ro_SeaLevel, INDEX_NONE,'Ro_SeaLevel =',
Legend:
Removed from v.1.41  
changed lines
  Added in v.1.69
  ViewVC Help
Powered by ViewVC 1.1.22