Philippe 16/09/2021: allow more flexibility for allocation of zero-size 3D and 4D arrays

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)