Philippe 05/04/2018: IO: store XHAT and YHAT when written by SURFEX (needed by PREP_PGD program)

src/MNH/write_surf_mnh.f90

@@ -257,6 +257,7 @@ USE MODE_TOOLS_ll
 USE MODE_WRITE_SURF_MNH_TOOLS
 !
 USE MODD_CONF_n,        ONLY: CSTORAGE_TYPE
+USE MODD_GRID_n,        ONLY: XXHAT, XYHAT
 USE MODD_IO_ll,         ONLY: TFILE_SURFEX
 USE MODD_IO_SURF_MNH,   ONLY :NMASK, CMASK,                          &
                               NIU, NJU, NIB, NJB, NIE, NJE,          &
@@ -408,6 +409,11 @@ IF (      (CSTORAGE_TYPE=='PG' .OR. CSTORAGE_TYPE=='SU')  &
   IF (HDIR=='H') THEN
     CALL FIND_FIELD_ID_FROM_MNHNAME('XHAT',IID,IRESP)
     CALL IO_WRITE_FIELD(TFILE_SURFEX,TFIELDLIST(IID),ZW1D(1+NHALO:IIU-NHALO),KRESP)
+    IF (.NOT. (ASSOCIATED(XXHAT))) THEN
+      !Store XXHAT if not yet done (necessary for PREP_PGD program when writing netCDF files)
+      ALLOCATE(XXHAT(IIU-2*NHALO))
+      XXHAT(:) = ZW1D(1+NHALO:IIU-NHALO)
+    END IF
   END IF
   DEALLOCATE(ZW1D)
 ELSE IF ( (CSTORAGE_TYPE=='PG' .OR. CSTORAGE_TYPE=='SU')  &
@@ -434,6 +440,11 @@ ELSE IF ( (CSTORAGE_TYPE=='PG' .OR. CSTORAGE_TYPE=='SU')  &
   IF (HDIR=='H') THEN
     CALL FIND_FIELD_ID_FROM_MNHNAME('YHAT',IID,IRESP)
     CALL IO_WRITE_FIELD(TFILE_SURFEX,TFIELDLIST(IID),ZW1D(1+NHALO:IJU-NHALO),KRESP)
+    IF (.NOT. (ASSOCIATED(XYHAT))) THEN
+      !Store XYHAT if not yet done (necessary for PREP_PGD program when writing netCDF files)
+      ALLOCATE(XYHAT(IJU-2*NHALO))
+      XYHAT(:) = ZW1D(1+NHALO:IJU-NHALO)
+    END IF
   END IF
   DEALLOCATE(ZW1D)
 ELSE IF (HDIR=='H') THEN