diff --git a/src/MNH/prep_pgd.f90 b/src/MNH/prep_pgd.f90
index be54707adea52b2ed84d1111017c3db405a2af41..c8f76c7c691cb870a5798b537799021ae62ef2d8 100644
--- a/src/MNH/prep_pgd.f90
+++ b/src/MNH/prep_pgd.f90
@@ -61,6 +61,7 @@
 !!    S.Bielli     23/04/2014     supress writing of LAt and LON in NETCDF case
 !!    S.Bielli     20/11/2014     add writing of LAt and LON in NETCDF case
 !!    J.Escobar : 15/09/2015 : WENO5 & JPHEXT <> 1 
+!!    J.Escobar : 05/10/2015 : missing JPHEXT for LAT/LON/ZS/ZSMT writing
 !----------------------------------------------------------------------------
 !
 !*    0.     DECLARATION
@@ -256,8 +257,8 @@ IF (LNETCDF.AND..NOT.LCARTESIAN) THEN
 !!!! WRITE LAT and LON
   CALL GET_DIM_PHYS_ll('B',IIMAX,IJMAX)
   ALLOCATE(ZWORK(IIMAX+NHALO*2,IJMAX+NHALO*2))
-  ALLOCATE(ZWORK_LAT(IIMAX+2,IJMAX+2))
-  ALLOCATE(ZWORK_LON(IIMAX+2,IJMAX+2))
+  ALLOCATE(ZWORK_LAT(IIMAX+2*JPHEXT,IJMAX+2*JPHEXT))
+  ALLOCATE(ZWORK_LON(IIMAX+2*JPHEXT,IJMAX+2*JPHEXT))
   ZWORK=RESHAPE(XLAT, (/ (IIMAX+NHALO*2),(IJMAX+NHALO*2) /) )
   ZWORK_LAT=ZWORK(NHALO:(IIMAX+NHALO+1),NHALO:(IJMAX+NHALO+1))
 !!
@@ -281,7 +282,7 @@ IF ( LNETCDF ) THEN
     CALL FMWRIT(COUTFMFILE,'LAT',CLUOUT0,'XY',ZWORK_LAT,1,21,'X_Y_latitude (degree)',IRESP)
     CALL FMWRIT(COUTFMFILE,'LON',CLUOUT0,'XY',ZWORK_LON,1,22,'X_Y_longitude (degree)',IRESP)
   END IF
-  ALLOCATE(ZZS(IIMAX+2,IJMAX+2))
+  ALLOCATE(ZZS(IIMAX+2*JPHEXT,IJMAX+2*JPHEXT))
 !!!!  writes smooth orography for SLEVE coordinate in netcdf
   YRECFM = 'ZS              '
   CALL FMREAD(COUTFMFILE,YRECFM,CLUOUT0,'XY',ZZS,IGRID,ILENCH,YCOMMENT,IRESP)
@@ -301,8 +302,8 @@ IF (.NOT.LCARTESIAN) THEN
 !!!! WRITE LAT and LON
    CALL GET_DIM_PHYS_ll('B',IIMAX,IJMAX)
    ALLOCATE(ZWORK(IIMAX+NHALO*2,IJMAX+NHALO*2))
-   ALLOCATE(ZWORK_LAT(IIMAX+2,IJMAX+2))
-   ALLOCATE(ZWORK_LON(IIMAX+2,IJMAX+2))
+   ALLOCATE(ZWORK_LAT(IIMAX+2*JPHEXT,IJMAX+2*JPHEXT))
+   ALLOCATE(ZWORK_LON(IIMAX+2*JPHEXT,IJMAX+2*JPHEXT))
    ZWORK=RESHAPE(XLAT, (/ (IIMAX+NHALO*2),(IJMAX+NHALO*2) /) )
    ZWORK_LAT=ZWORK(NHALO:(IIMAX+NHALO+1),NHALO:(IJMAX+NHALO+1))
    ZWORK=RESHAPE(XLON, (/ IIMAX+NHALO*2,IJMAX+NHALO*2 /) )