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

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

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


Revision 1.7 - (hide annotations) (download)
Thu May 21 18:25:48 1998 UTC (26 years ago) by cnh
Branch: MAIN
CVS Tags: checkpoint2
Changes since 1.6: +4 -2 lines
Added support for binary IO of model fields for restart and/or
postprocessing

1 cnh 1.7 C $Header: /u/gcmpack/models/MITgcmUV/model/src/config_summary.F,v 1.6 1998/04/27 04:24:22 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     CALL WRITE_1D_R8( f0,1, INDEX_NONE,'f0 =',
90     &' /* Reference coriolis parameter ( 1/s ) */')
91     CALL WRITE_1D_R8( beta,1, INDEX_NONE,'beta =',
92     &' /* Beta ( 1/(m.s) ) */')
93     WRITE(msgBuf,'(A)') '// '
94     CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)
95     WRITE(msgBuf,'(A)') '// Elliptic solver(s) paramters ( PARM02 in namelist ) '
96     CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)
97     WRITE(msgBuf,'(A)') '// '
98     CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)
99     CALL WRITE_1D_I( cg2dMaxIters,1, INDEX_NONE,'cg2dMaxIters =',
100     &' /* Upper limit on 2d con. grad iterations */')
101     CALL WRITE_1D_I( cg2dChkResFreq,1, INDEX_NONE,'cg2dChkResFreq =',
102     &' /* 2d con. grad convergence test frequency */')
103     CALL WRITE_1D_R8( cg2dTargetResidual,1, INDEX_NONE,'cg2dTargetResidual =',
104     &' /* 2d con. grad target residual */')
105    
106     WRITE(msgBuf,'(A)') '// '
107     CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)
108     WRITE(msgBuf,'(A)') '// Time stepping paramters ( PARM03 in namelist ) '
109     CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)
110     WRITE(msgBuf,'(A)') '// '
111     CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)
112     CALL WRITE_1D_I( nIter0,1, INDEX_NONE,'nIter0 =',
113     &' /* Base timestep number */')
114     CALL WRITE_1D_I( nTimeSteps,1, INDEX_NONE,'nTimeSteps =',
115     &' /* Number of timesteps */')
116     CALL WRITE_1D_R8( deltaTmom,1, INDEX_NONE,'deltatTmom =',
117     &' /* Momentum equation timestep ( s ) */')
118     CALL WRITE_1D_R8( deltaTtracer,1, INDEX_NONE,'deltatTtracer =',
119     &' /* Tracer equation timestep ( s ) */')
120     CALL WRITE_1D_R8( abEps,1, INDEX_NONE,'abEps =',
121     &' /* Adams-Bashforth stabilizing weight */')
122     CALL WRITE_1D_R8( tauCD,1, INDEX_NONE,'tauCD =',
123     &' /* CD coupling time-scale ( s ) */')
124     CALL WRITE_1D_R8( rCD,1, INDEX_NONE,'rCD =',
125     &' /* Normalised CD coupling parameter */')
126     CALL WRITE_1D_R8( startTime,1, INDEX_NONE,'startTime =',
127     &' /* Run start time ( s ). */')
128     CALL WRITE_1D_R8( endTime,1, INDEX_NONE,'endTime =',
129     &' /* Integration ending time ( s ). */')
130 cnh 1.7 CALL WRITE_1D_R8( pChkPtFreq,1, INDEX_NONE,'pChkPtFreq =',
131     &' /* Permanent restart/checkpoint file interval ( s ). */')
132 cnh 1.6 CALL WRITE_1D_R8( chkPtFreq,1, INDEX_NONE,'chkPtFreq =',
133 cnh 1.7 &' /* Rolling restart/checkpoint file interval ( s ). */')
134 cnh 1.6 CALL WRITE_1D_R8( dumpFreq,1, INDEX_NONE,'dumpFreq =',
135     &' /* Model state write out interval ( s ). */')
136    
137     WRITE(msgBuf,'(A)') '// '
138     CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)
139     WRITE(msgBuf,'(A)') '// Gridding paramters ( PARM04 in namelist ) '
140     CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)
141     WRITE(msgBuf,'(A)') '// '
142     CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)
143     CALL WRITE_1D_L( usingCartesianGrid,1, INDEX_NONE,'usingCartesianGrid =',
144     &' /* Cartesian coordinates flag ( True / False ) */')
145     CALL WRITE_1D_L( usingSphericalPolarGrid,1, INDEX_NONE,'usingSphericalPolarGrid =',
146     &' /* Spherical coordinates flag ( True / False ) */')
147     CALL WRITE_1D_R8( delZ,Nz, INDEX_K,'delZ = ',
148     &' /* W spacing ( m ) */')
149     CALL WRITE_1D_R8( delX, Nx, INDEX_I,'delX = ',
150     &' /* U spacing ( m - cartesian, degrees - spherical ) */')
151     CALL WRITE_1D_R8( delY, Ny, INDEX_J,'delY = ',
152     &' /* V spacing ( m - cartesian, degrees - spherical ) */')
153     CALL WRITE_1D_R8( phiMin, 1, INDEX_NONE,'phiMin = ',
154     &' /* Southern boundary ( ignored - cartesian, degrees - spherical ) */')
155     CALL WRITE_1D_R8( thetaMin, 1, INDEX_NONE,'thetaMin = ',
156     &' /* Western boundary ( ignored - cartesian, degrees - spherical ) */')
157     CALL WRITE_1D_R8( rSphere, 1, INDEX_NONE,'rSphere = ',
158     &' /* Radius ( ignored - cartesian, m - spherical ) */')
159     DO bi=1,nSx
160     DO I=1,sNx
161     xcoord((bi-1)*sNx+I) = xc(I,1,bi,1)
162     ENDDO
163     ENDDO
164     CALL WRITE_1D_R8( xcoord, Nx, INDEX_I,'xcoord = ',
165     &' /* P-point X coordinate ( m - cartesian, degrees - spherical ) */')
166     DO bj=1,nSy
167     DO J=1,sNy
168     ycoord((bj-1)*sNy+J) = yc(1,J,1,bj)
169     ENDDO
170     ENDDO
171     CALL WRITE_1D_R8( ycoord, Ny, INDEX_J,'ycoord = ',
172     &' /* P-point Y coordinate ( m - cartesian, degrees - spherical ) */')
173     DO K=1,Nz
174     zcoord(K) = zc(K)
175     ENDDO
176     CALL WRITE_1D_R8( zcoord, Nz, INDEX_K,'zcoord = ',
177     &' /* P-point Z coordinate ( m ) */')
178    
179 cnh 1.5
180    
181 cnh 1.1 WRITE(msgBuf,'(A)') ' '
182     CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
183     & SQUEEZE_RIGHT , 1)
184 cnh 1.5
185 cnh 1.1 _END_MASTER(myThid)
186     _BARRIER
187    
188    
189     RETURN
190     100 FORMAT(A,
191 cnh 1.4 &' '
192 cnh 1.1 &)
193     END
194    

  ViewVC Help
Powered by ViewVC 1.1.22