From 2760b0560765140e30935afada65e42515e53f90 Mon Sep 17 00:00:00 2001
From: Sibylle TECHENE <techenes@irene191.c-irene.tgcc.ccc.cea.fr>
Date: Sat, 19 Feb 2022 13:25:18 +0100
Subject: [PATCH] RK3: compute vertical speed for external mode
---
src/OCE/stp2d.F90 | 8 +++++---
1 file changed, 5 insertions(+), 3 deletions(-)
diff --git a/src/OCE/stp2d.F90 b/src/OCE/stp2d.F90
index 0cea844f1..6ca69903d 100644
--- a/src/OCE/stp2d.F90
+++ b/src/OCE/stp2d.F90
@@ -19,10 +19,8 @@ MODULE stp2d
!!----------------------------------------------------------------------
USE step_oce ! time stepping used modules
USE domqco ! quasi-eulerian coordinate (dom_qco_r3c routine)
- USE dynadv_cen2 ! centred flux form advection (dyn_adv_cen2 routine)
- USE dynadv_ubs ! UBS flux form advection (dyn_adv_ubs routine)
- USE dynkeg ! kinetic energy gradient (dyn_keg routine)
USE dynspg_ts ! 2D mode integration
+ USE sshwzv ! vertical speed
USE sbc_ice , ONLY : snwice_mass, snwice_mass_b
USE sbcapr ! surface boundary condition: atmospheric pressure
USE sbcwave, ONLY : bhd_wave
@@ -111,6 +109,10 @@ CONTAINS
!
! !* compute advection + coriolis *!
!
+ r3t(:,:,Kaa) = ssh(:,:,Kaa) * r1_ht_0(:,:) ! ratio at t-point at Kaa (n-1)
+ !
+ CALL wzv ( kt, Kbb, Kbb, Kaa , uu(:,:,:,Kbb), vv(:,:,:,Kbb), ww ) ! ww guess at Kbb (n)
+ !
CALL dyn_adv( kt, Kbb, Kbb , uu, vv, Krhs) !- vector form KEG+ZAD
! !- flux form ADV
CALL dyn_vor( kt, Kbb, uu, vv, Krhs ) !- vector form COR+RVO
--
GitLab