model/src/config_summary.F

C $Header: /u/gcmpack/models/MITgcmUV/model/src/config_summary.F,v 1.11 1998/06/10 17:05:59 cnh Exp $

#include "CPP_EEOPTIONS.h"

CStartOfInterface
      SUBROUTINE CONFIG_SUMMARY( myThid )
C     /==========================================================\
C     | SUBROUTINE CONFIG_SUMMARY                                |
C     | o Summarize model prognostic variables.                  |
C     |==========================================================|
C     | This routine writes a tabulated summary of the model     |
C     | configuration.                                           |
C     | Note                                                     |
C     |  1. Under multi-process parallelism the summary          |
C     |     is only given for the per-process data.              |
C     |  2. Under multi-threading the summary is produced by     |
C     |     the master thread. This threads reads data managed by|
C     |     other threads.                                       |
C     \==========================================================/

C     === Global variables ===
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "GRID.h"
#include "DYNVARS.h"

C     == Routine arguments ==
C     myThid -  Number of this instance of CONFIG_SUMMARY
      INTEGER myThid
CEndOfInterface

C     == Local variables ==
      CHARACTER*(MAX_LEN_MBUF) msgBuf
      INTEGER                  I,J,K
      INTEGER                  bi, bj
      REAL                     xcoord(Nx)
      REAL                     ycoord(Ny)
      REAL                     zcoord(Nz)


      _BARRIER
      _BEGIN_MASTER(myThid)

      WRITE(msgBuf,'(A)')
     &'// ======================================================='
      CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)
      WRITE(msgBuf,'(A)') '// Model configuration'
      CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)
      WRITE(msgBuf,'(A)')
     &'// ======================================================='
      CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
     &  SQUEEZE_RIGHT , 1)

      WRITE(msgBuf,'(A)') '//  '
      CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)
      WRITE(msgBuf,'(A)') '// "Physical" paramters ( PARM01 in namelist ) '
      CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)
      WRITE(msgBuf,'(A)') '//  '
      CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)
      CALL WRITE_1D_R8( tRef, Nz, INDEX_K,'tRef =',
     &'   /* Reference temperature profile ( oC or oK ) */')
      CALL WRITE_1D_R8( sRef, Nz, INDEX_K,'sRef =',
     &'   /* Reference salinity profile ( ppt ) */')
      CALL WRITE_1D_R8( viscAh, 1, INDEX_NONE,'viscAh =',
     &'   /* Lateral eddy viscosity ( m^2/s ) */')
      CALL WRITE_1D_R8( viscAz, 1, INDEX_NONE,'viscAz =',
     &'   /* Vertical eddy viscosity ( m^2/s ) */')
      CALL WRITE_1D_R8( diffKhT, 1, INDEX_NONE,'diffKhT =',
     &'   /* Laplacian diffusion of heat laterally ( m^2/s ) */')
      CALL WRITE_1D_R8( diffKzT, 1, INDEX_NONE,'diffKzT =',
     &'   /* Laplacian diffusion of heat vertically ( m^2/s ) */')
      CALL WRITE_1D_R8( diffK4T, 1, INDEX_NONE,'diffK4T =',
     &'   /* Bihaarmonic diffusion of heat laterally ( m^4/s ) */')
      CALL WRITE_1D_R8( diffKhS, 1, INDEX_NONE,'diffKhS =',
     &'   /* Laplacian diffusion of salt laterally ( m^2/s ) */')
      CALL WRITE_1D_R8( diffKzS, 1, INDEX_NONE,'diffKzS =',
     &'   /* Laplacian diffusion of salt vertically ( m^2/s ) */')
      CALL WRITE_1D_R8( diffK4S, 1, INDEX_NONE,'diffK4S =',
     &'   /* Bihaarmonic diffusion of salt laterally ( m^4/s ) */')
      CALL WRITE_1D_R8( tAlpha,1, INDEX_NONE,'tAlpha =',
     &'   /* Linear EOS thermal expansion coefficient ( 1/degree ) */')
      CALL WRITE_1D_R8( sBeta, 1, INDEX_NONE,'sBeta =',
     &'   /* Linear EOS haline contraction coefficient ( 1/ppt ) */')
      CALL WRITE_1D_R8( rhonil,1, INDEX_NONE,'rhonil =',
     &'   /* Reference density ( kg/m^3 ) */')
      CALL WRITE_1D_R8( gravity,1, INDEX_NONE,'gravity =',
     &'   /* Gravitational acceleration ( m/s^2 ) */')
      CALL WRITE_1D_R8( gBaro,1, INDEX_NONE,'gBaro =',
     &'   /* Barotropic gravity ( m/s^2 ) */')
      CALL WRITE_1D_R8( f0,1, INDEX_NONE,'f0 =',
     &'   /* Reference coriolis parameter ( 1/s ) */')
      CALL WRITE_1D_R8( beta,1, INDEX_NONE,'beta =',
     &'   /* Beta ( 1/(m.s) ) */')
      CALL WRITE_1D_R8( freeSurfFac,1, INDEX_NONE,'freeSurfFac =',
     &'   /* Implcit free surface factor */')
      CALL WRITE_1D_L( implicitFreeSurface,1, INDEX_NONE,
     &                 'implicitFreeSurface =',
     &'   /* Implicit free surface on/off flag */')
      CALL WRITE_1D_L( rigidLid,1, INDEX_NONE,
     &                 'rigidLid =',
     &'   /* Rigid lid on/off flag */')
      CALL WRITE_1D_L( momStepping,1, INDEX_NONE,
     & 'momStepping =', '  /* Momentum equation on/off flag */')
      CALL WRITE_1D_L( momAdvection,1, INDEX_NONE,
     & 'momAdvection =', '  /* Momentum advection on/off flag */')
      CALL WRITE_1D_L( momViscosity,1, INDEX_NONE,
     & 'momViscosity =', '  /* Momentum viscosity on/off flag */')
      CALL WRITE_1D_L( useCoriolis,1, INDEX_NONE,
     & 'useCoriolis =', '  /* Coriolis on/off flag */')
      CALL WRITE_1D_L( momForcing,1, INDEX_NONE,
     & 'momForcing =', '  /* Momentum forcing on/off flag */')
      CALL WRITE_1D_L( momPressureForcing,1, INDEX_NONE,
     & 'momPressureForcing =', '  /* Momentum pressure term on/off flag */')
      CALL WRITE_1D_L( tempStepping,1, INDEX_NONE,
     & 'tempStepping =', '  /* Temperature equation on/off flag */')
      CALL WRITE_1D_R8( GMMaxSlope,1, INDEX_NONE,'GMMaxSlope =',
     &'   /* Max. slope allowed in GM/Redi tensor */')
      CALL WRITE_1D_R8( GMLength,1, INDEX_NONE,'GMLength =',
     &'   /* Length to use in Visbeck et al. formula for K (m) */')
      CALL WRITE_1D_R8( GMAlpha,1, INDEX_NONE,'GMAlpha =',
     &'   /* alpha to use in Visbeck et al. formula for K */')
      CALL WRITE_1D_R8( GMdepth,1, INDEX_NONE,'GMdepth =',
     &'   /* Depth to integrate for Visbeck et. al Richardson # (m) */')
      CALL WRITE_1D_R8( GMkbackground,1, INDEX_NONE,'GMkbackground =',
     &'   /* background value of GM/Redi coefficient m^2/s */')
      WRITE(msgBuf,'(A)') '//  '
      CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)

      WRITE(msgBuf,'(A)') '// Elliptic solver(s) paramters ( PARM02 in namelist ) '
      CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)
      WRITE(msgBuf,'(A)') '//  '
      CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)
      CALL WRITE_1D_I( cg2dMaxIters,1, INDEX_NONE,'cg2dMaxIters =',
     &'   /* Upper limit on 2d con. grad iterations  */')
      CALL WRITE_1D_I( cg2dChkResFreq,1, INDEX_NONE,'cg2dChkResFreq =',
     &'   /* 2d con. grad convergence test frequency */')
      CALL WRITE_1D_R8( cg2dTargetResidual,1, INDEX_NONE,'cg2dTargetResidual =',
     &'   /* 2d con. grad target residual  */')

      WRITE(msgBuf,'(A)') '//  '
      CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)
      WRITE(msgBuf,'(A)') '// Time stepping paramters ( PARM03 in namelist ) '
      CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)
      WRITE(msgBuf,'(A)') '//  '
      CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)
      CALL WRITE_1D_I( nIter0,1, INDEX_NONE,'nIter0 =',
     &'   /* Base timestep number  */')
      CALL WRITE_1D_I( nTimeSteps,1, INDEX_NONE,'nTimeSteps =',
     &'   /* Number of timesteps */')
      CALL WRITE_1D_R8( deltaTmom,1, INDEX_NONE,'deltatTmom =',
     &'   /* Momentum equation timestep ( s ) */')
      CALL WRITE_1D_R8( deltaTtracer,1, INDEX_NONE,'deltatTtracer =',
     &'   /* Tracer equation timestep ( s ) */')
      CALL WRITE_1D_R8( deltaTClock ,1, INDEX_NONE,'deltatTClock  =',
     &'   /* Model clock timestep ( s ) */')
      CALL WRITE_1D_R8( cAdjFreq,1, INDEX_NONE,'cAdjFreq =',
     &'   /* Convective adjustment interval ( s ) */')
      CALL WRITE_1D_R8( abEps,1, INDEX_NONE,'abEps =',
     &'   /* Adams-Bashforth stabilizing weight */')
      CALL WRITE_1D_R8( tauCD,1, INDEX_NONE,'tauCD =',
     &'   /* CD coupling time-scale ( s ) */')
      CALL WRITE_1D_R8( rCD,1, INDEX_NONE,'rCD =',
     &'   /* Normalised CD coupling parameter */')
      CALL WRITE_1D_R8( startTime,1, INDEX_NONE,'startTime =',
     &'   /* Run start time ( s ). */')
      CALL WRITE_1D_R8( endTime,1, INDEX_NONE,'endTime =',
     &'   /* Integration ending time ( s ). */')
      CALL WRITE_1D_R8( pChkPtFreq,1, INDEX_NONE,'pChkPtFreq =',
     &'   /* Permanent restart/checkpoint file interval ( s ). */')
      CALL WRITE_1D_R8( chkPtFreq,1, INDEX_NONE,'chkPtFreq =',
     &'   /* Rolling restart/checkpoint file interval ( s ). */')
      CALL WRITE_1D_R8( dumpFreq,1, INDEX_NONE,'dumpFreq =',
     &'   /* Model state write out interval ( s ). */')

      WRITE(msgBuf,'(A)') '//  '
      CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)
      WRITE(msgBuf,'(A)') '// Gridding paramters ( PARM04 in namelist ) '
      CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)
      WRITE(msgBuf,'(A)') '//  '
      CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)
      CALL WRITE_1D_L( usingCartesianGrid,1, INDEX_NONE,'usingCartesianGrid =',
     &'   /* Cartesian coordinates flag ( True / False ) */')
      CALL WRITE_1D_L( usingSphericalPolarGrid,1, INDEX_NONE,'usingSphericalPolarGrid =',
     &'   /* Spherical coordinates flag ( True / False ) */')
      CALL WRITE_1D_R8( delZ,Nz, INDEX_K,'delZ = ',
     &'   /* W spacing ( m ) */')
      CALL WRITE_1D_R8( delX, Nx, INDEX_I,'delX = ',
     &'   /* U spacing ( m - cartesian, degrees - spherical ) */')
      CALL WRITE_1D_R8( delY, Ny, INDEX_J,'delY = ',
     &'   /* V spacing ( m - cartesian, degrees - spherical ) */')
      CALL WRITE_1D_R8( phiMin, 1, INDEX_NONE,'phiMin = ',
     &'   /* Southern boundary ( ignored - cartesian, degrees - spherical ) */')
      CALL WRITE_1D_R8( thetaMin, 1, INDEX_NONE,'thetaMin = ',
     &'   /* Western boundary ( ignored - cartesian, degrees - spherical ) */')
      CALL WRITE_1D_R8( rSphere, 1, INDEX_NONE,'rSphere = ',
     &'   /* Radius ( ignored - cartesian, m - spherical ) */')
      DO bi=1,nSx
       DO I=1,sNx
        xcoord((bi-1)*sNx+I) = xc(I,1,bi,1)
       ENDDO
      ENDDO
      CALL WRITE_1D_R8( xcoord, sNx*nSx, INDEX_I,'xcoord = ',
     &'   /* P-point X coordinate (  m - cartesian, degrees - spherical ) */')
      DO bj=1,nSy
       DO J=1,sNy
        ycoord((bj-1)*sNy+J) = yc(1,J,1,bj)
       ENDDO
      ENDDO
      CALL WRITE_1D_R8( ycoord, sNy*nSy, INDEX_J,'ycoord = ',
     &'   /* P-point Y coordinate (  m - cartesian, degrees - spherical ) */')
      DO K=1,Nz
       zcoord(K) = zc(K)
      ENDDO
      CALL WRITE_1D_R8( zcoord, Nz, INDEX_K,'zcoord = ',
     &'   /* P-point Z coordinate (  m ) */')


      WRITE(msgBuf,'(A)') ' '
      CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
     &  SQUEEZE_RIGHT , 1)

      _END_MASTER(myThid)
      _BARRIER


      RETURN
  100 FORMAT(A,
     &' '
     &)
      END

1	cnh	1.12	C $Header: /u/gcmpack/models/MITgcmUV/model/src/config_summary.F,v 1.11 1998/06/10 17:05:59 cnh Exp $
2	cnh	1.1
3			#include "CPP_EEOPTIONS.h"
4
5			CStartOfInterface
6			SUBROUTINE CONFIG_SUMMARY( myThid )
7			C /==========================================================\
8			C \| SUBROUTINE CONFIG_SUMMARY \|
9			C \| o Summarize model prognostic variables. \|
10			C \|==========================================================\|
11			C \| This routine writes a tabulated summary of the model \|
12			C \| configuration. \|
13			C \| Note \|
14			C \| 1. Under multi-process parallelism the summary \|
15			C \| is only given for the per-process data. \|
16			C \| 2. Under multi-threading the summary is produced by \|
17			C \| the master thread. This threads reads data managed by\|
18			C \| other threads. \|
19			C \==========================================================/
20
21			C === Global variables ===
22			#include "SIZE.h"
23			#include "EEPARAMS.h"
24			#include "PARAMS.h"
25			#include "GRID.h"
26			#include "DYNVARS.h"
27
28			C == Routine arguments ==
29			C myThid - Number of this instance of CONFIG_SUMMARY
30			INTEGER myThid
31			CEndOfInterface
32
33			C == Local variables ==
34			CHARACTER*(MAX_LEN_MBUF) msgBuf
35	cnh	1.5	INTEGER I,J,K
36	cnh	1.6	INTEGER bi, bj
37			REAL xcoord(Nx)
38			REAL ycoord(Ny)
39			REAL zcoord(Nz)
40	cnh	1.5
41	cnh	1.1
42			_BARRIER
43	cnh	1.5	_BEGIN_MASTER(myThid)
44	cnh	1.1
45			WRITE(msgBuf,'(A)')
46			&'// ======================================================='
47			CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)
48			WRITE(msgBuf,'(A)') '// Model configuration'
49			CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)
50			WRITE(msgBuf,'(A)')
51			&'// ======================================================='
52			CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
53			& SQUEEZE_RIGHT , 1)
54	cnh	1.5
55	cnh	1.6	WRITE(msgBuf,'(A)') '// '
56			CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)
57			WRITE(msgBuf,'(A)') '// "Physical" paramters ( PARM01 in namelist ) '
58			CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)
59			WRITE(msgBuf,'(A)') '// '
60			CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)
61	cnh	1.5	CALL WRITE_1D_R8( tRef, Nz, INDEX_K,'tRef =',
62			&' /* Reference temperature profile ( oC or oK ) */')
63	cnh	1.6	CALL WRITE_1D_R8( sRef, Nz, INDEX_K,'sRef =',
64			&' /* Reference salinity profile ( ppt ) */')
65	cnh	1.5	CALL WRITE_1D_R8( viscAh, 1, INDEX_NONE,'viscAh =',
66			&' /* Lateral eddy viscosity ( m^2/s ) */')
67			CALL WRITE_1D_R8( viscAz, 1, INDEX_NONE,'viscAz =',
68			&' /* Vertical eddy viscosity ( m^2/s ) */')
69			CALL WRITE_1D_R8( diffKhT, 1, INDEX_NONE,'diffKhT =',
70			&' /* Laplacian diffusion of heat laterally ( m^2/s ) */')
71			CALL WRITE_1D_R8( diffKzT, 1, INDEX_NONE,'diffKzT =',
72			&' /* Laplacian diffusion of heat vertically ( m^2/s ) */')
73			CALL WRITE_1D_R8( diffK4T, 1, INDEX_NONE,'diffK4T =',
74			&' /* Bihaarmonic diffusion of heat laterally ( m^4/s ) */')
75			CALL WRITE_1D_R8( diffKhS, 1, INDEX_NONE,'diffKhS =',
76			&' /* Laplacian diffusion of salt laterally ( m^2/s ) */')
77			CALL WRITE_1D_R8( diffKzS, 1, INDEX_NONE,'diffKzS =',
78			&' /* Laplacian diffusion of salt vertically ( m^2/s ) */')
79			CALL WRITE_1D_R8( diffK4S, 1, INDEX_NONE,'diffK4S =',
80			&' /* Bihaarmonic diffusion of salt laterally ( m^4/s ) */')
81	cnh	1.6	CALL WRITE_1D_R8( tAlpha,1, INDEX_NONE,'tAlpha =',
82			&' /* Linear EOS thermal expansion coefficient ( 1/degree ) */')
83			CALL WRITE_1D_R8( sBeta, 1, INDEX_NONE,'sBeta =',
84			&' /* Linear EOS haline contraction coefficient ( 1/ppt ) */')
85			CALL WRITE_1D_R8( rhonil,1, INDEX_NONE,'rhonil =',
86			&' /* Reference density ( kg/m^3 ) */')
87			CALL WRITE_1D_R8( gravity,1, INDEX_NONE,'gravity =',
88			&' /* Gravitational acceleration ( m/s^2 ) */')
89	cnh	1.8	CALL WRITE_1D_R8( gBaro,1, INDEX_NONE,'gBaro =',
90			&' /* Barotropic gravity ( m/s^2 ) */')
91	cnh	1.6	CALL WRITE_1D_R8( f0,1, INDEX_NONE,'f0 =',
92			&' /* Reference coriolis parameter ( 1/s ) */')
93			CALL WRITE_1D_R8( beta,1, INDEX_NONE,'beta =',
94			&' /* Beta ( 1/(m.s) ) */')
95	cnh	1.8	CALL WRITE_1D_R8( freeSurfFac,1, INDEX_NONE,'freeSurfFac =',
96			&' /* Implcit free surface factor */')
97			CALL WRITE_1D_L( implicitFreeSurface,1, INDEX_NONE,
98			& 'implicitFreeSurface =',
99			&' /* Implicit free surface on/off flag */')
100			CALL WRITE_1D_L( rigidLid,1, INDEX_NONE,
101			& 'rigidLid =',
102			&' /* Rigid lid on/off flag */')
103	cnh	1.10	CALL WRITE_1D_L( momStepping,1, INDEX_NONE,
104			& 'momStepping =', ' /* Momentum equation on/off flag */')
105	cnh	1.9	CALL WRITE_1D_L( momAdvection,1, INDEX_NONE,
106	cnh	1.10	& 'momAdvection =', ' /* Momentum advection on/off flag */')
107	cnh	1.9	CALL WRITE_1D_L( momViscosity,1, INDEX_NONE,
108			& 'momViscosity =', ' /* Momentum viscosity on/off flag */')
109			CALL WRITE_1D_L( useCoriolis,1, INDEX_NONE,
110			& 'useCoriolis =', ' /* Coriolis on/off flag */')
111			CALL WRITE_1D_L( momForcing,1, INDEX_NONE,
112			& 'momForcing =', ' /* Momentum forcing on/off flag */')
113			CALL WRITE_1D_L( momPressureForcing,1, INDEX_NONE,
114			& 'momPressureForcing =', ' /* Momentum pressure term on/off flag */')
115	cnh	1.10	CALL WRITE_1D_L( tempStepping,1, INDEX_NONE,
116			& 'tempStepping =', ' /* Temperature equation on/off flag */')
117	cnh	1.8	CALL WRITE_1D_R8( GMMaxSlope,1, INDEX_NONE,'GMMaxSlope =',
118			&' /* Max. slope allowed in GM/Redi tensor */')
119			CALL WRITE_1D_R8( GMLength,1, INDEX_NONE,'GMLength =',
120			&' /* Length to use in Visbeck et al. formula for K (m) */')
121			CALL WRITE_1D_R8( GMAlpha,1, INDEX_NONE,'GMAlpha =',
122			&' /* alpha to use in Visbeck et al. formula for K */')
123			CALL WRITE_1D_R8( GMdepth,1, INDEX_NONE,'GMdepth =',
124			&' /* Depth to integrate for Visbeck et. al Richardson # (m) */')
125			CALL WRITE_1D_R8( GMkbackground,1, INDEX_NONE,'GMkbackground =',
126			&' /* background value of GM/Redi coefficient m^2/s */')
127	cnh	1.6	WRITE(msgBuf,'(A)') '// '
128			CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)
129	cnh	1.9
130	cnh	1.6	WRITE(msgBuf,'(A)') '// Elliptic solver(s) paramters ( PARM02 in namelist ) '
131			CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)
132			WRITE(msgBuf,'(A)') '// '
133			CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)
134			CALL WRITE_1D_I( cg2dMaxIters,1, INDEX_NONE,'cg2dMaxIters =',
135			&' /* Upper limit on 2d con. grad iterations */')
136			CALL WRITE_1D_I( cg2dChkResFreq,1, INDEX_NONE,'cg2dChkResFreq =',
137			&' /* 2d con. grad convergence test frequency */')
138			CALL WRITE_1D_R8( cg2dTargetResidual,1, INDEX_NONE,'cg2dTargetResidual =',
139			&' /* 2d con. grad target residual */')
140
141			WRITE(msgBuf,'(A)') '// '
142			CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)
143			WRITE(msgBuf,'(A)') '// Time stepping paramters ( PARM03 in namelist ) '
144			CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)
145			WRITE(msgBuf,'(A)') '// '
146			CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)
147			CALL WRITE_1D_I( nIter0,1, INDEX_NONE,'nIter0 =',
148			&' /* Base timestep number */')
149			CALL WRITE_1D_I( nTimeSteps,1, INDEX_NONE,'nTimeSteps =',
150			&' /* Number of timesteps */')
151			CALL WRITE_1D_R8( deltaTmom,1, INDEX_NONE,'deltatTmom =',
152			&' /* Momentum equation timestep ( s ) */')
153			CALL WRITE_1D_R8( deltaTtracer,1, INDEX_NONE,'deltatTtracer =',
154			&' /* Tracer equation timestep ( s ) */')
155	cnh	1.12	CALL WRITE_1D_R8( deltaTClock ,1, INDEX_NONE,'deltatTClock =',
156			&' /* Model clock timestep ( s ) */')
157	cnh	1.9	CALL WRITE_1D_R8( cAdjFreq,1, INDEX_NONE,'cAdjFreq =',
158			&' /* Convective adjustment interval ( s ) */')
159	cnh	1.6	CALL WRITE_1D_R8( abEps,1, INDEX_NONE,'abEps =',
160			&' /* Adams-Bashforth stabilizing weight */')
161			CALL WRITE_1D_R8( tauCD,1, INDEX_NONE,'tauCD =',
162			&' /* CD coupling time-scale ( s ) */')
163			CALL WRITE_1D_R8( rCD,1, INDEX_NONE,'rCD =',
164			&' /* Normalised CD coupling parameter */')
165			CALL WRITE_1D_R8( startTime,1, INDEX_NONE,'startTime =',
166			&' /* Run start time ( s ). */')
167			CALL WRITE_1D_R8( endTime,1, INDEX_NONE,'endTime =',
168			&' /* Integration ending time ( s ). */')
169	cnh	1.7	CALL WRITE_1D_R8( pChkPtFreq,1, INDEX_NONE,'pChkPtFreq =',
170			&' /* Permanent restart/checkpoint file interval ( s ). */')
171	cnh	1.6	CALL WRITE_1D_R8( chkPtFreq,1, INDEX_NONE,'chkPtFreq =',
172	cnh	1.7	&' /* Rolling restart/checkpoint file interval ( s ). */')
173	cnh	1.6	CALL WRITE_1D_R8( dumpFreq,1, INDEX_NONE,'dumpFreq =',
174			&' /* Model state write out interval ( s ). */')
175
176			WRITE(msgBuf,'(A)') '// '
177			CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)
178			WRITE(msgBuf,'(A)') '// Gridding paramters ( PARM04 in namelist ) '
179			CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)
180			WRITE(msgBuf,'(A)') '// '
181			CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)
182			CALL WRITE_1D_L( usingCartesianGrid,1, INDEX_NONE,'usingCartesianGrid =',
183			&' /* Cartesian coordinates flag ( True / False ) */')
184			CALL WRITE_1D_L( usingSphericalPolarGrid,1, INDEX_NONE,'usingSphericalPolarGrid =',
185			&' /* Spherical coordinates flag ( True / False ) */')
186			CALL WRITE_1D_R8( delZ,Nz, INDEX_K,'delZ = ',
187			&' /* W spacing ( m ) */')
188			CALL WRITE_1D_R8( delX, Nx, INDEX_I,'delX = ',
189			&' /* U spacing ( m - cartesian, degrees - spherical ) */')
190			CALL WRITE_1D_R8( delY, Ny, INDEX_J,'delY = ',
191			&' /* V spacing ( m - cartesian, degrees - spherical ) */')
192			CALL WRITE_1D_R8( phiMin, 1, INDEX_NONE,'phiMin = ',
193			&' /* Southern boundary ( ignored - cartesian, degrees - spherical ) */')
194			CALL WRITE_1D_R8( thetaMin, 1, INDEX_NONE,'thetaMin = ',
195			&' /* Western boundary ( ignored - cartesian, degrees - spherical ) */')
196			CALL WRITE_1D_R8( rSphere, 1, INDEX_NONE,'rSphere = ',
197			&' /* Radius ( ignored - cartesian, m - spherical ) */')
198			DO bi=1,nSx
199			DO I=1,sNx
200			xcoord((bi-1)*sNx+I) = xc(I,1,bi,1)
201			ENDDO
202			ENDDO
203	cnh	1.11	CALL WRITE_1D_R8( xcoord, sNx*nSx, INDEX_I,'xcoord = ',
204	cnh	1.6	&' /* P-point X coordinate ( m - cartesian, degrees - spherical ) */')
205			DO bj=1,nSy
206			DO J=1,sNy
207			ycoord((bj-1)*sNy+J) = yc(1,J,1,bj)
208			ENDDO
209			ENDDO
210	cnh	1.11	CALL WRITE_1D_R8( ycoord, sNy*nSy, INDEX_J,'ycoord = ',
211	cnh	1.6	&' /* P-point Y coordinate ( m - cartesian, degrees - spherical ) */')
212			DO K=1,Nz
213			zcoord(K) = zc(K)
214			ENDDO
215			CALL WRITE_1D_R8( zcoord, Nz, INDEX_K,'zcoord = ',
216			&' /* P-point Z coordinate ( m ) */')
217
218	cnh	1.5
219
220	cnh	1.1	WRITE(msgBuf,'(A)') ' '
221			CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
222			& SQUEEZE_RIGHT , 1)
223	cnh	1.5
224	cnh	1.1	_END_MASTER(myThid)
225			_BARRIER
226
227
228			RETURN
229			100 FORMAT(A,
230	cnh	1.4	&' '
231	cnh	1.1	&)
232			END
233