diff --git a/src/array_mod.F90 b/src/array_mod.F90
index 4f4096df0c4114c554a8ed9da3abec198ce6ce0e..213c108be9650f8109646c20165a7681ab7fdcf0 100644
--- a/src/array_mod.F90
+++ b/src/array_mod.F90
@@ -14,10 +14,10 @@ MODULE array
   ! Derivatives and interpolated moments
   COMPLEX(dp), DIMENSION(:,:,:,:,:), ALLOCATABLE :: ddz_nepj
   COMPLEX(dp), DIMENSION(:,:,:,:,:), ALLOCATABLE :: interp_nepj
-  COMPLEX(dp), DIMENSION(:,:,:,:,:), ALLOCATABLE :: ddzND_nepj
+  COMPLEX(dp), DIMENSION(:,:,:,:,:), ALLOCATABLE :: ddzND_Nepj
   COMPLEX(dp), DIMENSION(:,:,:,:,:), ALLOCATABLE :: ddz_nipj
   COMPLEX(dp), DIMENSION(:,:,:,:,:), ALLOCATABLE :: interp_nipj
-  COMPLEX(dp), DIMENSION(:,:,:,:,:), ALLOCATABLE :: ddzND_nipj
+  COMPLEX(dp), DIMENSION(:,:,:,:,:), ALLOCATABLE :: ddzND_Nipj
 
   ! Arrays to store special initial modes (semi linear simulation)
   ! Zonal ones (ky=0)
@@ -46,6 +46,9 @@ MODULE array
   ! dnjs coefficient storage (in, ij, is)
   COMPLEX(dp), DIMENSION(:,:,:), ALLOCATABLE :: dnjs
 
+  ! Hermite fourth derivative coeff storage 4*sqrt(p!/(p-4)!)
+  COMPLEX(dp), DIMENSION(:), ALLOCATABLE :: dv4_Hp_coeff
+
   ! lin rhs p,j coefficient storage (ip,ij)
   REAL(dp), DIMENSION(:,:), ALLOCATABLE :: xnepj,xnipj
   REAL(dp), DIMENSION(:),   ALLOCATABLE :: xnepp1j, xnepm1j,   xnepp2j,   xnepm2j, xnepjp1, xnepjm1