diff --git a/src/MNH/mean_field.f90 b/src/MNH/mean_field.f90
index 28153fcdc72116cdf0e503407e88e1fa4fc7a5d3..2efc04c8089bbc5dd6fc818b80ba05701c4a2bda 100644
--- a/src/MNH/mean_field.f90
+++ b/src/MNH/mean_field.f90
@@ -81,6 +81,7 @@ USE MODD_EOL_ALM, ONLY: XFAERO_RE_SUM, XFAERO_RE_GLB
 USE MODI_SHUMAN,  ONLY: MXF, MYF, MZF
 !
 USE MODE_MODELN_HANDLER
+USE MODE_THERMO
 USE MODI_UPDATE_WELFORD
 ! 
 USE MODI_SHUMAN
@@ -126,9 +127,6 @@ REAL, DIMENSION(:,:),   ALLOCATABLE :: ZRH_I_MAXCOL
 REAL, DIMENSION(:,:),   ALLOCATABLE :: ZRH_P_MAXCOL
 INTEGER           :: JI,JJ,JK, JLOOP   ! Loop indexes
 REAL, DIMENSION(SIZE(PUT,1),SIZE(PUT,2),SIZE(PUT,3)) ::  ZTEMPT
-REAL, DIMENSION(SIZE(PUT,1),SIZE(PUT,2),SIZE(PUT,3)) :: ZUMEAN_OLD
-REAL, DIMENSION(SIZE(PUT,1),SIZE(PUT,2),SIZE(PUT,3)) :: ZWMEAN_OLD
-REAL, DIMENSION(SIZE(PUT,1),SIZE(PUT,2),SIZE(PUT,3)) :: ZTHMEAN_OLD
 REAL, DIMENSION(SIZE(PUT,1),SIZE(PUT,2),SIZE(PUT,3)) :: ZTHV ! 
 REAL, DIMENSION(SIZE(PUT,1),SIZE(PUT,2),SIZE(PUT,3)) ::  ZVOX, ZVOY, ZVOZ ! 
 REAL, DIMENSION(SIZE(PUT,1),SIZE(PUT,2),SIZE(PUT,3)) ::  ZWM !