Skip to content
Snippets Groups Projects
Commit e1d42556 authored by WAUTELET Philippe's avatar WAUTELET Philippe
Browse files

Philippe 06/02/2023: WRITE_SURFX2_MNH: treat 2D case separately

parent f41638c8
No related branches found
No related tags found
No related merge requests found
...@@ -758,6 +758,7 @@ END SUBROUTINE WRITE_SURFX2COV_MNH ...@@ -758,6 +758,7 @@ END SUBROUTINE WRITE_SURFX2COV_MNH
!! original 01/08/03 !! original 01/08/03
!! G.TANGUY 03/2009 add replace ZUNDEF(surfex) by XUNDEF(MNH) !! G.TANGUY 03/2009 add replace ZUNDEF(surfex) by XUNDEF(MNH)
!! Philippe Wautelet: 05/2016-04/2018: new data structures and calls for I/O !! Philippe Wautelet: 05/2016-04/2018: new data structures and calls for I/O
! P. Wautelet 07/02/2023: treat 2D case separately
!---------------------------------------------------------------------------- !----------------------------------------------------------------------------
! !
!* 0. DECLARATIONS !* 0. DECLARATIONS
...@@ -843,12 +844,24 @@ IF (HDIR=='H' .OR. HDIR=='A') THEN ...@@ -843,12 +844,24 @@ IF (HDIR=='H' .OR. HDIR=='A') THEN
WHERE (ZWORK==ZUNDEF) ZWORK=XUNDEF WHERE (ZWORK==ZUNDEF) ZWORK=XUNDEF
! !
IF (HDIR=='H') THEN IF (HDIR=='H') THEN
CALL PREPARE_METADATA_WRITE_SURF(HREC,'XY',HCOMMENT,4,TYPEREAL,3,'WRITE_SURFX2_MNH',TZFIELD) IF ( SIZE(PFIELD,2) == 1 ) THEN
CALL IO_Field_write(TFILE_SURFEX,TZFIELD,ZWORK(1+NHALO:IIU-NHALO,1+NHALO:IJU-NHALO,:),KRESP) ! 2D case
CALL PREPARE_METADATA_WRITE_SURF(HREC,'XY',HCOMMENT,4,TYPEREAL,2,'WRITE_SURFX2_MNH',TZFIELD)
CALL IO_Field_write(TFILE_SURFEX,TZFIELD,ZWORK(1+NHALO:IIU-NHALO,1+NHALO:IJU-NHALO,1),KRESP)
ELSE
CALL PREPARE_METADATA_WRITE_SURF(HREC,'XY',HCOMMENT,4,TYPEREAL,3,'WRITE_SURFX2_MNH',TZFIELD)
CALL IO_Field_write(TFILE_SURFEX,TZFIELD,ZWORK(1+NHALO:IIU-NHALO,1+NHALO:IJU-NHALO,:),KRESP)
END IF
END IF END IF
IF (HDIR=='A') THEN IF (HDIR=='A') THEN
CALL PREPARE_METADATA_WRITE_SURF(HREC,'--',HCOMMENT,4,TYPEREAL,3,'WRITE_SURFX2_MNH',TZFIELD) IF ( SIZE(PFIELD,2) == 1 ) THEN
CALL IO_Field_write(TFILE_SURFEX,TZFIELD,ZWORK(:,:,:),KRESP) ! 2D case
CALL PREPARE_METADATA_WRITE_SURF(HREC,'--',HCOMMENT,4,TYPEREAL,2,'WRITE_SURFX2_MNH',TZFIELD)
CALL IO_Field_write(TFILE_SURFEX,TZFIELD,ZWORK(:,:,1),KRESP)
ELSE
CALL PREPARE_METADATA_WRITE_SURF(HREC,'--',HCOMMENT,4,TYPEREAL,3,'WRITE_SURFX2_MNH',TZFIELD)
CALL IO_Field_write(TFILE_SURFEX,TZFIELD,ZWORK(:,:,:),KRESP)
END IF
END IF END IF
! !
DEALLOCATE(ZWORK) DEALLOCATE(ZWORK)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment