/[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.34 by mlosch, Wed Aug 7 16:55:52 2002 UTC revision 1.44 by jmc, Wed Oct 29 00:19:22 2003 UTC
# Line 82  CEOP Line 82  CEOP
82        WRITE(msgBuf,'(A)') '//  '        WRITE(msgBuf,'(A)') '//  '
83        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
84       &                    SQUEEZE_RIGHT , 1)       &                    SQUEEZE_RIGHT , 1)
85          WRITE(msgBuf,'(A,A40)') 'buoyancyRelation = ', buoyancyRelation
86          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
87         &                    SQUEEZE_RIGHT , 1)
88        CALL WRITE_1D_R8( tRef, Nr, INDEX_K,'tRef =',        CALL WRITE_1D_R8( tRef, Nr, INDEX_K,'tRef =',
89       &'   /* Reference temperature profile ( oC or oK ) */')       &'   /* Reference temperature profile ( oC or oK ) */')
90        CALL WRITE_1D_R8( sRef, Nr, INDEX_K,'sRef =',        CALL WRITE_1D_R8( sRef, Nr, INDEX_K,'sRef =',
# Line 92  CEOP Line 95  CEOP
95       &'   /* Lateral biharmonic viscosity ( m^4/s ) */')       &'   /* Lateral biharmonic viscosity ( m^4/s ) */')
96        CALL WRITE_0D_L( no_slip_sides, INDEX_NONE,        CALL WRITE_0D_L( no_slip_sides, INDEX_NONE,
97       & 'no_slip_sides =', '  /* Viscous BCs: No-slip sides */')       & 'no_slip_sides =', '  /* Viscous BCs: No-slip sides */')
98        IF ( viscAz .NE. UNSET_RL ) THEN  c     IF ( viscAz .NE. UNSET_RL ) THEN
99         CALL WRITE_0D_R8( viscAz, INDEX_NONE,'viscAz =',  c      CALL WRITE_0D_R8( viscAz, INDEX_NONE,'viscAz =',
100       & '   /* Vertical eddy viscosity ( m^2/s ) */')  c    & '   /* Vertical eddy viscosity ( m^2/s ) */')
101        ENDIF  c     ENDIF
102        IF ( viscAp .NE. UNSET_RL ) THEN  c     IF ( viscAp .NE. UNSET_RL ) THEN
103         CALL WRITE_0D_R8( viscAp, INDEX_NONE,'viscAp =',  c      CALL WRITE_0D_R8( viscAp, INDEX_NONE,'viscAp =',
104       & '   /* Vertical eddy viscosity ( Pa^2/s ) */')  c    & '   /* Vertical eddy viscosity ( Pa^2/s ) */')
105        ENDIF  c     ENDIF
106        CALL WRITE_0D_R8( viscAr,  INDEX_NONE,'viscAr =',        CALL WRITE_0D_R8( viscAr,  INDEX_NONE,'viscAr =',
107       &'   /* Vertical eddy viscosity ( units of r^2/s ) */')       &'   /* Vertical eddy viscosity ( units of r^2/s ) */')
108        CALL WRITE_0D_R8( diffKhT, INDEX_NONE,'diffKhT =',        CALL WRITE_0D_R8( diffKhT, INDEX_NONE,'diffKhT =',
109       &'   /* Laplacian diffusion of heat laterally ( m^2/s ) */')       &'   /* Laplacian diffusion of heat laterally ( m^2/s ) */')
110        CALL WRITE_0D_R8( diffK4T, INDEX_NONE,'diffK4T =',        CALL WRITE_0D_R8( diffK4T, INDEX_NONE,'diffK4T =',
111       &'   /* Bihaarmonic diffusion of heat laterally ( m^4/s ) */')       &'   /* Bihaarmonic diffusion of heat laterally ( m^4/s ) */')
112        CALL WRITE_0D_R8( diffKzT, INDEX_NONE,'diffKzT =',  c     CALL WRITE_0D_R8( diffKzT, INDEX_NONE,'diffKzT =',
113       &'   /* Laplacian diffusion of heat vertically ( m^2/s ) */')  c    &'   /* Laplacian diffusion of heat vertically ( m^2/s ) */')
114        CALL WRITE_0D_R8( diffKrT, INDEX_NONE,'diffKrT =',        CALL WRITE_0D_R8( diffKrT, INDEX_NONE,'diffKrT =',
115       &'   /* Laplacian diffusion of heat vertically ( m^2/s ) */')       &'   /* Laplacian diffusion of heat vertically ( m^2/s ) */')
116        CALL WRITE_0D_R8( diffKhS, INDEX_NONE,'diffKhS =',        CALL WRITE_0D_R8( diffKhS, INDEX_NONE,'diffKhS =',
117       &'   /* Laplacian diffusion of salt laterally ( m^2/s ) */')       &'   /* Laplacian diffusion of salt laterally ( m^2/s ) */')
118        CALL WRITE_0D_R8( diffK4S, INDEX_NONE,'diffK4S =',        CALL WRITE_0D_R8( diffK4S, INDEX_NONE,'diffK4S =',
119       &'   /* Bihaarmonic diffusion of salt laterally ( m^4/s ) */')       &'   /* Bihaarmonic diffusion of salt laterally ( m^4/s ) */')
120        CALL WRITE_0D_R8( diffKzS, INDEX_NONE,'diffKzS =',  c     CALL WRITE_0D_R8( diffKzS, INDEX_NONE,'diffKzS =',
121       &'   /* Laplacian diffusion of salt vertically ( m^2/s ) */')  c    &'   /* Laplacian diffusion of salt vertically ( m^2/s ) */')
122        CALL WRITE_0D_R8( diffKrS, INDEX_NONE,'diffKrS =',        CALL WRITE_0D_R8( diffKrS, INDEX_NONE,'diffKrS =',
123       &'   /* Laplacian diffusion of salt vertically ( m^2/s ) */')       &'   /* Laplacian diffusion of salt vertically ( m^2/s ) */')
124          WRITE(msgBuf,'(2A)') ' Equation of State : eosType = ', eosType
125          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
126         &                    SQUEEZE_RIGHT , 1)
127        CALL WRITE_0D_R8( tAlpha,  INDEX_NONE,'tAlpha =',        CALL WRITE_0D_R8( tAlpha,  INDEX_NONE,'tAlpha =',
128       &'   /* Linear EOS thermal expansion coefficient ( 1/degree ) */')       &'   /* Linear EOS thermal expansion coefficient ( 1/degree ) */')
129        CALL WRITE_0D_R8( sBeta,   INDEX_NONE,'sBeta =',        CALL WRITE_0D_R8( sBeta,   INDEX_NONE,'sBeta =',
# Line 132  CEOP Line 138  CEOP
138       &                    SQUEEZE_RIGHT , 1)       &                    SQUEEZE_RIGHT , 1)
139          ENDDO          ENDDO
140        ENDIF        ENDIF
141          IF ( buoyancyRelation .eq. 'ATMOSPHERIC' ) THEN
142           CALL WRITE_0D_R8( atm_Rd, INDEX_NONE, 'atm_Rd =',
143         & '  /* gas constant for dry air ( J/kg/K ) */')
144          CALL WRITE_0D_R8( atm_Cp, INDEX_NONE, 'atm_Cp =',
145         & '  /* specific heat (Cp) of dry air ( J/kg/K ) */')
146          CALL WRITE_0D_R8( atm_kappa, INDEX_NONE, 'atm_kappa =',
147         & '  /* kappa (=Rd/Cp ) of dry air */')
148           CALL WRITE_0D_R8( atm_Rq, INDEX_NONE, 'atm_Rq =',
149         &  ' /* water vap. specific vol. anomaly relative to dry air */')
150          CALL WRITE_0D_R8( atm_Po, INDEX_NONE, 'atm_Po =',
151         & '  /* standard reference pressure ( Pa ) */')
152          CALL WRITE_0D_I( integr_GeoPot, INDEX_NONE, 'integr_GeoPot =',
153         & '  /* select how the geopotential is integrated */')
154          CALL WRITE_0D_I( selectFindRoSurf, INDEX_NONE,
155         & 'selectFindRoSurf=',
156         & '  /* select how Surf.Ref. pressure is defined */')
157          ENDIF
158        CALL WRITE_0D_R8( rhonil,  INDEX_NONE,'rhonil =',        CALL WRITE_0D_R8( rhonil,  INDEX_NONE,'rhonil =',
159       &'   /* Reference density ( kg/m^3 ) */')       &'   /* Reference density ( kg/m^3 ) */')
160        CALL WRITE_0D_R8( rhoConst, INDEX_NONE,'rhoConst =',        CALL WRITE_0D_R8( rhoConst, INDEX_NONE,'rhoConst =',
161       &'   /* Reference density ( kg/m^3 ) */')       &'   /* Reference density ( kg/m^3 ) */')
162          CALL WRITE_0D_R8( rhoConstFresh, INDEX_NONE,'rhoConstFresh =',
163         &'   /* Reference density ( kg/m^3 ) */')
164        CALL WRITE_0D_R8( gravity, INDEX_NONE,'gravity =',        CALL WRITE_0D_R8( gravity, INDEX_NONE,'gravity =',
165       &'   /* Gravitational acceleration ( m/s^2 ) */')       &'   /* Gravitational acceleration ( m/s^2 ) */')
166        CALL WRITE_0D_R8( gBaro,   INDEX_NONE,'gBaro =',        CALL WRITE_0D_R8( gBaro,   INDEX_NONE,'gBaro =',
167       &'   /* Barotropic gravity ( m/s^2 ) */')       &'   /* Barotropic gravity ( m/s^2 ) */')
168          CALL WRITE_0D_R8(rotationPeriod,INDEX_NONE,'rotationPeriod =',
169         &'   /* Rotation Period ( s ) */')
170          CALL WRITE_0D_R8( omega,   INDEX_NONE,'omega =',
171         &'   /* Angular velocity ( rad/s ) */')
172        CALL WRITE_0D_R8( f0,      INDEX_NONE,'f0 =',        CALL WRITE_0D_R8( f0,      INDEX_NONE,'f0 =',
173       &'   /* Reference coriolis parameter ( 1/s ) */')       &'   /* Reference coriolis parameter ( 1/s ) */')
174        CALL WRITE_0D_R8( beta,    INDEX_NONE,'beta =',        CALL WRITE_0D_R8( beta,    INDEX_NONE,'beta =',
# Line 178  CEOP Line 207  CEOP
207        CALL WRITE_0D_R8( hFacSup, INDEX_NONE,        CALL WRITE_0D_R8( hFacSup, INDEX_NONE,
208       &'hFacSup =',       &'hFacSup =',
209       &'   /* upper threshold for hFac (nonlinFreeSurf only)*/')       &'   /* upper threshold for hFac (nonlinFreeSurf only)*/')
210          CALL WRITE_0D_I( select_rStar, INDEX_NONE,
211         &'select_rStar =',
212         &'   /* r* Coordinate options (not yet implemented)*/')
213        CALL WRITE_0D_L( useRealFreshWaterFlux, INDEX_NONE,        CALL WRITE_0D_L( useRealFreshWaterFlux, INDEX_NONE,
214       &'useRealFreshWaterFlux =',       &'useRealFreshWaterFlux =',
215       &'   /* Real Fresh Water Flux on/off flag*/')       &'   /* Real Fresh Water Flux on/off flag*/')
# Line 230  CEOP Line 262  CEOP
262       & 'saltForcing  =', '  /* Salinity forcing on/off flag */')       & 'saltForcing  =', '  /* Salinity forcing on/off flag */')
263        CALL WRITE_0D_L( nonHydrostatic, INDEX_NONE,        CALL WRITE_0D_L( nonHydrostatic, INDEX_NONE,
264       & 'nonHydrostatic =', '  /* Non-Hydrostatic on/off flag */')       & 'nonHydrostatic =', '  /* Non-Hydrostatic on/off flag */')
265          CALL WRITE_0D_L( useCDscheme,  INDEX_NONE,
266         & 'useCDscheme =', '  /* CD scheme on/off flag */')
267          CALL WRITE_0D_L( useJamartWetPoints,  INDEX_NONE,
268         & 'useJamartWetPoints=',' /* Coriolis WetPoints method flag */')
269        WRITE(msgBuf,'(A)') '//  '        WRITE(msgBuf,'(A)') '//  '
270        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
271       &                    SQUEEZE_RIGHT , 1)       &                    SQUEEZE_RIGHT , 1)
# Line 265  CEOP Line 301  CEOP
301       &'   /* Number of timesteps */')       &'   /* Number of timesteps */')
302        CALL WRITE_0D_R8( deltaTmom, INDEX_NONE,'deltatTmom =',        CALL WRITE_0D_R8( deltaTmom, INDEX_NONE,'deltatTmom =',
303       &'   /* Momentum equation timestep ( s ) */')       &'   /* Momentum equation timestep ( s ) */')
304          CALL WRITE_0D_R8( deltaTfreesurf,INDEX_NONE,'deltaTfreesurf =',
305         &'   /* FreeSurface equation timestep ( s ) */')
306        CALL WRITE_0D_R8( deltaTtracer, INDEX_NONE,'deltatTtracer =',        CALL WRITE_0D_R8( deltaTtracer, INDEX_NONE,'deltatTtracer =',
307       &'   /* Tracer equation timestep ( s ) */')       &'   /* Tracer equation timestep ( s ) */')
308        CALL WRITE_0D_R8( deltaTClock, INDEX_NONE,'deltatTClock  =',        CALL WRITE_0D_R8( deltaTClock, INDEX_NONE,'deltatTClock  =',
# Line 275  CEOP Line 313  CEOP
313       &'   /* put T,S Forcing in Adams-Bash. stepping */')       &'   /* put T,S Forcing in Adams-Bash. stepping */')
314        CALL WRITE_0D_R8( abeps, INDEX_NONE,'abeps =',        CALL WRITE_0D_R8( abeps, INDEX_NONE,'abeps =',
315       &'   /* Adams-Bashforth stabilizing weight */')       &'   /* Adams-Bashforth stabilizing weight */')
316          IF (useCDscheme) THEN
317        CALL WRITE_0D_R8( tauCD, INDEX_NONE,'tauCD =',        CALL WRITE_0D_R8( tauCD, INDEX_NONE,'tauCD =',
318       &'   /* CD coupling time-scale ( s ) */')       &'   /* CD coupling time-scale ( s ) */')
319        CALL WRITE_0D_R8( rCD, INDEX_NONE,'rCD =',        CALL WRITE_0D_R8( rCD, INDEX_NONE,'rCD =',
320       &'   /* Normalised CD coupling parameter */')       &'   /* Normalised CD coupling parameter */')
321          ENDIF
322        CALL WRITE_0D_R8( startTime, INDEX_NONE,'startTime =',        CALL WRITE_0D_R8( startTime, INDEX_NONE,'startTime =',
323       &'   /* Run start time ( s ). */')       &'   /* Run start time ( s ). */')
324        CALL WRITE_0D_R8( endTime, INDEX_NONE,'endTime =',        CALL WRITE_0D_R8( endTime, INDEX_NONE,'endTime =',
# Line 289  CEOP Line 329  CEOP
329       &'   /* Rolling restart/checkpoint file interval ( s ). */')       &'   /* Rolling restart/checkpoint file interval ( s ). */')
330        CALL WRITE_0D_R8( dumpFreq, INDEX_NONE,'dumpFreq =',        CALL WRITE_0D_R8( dumpFreq, INDEX_NONE,'dumpFreq =',
331       &'   /* Model state write out interval ( s ). */')       &'   /* Model state write out interval ( s ). */')
332          CALL WRITE_0D_R8( externForcingPeriod, INDEX_NONE,
333         &   'externForcingPeriod =', '   /* forcing period (s) */')
334          CALL WRITE_0D_R8( externForcingCycle, INDEX_NONE,
335         &   'externForcingCycle =', '   /* period of the cyle (s). */')
336          CALL WRITE_0D_R8( tauThetaClimRelax, INDEX_NONE,
337         &   'tauThetaClimRelax =', '   /* relaxation time scale (s) */')
338          CALL WRITE_0D_R8( tauSaltClimRelax, INDEX_NONE,
339         &   'tauSaltClimRelax =',  '   /* relaxation time scale (s) */')
340          CALL WRITE_0D_R8( latBandClimRelax, INDEX_NONE,
341         &   'latBandClimRelax =', '   /* max. Lat. where relaxation */')
342        WRITE(msgBuf,'(A)') '//  '        WRITE(msgBuf,'(A)') '//  '
343        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
344       &                    SQUEEZE_RIGHT , 1)       &                    SQUEEZE_RIGHT , 1)
# Line 308  CEOP Line 357  CEOP
357       &'   /* Spherical coordinates flag ( True / False ) */')       &'   /* Spherical coordinates flag ( True / False ) */')
358        CALL WRITE_0D_L( groundAtK1, INDEX_NONE, 'groundAtK1 =',        CALL WRITE_0D_L( groundAtK1, INDEX_NONE, 'groundAtK1 =',
359       &'   /* Lower Boundary (ground) at the surface(k=1) ( T / F ) */')       &'   /* Lower Boundary (ground) at the surface(k=1) ( T / F ) */')
360        CALL WRITE_1D_R8( Ro_SeaLevel,1, INDEX_NONE,'Ro_SeaLevel =',        CALL WRITE_0D_R8( Ro_SeaLevel, INDEX_NONE,'Ro_SeaLevel =',
361       &'   /* r(1) ( units of r ) */')       &'   /* r(1) ( units of r ) */')
362        CALL WRITE_1D_R8( rkFac,1, INDEX_NONE,'rkFac =',        CALL WRITE_0D_R8( rkFac, INDEX_NONE,'rkFac =',
363       &'   /* minus Vertical index orientation  */')       &'   /* minus Vertical index orientation  */')
364        CALL WRITE_1D_R8( horiVertRatio,1, INDEX_NONE,'horiVertRatio =',        CALL WRITE_0D_R8( horiVertRatio, INDEX_NONE,'horiVertRatio =',
365       &'   /* Ratio on units : Horiz - Vertical */')       &'   /* Ratio on units : Horiz - Vertical */')
366  c     CALL WRITE_1D_R8( delZ,Nr, INDEX_K,'delZ = ',  c     CALL WRITE_1D_R8( delZ,Nr, INDEX_K,'delZ = ',
367  c    &'   /* W spacing ( m ) */')  c    &'   /* W spacing ( m ) */')
Legend:
Removed from v.1.34  
changed lines
  Added in v.1.44
  ViewVC Help
Powered by ViewVC 1.1.22