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

Legend:
Removed from v.1.2  
changed lines
  Added in v.1.20

  ViewVC Help
Powered by ViewVC 1.1.22