124 |
& gravity, gBaro, rhonil, tAlpha, sBeta, |
& gravity, gBaro, rhonil, tAlpha, sBeta, |
125 |
& f0, beta, omega, rotationPeriod, |
& f0, beta, omega, rotationPeriod, |
126 |
& viscAh, viscAhW, viscAhMax, viscAhGrid, viscC2leith, |
& viscAh, viscAhW, viscAhMax, viscAhGrid, viscC2leith, |
127 |
|
& viscAhD, viscAhZ, viscA4D, viscA4Z, |
128 |
& viscA4, viscA4Max, viscA4Grid, viscC4leith, |
& viscA4, viscA4Max, viscA4Grid, viscC4leith, |
129 |
& viscAz, cosPower, viscAstrain, viscAtension, |
& viscAz, cosPower, viscAstrain, viscAtension, |
130 |
& diffKhT, diffKzT, diffK4T, |
& diffKhT, diffKzT, diffK4T, |
288 |
|
|
289 |
C-- Set default "physical" parameters |
C-- Set default "physical" parameters |
290 |
viscAhW = UNSET_RL |
viscAhW = UNSET_RL |
291 |
|
viscAhD = UNSET_RL |
292 |
|
viscAhZ = UNSET_RL |
293 |
|
viscA4D = UNSET_RL |
294 |
|
viscA4Z = UNSET_RL |
295 |
viscAz = UNSET_RL |
viscAz = UNSET_RL |
296 |
viscAr = UNSET_RL |
viscAr = UNSET_RL |
297 |
viscAp = UNSET_RL |
viscAp = UNSET_RL |
449 |
IF (saltVertAdvScheme.EQ.0) saltVertAdvScheme = saltAdvScheme |
IF (saltVertAdvScheme.EQ.0) saltVertAdvScheme = saltAdvScheme |
450 |
C-- horizontal viscosity for vertical moments |
C-- horizontal viscosity for vertical moments |
451 |
IF ( viscAhW .EQ. UNSET_RL ) viscAhW = viscAh |
IF ( viscAhW .EQ. UNSET_RL ) viscAhW = viscAh |
452 |
|
C-- horizontal viscosity (acting on Divergence or Vorticity) |
453 |
|
IF ( viscAhD .EQ. UNSET_RL ) viscAhD = viscAh |
454 |
|
IF ( viscAhZ .EQ. UNSET_RL ) viscAhZ = viscAh |
455 |
|
IF ( viscA4D .EQ. UNSET_RL ) viscA4D = viscA4 |
456 |
|
IF ( viscA4Z .EQ. UNSET_RL ) viscA4Z = viscA4 |
457 |
C-- z,p,r coord input switching. |
C-- z,p,r coord input switching. |
458 |
IF ( viscAz .NE. UNSET_RL ) zCoordInputData = .TRUE. |
IF ( viscAz .NE. UNSET_RL ) zCoordInputData = .TRUE. |
459 |
IF ( viscAp .NE. UNSET_RL ) pCoordInputData = .TRUE. |
IF ( viscAp .NE. UNSET_RL ) pCoordInputData = .TRUE. |