/[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.1 by cnh, Wed Apr 22 19:15:30 1998 UTC revision 1.16 by cnh, Sat Sep 5 17:52:13 1998 UTC
# Line 1  Line 1 
1  C $Id$  C $Header$
2    
3  #include "CPP_EEOPTIONS.h"  #include "CPP_EEOPTIONS.h"
4    
# Line 32  CEndOfInterface Line 32  CEndOfInterface
32    
33  C     == Local variables ==  C     == Local variables ==
34        CHARACTER*(MAX_LEN_MBUF) msgBuf        CHARACTER*(MAX_LEN_MBUF) msgBuf
35          INTEGER                  I,J,K
36          INTEGER                  bi, bj
37          REAL                     xcoord(Nx)
38          REAL                     ycoord(Ny)
39          REAL                     rcoord(Nr)
40    
       _BARRIER  
41    
42          _BARRIER
43        _BEGIN_MASTER(myThid)        _BEGIN_MASTER(myThid)
44        WRITE(msgBuf,100) '// '  
       CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)  
45        WRITE(msgBuf,'(A)')        WRITE(msgBuf,'(A)')
46       &'// ======================================================='       &'// ======================================================='
47        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)
# Line 47  C     == Local variables == Line 51  C     == Local variables ==
51       &'// ======================================================='       &'// ======================================================='
52        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
53       &  SQUEEZE_RIGHT , 1)       &  SQUEEZE_RIGHT , 1)
54    
55          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          CALL WRITE_1D_R8( tRef, Nr, INDEX_K,'tRef =',
62         &'   /* Reference temperature profile ( oC or oK ) */')
63          CALL WRITE_1D_R8( sRef, Nr, INDEX_K,'sRef =',
64         &'   /* Reference salinity profile ( ppt ) */')
65          CALL WRITE_1D_R8( viscAh, 1, INDEX_NONE,'viscAh =',
66         &'   /* Lateral eddy viscosity ( m^2/s ) */')
67          CALL WRITE_1D_R8( viscA4, 1, INDEX_NONE,'viscAh =',
68         &'   /* Lateral biharmonic viscosity ( m^4/s ) */')
69          IF ( viscAz .NE. UNSET_RL ) THEN
70           CALL WRITE_1D_R8( viscAz, 1, INDEX_NONE,'viscAz =',
71         & '   /* Vertical eddy viscosity ( m^2/s ) */')
72          ENDIF
73          IF ( viscAp .NE. UNSET_RL ) THEN
74           CALL WRITE_1D_R8( viscAp, 1, INDEX_NONE,'viscAp =',
75         & '   /* Vertical eddy viscosity ( Pa^2/s ) */')
76          ENDIF
77          CALL WRITE_1D_R8( viscAr, 1, INDEX_NONE,'viscAr =',
78         &'   /* Vertical eddy viscosity ( units of r^2/s ) */')
79          CALL WRITE_1D_R8( diffKhT, 1, INDEX_NONE,'diffKhT =',
80         &'   /* Laplacian diffusion of heat laterally ( m^2/s ) */')
81          CALL WRITE_1D_R8( diffKzT, 1, INDEX_NONE,'diffKzT =',
82         &'   /* Laplacian diffusion of heat vertically ( m^2/s ) */')
83          CALL WRITE_1D_R8( diffK4T, 1, INDEX_NONE,'diffK4T =',
84         &'   /* Bihaarmonic diffusion of heat laterally ( m^4/s ) */')
85          CALL WRITE_1D_R8( diffKhS, 1, INDEX_NONE,'diffKhS =',
86         &'   /* Laplacian diffusion of salt laterally ( m^2/s ) */')
87          CALL WRITE_1D_R8( diffKzS, 1, INDEX_NONE,'diffKzS =',
88         &'   /* Laplacian diffusion of salt vertically ( m^2/s ) */')
89          CALL WRITE_1D_R8( diffK4S, 1, INDEX_NONE,'diffK4S =',
90         &'   /* Bihaarmonic diffusion of salt laterally ( m^4/s ) */')
91          CALL WRITE_1D_R8( tAlpha,1, INDEX_NONE,'tAlpha =',
92         &'   /* Linear EOS thermal expansion coefficient ( 1/degree ) */')
93          CALL WRITE_1D_R8( sBeta, 1, INDEX_NONE,'sBeta =',
94         &'   /* Linear EOS haline contraction coefficient ( 1/ppt ) */')
95          IF ( eosType .EQ. 'POLY3' ) THEN
96            WRITE(msgBuf,'(A)') '// Polynomial EQS parameters ( from POLY3.COEFFS ) '
97            DO K = 1, Nr
98             WRITE(msgBuf,'(I3,13F8.3)')
99         &   K,eosRefT(K),eosRefS(K),eosSig0(K), (eosC(I,K),I=1,9)
100             CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)
101            ENDDO
102          ENDIF
103          CALL WRITE_1D_R8( rhonil,1, INDEX_NONE,'rhonil =',
104         &'   /* Reference density ( kg/m^3 ) */')
105          CALL WRITE_1D_R8( gravity,1, INDEX_NONE,'gravity =',
106         &'   /* Gravitational acceleration ( m/s^2 ) */')
107          CALL WRITE_1D_R8( gBaro,1, INDEX_NONE,'gBaro =',
108         &'   /* Barotropic gravity ( m/s^2 ) */')
109          CALL WRITE_1D_R8( f0,1, INDEX_NONE,'f0 =',
110         &'   /* Reference coriolis parameter ( 1/s ) */')
111          CALL WRITE_1D_R8( beta,1, INDEX_NONE,'beta =',
112         &'   /* Beta ( 1/(m.s) ) */')
113          CALL WRITE_1D_R8( freeSurfFac,1, INDEX_NONE,'freeSurfFac =',
114         &'   /* Implcit free surface factor */')
115          CALL WRITE_1D_L( implicitFreeSurface,1, INDEX_NONE,
116         &                 'implicitFreeSurface =',
117         &'   /* Implicit free surface on/off flag */')
118          CALL WRITE_1D_L( rigidLid,1, INDEX_NONE,
119         &                 'rigidLid =',
120         &'   /* Rigid lid on/off flag */')
121          CALL WRITE_1D_L( momStepping,1, INDEX_NONE,
122         & 'momStepping =', '  /* Momentum equation on/off flag */')
123          CALL WRITE_1D_L( momAdvection,1, INDEX_NONE,
124         & 'momAdvection =', '  /* Momentum advection on/off flag */')
125          CALL WRITE_1D_L( momViscosity,1, INDEX_NONE,
126         & 'momViscosity =', '  /* Momentum viscosity on/off flag */')
127          CALL WRITE_1D_L( useCoriolis,1, INDEX_NONE,
128         & 'useCoriolis =', '  /* Coriolis on/off flag */')
129          CALL WRITE_1D_L( momForcing,1, INDEX_NONE,
130         & 'momForcing =', '  /* Momentum forcing on/off flag */')
131          CALL WRITE_1D_L( momPressureForcing,1, INDEX_NONE,
132         & 'momPressureForcing =', '  /* Momentum pressure term on/off flag */')
133          CALL WRITE_1D_L( tempStepping,1, INDEX_NONE,
134         & 'tempStepping =', '  /* Temperature equation on/off flag */')
135          CALL WRITE_1D_R8( GMMaxSlope,1, INDEX_NONE,'GMMaxSlope =',
136         &'   /* Max. slope allowed in GM/Redi tensor */')
137          CALL WRITE_1D_R8( GMLength,1, INDEX_NONE,'GMLength =',
138         &'   /* Length to use in Visbeck et al. formula for K (m) */')
139          CALL WRITE_1D_R8( GMAlpha,1, INDEX_NONE,'GMAlpha =',
140         &'   /* alpha to use in Visbeck et al. formula for K */')
141          CALL WRITE_1D_R8( GMdepth,1, INDEX_NONE,'GMdepth =',
142         &'   /* Depth to integrate for Visbeck et. al Richardson # (m) */')
143          CALL WRITE_1D_R8( GMkbackground,1, INDEX_NONE,'GMkbackground =',
144         &'   /* background value of GM/Redi coefficient m^2/s */')
145          WRITE(msgBuf,'(A)') '//  '
146          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)
147    
148          WRITE(msgBuf,'(A)') '// Elliptic solver(s) paramters ( PARM02 in namelist ) '
149          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)
150          WRITE(msgBuf,'(A)') '//  '
151          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)
152          CALL WRITE_1D_I( cg2dMaxIters,1, INDEX_NONE,'cg2dMaxIters =',
153         &'   /* Upper limit on 2d con. grad iterations  */')
154          CALL WRITE_1D_I( cg2dChkResFreq,1, INDEX_NONE,'cg2dChkResFreq =',
155         &'   /* 2d con. grad convergence test frequency */')
156          CALL WRITE_1D_R8( cg2dTargetResidual,1, INDEX_NONE,'cg2dTargetResidual =',
157         &'   /* 2d con. grad target residual  */')
158    
159          WRITE(msgBuf,'(A)') '//  '
160          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)
161          WRITE(msgBuf,'(A)') '// Time stepping paramters ( PARM03 in namelist ) '
162          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)
163          WRITE(msgBuf,'(A)') '//  '
164          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)
165          CALL WRITE_1D_I( nIter0,1, INDEX_NONE,'nIter0 =',
166         &'   /* Base timestep number  */')
167          CALL WRITE_1D_I( nTimeSteps,1, INDEX_NONE,'nTimeSteps =',
168         &'   /* Number of timesteps */')
169          CALL WRITE_1D_R8( deltaTmom,1, INDEX_NONE,'deltatTmom =',
170         &'   /* Momentum equation timestep ( s ) */')
171          CALL WRITE_1D_R8( deltaTtracer,1, INDEX_NONE,'deltatTtracer =',
172         &'   /* Tracer equation timestep ( s ) */')
173          CALL WRITE_1D_R8( deltaTClock ,1, INDEX_NONE,'deltatTClock  =',
174         &'   /* Model clock timestep ( s ) */')
175          CALL WRITE_1D_R8( cAdjFreq,1, INDEX_NONE,'cAdjFreq =',
176         &'   /* Convective adjustment interval ( s ) */')
177          CALL WRITE_1D_R8( abEps,1, INDEX_NONE,'abEps =',
178         &'   /* Adams-Bashforth stabilizing weight */')
179          CALL WRITE_1D_R8( tauCD,1, INDEX_NONE,'tauCD =',
180         &'   /* CD coupling time-scale ( s ) */')
181          CALL WRITE_1D_R8( rCD,1, INDEX_NONE,'rCD =',
182         &'   /* Normalised CD coupling parameter */')
183          CALL WRITE_1D_R8( startTime,1, INDEX_NONE,'startTime =',
184         &'   /* Run start time ( s ). */')
185          CALL WRITE_1D_R8( endTime,1, INDEX_NONE,'endTime =',
186         &'   /* Integration ending time ( s ). */')
187          CALL WRITE_1D_R8( pChkPtFreq,1, INDEX_NONE,'pChkPtFreq =',
188         &'   /* Permanent restart/checkpoint file interval ( s ). */')
189          CALL WRITE_1D_R8( chkPtFreq,1, INDEX_NONE,'chkPtFreq =',
190         &'   /* Rolling restart/checkpoint file interval ( s ). */')
191          CALL WRITE_1D_R8( dumpFreq,1, INDEX_NONE,'dumpFreq =',
192         &'   /* Model state write out interval ( s ). */')
193    
194          WRITE(msgBuf,'(A)') '//  '
195          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)
196          WRITE(msgBuf,'(A)') '// Gridding paramters ( PARM04 in namelist ) '
197          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)
198          WRITE(msgBuf,'(A)') '//  '
199          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)
200          CALL WRITE_1D_L( usingCartesianGrid,1, INDEX_NONE,'usingCartesianGrid =',
201         &'   /* Cartesian coordinates flag ( True / False ) */')
202          CALL WRITE_1D_L( usingSphericalPolarGrid,1, INDEX_NONE,'usingSphericalPolarGrid =',
203         &'   /* Spherical coordinates flag ( True / False ) */')
204          CALL WRITE_1D_R8( delZ,Nr, INDEX_K,'delZ = ',
205         &'   /* W spacing ( m ) */')
206          CALL WRITE_1D_R8( delP,Nr, INDEX_K,'delP = ',
207         &'   /* W spacing ( Pa ) */')
208          CALL WRITE_1D_R8( delR,Nr, INDEX_K,'delR = ',
209         &'   /* W spacing ( units of r ) */')
210          CALL WRITE_1D_R8( delX, Nx, INDEX_I,'delX = ',
211         &'   /* U spacing ( m - cartesian, degrees - spherical ) */')
212          CALL WRITE_1D_R8( delY, Ny, INDEX_J,'delY = ',
213         &'   /* V spacing ( m - cartesian, degrees - spherical ) */')
214          CALL WRITE_1D_R8( phiMin, 1, INDEX_NONE,'phiMin = ',
215         &'   /* Southern boundary ( ignored - cartesian, degrees - spherical ) */')
216          CALL WRITE_1D_R8( thetaMin, 1, INDEX_NONE,'thetaMin = ',
217         &'   /* Western boundary ( ignored - cartesian, degrees - spherical ) */')
218          CALL WRITE_1D_R8( rSphere, 1, INDEX_NONE,'rSphere = ',
219         &'   /* Radius ( ignored - cartesian, m - spherical ) */')
220          DO bi=1,nSx
221           DO I=1,sNx
222            xcoord((bi-1)*sNx+I) = xc(I,1,bi,1)
223           ENDDO
224          ENDDO
225          CALL WRITE_1D_R8( xcoord, sNx*nSx, INDEX_I,'xcoord = ',
226         &'   /* P-point X coordinate (  m - cartesian, degrees - spherical ) */')
227          DO bj=1,nSy
228           DO J=1,sNy
229            ycoord((bj-1)*sNy+J) = yc(1,J,1,bj)
230           ENDDO
231          ENDDO
232          CALL WRITE_1D_R8( ycoord, sNy*nSy, INDEX_J,'ycoord = ',
233         &'   /* P-point Y coordinate (  m - cartesian, degrees - spherical ) */')
234          DO K=1,Nr
235           rcoord(K) = rc(K)
236          ENDDO
237          CALL WRITE_1D_R8( rcoord, Nr, INDEX_K,'rcoord = ',
238         &'   /* P-point R coordinate (  units of r ) */')
239    
240    
241    
242        WRITE(msgBuf,'(A)') ' '        WRITE(msgBuf,'(A)') ' '
243        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
244       &  SQUEEZE_RIGHT , 1)       &  SQUEEZE_RIGHT , 1)
       _END_MASTER(myThid)  
245    
246          _END_MASTER(myThid)
247        _BARRIER        _BARRIER
248    
249    
250        RETURN        RETURN
251    100 FORMAT(A,    100 FORMAT(A,
252       &'$Id$'       &' '
253       &)       &)
254        END        END
255    
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.16
  ViewVC Help
Powered by ViewVC 1.1.22