New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
#2745 (zdfric not compatible with p_sh2 calculation) – NEMO

Opened 6 months ago

Last modified 5 months ago

#2745 assigned Bug

zdfric not compatible with p_sh2 calculation

Reported by: amoulin Owned by: systeam
Priority: low Milestone:
Component: ZDF Version: trunk
Severity: minor Keywords:
Cc: emanuela.clementi@…

Description

Context

...

Analysis

Richardson number vertical mixing is not compatible with the shear production term (p_sh2) calculated in zdfsh2 (using the before and after velocities)

(Simulation in the MEDSEA gave really high RMSE in the temperature and salinity fields)

Fix

We propose to reintroduce the calculation of the shear production term as it was in the past version (using only before velocities)

zdku = 0.5 / e3uw(ji,jj,jk,Kbb) * ( uu(ji-1,jj,jk-1,Kbb) + uu(ji,jj,jk-1,Kbb) &

& - uu(ji-1,jj,jk,Kbb) - uu(ji,jj,jk,Kbb) )* wumask(ji,jj,jk)

zdkv = 0.5 / e3vw(ji,jj,jk,Kbb) * ( vv(ji,jj-1,jk-1,Kbb) + vv(ji,jj,jk-1,Kbb) &

& - vv(ji,jj-1,jk,Kbb) - vv(ji,jj,jk,Kbb) )* wvmask(ji,jj,jk)

zwx = zdku*zdku + zdkv*zdkv


zcfRi = 1._wp / ( 1._wp + rn_alp * MAX( 0._wp , rn2(ji,jj,jk) / ( zwx + 1.e-20 ) ) )

Commit History (0)

(No commits)

Change History (1)

comment:1 Changed 5 months ago by nemo

  • Owner set to systeam
  • Status changed from new to assigned
Note: See TracTickets for help on using tickets.