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