Juan 05/10/2015: missing JPHEXT for LAT/LON/ZS/ZSMT writing
authorJuan Escobar <juan.escobar@aero.obs-mip.fr>
Mon, 5 Oct 2015 15:42:08 +0000 (15:42 +0000)
committerPhilippe WAUTELET <philippe.wautelet@aero.obs-mip.fr>
Thu, 19 May 2016 14:44:49 +0000 (16:44 +0200)
src/MNH/prep_pgd.f90

index be54707..c8f76c7 100644 (file)
@@ -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 /) )