66 |
&' /* Lateral eddy viscosity ( m^2/s ) */') |
&' /* Lateral eddy viscosity ( m^2/s ) */') |
67 |
CALL WRITE_1D_R8( viscA4, 1, INDEX_NONE,'viscAh =', |
CALL WRITE_1D_R8( viscA4, 1, INDEX_NONE,'viscAh =', |
68 |
&' /* Lateral biharmonic viscosity ( m^4/s ) */') |
&' /* Lateral biharmonic viscosity ( m^4/s ) */') |
69 |
CALL WRITE_1D_R8( viscAz, 1, INDEX_NONE,'viscAz =', |
IF ( viscAz .NE. UNSET_RL ) THEN |
70 |
&' /* Vertical eddy viscosity ( m^2/s ) */') |
CALL WRITE_1D_R8( viscAz, 1, INDEX_NONE,'viscAz =', |
71 |
|
& ' /* Vertical eddy viscosity ( m^2/s ) */') |
72 |
|
ENDIF |
73 |
|
IF ( viscAp .NE. UNSET_RL ) THEN |
74 |
|
CALL WRITE_1D_R8( viscAp, 1, INDEX_NONE,'viscAp =', |
75 |
|
& ' /* Vertical eddy viscosity ( Pa^2/s ) */') |
76 |
|
ENDIF |
77 |
|
CALL WRITE_1D_R8( viscAr, 1, INDEX_NONE,'viscAr =', |
78 |
|
&' /* Vertical eddy viscosity ( units of r^2/s ) */') |
79 |
CALL WRITE_1D_R8( diffKhT, 1, INDEX_NONE,'diffKhT =', |
CALL WRITE_1D_R8( diffKhT, 1, INDEX_NONE,'diffKhT =', |
80 |
&' /* Laplacian diffusion of heat laterally ( m^2/s ) */') |
&' /* Laplacian diffusion of heat laterally ( m^2/s ) */') |
81 |
CALL WRITE_1D_R8( diffKzT, 1, INDEX_NONE,'diffKzT =', |
CALL WRITE_1D_R8( diffKzT, 1, INDEX_NONE,'diffKzT =', |
92 |
&' /* Linear EOS thermal expansion coefficient ( 1/degree ) */') |
&' /* Linear EOS thermal expansion coefficient ( 1/degree ) */') |
93 |
CALL WRITE_1D_R8( sBeta, 1, INDEX_NONE,'sBeta =', |
CALL WRITE_1D_R8( sBeta, 1, INDEX_NONE,'sBeta =', |
94 |
&' /* Linear EOS haline contraction coefficient ( 1/ppt ) */') |
&' /* Linear EOS haline contraction coefficient ( 1/ppt ) */') |
95 |
|
IF ( eosType .EQ. 'POLY3' ) THEN |
96 |
|
WRITE(msgBuf,'(A)') '// Polynomial EQS parameters ( from POLY3.COEFFS ) ' |
97 |
|
DO K = 1, Nr |
98 |
|
WRITE(msgBuf,'(I3,13F8.3)') |
99 |
|
& K,eosRefT(K),eosRefS(K),eosSig0(K), (eosC(I,K),I=1,9) |
100 |
|
CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1) |
101 |
|
ENDDO |
102 |
|
ENDIF |
103 |
CALL WRITE_1D_R8( rhonil,1, INDEX_NONE,'rhonil =', |
CALL WRITE_1D_R8( rhonil,1, INDEX_NONE,'rhonil =', |
104 |
&' /* Reference density ( kg/m^3 ) */') |
&' /* Reference density ( kg/m^3 ) */') |
105 |
CALL WRITE_1D_R8( gravity,1, INDEX_NONE,'gravity =', |
CALL WRITE_1D_R8( gravity,1, INDEX_NONE,'gravity =', |