From acc4a1fb8e23a9c3962b94c00edd3781fa29091b Mon Sep 17 00:00:00 2001
From: Guillaume Samson <guillaume.samson@mercator-ocean.fr>
Date: Thu, 12 May 2022 16:15:00 +0000
Subject: [PATCH] Resolve "ah_wslp2 variable not initialised when using
ln_trcldf_OFF = .true. in TOP"
---
src/OCE/LDF/ldfslp.F90 | 5 +++++
src/OCE/TRA/traldf_iso.F90 | 5 -----
src/OCE/TRA/traldf_triad.F90 | 5 -----
3 files changed, 5 insertions(+), 10 deletions(-)
diff --git a/src/OCE/LDF/ldfslp.F90 b/src/OCE/LDF/ldfslp.F90
index 4a85f6991..d48046735 100644
--- a/src/OCE/LDF/ldfslp.F90
+++ b/src/OCE/LDF/ldfslp.F90
@@ -684,6 +684,11 @@ CONTAINS
ALLOCATE( ah_wslp2(jpi,jpj,jpk) , akz(jpi,jpj,jpk) , STAT=ierr )
IF( ierr > 0 ) CALL ctl_stop( 'STOP', 'ldf_slp_init : unable to allocate ah_slp2 or akz' )
!
+ DO_3D( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1, 1, jpk )
+ akz (ji,jj,jk) = 0._wp
+ ah_wslp2(ji,jj,jk) = 0._wp
+ END_3D
+ !
IF( ln_traldf_triad ) THEN ! Griffies operator : triad of slopes
IF(lwp) WRITE(numout,*) ' ==>>> triad) operator (Griffies)'
ALLOCATE( triadi_g(jpi,jpj,jpk,0:1,0:1) , triadj_g(jpi,jpj,jpk,0:1,0:1) , &
diff --git a/src/OCE/TRA/traldf_iso.F90 b/src/OCE/TRA/traldf_iso.F90
index f93bf2dcd..8d84e7689 100644
--- a/src/OCE/TRA/traldf_iso.F90
+++ b/src/OCE/TRA/traldf_iso.F90
@@ -145,11 +145,6 @@ CONTAINS
IF(lwp) WRITE(numout,*) 'tra_ldf_iso : rotated laplacian diffusion operator on ', cdtype
IF(lwp) WRITE(numout,*) '~~~~~~~~~~~'
ENDIF
- !
- DO_3D_OVR( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1, 1, jpk )
- akz (ji,jj,jk) = 0._wp
- ah_wslp2(ji,jj,jk) = 0._wp
- END_3D
ENDIF
!
IF( .NOT. l_istiled .OR. ntile == 1 ) THEN ! Do only on the first tile
diff --git a/src/OCE/TRA/traldf_triad.F90 b/src/OCE/TRA/traldf_triad.F90
index 70eddec54..19039b88c 100644
--- a/src/OCE/TRA/traldf_triad.F90
+++ b/src/OCE/TRA/traldf_triad.F90
@@ -147,11 +147,6 @@ CONTAINS
!!----------------------------------------------------------------------
!
IF( kpass == 1 ) THEN !== first pass only and whatever the tracer is ==!
- !
- DO_3D_OVR( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1, 1, jpk )
- akz (ji,jj,jk) = 0._wp
- ah_wslp2(ji,jj,jk) = 0._wp
- END_3D
!
DO kp = 0, 1 ! i-k triads
DO_3D_OVR( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1, 1, jpkm1 )
--
GitLab