71 |
&' /* Reference salinity profile ( ppt ) */') |
&' /* Reference salinity profile ( ppt ) */') |
72 |
CALL WRITE_1D_R8( viscAh, 1, INDEX_NONE,'viscAh =', |
CALL WRITE_1D_R8( viscAh, 1, 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,'viscAh =', |
CALL WRITE_1D_R8( viscA4, 1, INDEX_NONE,'viscA4 =', |
75 |
&' /* Lateral biharmonic viscosity ( m^4/s ) */') |
&' /* Lateral biharmonic viscosity ( m^4/s ) */') |
76 |
|
CALL WRITE_1D_L( no_slip_sides,1, INDEX_NONE, |
77 |
|
& '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_1D_R8( viscAz, 1, INDEX_NONE,'viscAz =', |
80 |
& ' /* Vertical eddy viscosity ( m^2/s ) */') |
& ' /* Vertical eddy viscosity ( m^2/s ) */') |
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_1D_R8( diffKhT, 1, INDEX_NONE,'diffKhT =', |
89 |
&' /* Laplacian diffusion of heat laterally ( m^2/s ) */') |
&' /* Laplacian diffusion of heat laterally ( m^2/s ) */') |
|
CALL WRITE_1D_R8( diffKzT, 1, INDEX_NONE,'diffKzT =', |
|
|
&' /* Laplacian diffusion of heat vertically ( m^2/s ) */') |
|
90 |
CALL WRITE_1D_R8( diffK4T, 1, INDEX_NONE,'diffK4T =', |
CALL WRITE_1D_R8( diffK4T, 1, 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 =', |
93 |
|
&' /* Laplacian diffusion of heat vertically ( m^2/s ) */') |
94 |
|
CALL WRITE_1D_R8( diffKrT, 1, INDEX_NONE,'diffKrT =', |
95 |
|
&' /* Laplacian diffusion of heat vertically ( m^2/s ) */') |
96 |
CALL WRITE_1D_R8( diffKhS, 1, INDEX_NONE,'diffKhS =', |
CALL WRITE_1D_R8( diffKhS, 1, INDEX_NONE,'diffKhS =', |
97 |
&' /* Laplacian diffusion of salt laterally ( m^2/s ) */') |
&' /* Laplacian diffusion of salt laterally ( m^2/s ) */') |
|
CALL WRITE_1D_R8( diffKzS, 1, INDEX_NONE,'diffKzS =', |
|
|
&' /* Laplacian diffusion of salt vertically ( m^2/s ) */') |
|
98 |
CALL WRITE_1D_R8( diffK4S, 1, INDEX_NONE,'diffK4S =', |
CALL WRITE_1D_R8( diffK4S, 1, 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 =', |
101 |
|
&' /* Laplacian diffusion of salt vertically ( m^2/s ) */') |
102 |
|
CALL WRITE_1D_R8( diffKrS, 1, INDEX_NONE,'diffKrS =', |
103 |
|
&' /* Laplacian diffusion of salt vertically ( m^2/s ) */') |
104 |
CALL WRITE_1D_R8( tAlpha,1, INDEX_NONE,'tAlpha =', |
CALL WRITE_1D_R8( tAlpha,1, 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_1D_R8( sBeta, 1, INDEX_NONE,'sBeta =', |
117 |
ENDIF |
ENDIF |
118 |
CALL WRITE_1D_R8( rhonil,1, INDEX_NONE,'rhonil =', |
CALL WRITE_1D_R8( rhonil,1, INDEX_NONE,'rhonil =', |
119 |
&' /* Reference density ( kg/m^3 ) */') |
&' /* Reference density ( kg/m^3 ) */') |
120 |
|
CALL WRITE_1D_R8( rhoConst,1, INDEX_NONE,'rhoConst =', |
121 |
|
&' /* Reference density ( kg/m^3 ) */') |
122 |
CALL WRITE_1D_R8( gravity,1, INDEX_NONE,'gravity =', |
CALL WRITE_1D_R8( gravity,1, 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_1D_R8( gBaro,1, INDEX_NONE,'gBaro =', |
150 |
& ' /* Momentum pressure term on/off flag */') |
& ' /* Momentum pressure term on/off flag */') |
151 |
CALL WRITE_1D_L( tempStepping,1, INDEX_NONE, |
CALL WRITE_1D_L( tempStepping,1, INDEX_NONE, |
152 |
& 'tempStepping =', ' /* Temperature equation on/off flag */') |
& 'tempStepping =', ' /* Temperature equation on/off flag */') |
153 |
|
CALL WRITE_1D_L( openBoundaries,1, INDEX_NONE, |
154 |
|
& 'openBoundaries =', ' /* OpenBoundaries on/off flag */') |
155 |
|
CALL WRITE_1D_L( nonHydrostatic,1, INDEX_NONE, |
156 |
|
& 'nonHydrostatic =', ' /* Non-Hydrostatic on/off flag */') |
157 |
CALL WRITE_1D_R8( GMMaxSlope,1, INDEX_NONE,'GMMaxSlope =', |
CALL WRITE_1D_R8( GMMaxSlope,1, 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_1D_R8( GMLength,1, INDEX_NONE,'GMLength =', |
197 |
&' /* Base timestep number */') |
&' /* Base timestep number */') |
198 |
CALL WRITE_1D_I( nTimeSteps,1, INDEX_NONE,'nTimeSteps =', |
CALL WRITE_1D_I( nTimeSteps,1, INDEX_NONE,'nTimeSteps =', |
199 |
&' /* Number of timesteps */') |
&' /* Number of timesteps */') |
200 |
CALL WRITE_1D_R8( deltaTmom,1, INDEX_NONE,'deltatTmom =', |
CALL WRITE_1D_R8( deltaTmom,1, INDEX_NONE,'deltaTmom =', |
201 |
&' /* Momentum equation timestep ( s ) */') |
&' /* Momentum equation timestep ( s ) */') |
202 |
CALL WRITE_1D_R8( deltaTtracer,1, INDEX_NONE,'deltatTtracer =', |
CALL WRITE_1D_R8( deltaTtracer,1, INDEX_NONE,'deltaTtracer =', |
203 |
&' /* Tracer equation timestep ( s ) */') |
&' /* Tracer equation timestep ( s ) */') |
204 |
CALL WRITE_1D_R8( deltaTClock ,1, INDEX_NONE,'deltatTClock =', |
CALL WRITE_1D_R8( deltaTClock ,1, INDEX_NONE,'deltaTClock =', |
205 |
&' /* Model clock timestep ( s ) */') |
&' /* Model clock timestep ( s ) */') |
206 |
CALL WRITE_1D_R8( cAdjFreq,1, INDEX_NONE,'cAdjFreq =', |
CALL WRITE_1D_R8( cAdjFreq,1, INDEX_NONE,'cAdjFreq =', |
207 |
&' /* Convective adjustment interval ( s ) */') |
&' /* Convective adjustment interval ( s ) */') |