28 |
#include "SIZE.h" |
#include "SIZE.h" |
29 |
#include "EEPARAMS.h" |
#include "EEPARAMS.h" |
30 |
#include "PARAMS.h" |
#include "PARAMS.h" |
31 |
|
#include "EOS.h" |
32 |
#include "GRID.h" |
#include "GRID.h" |
33 |
#include "DYNVARS.h" |
#include "DYNVARS.h" |
34 |
|
|
136 |
&' /* Reference density ( kg/m^3 ) */') |
&' /* Reference density ( kg/m^3 ) */') |
137 |
CALL WRITE_0D_R8( rhoConst, INDEX_NONE,'rhoConst =', |
CALL WRITE_0D_R8( rhoConst, INDEX_NONE,'rhoConst =', |
138 |
&' /* Reference density ( kg/m^3 ) */') |
&' /* Reference density ( kg/m^3 ) */') |
139 |
|
CALL WRITE_0D_R8( rhoConstFresh, INDEX_NONE,'rhoConstFresh =', |
140 |
|
&' /* Reference density ( kg/m^3 ) */') |
141 |
CALL WRITE_0D_R8( gravity, INDEX_NONE,'gravity =', |
CALL WRITE_0D_R8( gravity, INDEX_NONE,'gravity =', |
142 |
&' /* Gravitational acceleration ( m/s^2 ) */') |
&' /* Gravitational acceleration ( m/s^2 ) */') |
143 |
CALL WRITE_0D_R8( gBaro, INDEX_NONE,'gBaro =', |
CALL WRITE_0D_R8( gBaro, INDEX_NONE,'gBaro =', |
199 |
&' /* convert F.W. Flux to Salt Flux (-1=use local S)(ppt)*/') |
&' /* convert F.W. Flux to Salt Flux (-1=use local S)(ppt)*/') |
200 |
ENDIF |
ENDIF |
201 |
|
|
202 |
|
CALL WRITE_0D_L( multiDimAdvection, INDEX_NONE, |
203 |
|
& 'multiDimAdvection =', |
204 |
|
&' /* enable/disable Multi-Dim Advection */') |
205 |
CALL WRITE_0D_L( staggerTimeStep, INDEX_NONE, |
CALL WRITE_0D_L( staggerTimeStep, INDEX_NONE, |
206 |
& 'staggerTimeStep =', |
& 'staggerTimeStep =', |
207 |
&' /* Stagger time stepping on/off flag */') |
&' /* Stagger time stepping on/off flag */') |
220 |
& ' /* Momentum pressure term on/off flag */') |
& ' /* Momentum pressure term on/off flag */') |
221 |
CALL WRITE_0D_L( tempStepping, INDEX_NONE, |
CALL WRITE_0D_L( tempStepping, INDEX_NONE, |
222 |
& 'tempStepping =', ' /* Temperature equation on/off flag */') |
& 'tempStepping =', ' /* Temperature equation on/off flag */') |
223 |
|
CALL WRITE_0D_L( tempAdvection, INDEX_NONE, |
224 |
|
& 'tempAdvection=', ' /* Temperature advection on/off flag */') |
225 |
|
CALL WRITE_0D_L( tempForcing, INDEX_NONE, |
226 |
|
& 'tempForcing =', ' /* Temperature forcing on/off flag */') |
227 |
|
CALL WRITE_0D_L( saltStepping, INDEX_NONE, |
228 |
|
& 'saltStepping =', ' /* Salinity equation on/off flag */') |
229 |
|
CALL WRITE_0D_L( saltAdvection, INDEX_NONE, |
230 |
|
& 'saltAdvection=', ' /* Salinity advection on/off flag */') |
231 |
|
CALL WRITE_0D_L( saltForcing, INDEX_NONE, |
232 |
|
& 'saltForcing =', ' /* Salinity forcing on/off flag */') |
233 |
CALL WRITE_0D_L( nonHydrostatic, INDEX_NONE, |
CALL WRITE_0D_L( nonHydrostatic, INDEX_NONE, |
234 |
& 'nonHydrostatic =', ' /* Non-Hydrostatic on/off flag */') |
& 'nonHydrostatic =', ' /* Non-Hydrostatic on/off flag */') |
235 |
WRITE(msgBuf,'(A)') '// ' |
WRITE(msgBuf,'(A)') '// ' |
273 |
&' /* Model clock timestep ( s ) */') |
&' /* Model clock timestep ( s ) */') |
274 |
CALL WRITE_0D_R8( cAdjFreq, INDEX_NONE,'cAdjFreq =', |
CALL WRITE_0D_R8( cAdjFreq, INDEX_NONE,'cAdjFreq =', |
275 |
&' /* Convective adjustment interval ( s ) */') |
&' /* Convective adjustment interval ( s ) */') |
276 |
|
CALL WRITE_0D_L( forcing_In_AB,INDEX_NONE,'forcing_In_AB =', |
277 |
|
&' /* put T,S Forcing in Adams-Bash. stepping */') |
278 |
CALL WRITE_0D_R8( abeps, INDEX_NONE,'abeps =', |
CALL WRITE_0D_R8( abeps, INDEX_NONE,'abeps =', |
279 |
&' /* Adams-Bashforth stabilizing weight */') |
&' /* Adams-Bashforth stabilizing weight */') |
280 |
CALL WRITE_0D_R8( tauCD, INDEX_NONE,'tauCD =', |
CALL WRITE_0D_R8( tauCD, INDEX_NONE,'tauCD =', |
310 |
&' /* Spherical coordinates flag ( True / False ) */') |
&' /* Spherical coordinates flag ( True / False ) */') |
311 |
CALL WRITE_0D_L( groundAtK1, INDEX_NONE, 'groundAtK1 =', |
CALL WRITE_0D_L( groundAtK1, INDEX_NONE, 'groundAtK1 =', |
312 |
&' /* Lower Boundary (ground) at the surface(k=1) ( T / F ) */') |
&' /* Lower Boundary (ground) at the surface(k=1) ( T / F ) */') |
313 |
CALL WRITE_1D_R8( Ro_SeaLevel,1, INDEX_NONE,'Ro_SeaLevel =', |
CALL WRITE_0D_R8( Ro_SeaLevel, INDEX_NONE,'Ro_SeaLevel =', |
314 |
&' /* r(1) ( units of r ) */') |
&' /* r(1) ( units of r ) */') |
315 |
CALL WRITE_1D_R8( rkFac,1, INDEX_NONE,'rkFac =', |
CALL WRITE_0D_R8( rkFac, INDEX_NONE,'rkFac =', |
316 |
&' /* minus Vertical index orientation */') |
&' /* minus Vertical index orientation */') |
317 |
CALL WRITE_1D_R8( horiVertRatio,1, INDEX_NONE,'horiVertRatio =', |
CALL WRITE_0D_R8( horiVertRatio, INDEX_NONE,'horiVertRatio =', |
318 |
&' /* Ratio on units : Horiz - Vertical */') |
&' /* Ratio on units : Horiz - Vertical */') |
319 |
c CALL WRITE_1D_R8( delZ,Nr, INDEX_K,'delZ = ', |
c CALL WRITE_1D_R8( delZ,Nr, INDEX_K,'delZ = ', |
320 |
c &' /* W spacing ( m ) */') |
c &' /* W spacing ( m ) */') |