From c2b6c44e8a0259a3b1d01edeee96f46b81f2ce2c Mon Sep 17 00:00:00 2001
From: sebastien masson <smasson@amphitrite.locean-ipsl.upmc.fr>
Date: Thu, 22 Sep 2022 11:19:10 +0200
Subject: [PATCH] minor bugfix in icesbc.F90
---
src/ICE/icesbc.F90 | 2 +-
tests/BENCH/MY_SRC/usrdef_sbc.F90 | 7 +++++++
2 files changed, 8 insertions(+), 1 deletion(-)
diff --git a/src/ICE/icesbc.F90 b/src/ICE/icesbc.F90
index 03d5b83a6..d9406ce66 100644
--- a/src/ICE/icesbc.F90
+++ b/src/ICE/icesbc.F90
@@ -100,7 +100,7 @@ CONTAINS
!
ENDIF
!
- CALL lbc_lnk( 'icesbc', utau_ice, 'U', -1.0_wp, vtau_ice, 'V', -1.0_wp )
+ CALL lbc_lnk( 'icesbc', utau_ice, 'T', -1.0_wp, vtau_ice, 'T', -1.0_wp )
!
IF( ln_timing ) CALL timing_stop('icesbc')
!
diff --git a/tests/BENCH/MY_SRC/usrdef_sbc.F90 b/tests/BENCH/MY_SRC/usrdef_sbc.F90
index 7ad3b7064..9d7d2808d 100644
--- a/tests/BENCH/MY_SRC/usrdef_sbc.F90
+++ b/tests/BENCH/MY_SRC/usrdef_sbc.F90
@@ -109,6 +109,13 @@ CONTAINS
vtau_ice(ji,jj) = 0.1_wp + zztmp
END_2D
+ IF( l_NFold .AND. c_NFtype == 'T' ) THEN ! force 0 at the folding points
+ utau_ice(mi0(jpiglo/2+1,nn_hls):mi1(jpiglo/2+1,nn_hls),mj0(jpjglo-nn_hls,nn_hls):mj1(jpjglo-nn_hls,nn_hls)) = 0._wp
+ vtau_ice(mi0(jpiglo/2+1,nn_hls):mi1(jpiglo/2+1,nn_hls),mj0(jpjglo-nn_hls,nn_hls):mj1(jpjglo-nn_hls,nn_hls)) = 0._wp
+ utau_ice(mi0( nn_hls+1,nn_hls):mi1( nn_hls+1,nn_hls),mj0(jpjglo-nn_hls,nn_hls):mj1(jpjglo-nn_hls,nn_hls)) = 0._wp
+ vtau_ice(mi0( nn_hls+1,nn_hls):mi1( nn_hls+1,nn_hls),mj0(jpjglo-nn_hls,nn_hls):mj1(jpjglo-nn_hls,nn_hls)) = 0._wp
+ ENDIF
+
CALL lbc_lnk( 'usrdef_sbc', utau_ice, 'T', -1., vtau_ice, 'T', -1. )
#endif
!
--
GitLab