Juan 01/02/2023:ZSOLVER/tridz.f90, BUG correction, Bypass SegFault for...

Juan 01/02/2023:ZSOLVER/tridz.f90, BUG correction, Bypass SegFault for unnecessary computation for BSPLITTING ( of YSPLITTING)

src/ZSOLVER/tridz.f90

@@ -672,6 +672,8 @@ IF (L2D) THEN
                                ZDZM(IKE+1) **2
   !
 ELSE
+  IF ( CSPLIT /= 'BSPLITTING' ) THEN
+  print*,"WARNING CSPLIT /= 'BSPLITTING =", CSPLIT     
   DO JK= IKB,IKE
     DO JJ= 1, IYMODEY_ll
       DO JI= 1, IXMODEY_ll
@@ -689,6 +691,7 @@ ELSE
   !
   PBFY(1:IYMODEY_ll,1:IXMODEY_ll,IKE+1) =  0.5 * ( PRHOM(IKE)   + PRHOM(IKE+1) ) /  &
                                ZDZM(IKE+1) **2
+  END IF
   !
 
 !JUAN Z_SPLITTING