From 69ff4fa3ce2df6e6f12b20da6734266ff4deee82 Mon Sep 17 00:00:00 2001
From: Christian Ethe <christian.ethe@ipsl.fr>
Date: Mon, 27 Jun 2022 14:58:32 +0000
Subject: [PATCH] Resolve "Missing boundaries conditions in NEMO-TOP 1D"
---
src/TOP/TRP/trctrp.F90 | 3 +++
1 file changed, 3 insertions(+)
diff --git a/src/TOP/TRP/trctrp.F90 b/src/TOP/TRP/trctrp.F90
index f164444ed..135466dc8 100644
--- a/src/TOP/TRP/trctrp.F90
+++ b/src/TOP/TRP/trctrp.F90
@@ -116,6 +116,9 @@ CONTAINS
!
ELSE ! 1D vertical configuration
CALL trc_sbc( kt, Kmm, tr, Krhs ) ! surface boundary condition
+ IF( ln_trcbc .AND. lltrcbc .AND. kt /= nit000 ) &
+ CALL trc_bc ( kt, Kmm, tr, Krhs ) ! tracers: surface and lateral Boundary Conditions
+ IF( ln_trcais ) CALL trc_ais ( kt, Kmm, tr, Krhs ) ! tracers from Antarctic Ice Sheet (icb, isf)
IF( ln_trcdmp ) CALL trc_dmp( kt, Kbb, Kmm, tr, Krhs ) ! internal damping trends
CALL trc_zdf( kt, Kbb, Kmm, Krhs, tr, Kaa ) ! vert. mixing & after tracer ==> after
CALL trc_atf( kt, Kbb, Kmm, Kaa , tr ) ! time filtering of "now" tracer fields
--
GitLab