/[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.21 by adcroft, Fri Jul 30 15:42:58 1999 UTC revision 1.22 by heimbach, Fri Jun 9 02:45:04 2000 UTC
# Line 4  C $Header$ Line 4  C $Header$
4    
5  CStartOfInterface  CStartOfInterface
6        SUBROUTINE CONFIG_SUMMARY( myThid )        SUBROUTINE CONFIG_SUMMARY( myThid )
7  C     /==========================================================\  C     /==========================================================
8  C     | SUBROUTINE CONFIG_SUMMARY                                |  C     | SUBROUTINE CONFIG_SUMMARY                                |
9  C     | o Summarize model prognostic variables.                  |  C     | o Summarize model prognostic variables.                  |
10  C     |==========================================================|  C     |==========================================================|
# Line 35  C     == Local variables == Line 35  C     == Local variables ==
35        CHARACTER*(MAX_LEN_MBUF) msgBuf        CHARACTER*(MAX_LEN_MBUF) msgBuf
36        INTEGER                  I,J,K        INTEGER                  I,J,K
37        INTEGER                  bi, bj        INTEGER                  bi, bj
38        _RL                      xcoord(Nx)        _RL                     xcoord(Nx)
39        _RL                      ycoord(Ny)        _RL                     ycoord(Ny)
40        _RL                      rcoord(Nr)        _RL                     rcoord(Nr)
41    
42    
43        _BARRIER        _BARRIER
# Line 69  C     == Local variables == Line 69  C     == Local variables ==
69       &'   /* Reference temperature profile ( oC or oK ) */')       &'   /* Reference temperature profile ( oC or oK ) */')
70        CALL WRITE_1D_R8( sRef, Nr, INDEX_K,'sRef =',        CALL WRITE_1D_R8( sRef, Nr, INDEX_K,'sRef =',
71       &'   /* Reference salinity profile ( ppt ) */')       &'   /* Reference salinity profile ( ppt ) */')
72        CALL WRITE_1D_R8( viscAh, 1, INDEX_NONE,'viscAh =',        CALL WRITE_0D_R8( viscAh, INDEX_NONE,'viscAh =',
73       &'   /* Lateral eddy viscosity ( m^2/s ) */')       &'   /* Lateral eddy viscosity ( m^2/s ) */')
74        CALL WRITE_1D_R8( viscA4, 1, INDEX_NONE,'viscA4 =',        CALL WRITE_0D_R8( viscA4, INDEX_NONE,'viscAh =',
75       &'   /* Lateral biharmonic viscosity ( m^4/s ) */')       &'   /* Lateral biharmonic viscosity ( m^4/s ) */')
76        CALL WRITE_1D_L( no_slip_sides,1, INDEX_NONE,        CALL WRITE_0D_L( no_slip_sides, INDEX_NONE,
77       & 'no_slip_sides =', '  /* Viscous BCs: No-slip sides */')       & 'no_slip_sides =', '  /* Viscous BCs: No-slip sides */')
78        IF ( viscAz .NE. UNSET_RL ) THEN        IF ( viscAz .NE. UNSET_RL ) THEN
79         CALL WRITE_1D_R8( viscAz, 1, INDEX_NONE,'viscAz =',         CALL WRITE_0D_R8( viscAz, INDEX_NONE,'viscAz =',
80       & '   /* Vertical eddy viscosity ( m^2/s ) */')       & '   /* Vertical eddy viscosity ( m^2/s ) */')
81        ENDIF        ENDIF
82        IF ( viscAp .NE. UNSET_RL ) THEN        IF ( viscAp .NE. UNSET_RL ) THEN
83         CALL WRITE_1D_R8( viscAp, 1, INDEX_NONE,'viscAp =',         CALL WRITE_0D_R8( viscAp, INDEX_NONE,'viscAp =',
84       & '   /* Vertical eddy viscosity ( Pa^2/s ) */')       & '   /* Vertical eddy viscosity ( Pa^2/s ) */')
85        ENDIF        ENDIF
86        CALL WRITE_1D_R8( viscAr, 1, INDEX_NONE,'viscAr =',        CALL WRITE_0D_R8( viscAr,  INDEX_NONE,'viscAr =',
87       &'   /* Vertical eddy viscosity ( units of r^2/s ) */')       &'   /* Vertical eddy viscosity ( units of r^2/s ) */')
88        CALL WRITE_1D_R8( diffKhT, 1, INDEX_NONE,'diffKhT =',        CALL WRITE_0D_R8( diffKhT, INDEX_NONE,'diffKhT =',
89       &'   /* Laplacian diffusion of heat laterally ( m^2/s ) */')       &'   /* Laplacian diffusion of heat laterally ( m^2/s ) */')
90        CALL WRITE_1D_R8( diffK4T, 1, INDEX_NONE,'diffK4T =',        CALL WRITE_0D_R8( diffK4T, INDEX_NONE,'diffK4T =',
91       &'   /* Bihaarmonic diffusion of heat laterally ( m^4/s ) */')       &'   /* Bihaarmonic diffusion of heat laterally ( m^4/s ) */')
92        CALL WRITE_1D_R8( diffKzT, 1, INDEX_NONE,'diffKzT =',        CALL WRITE_0D_R8( diffKzT, INDEX_NONE,'diffKzT =',
93       &'   /* Laplacian diffusion of heat vertically ( m^2/s ) */')       &'   /* Laplacian diffusion of heat vertically ( m^2/s ) */')
94        CALL WRITE_1D_R8( diffKrT, 1, INDEX_NONE,'diffKrT =',        CALL WRITE_0D_R8( diffKrT, INDEX_NONE,'diffKrT =',
95       &'   /* Laplacian diffusion of heat vertically ( m^2/s ) */')       &'   /* Laplacian diffusion of heat vertically ( m^2/s ) */')
96        CALL WRITE_1D_R8( diffKhS, 1, INDEX_NONE,'diffKhS =',        CALL WRITE_0D_R8( diffKhS, INDEX_NONE,'diffKhS =',
97       &'   /* Laplacian diffusion of salt laterally ( m^2/s ) */')       &'   /* Laplacian diffusion of salt laterally ( m^2/s ) */')
98        CALL WRITE_1D_R8( diffK4S, 1, INDEX_NONE,'diffK4S =',        CALL WRITE_0D_R8( diffK4S, INDEX_NONE,'diffK4S =',
99       &'   /* Bihaarmonic diffusion of salt laterally ( m^4/s ) */')       &'   /* Bihaarmonic diffusion of salt laterally ( m^4/s ) */')
100        CALL WRITE_1D_R8( diffKzS, 1, INDEX_NONE,'diffKzS =',        CALL WRITE_0D_R8( diffKzS, INDEX_NONE,'diffKzS =',
101       &'   /* Laplacian diffusion of salt vertically ( m^2/s ) */')       &'   /* Laplacian diffusion of salt vertically ( m^2/s ) */')
102        CALL WRITE_1D_R8( diffKrS, 1, INDEX_NONE,'diffKrS =',        CALL WRITE_0D_R8( diffKrS, INDEX_NONE,'diffKrS =',
103       &'   /* Laplacian diffusion of salt vertically ( m^2/s ) */')       &'   /* Laplacian diffusion of salt vertically ( m^2/s ) */')
104        CALL WRITE_1D_R8( tAlpha,1, INDEX_NONE,'tAlpha =',        CALL WRITE_0D_R8( tAlpha,  INDEX_NONE,'tAlpha =',
105       &'   /* Linear EOS thermal expansion coefficient ( 1/degree ) */')       &'   /* Linear EOS thermal expansion coefficient ( 1/degree ) */')
106        CALL WRITE_1D_R8( sBeta, 1, INDEX_NONE,'sBeta =',        CALL WRITE_0D_R8( sBeta,   INDEX_NONE,'sBeta =',
107       &'   /* Linear EOS haline contraction coefficient ( 1/ppt ) */')       &'   /* Linear EOS haline contraction coefficient ( 1/ppt ) */')
108        IF ( eosType .EQ. 'POLY3' ) THEN        IF ( eosType .EQ. 'POLY3' ) THEN
109          WRITE(msgBuf,'(A)')          WRITE(msgBuf,'(A)')
# Line 115  C     == Local variables == Line 115  C     == Local variables ==
115       &                    SQUEEZE_RIGHT , 1)       &                    SQUEEZE_RIGHT , 1)
116          ENDDO          ENDDO
117        ENDIF        ENDIF
118        CALL WRITE_1D_R8( rhonil,1, INDEX_NONE,'rhonil =',        CALL WRITE_0D_R8( rhonil,  INDEX_NONE,'rhonil =',
119       &'   /* Reference density ( kg/m^3 ) */')       &'   /* Reference density ( kg/m^3 ) */')
120        CALL WRITE_1D_R8( rhoConst,1, INDEX_NONE,'rhoConst =',        CALL WRITE_0D_R8( rhoConst, INDEX_NONE,'rhoConst =',
121       &'   /* Reference density ( kg/m^3 ) */')       &'   /* Reference density ( kg/m^3 ) */')
122        CALL WRITE_1D_R8( gravity,1, INDEX_NONE,'gravity =',        CALL WRITE_0D_R8( gravity, INDEX_NONE,'gravity =',
123       &'   /* Gravitational acceleration ( m/s^2 ) */')       &'   /* Gravitational acceleration ( m/s^2 ) */')
124        CALL WRITE_1D_R8( gBaro,1, INDEX_NONE,'gBaro =',        CALL WRITE_0D_R8( gBaro,   INDEX_NONE,'gBaro =',
125       &'   /* Barotropic gravity ( m/s^2 ) */')       &'   /* Barotropic gravity ( m/s^2 ) */')
126        CALL WRITE_1D_R8( f0,1, INDEX_NONE,'f0 =',        CALL WRITE_0D_R8( f0,      INDEX_NONE,'f0 =',
127       &'   /* Reference coriolis parameter ( 1/s ) */')       &'   /* Reference coriolis parameter ( 1/s ) */')
128        CALL WRITE_1D_R8( beta,1, INDEX_NONE,'beta =',        CALL WRITE_0D_R8( beta,    INDEX_NONE,'beta =',
129       &'   /* Beta ( 1/(m.s) ) */')       &'   /* Beta ( 1/(m.s) ) */')
130        CALL WRITE_1D_R8( freeSurfFac,1, INDEX_NONE,'freeSurfFac =',        CALL WRITE_0D_R8( freeSurfFac, INDEX_NONE,'freeSurfFac =',
131       &'   /* Implcit free surface factor */')       &'   /* Implcit free surface factor */')
132        CALL WRITE_1D_L( implicitFreeSurface,1, INDEX_NONE,        CALL WRITE_0D_L( implicitFreeSurface, INDEX_NONE,
133       &                 'implicitFreeSurface =',       &                 'implicitFreeSurface =',
134       &'   /* Implicit free surface on/off flag */')       &'   /* Implicit free surface on/off flag */')
135        CALL WRITE_1D_L( rigidLid,1, INDEX_NONE,        CALL WRITE_0D_L( rigidLid, INDEX_NONE,
136       &                 'rigidLid =',       &                 'rigidLid =',
137       &'   /* Rigid lid on/off flag */')       &'   /* Rigid lid on/off flag */')
138        CALL WRITE_1D_L( momStepping,1, INDEX_NONE,        CALL WRITE_0D_L( momStepping,  INDEX_NONE,
139       & 'momStepping =', '  /* Momentum equation on/off flag */')       & 'momStepping =', '  /* Momentum equation on/off flag */')
140        CALL WRITE_1D_L( momAdvection,1, INDEX_NONE,        CALL WRITE_0D_L( momAdvection, INDEX_NONE,
141       & 'momAdvection =', '  /* Momentum advection on/off flag */')       & 'momAdvection =', '  /* Momentum advection on/off flag */')
142        CALL WRITE_1D_L( momViscosity,1, INDEX_NONE,        CALL WRITE_0D_L( momViscosity, INDEX_NONE,
143       & 'momViscosity =', '  /* Momentum viscosity on/off flag */')       & 'momViscosity =', '  /* Momentum viscosity on/off flag */')
144        CALL WRITE_1D_L( useCoriolis,1, INDEX_NONE,        CALL WRITE_0D_L( useCoriolis,  INDEX_NONE,
145       & 'useCoriolis =', '  /* Coriolis on/off flag */')       & 'useCoriolis =', '  /* Coriolis on/off flag */')
146        CALL WRITE_1D_L( momForcing,1, INDEX_NONE,        CALL WRITE_0D_L( momForcing,   INDEX_NONE,
147       & 'momForcing =', '  /* Momentum forcing on/off flag */')       & 'momForcing =', '  /* Momentum forcing on/off flag */')
148        CALL WRITE_1D_L( momPressureForcing,1, INDEX_NONE,        CALL WRITE_0D_L( momPressureForcing, INDEX_NONE,
149       & 'momPressureForcing =',         & 'momPressureForcing =',  
150       & '  /* Momentum pressure term on/off flag */')       & '  /* Momentum pressure term on/off flag */')
151        CALL WRITE_1D_L( tempStepping,1, INDEX_NONE,        CALL WRITE_0D_L( tempStepping,  INDEX_NONE,
152       & 'tempStepping =', '  /* Temperature equation on/off flag */')       & 'tempStepping =', '  /* Temperature equation on/off flag */')
153        CALL WRITE_1D_L( openBoundaries,1, INDEX_NONE,        CALL WRITE_0D_L( openBoundaries, INDEX_NONE,
154       & 'openBoundaries =', '  /* OpenBoundaries on/off flag */')       & 'openBoundaries =', '  /* OpenBoundaries on/off flag */')
155        CALL WRITE_1D_L( nonHydrostatic,1, INDEX_NONE,        CALL WRITE_0D_L( nonHydrostatic, INDEX_NONE,
156       & 'nonHydrostatic =', '  /* Non-Hydrostatic on/off flag */')       & 'nonHydrostatic =', '  /* Non-Hydrostatic on/off flag */')
157        CALL WRITE_1D_R8( GMMaxSlope,1, INDEX_NONE,'GMMaxSlope =',        CALL WRITE_0D_R8( GMmaxslope,   INDEX_NONE,'GMMaxslope =',
158       &'   /* Max. slope allowed in GM/Redi tensor */')       &'   /* Max. slope allowed in GM/Redi tensor */')
159        CALL WRITE_1D_R8( GMLength,1, INDEX_NONE,'GMLength =',        CALL WRITE_0D_R8( GMlength,     INDEX_NONE,'GMlength =',
160       &'   /* Length to use in Visbeck et al. formula for K (m) */')       &'   /* Length to use in Visbeck et al. formula for K (m) */')
161        CALL WRITE_1D_R8( GMAlpha,1, INDEX_NONE,'GMAlpha =',        CALL WRITE_0D_R8( GMalpha,      INDEX_NONE,'GMalpha =',
162       &'   /* alpha to use in Visbeck et al. formula for K */')       &'   /* alpha to use in Visbeck et al. formula for K */')
163        CALL WRITE_1D_R8( GMdepth,1, INDEX_NONE,'GMdepth =',        CALL WRITE_0D_R8( GMdepth,      INDEX_NONE,'GMdepth =',
164       &'   /* Depth to integrate for Visbeck et. al Richardson # (m) */')       &'   /* Depth to integrate for Visbeck et. al Richardson # (m) */')
165        CALL WRITE_1D_R8( GMkbackground,1, INDEX_NONE,'GMkbackground =',        CALL WRITE_0D_R8( GMkbackground,INDEX_NONE,'GMkbackground =',
166       &'   /* background value of GM/Redi coefficient m^2/s */')       &'   /* background value of GM/Redi coefficient m^2/s */')
167        WRITE(msgBuf,'(A)') '//  '        WRITE(msgBuf,'(A)') '//  '
168        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
# Line 175  C     == Local variables == Line 175  C     == Local variables ==
175        WRITE(msgBuf,'(A)') '//  '        WRITE(msgBuf,'(A)') '//  '
176        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
177       &                    SQUEEZE_RIGHT , 1)       &                    SQUEEZE_RIGHT , 1)
178        CALL WRITE_1D_I( cg2dMaxIters,1, INDEX_NONE,'cg2dMaxIters =',        CALL WRITE_0D_I( cg2dMaxIters,   INDEX_NONE,'cg2dMaxIters =',
179       &'   /* Upper limit on 2d con. grad iterations  */')       &'   /* Upper limit on 2d con. grad iterations  */')
180        CALL WRITE_1D_I( cg2dChkResFreq,1, INDEX_NONE,'cg2dChkResFreq =',        CALL WRITE_0D_I( cg2dChkResFreq, INDEX_NONE,'cg2dChkResFreq =',
181       &'   /* 2d con. grad convergence test frequency */')       &'   /* 2d con. grad convergence test frequency */')
182        CALL WRITE_1D_R8( cg2dTargetResidual,1, INDEX_NONE,        CALL WRITE_0D_R8( cg2dTargetResidual, INDEX_NONE,
183       & 'cg2dTargetResidual =',       & 'cg2dTargetResidual =',
184       &'   /* 2d con. grad target residual  */')       &'   /* 2d con. grad target residual  */')
185    
# Line 193  C     == Local variables == Line 193  C     == Local variables ==
193        WRITE(msgBuf,'(A)') '//  '        WRITE(msgBuf,'(A)') '//  '
194        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
195       &                    SQUEEZE_RIGHT , 1)       &                    SQUEEZE_RIGHT , 1)
196        CALL WRITE_1D_I( nIter0,1, INDEX_NONE,'nIter0 =',        CALL WRITE_0D_I( nIter0, INDEX_NONE,'nIter0 =',
197       &'   /* Base timestep number  */')       &'   /* Base timestep number  */')
198        CALL WRITE_1D_I( nTimeSteps,1, INDEX_NONE,'nTimeSteps =',        CALL WRITE_0D_I( nTimeSteps, INDEX_NONE,'nTimeSteps =',
199       &'   /* Number of timesteps */')       &'   /* Number of timesteps */')
200        CALL WRITE_1D_R8( deltaTmom,1, INDEX_NONE,'deltaTmom =',        CALL WRITE_0D_R8( deltaTmom, INDEX_NONE,'deltatTmom =',
201       &'   /* Momentum equation timestep ( s ) */')       &'   /* Momentum equation timestep ( s ) */')
202        CALL WRITE_1D_R8( deltaTtracer,1, INDEX_NONE,'deltaTtracer =',        CALL WRITE_0D_R8( deltaTtracer, INDEX_NONE,'deltatTtracer =',
203       &'   /* Tracer equation timestep ( s ) */')       &'   /* Tracer equation timestep ( s ) */')
204        CALL WRITE_1D_R8( deltaTClock ,1, INDEX_NONE,'deltaTClock  =',        CALL WRITE_0D_R8( deltaTClock, INDEX_NONE,'deltatTClock  =',
205       &'   /* Model clock timestep ( s ) */')       &'   /* Model clock timestep ( s ) */')
206        CALL WRITE_1D_R8( cAdjFreq,1, INDEX_NONE,'cAdjFreq =',        CALL WRITE_0D_R8( cAdjFreq, INDEX_NONE,'cAdjFreq =',
207       &'   /* Convective adjustment interval ( s ) */')       &'   /* Convective adjustment interval ( s ) */')
208        CALL WRITE_1D_R8( abEps,1, INDEX_NONE,'abEps =',        CALL WRITE_0D_R8( abeps, INDEX_NONE,'abeps =',
209       &'   /* Adams-Bashforth stabilizing weight */')       &'   /* Adams-Bashforth stabilizing weight */')
210        CALL WRITE_1D_R8( tauCD,1, INDEX_NONE,'tauCD =',        CALL WRITE_0D_R8( tauCD, INDEX_NONE,'tauCD =',
211       &'   /* CD coupling time-scale ( s ) */')       &'   /* CD coupling time-scale ( s ) */')
212        CALL WRITE_1D_R8( rCD,1, INDEX_NONE,'rCD =',        CALL WRITE_0D_R8( rCD, INDEX_NONE,'rCD =',
213       &'   /* Normalised CD coupling parameter */')       &'   /* Normalised CD coupling parameter */')
214        CALL WRITE_1D_R8( startTime,1, INDEX_NONE,'startTime =',        CALL WRITE_0D_R8( startTime, INDEX_NONE,'startTime =',
215       &'   /* Run start time ( s ). */')       &'   /* Run start time ( s ). */')
216        CALL WRITE_1D_R8( endTime,1, INDEX_NONE,'endTime =',        CALL WRITE_0D_R8( endTime, INDEX_NONE,'endTime =',
217       &'   /* Integration ending time ( s ). */')       &'   /* Integration ending time ( s ). */')
218        CALL WRITE_1D_R8( pChkPtFreq,1, INDEX_NONE,'pChkPtFreq =',        CALL WRITE_0D_R8( pChkPtFreq, INDEX_NONE,'pChkPtFreq =',
219       &'   /* Permanent restart/checkpoint file interval ( s ). */')       &'   /* Permanent restart/checkpoint file interval ( s ). */')
220        CALL WRITE_1D_R8( chkPtFreq,1, INDEX_NONE,'chkPtFreq =',        CALL WRITE_0D_R8( chkPtFreq, INDEX_NONE,'chkPtFreq =',
221       &'   /* Rolling restart/checkpoint file interval ( s ). */')       &'   /* Rolling restart/checkpoint file interval ( s ). */')
222        CALL WRITE_1D_R8( dumpFreq,1, INDEX_NONE,'dumpFreq =',        CALL WRITE_0D_R8( dumpFreq, INDEX_NONE,'dumpFreq =',
223       &'   /* Model state write out interval ( s ). */')       &'   /* Model state write out interval ( s ). */')
224    
225        WRITE(msgBuf,'(A)') '//  '        WRITE(msgBuf,'(A)') '//  '
# Line 232  C     == Local variables == Line 232  C     == Local variables ==
232        WRITE(msgBuf,'(A)') '//  '        WRITE(msgBuf,'(A)') '//  '
233        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
234       &                    SQUEEZE_RIGHT , 1)       &                    SQUEEZE_RIGHT , 1)
235        CALL WRITE_1D_L( usingCartesianGrid,1, INDEX_NONE,        CALL WRITE_0D_L( usingCartesianGrid, INDEX_NONE,
236       & 'usingCartesianGrid =',       & 'usingCartesianGrid =',
237       &'   /* Cartesian coordinates flag ( True / False ) */')       &'   /* Cartesian coordinates flag ( True / False ) */')
238        CALL WRITE_1D_L( usingSphericalPolarGrid,1, INDEX_NONE,        CALL WRITE_0D_L( usingSphericalPolarGrid, INDEX_NONE,
239       & 'usingSphericalPolarGrid =',       & 'usingSphericalPolarGrid =',
240       &'   /* Spherical coordinates flag ( True / False ) */')       &'   /* Spherical coordinates flag ( True / False ) */')
241        CALL WRITE_1D_R8( delZ,Nr, INDEX_K,'delZ = ',        CALL WRITE_1D_R8( delZ,Nr, INDEX_K,'delZ = ',
# Line 248  C     == Local variables == Line 248  C     == Local variables ==
248       &'   /* U spacing ( m - cartesian, degrees - spherical ) */')       &'   /* U spacing ( m - cartesian, degrees - spherical ) */')
249        CALL WRITE_1D_R8( delY, Ny, INDEX_J,'delY = ',        CALL WRITE_1D_R8( delY, Ny, INDEX_J,'delY = ',
250       &'   /* V spacing ( m - cartesian, degrees - spherical ) */')       &'   /* V spacing ( m - cartesian, degrees - spherical ) */')
251        CALL WRITE_1D_R8( phiMin, 1, INDEX_NONE,'phiMin = ',        CALL WRITE_0D_R8( phiMin, INDEX_NONE,'phiMin = ',
252       &'   /* South edge (ignored - cartesian, degrees - spherical ) */')       &'   /* South edge (ignored - cartesian, degrees - spherical ) */')
253        CALL WRITE_1D_R8( thetaMin, 1, INDEX_NONE,'thetaMin = ',        CALL WRITE_0D_R8( thetaMin, INDEX_NONE,'thetaMin = ',
254       &'   /* West edge ( ignored - cartesian, degrees - spherical ) */')       &'   /* West edge ( ignored - cartesian, degrees - spherical ) */')
255        CALL WRITE_1D_R8( rSphere, 1, INDEX_NONE,'rSphere = ',        CALL WRITE_0D_R8( rSphere, INDEX_NONE,'rSphere = ',
256       &'   /* Radius ( ignored - cartesian, m - spherical ) */')       &'   /* Radius ( ignored - cartesian, m - spherical ) */')
257        DO bi=1,nSx        DO bi=1,nSx
258         DO I=1,sNx         DO I=1,sNx
259          xcoord((bi-1)*sNx+I) = xc(I,1,bi,1)          xcoord((bi-1)*sNx+I) = xC(I,1,bi,1)
260         ENDDO         ENDDO
261        ENDDO        ENDDO
262        CALL WRITE_1D_R8( xcoord, sNx*nSx, INDEX_I,'xcoord = ',        CALL WRITE_1D_R8( xcoord, sNx*nSx, INDEX_I,'xcoord = ',
263       &'   /* P-point X coord ( m - cartesian, degrees - spherical ) */')       &'   /* P-point X coord ( m - cartesian, degrees - spherical ) */')
264        DO bj=1,nSy        DO bj=1,nSy
265         DO J=1,sNy         DO J=1,sNy
266          ycoord((bj-1)*sNy+J) = yc(1,J,1,bj)          ycoord((bj-1)*sNy+J) = yC(1,J,1,bj)
267         ENDDO         ENDDO
268        ENDDO        ENDDO
269        CALL WRITE_1D_R8( ycoord, sNy*nSy, INDEX_J,'ycoord = ',        CALL WRITE_1D_R8( ycoord, sNy*nSy, INDEX_J,'ycoord = ',
270       &'   /* P-point Y coord ( m - cartesian, degrees - spherical ) */')       &'   /* P-point Y coord ( m - cartesian, degrees - spherical ) */')
271        DO K=1,Nr        DO K=1,Nr
272         rcoord(K) = rc(K)         rcoord(K) = rC(K)
273        ENDDO        ENDDO
274        CALL WRITE_1D_R8( rcoord, Nr, INDEX_K,'rcoord = ',        CALL WRITE_1D_R8( rcoord, Nr, INDEX_K,'rcoord = ',
275       &'   /* P-point R coordinate (  units of r ) */')       &'   /* P-point R coordinate (  units of r ) */')
Legend:
Removed from v.1.21  
changed lines
  Added in v.1.22
  ViewVC Help
Powered by ViewVC 1.1.22