Tidal mixing energy in zdt_tmx with VVL

[timgraham]

In zdf_tmx the vertical structure of the tidal mixing energy (az_tmx) is calculated in zdf_tmx_init at the start of each model run and is dependent on the model depths. Ticket #1423 changed this to use gdepw_0 instead of fsdepw so that the model is restartable when using VVL.

In order to work correctly should az_tmx be calculated at each time step in the main zdf_tmx routine using fsdepw? I think this would involve moving the 2 loops below into the main routine.

      DO jj = 1, jpj                ! part independent of the level
         DO ji = 1, jpi
            zhdep(ji,jj) = gdepw_0(ji,jj,mbkt(ji,jj)+1)       ! depth of the ocean
            zfact(ji,jj) = rau0 * rn_htmx * ( 1. - EXP( -zhdep(ji,jj) / rn_htmx ) )
            IF( zfact(ji,jj) /= 0 )   zfact(ji,jj) = en_tmx(ji,jj) / zfact(ji,jj)
         END DO
      END DO
      DO jk= 1, jpk                 ! complete with the level-dependent part
         DO jj = 1, jpj
            DO ji = 1, jpi
               az_tmx(ji,jj,jk) = zfact(ji,jj) * EXP( -( zhdep(ji,jj)-gdepw_0(ji,jj,jk) ) / rn_htmx ) * tmask(ji,jj,jk)
            END DO
         END DO
      END DO

[clevy]

Suppressing tickets opened on the trunk since more than 2 years... cleaning started!