/[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.16 by cnh, Sat Sep 5 17:52:13 1998 UTC revision 1.17 by cnh, Wed Oct 28 03:11:36 1998 UTC
# Line 44  C     == Local variables == Line 44  C     == Local variables ==
44    
45        WRITE(msgBuf,'(A)')        WRITE(msgBuf,'(A)')
46       &'// ======================================================='       &'// ======================================================='
47        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
48         &                    SQUEEZE_RIGHT , 1)
49        WRITE(msgBuf,'(A)') '// Model configuration'        WRITE(msgBuf,'(A)') '// Model configuration'
50        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
51         &                    SQUEEZE_RIGHT , 1)
52        WRITE(msgBuf,'(A)')        WRITE(msgBuf,'(A)')
53       &'// ======================================================='       &'// ======================================================='
54        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
55       &  SQUEEZE_RIGHT , 1)       &  SQUEEZE_RIGHT , 1)
56    
57        WRITE(msgBuf,'(A)') '//  '        WRITE(msgBuf,'(A)') '//  '
58        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
59        WRITE(msgBuf,'(A)') '// "Physical" paramters ( PARM01 in namelist ) '       &                    SQUEEZE_RIGHT , 1)
60        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)        WRITE(msgBuf,'(A)')
61         & '// "Physical" paramters ( PARM01 in namelist ) '
62          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
63         &                    SQUEEZE_RIGHT , 1)
64        WRITE(msgBuf,'(A)') '//  '        WRITE(msgBuf,'(A)') '//  '
65        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
66         &                    SQUEEZE_RIGHT , 1)
67        CALL WRITE_1D_R8( tRef, Nr, INDEX_K,'tRef =',        CALL WRITE_1D_R8( tRef, Nr, INDEX_K,'tRef =',
68       &'   /* Reference temperature profile ( oC or oK ) */')       &'   /* Reference temperature profile ( oC or oK ) */')
69        CALL WRITE_1D_R8( sRef, Nr, INDEX_K,'sRef =',        CALL WRITE_1D_R8( sRef, Nr, INDEX_K,'sRef =',
# Line 93  C     == Local variables == Line 99  C     == Local variables ==
99        CALL WRITE_1D_R8( sBeta, 1, INDEX_NONE,'sBeta =',        CALL WRITE_1D_R8( sBeta, 1, INDEX_NONE,'sBeta =',
100       &'   /* Linear EOS haline contraction coefficient ( 1/ppt ) */')       &'   /* Linear EOS haline contraction coefficient ( 1/ppt ) */')
101        IF ( eosType .EQ. 'POLY3' ) THEN        IF ( eosType .EQ. 'POLY3' ) THEN
102          WRITE(msgBuf,'(A)') '// Polynomial EQS parameters ( from POLY3.COEFFS ) '          WRITE(msgBuf,'(A)')
103         &   '// Polynomial EQS parameters ( from POLY3.COEFFS ) '
104          DO K = 1, Nr          DO K = 1, Nr
105           WRITE(msgBuf,'(I3,13F8.3)')           WRITE(msgBuf,'(I3,13F8.3)')
106       &   K,eosRefT(K),eosRefS(K),eosSig0(K), (eosC(I,K),I=1,9)       &   K,eosRefT(K),eosRefS(K),eosSig0(K), (eosC(I,K),I=1,9)
107           CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)           CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
108         &                    SQUEEZE_RIGHT , 1)
109          ENDDO          ENDDO
110        ENDIF        ENDIF
111        CALL WRITE_1D_R8( rhonil,1, INDEX_NONE,'rhonil =',        CALL WRITE_1D_R8( rhonil,1, INDEX_NONE,'rhonil =',
# Line 129  C     == Local variables == Line 137  C     == Local variables ==
137        CALL WRITE_1D_L( momForcing,1, INDEX_NONE,        CALL WRITE_1D_L( momForcing,1, INDEX_NONE,
138       & 'momForcing =', '  /* Momentum forcing on/off flag */')       & 'momForcing =', '  /* Momentum forcing on/off flag */')
139        CALL WRITE_1D_L( momPressureForcing,1, INDEX_NONE,        CALL WRITE_1D_L( momPressureForcing,1, INDEX_NONE,
140       & 'momPressureForcing =', '  /* Momentum pressure term on/off flag */')       & 'momPressureForcing =',  
141         & '  /* Momentum pressure term on/off flag */')
142        CALL WRITE_1D_L( tempStepping,1, INDEX_NONE,        CALL WRITE_1D_L( tempStepping,1, INDEX_NONE,
143       & 'tempStepping =', '  /* Temperature equation on/off flag */')       & 'tempStepping =', '  /* Temperature equation on/off flag */')
144        CALL WRITE_1D_R8( GMMaxSlope,1, INDEX_NONE,'GMMaxSlope =',        CALL WRITE_1D_R8( GMMaxSlope,1, INDEX_NONE,'GMMaxSlope =',
# Line 143  C     == Local variables == Line 152  C     == Local variables ==
152        CALL WRITE_1D_R8( GMkbackground,1, INDEX_NONE,'GMkbackground =',        CALL WRITE_1D_R8( GMkbackground,1, INDEX_NONE,'GMkbackground =',
153       &'   /* background value of GM/Redi coefficient m^2/s */')       &'   /* background value of GM/Redi coefficient m^2/s */')
154        WRITE(msgBuf,'(A)') '//  '        WRITE(msgBuf,'(A)') '//  '
155        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
156         &                    SQUEEZE_RIGHT , 1)
157    
158        WRITE(msgBuf,'(A)') '// Elliptic solver(s) paramters ( PARM02 in namelist ) '        WRITE(msgBuf,'(A)')
159        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)       & '// Elliptic solver(s) paramters ( PARM02 in namelist ) '
160          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
161         &                    SQUEEZE_RIGHT , 1)
162        WRITE(msgBuf,'(A)') '//  '        WRITE(msgBuf,'(A)') '//  '
163        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
164         &                    SQUEEZE_RIGHT , 1)
165        CALL WRITE_1D_I( cg2dMaxIters,1, INDEX_NONE,'cg2dMaxIters =',        CALL WRITE_1D_I( cg2dMaxIters,1, INDEX_NONE,'cg2dMaxIters =',
166       &'   /* Upper limit on 2d con. grad iterations  */')       &'   /* Upper limit on 2d con. grad iterations  */')
167        CALL WRITE_1D_I( cg2dChkResFreq,1, INDEX_NONE,'cg2dChkResFreq =',        CALL WRITE_1D_I( cg2dChkResFreq,1, INDEX_NONE,'cg2dChkResFreq =',
168       &'   /* 2d con. grad convergence test frequency */')       &'   /* 2d con. grad convergence test frequency */')
169        CALL WRITE_1D_R8( cg2dTargetResidual,1, INDEX_NONE,'cg2dTargetResidual =',        CALL WRITE_1D_R8( cg2dTargetResidual,1, INDEX_NONE,
170         & 'cg2dTargetResidual =',
171       &'   /* 2d con. grad target residual  */')       &'   /* 2d con. grad target residual  */')
172    
173        WRITE(msgBuf,'(A)') '//  '        WRITE(msgBuf,'(A)') '//  '
174        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
175        WRITE(msgBuf,'(A)') '// Time stepping paramters ( PARM03 in namelist ) '       &                    SQUEEZE_RIGHT , 1)
176        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)        WRITE(msgBuf,'(A)')
177         & '// Time stepping paramters ( PARM03 in namelist ) '
178          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
179         &                    SQUEEZE_RIGHT , 1)
180        WRITE(msgBuf,'(A)') '//  '        WRITE(msgBuf,'(A)') '//  '
181        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
182         &                    SQUEEZE_RIGHT , 1)
183        CALL WRITE_1D_I( nIter0,1, INDEX_NONE,'nIter0 =',        CALL WRITE_1D_I( nIter0,1, INDEX_NONE,'nIter0 =',
184       &'   /* Base timestep number  */')       &'   /* Base timestep number  */')
185        CALL WRITE_1D_I( nTimeSteps,1, INDEX_NONE,'nTimeSteps =',        CALL WRITE_1D_I( nTimeSteps,1, INDEX_NONE,'nTimeSteps =',
# Line 192  C     == Local variables == Line 210  C     == Local variables ==
210       &'   /* Model state write out interval ( s ). */')       &'   /* Model state write out interval ( s ). */')
211    
212        WRITE(msgBuf,'(A)') '//  '        WRITE(msgBuf,'(A)') '//  '
213        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
214        WRITE(msgBuf,'(A)') '// Gridding paramters ( PARM04 in namelist ) '       &                    SQUEEZE_RIGHT , 1)
215        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)        WRITE(msgBuf,'(A)')
216        WRITE(msgBuf,'(A)') '//  '       & '// Gridding paramters ( PARM04 in namelist ) '
217        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
218        CALL WRITE_1D_L( usingCartesianGrid,1, INDEX_NONE,'usingCartesianGrid =',       &                    SQUEEZE_RIGHT , 1)
219          WRITE(msgBuf,'(A)') '//  '
220          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
221         &                    SQUEEZE_RIGHT , 1)
222          CALL WRITE_1D_L( usingCartesianGrid,1, INDEX_NONE,
223         & 'usingCartesianGrid =',
224       &'   /* Cartesian coordinates flag ( True / False ) */')       &'   /* Cartesian coordinates flag ( True / False ) */')
225        CALL WRITE_1D_L( usingSphericalPolarGrid,1, INDEX_NONE,'usingSphericalPolarGrid =',        CALL WRITE_1D_L( usingSphericalPolarGrid,1, INDEX_NONE,
226         & 'usingSphericalPolarGrid =',
227       &'   /* Spherical coordinates flag ( True / False ) */')       &'   /* Spherical coordinates flag ( True / False ) */')
228        CALL WRITE_1D_R8( delZ,Nr, INDEX_K,'delZ = ',        CALL WRITE_1D_R8( delZ,Nr, INDEX_K,'delZ = ',
229       &'   /* W spacing ( m ) */')       &'   /* W spacing ( m ) */')
# Line 212  C     == Local variables == Line 236  C     == Local variables ==
236        CALL WRITE_1D_R8( delY, Ny, INDEX_J,'delY = ',        CALL WRITE_1D_R8( delY, Ny, INDEX_J,'delY = ',
237       &'   /* V spacing ( m - cartesian, degrees - spherical ) */')       &'   /* V spacing ( m - cartesian, degrees - spherical ) */')
238        CALL WRITE_1D_R8( phiMin, 1, INDEX_NONE,'phiMin = ',        CALL WRITE_1D_R8( phiMin, 1, INDEX_NONE,'phiMin = ',
239       &'   /* Southern boundary ( ignored - cartesian, degrees - spherical ) */')       &'   /* South edge (ignored - cartesian, degrees - spherical ) */')
240        CALL WRITE_1D_R8( thetaMin, 1, INDEX_NONE,'thetaMin = ',        CALL WRITE_1D_R8( thetaMin, 1, INDEX_NONE,'thetaMin = ',
241       &'   /* Western boundary ( ignored - cartesian, degrees - spherical ) */')       &'   /* West edge ( ignored - cartesian, degrees - spherical ) */')
242        CALL WRITE_1D_R8( rSphere, 1, INDEX_NONE,'rSphere = ',        CALL WRITE_1D_R8( rSphere, 1, INDEX_NONE,'rSphere = ',
243       &'   /* Radius ( ignored - cartesian, m - spherical ) */')       &'   /* Radius ( ignored - cartesian, m - spherical ) */')
244        DO bi=1,nSx        DO bi=1,nSx
# Line 223  C     == Local variables == Line 247  C     == Local variables ==
247         ENDDO         ENDDO
248        ENDDO        ENDDO
249        CALL WRITE_1D_R8( xcoord, sNx*nSx, INDEX_I,'xcoord = ',        CALL WRITE_1D_R8( xcoord, sNx*nSx, INDEX_I,'xcoord = ',
250       &'   /* P-point X coordinate (  m - cartesian, degrees - spherical ) */')       &'   /* P-point X coord ( m - cartesian, degrees - spherical ) */')
251        DO bj=1,nSy        DO bj=1,nSy
252         DO J=1,sNy         DO J=1,sNy
253          ycoord((bj-1)*sNy+J) = yc(1,J,1,bj)          ycoord((bj-1)*sNy+J) = yc(1,J,1,bj)
254         ENDDO         ENDDO
255        ENDDO        ENDDO
256        CALL WRITE_1D_R8( ycoord, sNy*nSy, INDEX_J,'ycoord = ',        CALL WRITE_1D_R8( ycoord, sNy*nSy, INDEX_J,'ycoord = ',
257       &'   /* P-point Y coordinate (  m - cartesian, degrees - spherical ) */')       &'   /* P-point Y coord ( m - cartesian, degrees - spherical ) */')
258        DO K=1,Nr        DO K=1,Nr
259         rcoord(K) = rc(K)         rcoord(K) = rc(K)
260        ENDDO        ENDDO
Legend:
Removed from v.1.16  
changed lines
  Added in v.1.17
  ViewVC Help
Powered by ViewVC 1.1.22