diff --git a/src/MNH/mode_mnh_zwork.f90 b/src/MNH/mode_mnh_zwork.f90 index 95fc395d786492dc2b7a6cac48345b07d05ef94c..7e2e6bdcc2c90bfd6c4c4ee0b88a7760ced24cc0 100644 --- a/src/MNH/mode_mnh_zwork.f90 +++ b/src/MNH/mode_mnh_zwork.f90 @@ -406,8 +406,10 @@ CONTAINS INTEGER :: KINDEX !local - - IF ( (KI .EQ. IIU) .AND. (KJ .EQ. IJU) ) THEN + IF ( KI == 0 .OR. KJ == 0 ) THEN + CALL MNH_GET_ZT3D_N0(KINDEX) + PTAB => ZT3D_OSIZE + ELSE IF ( (KI .EQ. IIU) .AND. (KJ .EQ. IJU) ) THEN CALL MNH_GET_ZT3D_N0(KINDEX) IF (KK .LE. IKU) THEN PTAB => ZT3D(:,:,1:KK,KINDEX) @@ -432,7 +434,10 @@ CONTAINS INTEGER :: KINDEX_END - IF ( (KI .EQ. IIU) .AND. (KJ .EQ. IJU) .AND. (KK .EQ. IKU) ) THEN + IF ( KI == 0 .OR. KJ == 0 .OR. KK == 0 ) THEN + PTAB => ZT4D_OSIZE + KINDEX_BEG = 0 + ELSE IF ( (KI .EQ. IIU) .AND. (KJ .EQ. IJU) .AND. (KK .EQ. IKU) ) THEN IF ( KL .GE. 1 ) THEN CALL MNH_GET_ZT4D(KL,KINDEX_BEG,KINDEX_END) PTAB => ZT3D(:,:,:,KINDEX_BEG:KINDEX_END)