Opened 2 years ago
Last modified 2 years 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 2 years ago by nemo
- Owner set to systeam
- Status changed from new to assigned
