110 |
&'useFullLeith =', |
&'useFullLeith =', |
111 |
&' /* Use Full Form of Leith Viscosity on/off flag*/') |
&' /* Use Full Form of Leith Viscosity on/off flag*/') |
112 |
CALL WRITE_0D_R8( viscC2leith, INDEX_NONE,'viscC2leith =', |
CALL WRITE_0D_R8( viscC2leith, INDEX_NONE,'viscC2leith =', |
113 |
&' /* Leith harmonic viscosity factor ( on grad(vort), non-dom. ) */') |
&' /* Leith harmonic viscosity factor ( on grad(vort), non-dim. ) */') |
114 |
CALL WRITE_0D_R8( viscC2leithD, INDEX_NONE,'viscC2leithD =', |
CALL WRITE_0D_R8( viscC2leithD, INDEX_NONE,'viscC2leithD =', |
115 |
&' /* Leith harmonic viscosity factor ( on grad(div), non-dom. ) */') |
&' /* Leith harmonic viscosity factor ( on grad(div), non-dim. ) */') |
116 |
|
CALL WRITE_0D_R8( viscC2smag, INDEX_NONE,'viscC2smag =', |
117 |
|
&' /* Smagorinsky harmonic viscosity factor (non-dim.) */') |
118 |
CALL WRITE_0D_R8( viscA4, INDEX_NONE,'viscA4 =', |
CALL WRITE_0D_R8( viscA4, INDEX_NONE,'viscA4 =', |
119 |
&' /* Lateral biharmonic viscosity ( m^4/s ) */') |
&' /* Lateral biharmonic viscosity ( m^4/s ) */') |
120 |
IF ( viscA4D.NE.viscA4 ) |
IF ( viscA4D.NE.viscA4 ) |