Philippe 28/11/2022: FFT: fix for TRIDZ

src/MNH/tridz.f90

@@ -179,6 +179,7 @@ USE MODD_CONF
 USE MODD_LUNIT_n, ONLY: TLUOUT
 USE MODD_PARAMETERS
 !
+USE MODE_FFT,            ONLY: SET99
 USE MODE_ll
 USE MODE_MSG
 !JUAN P1/P2 SPLITTING
@@ -630,11 +631,9 @@ ELSE
 ! at the upper and lower levels PBFB is computed using the Neumann
 ! condition
 !
-  PBFB(IIB:IIE,IJB:IJE,IKB-1) = -0.5 * ( PRHOM(IKB-1) + PRHOM(IKB)   ) /  &
-                               ZDZM(IKB)   **2
+  PBFB(IIB:IIE,IJB:IJE,IKB-1) = - 0.5 * ( PRHOM(IKB-1) + PRHOM(IKB)   ) / ZDZM(IKB)   **2
   !
-  PBFB(IIB:IIE,IJB:IJE,IKE+1) =  0.5 * ( PRHOM(IKE)   + PRHOM(IKE+1) ) /  &
-                               ZDZM(IKE+1) **2
+  PBFB(IIB:IIE,IJB:IJE,IKE+1) = + 0.5 * ( PRHOM(IKE)   + PRHOM(IKE+1) ) / ZDZM(IKE+1) **2
 !
 IF (HLBCX(1) == 'CYCL' .AND. .NOT.(L2D) ) THEN
    !JUAN
@@ -695,7 +694,10 @@ DEALLOCATE(ZEIGEN_ll)
 !              -------------------------------------------------------
 !
 !        8.1 x lateral boundary conditions
-!                  
+!
+!Initialise PTRIGSX and KIFAXX to allow comparisons (they are not fully set in SET99)
+PTRIGSX(:) = XUNDEF
+KIFAXX = NUNDEF
 CALL SET99(PTRIGSX,KIFAXX,IIMAX_ll)
 !
 ! test on the value of KIMAX: KIMAX must be factorizable as a product
@@ -728,6 +730,9 @@ ENDIF
 !       8.2 y lateral boundary conditions
 !
 IF (.NOT. L2D) THEN 
+  !Initialise PTRIGSY and KIFAXY to allow comparisons (they are not fully set in SET99)
+  PTRIGSY(:) = XUNDEF
+  KIFAXY = NUNDEF
   CALL SET99(PTRIGSY,KIFAXY,IJMAX_ll)
   !
   ! test on the value of KJMAX: KJMAX must be factorizable as a product