/[MITgcm]/MITgcm_contrib/osse/codemod/config_summary.F
ViewVC logotype

Diff of /MITgcm_contrib/osse/codemod/config_summary.F

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

revision 1.1 by afe, Tue Jun 22 19:44:40 2004 UTC revision 1.3 by afe, Thu Jun 24 17:52:38 2004 UTC
# Line 91  CEOP Line 91  CEOP
91       &'   /* Reference salinity profile ( ppt ) */')       &'   /* Reference salinity profile ( ppt ) */')
92        CALL WRITE_0D_R8( viscAh, INDEX_NONE,'viscAh =',        CALL WRITE_0D_R8( viscAh, INDEX_NONE,'viscAh =',
93       &'   /* Lateral eddy viscosity ( m^2/s ) */')       &'   /* Lateral eddy viscosity ( m^2/s ) */')
94        CALL WRITE_0D_R8( viscA4, INDEX_NONE,'viscAh =',        CALL WRITE_0D_R8( viscAhMax, INDEX_NONE,'viscAhMax =',
95         &'   /* Maximum lateral eddy viscosity ( m^2/s ) */')
96          CALL WRITE_0D_R8( viscAhGrid, INDEX_NONE,'viscAhGrid =',
97         &'   /* Grid dependent lateral eddy viscosity ( non-dim. ) */')
98          CALL WRITE_0D_R8( viscC2leith, INDEX_NONE,'viscC2leith =',
99         &'   /* Leith harmonic viscosity factor ( non-dom. ) */')
100          CALL WRITE_0D_R8( viscA4, INDEX_NONE,'viscA4 =',
101       &'   /* Lateral biharmonic viscosity ( m^4/s ) */')       &'   /* Lateral biharmonic viscosity ( m^4/s ) */')
102          CALL WRITE_0D_R8( viscA4Max, INDEX_NONE,'viscA4Max =',
103         &'   /* Maximum biharmonic viscosity ( m^2/s ) */')
104          CALL WRITE_0D_R8( viscA4Grid, INDEX_NONE,'viscA4Grid =',
105         &'   /* Grid dependent biharmonic viscosity ( non-dim. ) */')
106          CALL WRITE_0D_R8( viscC4leith, INDEX_NONE,'viscC4leith =',
107         &'   /* Leith biharmonic viscosity factor ( non-dom. ) */')
108        CALL WRITE_0D_L( no_slip_sides, INDEX_NONE,        CALL WRITE_0D_L( no_slip_sides, INDEX_NONE,
109       & 'no_slip_sides =', '  /* Viscous BCs: No-slip sides */')       & 'no_slip_sides =', '  /* Viscous BCs: No-slip sides */')
110        IF ( viscAz .NE. UNSET_RL ) THEN  c     IF ( viscAz .NE. UNSET_RL ) THEN
111         CALL WRITE_0D_R8( viscAz, INDEX_NONE,'viscAz =',  c      CALL WRITE_0D_R8( viscAz, INDEX_NONE,'viscAz =',
112       & '   /* Vertical eddy viscosity ( m^2/s ) */')  c    & '   /* Vertical eddy viscosity ( m^2/s ) */')
113        ENDIF  c     ENDIF
114        IF ( viscAp .NE. UNSET_RL ) THEN  c     IF ( viscAp .NE. UNSET_RL ) THEN
115         CALL WRITE_0D_R8( viscAp, INDEX_NONE,'viscAp =',  c      CALL WRITE_0D_R8( viscAp, INDEX_NONE,'viscAp =',
116       & '   /* Vertical eddy viscosity ( Pa^2/s ) */')  c    & '   /* Vertical eddy viscosity ( Pa^2/s ) */')
117        ENDIF  c     ENDIF
118        CALL WRITE_0D_R8( viscAr,  INDEX_NONE,'viscAr =',        CALL WRITE_0D_R8( viscAr,  INDEX_NONE,'viscAr =',
119       &'   /* Vertical eddy viscosity ( units of r^2/s ) */')       &'   /* Vertical eddy viscosity ( units of r^2/s ) */')
120        CALL WRITE_0D_R8( diffKhT, INDEX_NONE,'diffKhT =',        CALL WRITE_0D_R8( diffKhT, INDEX_NONE,'diffKhT =',
121       &'   /* Laplacian diffusion of heat laterally ( m^2/s ) */')       &'   /* Laplacian diffusion of heat laterally ( m^2/s ) */')
122        CALL WRITE_0D_R8( diffK4T, INDEX_NONE,'diffK4T =',        CALL WRITE_0D_R8( diffK4T, INDEX_NONE,'diffK4T =',
123       &'   /* Bihaarmonic diffusion of heat laterally ( m^4/s ) */')       &'   /* Bihaarmonic diffusion of heat laterally ( m^4/s ) */')
124        CALL WRITE_0D_R8( diffKzT, INDEX_NONE,'diffKzT =',  c     CALL WRITE_0D_R8( diffKzT, INDEX_NONE,'diffKzT =',
125       &'   /* Laplacian diffusion of heat vertically ( m^2/s ) */')  c    &'   /* Laplacian diffusion of heat vertically ( m^2/s ) */')
126        CALL WRITE_0D_R8( diffKrT, INDEX_NONE,'diffKrT =',        CALL WRITE_0D_R8( diffKrT, INDEX_NONE,'diffKrT =',
127       &'   /* Laplacian diffusion of heat vertically ( m^2/s ) */')       &'   /* Laplacian diffusion of heat vertically ( m^2/s ) */')
128        CALL WRITE_0D_R8( diffKhS, INDEX_NONE,'diffKhS =',        CALL WRITE_0D_R8( diffKhS, INDEX_NONE,'diffKhS =',
129       &'   /* Laplacian diffusion of salt laterally ( m^2/s ) */')       &'   /* Laplacian diffusion of salt laterally ( m^2/s ) */')
130        CALL WRITE_0D_R8( diffK4S, INDEX_NONE,'diffK4S =',        CALL WRITE_0D_R8( diffK4S, INDEX_NONE,'diffK4S =',
131       &'   /* Bihaarmonic diffusion of salt laterally ( m^4/s ) */')       &'   /* Bihaarmonic diffusion of salt laterally ( m^4/s ) */')
132        CALL WRITE_0D_R8( diffKzS, INDEX_NONE,'diffKzS =',  c     CALL WRITE_0D_R8( diffKzS, INDEX_NONE,'diffKzS =',
133       &'   /* Laplacian diffusion of salt vertically ( m^2/s ) */')  c    &'   /* Laplacian diffusion of salt vertically ( m^2/s ) */')
134        CALL WRITE_0D_R8( diffKrS, INDEX_NONE,'diffKrS =',        CALL WRITE_0D_R8( diffKrS, INDEX_NONE,'diffKrS =',
135       &'   /* Laplacian diffusion of salt vertically ( m^2/s ) */')       &'   /* Laplacian diffusion of salt vertically ( m^2/s ) */')
136          CALL WRITE_0D_R8( diffKrBL79surf, INDEX_NONE,'diffKrBL79surf =',
137         &'   /* Surface diffusion for Bryan and Lewis 1979 ( m^2/s ) */')
138          CALL WRITE_0D_R8( diffKrBL79deep, INDEX_NONE,'diffKrBL79deep =',
139         &'   /* Deep diffusion for Bryan and Lewis 1979 ( m^2/s ) */')
140          CALL WRITE_0D_R8( diffKrBL79scl, INDEX_NONE,'diffKrBL79scl =',
141         &'   /* Depth scale for Bryan and Lewis 1979 ( m ) */')
142          CALL WRITE_0D_R8( diffKrBL79Ho, INDEX_NONE,'diffKrBL79Ho =',
143         &'   /* Turning depth for Bryan and Lewis 1979 ( m ) */')
144        WRITE(msgBuf,'(2A)') ' Equation of State : eosType = ', eosType        WRITE(msgBuf,'(2A)') ' Equation of State : eosType = ', eosType
145        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
146       &                    SQUEEZE_RIGHT , 1)       &                    SQUEEZE_RIGHT , 1)
# Line 145  CEOP Line 165  CEOP
165       & '  /* specific heat (Cp) of dry air ( J/kg/K ) */')       & '  /* specific heat (Cp) of dry air ( J/kg/K ) */')
166        CALL WRITE_0D_R8( atm_kappa, INDEX_NONE, 'atm_kappa =',        CALL WRITE_0D_R8( atm_kappa, INDEX_NONE, 'atm_kappa =',
167       & '  /* kappa (=Rd/Cp ) of dry air */')       & '  /* kappa (=Rd/Cp ) of dry air */')
168           CALL WRITE_0D_R8( atm_Rq, INDEX_NONE, 'atm_Rq =',
169         &  ' /* water vap. specific vol. anomaly relative to dry air */')
170        CALL WRITE_0D_R8( atm_Po, INDEX_NONE, 'atm_Po =',        CALL WRITE_0D_R8( atm_Po, INDEX_NONE, 'atm_Po =',
171       & '  /* standard reference pressure ( Pa ) */')       & '  /* standard reference pressure ( Pa ) */')
172        CALL WRITE_0D_I( integr_GeoPot, INDEX_NONE, 'integr_GeoPot =',        CALL WRITE_0D_I( integr_GeoPot, INDEX_NONE, 'integr_GeoPot =',
# Line 163  CEOP Line 185  CEOP
185       &'   /* Gravitational acceleration ( m/s^2 ) */')       &'   /* Gravitational acceleration ( m/s^2 ) */')
186        CALL WRITE_0D_R8( gBaro,   INDEX_NONE,'gBaro =',        CALL WRITE_0D_R8( gBaro,   INDEX_NONE,'gBaro =',
187       &'   /* Barotropic gravity ( m/s^2 ) */')       &'   /* Barotropic gravity ( m/s^2 ) */')
188          CALL WRITE_0D_R8(rotationPeriod,INDEX_NONE,'rotationPeriod =',
189         &'   /* Rotation Period ( s ) */')
190          CALL WRITE_0D_R8( omega,   INDEX_NONE,'omega =',
191         &'   /* Angular velocity ( rad/s ) */')
192        CALL WRITE_0D_R8( f0,      INDEX_NONE,'f0 =',        CALL WRITE_0D_R8( f0,      INDEX_NONE,'f0 =',
193       &'   /* Reference coriolis parameter ( 1/s ) */')       &'   /* Reference coriolis parameter ( 1/s ) */')
194        CALL WRITE_0D_R8( beta,    INDEX_NONE,'beta =',        CALL WRITE_0D_R8( beta,    INDEX_NONE,'beta =',
# Line 201  CEOP Line 227  CEOP
227        CALL WRITE_0D_R8( hFacSup, INDEX_NONE,        CALL WRITE_0D_R8( hFacSup, INDEX_NONE,
228       &'hFacSup =',       &'hFacSup =',
229       &'   /* upper threshold for hFac (nonlinFreeSurf only)*/')       &'   /* upper threshold for hFac (nonlinFreeSurf only)*/')
230          CALL WRITE_0D_I( select_rStar, INDEX_NONE,
231         &'select_rStar =',
232         &'   /* r* Coordinate options (not yet implemented)*/')
233        CALL WRITE_0D_L( useRealFreshWaterFlux, INDEX_NONE,        CALL WRITE_0D_L( useRealFreshWaterFlux, INDEX_NONE,
234       &'useRealFreshWaterFlux =',       &'useRealFreshWaterFlux =',
235       &'   /* Real Fresh Water Flux on/off flag*/')       &'   /* Real Fresh Water Flux on/off flag*/')
# Line 220  CEOP Line 249  CEOP
249       &' /* convert F.W. Flux to Salt Flux (-1=use local S)(ppt)*/')       &' /* convert F.W. Flux to Salt Flux (-1=use local S)(ppt)*/')
250        ENDIF        ENDIF
251    
252        CALL WRITE_0D_L( multiDimAdvection, INDEX_NONE,        CALL WRITE_0D_L( nonHydrostatic, INDEX_NONE,
253       & 'multiDimAdvection =',       & 'nonHydrostatic =', '  /* Non-Hydrostatic on/off flag */')
      &'   /* enable/disable Multi-Dim Advection */')  
       CALL WRITE_0D_L( staggerTimeStep, INDEX_NONE,  
      &                 'staggerTimeStep =',  
      &'   /* Stagger time stepping on/off flag */')  
254        CALL WRITE_0D_L( momStepping,  INDEX_NONE,        CALL WRITE_0D_L( momStepping,  INDEX_NONE,
255       & 'momStepping =', '  /* Momentum equation on/off flag */')       & 'momStepping =', '  /* Momentum equation on/off flag */')
256        CALL WRITE_0D_L( momAdvection, INDEX_NONE,        CALL WRITE_0D_L( momAdvection, INDEX_NONE,
257       & 'momAdvection =', '  /* Momentum advection on/off flag */')       & 'momAdvection =', '  /* Momentum advection on/off flag */')
258        CALL WRITE_0D_L( momViscosity, INDEX_NONE,        CALL WRITE_0D_L( momViscosity, INDEX_NONE,
259       & 'momViscosity =', '  /* Momentum viscosity on/off flag */')       & 'momViscosity =', '  /* Momentum viscosity on/off flag */')
260          CALL WRITE_0D_L( momImplVertAdv, INDEX_NONE, 'momImplVertAdv =',
261         &                '/* Momentum implicit vert. advection on/off*/')
262          CALL WRITE_0D_L( implicitViscosity, INDEX_NONE,
263         & 'implicitViscosity =', ' /* Implicit viscosity on/off flag */')
264        CALL WRITE_0D_L( useCoriolis,  INDEX_NONE,        CALL WRITE_0D_L( useCoriolis,  INDEX_NONE,
265       & 'useCoriolis =', '  /* Coriolis on/off flag */')       & 'useCoriolis =', '  /* Coriolis on/off flag */')
266          CALL WRITE_0D_L( useCDscheme,  INDEX_NONE,
267         & 'useCDscheme =', '  /* CD scheme on/off flag */')
268          CALL WRITE_0D_L( useJamartWetPoints,  INDEX_NONE,
269         & 'useJamartWetPoints=',' /* Coriolis WetPoints method flag */')
270          CALL WRITE_0D_L( useJamartMomAdv,  INDEX_NONE,
271         & 'useJamartMomAdv=',' /* V.I. Non-linear terms Jamart flag */')
272          CALL WRITE_0D_L( SadournyCoriolis,  INDEX_NONE,
273         & 'SadournyCoriolis=',' /* Sadourny Coriolis discr. flag */')
274          CALL WRITE_0D_L( upwindVorticity,  INDEX_NONE,
275         & 'upwindVorticity=',' /* Upwind bias vorticity flag */')
276          CALL WRITE_0D_L( useAbsVorticity,  INDEX_NONE,
277         & 'useAbsVorticity=',' /* Work with f+zeta in Coriolis */')
278          CALL WRITE_0D_L( highOrderVorticity,  INDEX_NONE,
279         & 'highOrderVorticity=',' /* High order interp. of vort. flag */')
280        CALL WRITE_0D_L( momForcing,   INDEX_NONE,        CALL WRITE_0D_L( momForcing,   INDEX_NONE,
281       & 'momForcing =', '  /* Momentum forcing on/off flag */')       & 'momForcing =', '  /* Momentum forcing on/off flag */')
282        CALL WRITE_0D_L( momPressureForcing, INDEX_NONE,        CALL WRITE_0D_L( momPressureForcing, INDEX_NONE,
283       & 'momPressureForcing =',         & 'momPressureForcing =',  
284       & '  /* Momentum pressure term on/off flag */')       & '  /* Momentum pressure term on/off flag */')
285          CALL WRITE_0D_L( staggerTimeStep, INDEX_NONE,
286         &                 'staggerTimeStep =',
287         &'   /* Stagger time stepping on/off flag */')
288          CALL WRITE_0D_L( multiDimAdvection, INDEX_NONE,
289         & 'multiDimAdvection =',
290         &'   /* enable/disable Multi-Dim Advection */')
291          CALL WRITE_0D_L( implicitDiffusion, INDEX_NONE,
292         & 'implicitDiffusion =','/* Implicit Diffusion on/off flag */')
293        CALL WRITE_0D_L( tempStepping,  INDEX_NONE,        CALL WRITE_0D_L( tempStepping,  INDEX_NONE,
294       & 'tempStepping =', '  /* Temperature equation on/off flag */')       & 'tempStepping =', '  /* Temperature equation on/off flag */')
295        CALL WRITE_0D_L( tempAdvection,  INDEX_NONE,        CALL WRITE_0D_L( tempAdvection,  INDEX_NONE,
296       & 'tempAdvection=', '  /* Temperature advection on/off flag */')       & 'tempAdvection=', '  /* Temperature advection on/off flag */')
297          CALL WRITE_0D_L( tempImplVertAdv,INDEX_NONE,'tempImplVertAdv =',
298         &                '/* Temp. implicit vert. advection on/off */')
299        CALL WRITE_0D_L( tempForcing,  INDEX_NONE,        CALL WRITE_0D_L( tempForcing,  INDEX_NONE,
300       & 'tempForcing  =', '  /* Temperature forcing on/off flag */')       & 'tempForcing  =', '  /* Temperature forcing on/off flag */')
301        CALL WRITE_0D_L( saltStepping,  INDEX_NONE,        CALL WRITE_0D_L( saltStepping,  INDEX_NONE,
302       & 'saltStepping =', '  /* Salinity equation on/off flag */')       & 'saltStepping =', '  /* Salinity equation on/off flag */')
303        CALL WRITE_0D_L( saltAdvection,  INDEX_NONE,        CALL WRITE_0D_L( saltAdvection,  INDEX_NONE,
304       & 'saltAdvection=', '  /* Salinity advection on/off flag */')       & 'saltAdvection=', '  /* Salinity advection on/off flag */')
305          CALL WRITE_0D_L( saltImplVertAdv,INDEX_NONE,'saltImplVertAdv =',
306         &                '/* Sali. implicit vert. advection on/off */')
307        CALL WRITE_0D_L( saltForcing,  INDEX_NONE,        CALL WRITE_0D_L( saltForcing,  INDEX_NONE,
308       & 'saltForcing  =', '  /* Salinity forcing on/off flag */')       & 'saltForcing  =', '  /* Salinity forcing on/off flag */')
       CALL WRITE_0D_L( nonHydrostatic, INDEX_NONE,  
      & 'nonHydrostatic =', '  /* Non-Hydrostatic on/off flag */')  
309        WRITE(msgBuf,'(A)') '//  '        WRITE(msgBuf,'(A)') '//  '
310        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
311       &                    SQUEEZE_RIGHT , 1)       &                    SQUEEZE_RIGHT , 1)
# Line 300  CEOP Line 353  CEOP
353       &'   /* put T,S Forcing in Adams-Bash. stepping */')       &'   /* put T,S Forcing in Adams-Bash. stepping */')
354        CALL WRITE_0D_R8( abeps, INDEX_NONE,'abeps =',        CALL WRITE_0D_R8( abeps, INDEX_NONE,'abeps =',
355       &'   /* Adams-Bashforth stabilizing weight */')       &'   /* Adams-Bashforth stabilizing weight */')
356          IF (useCDscheme) THEN
357        CALL WRITE_0D_R8( tauCD, INDEX_NONE,'tauCD =',        CALL WRITE_0D_R8( tauCD, INDEX_NONE,'tauCD =',
358       &'   /* CD coupling time-scale ( s ) */')       &'   /* CD coupling time-scale ( s ) */')
359        CALL WRITE_0D_R8( rCD, INDEX_NONE,'rCD =',        CALL WRITE_0D_R8( rCD, INDEX_NONE,'rCD =',
360       &'   /* Normalised CD coupling parameter */')       &'   /* Normalised CD coupling parameter */')
361          ENDIF
362        CALL WRITE_0D_R8( startTime, INDEX_NONE,'startTime =',        CALL WRITE_0D_R8( startTime, INDEX_NONE,'startTime =',
363       &'   /* Run start time ( s ). */')       &'   /* Run start time ( s ). */')
364        CALL WRITE_0D_R8( endTime, INDEX_NONE,'endTime =',        CALL WRITE_0D_R8( endTime, INDEX_NONE,'endTime =',
# Line 314  CEOP Line 369  CEOP
369       &'   /* Rolling restart/checkpoint file interval ( s ). */')       &'   /* Rolling restart/checkpoint file interval ( s ). */')
370        CALL WRITE_0D_R8( dumpFreq, INDEX_NONE,'dumpFreq =',        CALL WRITE_0D_R8( dumpFreq, INDEX_NONE,'dumpFreq =',
371       &'   /* Model state write out interval ( s ). */')       &'   /* Model state write out interval ( s ). */')
372          CALL WRITE_0D_R8( externForcingPeriod, INDEX_NONE,
373         &   'externForcingPeriod =', '   /* forcing period (s) */')
374          CALL WRITE_0D_R8( externForcingCycle, INDEX_NONE,
375         &   'externForcingCycle =', '   /* period of the cyle (s). */')
376          CALL WRITE_0D_R8( tauThetaClimRelax, INDEX_NONE,
377         &   'tauThetaClimRelax =', '   /* relaxation time scale (s) */')
378          CALL WRITE_0D_R8( tauSaltClimRelax, INDEX_NONE,
379         &   'tauSaltClimRelax =',  '   /* relaxation time scale (s) */')
380          CALL WRITE_0D_R8( latBandClimRelax, INDEX_NONE,
381         &   'latBandClimRelax =', '   /* max. Lat. where relaxation */')
382        WRITE(msgBuf,'(A)') '//  '        WRITE(msgBuf,'(A)') '//  '
383        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
384       &                    SQUEEZE_RIGHT , 1)       &                    SQUEEZE_RIGHT , 1)
# Line 331  CEOP Line 395  CEOP
395        CALL WRITE_0D_L( usingSphericalPolarGrid, INDEX_NONE,        CALL WRITE_0D_L( usingSphericalPolarGrid, INDEX_NONE,
396       & 'usingSphericalPolarGrid =',       & 'usingSphericalPolarGrid =',
397       &'   /* Spherical coordinates flag ( True / False ) */')       &'   /* Spherical coordinates flag ( True / False ) */')
398        CALL WRITE_0D_L( bUseCylindricalGrid, INDEX_NONE,        CALL WRITE_0D_L( usingCylindricalGrid, INDEX_NONE,
399       & 'usingCylindricalGrid =',       & 'usingCylindricalGrid =',
400       &'   /* Spherical coordinates flag ( True / False ) */')       &'   /* Spherical coordinates flag ( True / False ) */')
401        CALL WRITE_0D_L( groundAtK1, INDEX_NONE, 'groundAtK1 =',        CALL WRITE_0D_L( groundAtK1, INDEX_NONE, 'groundAtK1 =',
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22