From fd42a0e7860fd368627e6d59b5de6ec5a08c60f0 Mon Sep 17 00:00:00 2001
From: Sebastien Masson <sebastien.masson@locean.ipsl.fr>
Date: Mon, 9 Jan 2023 12:50:33 +0100
Subject: [PATCH] minor fix, #150
---
src/OCE/SBC/sbcssr.F90 | 11 +++++------
1 file changed, 5 insertions(+), 6 deletions(-)
diff --git a/src/OCE/SBC/sbcssr.F90 b/src/OCE/SBC/sbcssr.F90
index 2b607dbb..1cfc838b 100644
--- a/src/OCE/SBC/sbcssr.F90
+++ b/src/OCE/SBC/sbcssr.F90
@@ -93,15 +93,14 @@ CONTAINS
IF( MOD( kt-1, nn_fsbc ) == 0 ) THEN ! Add restoring term !
! ! ========================= !
!
- qrp(:,:) = 0._wp ! necessary init
- erp(:,:) = 0._wp
- !
- IF( nn_sstr == 1 ) THEN !* Temperature restoring term
+ IF( nn_sstr == 1 ) THEN !* Temperature restoring term
DO_2D( 0, 0, 0, 0 )
zqrp = rn_dqdt * ( sst_m(ji,jj) - sf_sst(1)%fnow(ji,jj,1) ) * smask0(ji,jj)
qns(ji,jj) = qns(ji,jj) + zqrp
qrp(ji,jj) = zqrp
END_2D
+ ELSEIF( nn_sssr == 2 ) THEN
+ qrp(:,:) = 0._wp ! necessary init, see bellow: qrp(ji,jj) = qrp(ji,jj) - ...
ENDIF
!
IF( nn_sssr /= 0 .AND. nn_sssr_ice /= 1 ) THEN
@@ -233,8 +232,8 @@ CONTAINS
!
coefice(:,:) = 1._wp ! Initialise coefice to 1._wp ; will not need to be changed if nn_sssr_ice=1
! !* Initialize qrp and erp if no restoring
- IF( nn_sstr /= 1 ) qrp(:,:) = 0._wp
- IF( nn_sssr /= 1 .OR. nn_sssr /= 2 ) erp(:,:) = 0._wp
+ IF( nn_sstr /= 1 .AND. nn_sssr /= 2 ) qrp(:,:) = 0._wp
+ IF( nn_sssr /= 1 .AND. nn_sssr /= 2 ) erp(:,:) = 0._wp
!
END SUBROUTINE sbc_ssr_init
--
GitLab