write_diag_seb_isban.F90

Go to the documentation of this file.

!SFX_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
!SFX_LIC This is part of the SURFEX software governed by the CeCILL-C licence
!SFX_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt  
!SFX_LIC for details. version 1.
!     #########
      SUBROUTINE write_diag_seb_isba_n ( DTCO, DUO, U, NCHI, CHI, ID, NDST, GB, &
                                         IO, S, NP, NPE, HPROGRAM)
!     #################################
!
!!****  *WRITE_DIAG_SEB_ISBA* - writes the ISBA diagnostic fields
!!
!!    PURPOSE
!!    -------
!!
!!
!!**  METHOD
!!    ------
!!
!!    REFERENCE
!!    ---------
!!
!!
!!    AUTHOR
!!    ------
!!      V. Masson   *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    01/2004
!!      B. Decharme 06/2009  key to write (or not) patch result
!!      B. Decharme 08/2009  cumulative radiative budget
!!      B. Decharme  09/2012 : Bug in local variables declaration in PROVAR_TO_DIAG
!!      B. Decharme 09/2012  New diag :
!!                           carbon fluxes and reservoirs
!!                           soil liquid and ice water content in kg/m2 and m3/m3
!!      B. Decharme  06/13   Add diags (sublimation, lateral drainage)
!!                           All snow outputs noted SN
!!                           delete NWG_SIZE
!!      S. Belamari 06/2014 : Introduce GRESET to avoid errors due to NBLOCK=0
!!                            when coupled with ARPEGE/ALADIN/AROME
!!      P. Samuelsson 10/2014 MEB
!!      B. Decharme 02/2016 : NBLOCK instead of LCOUNTW for compilation in AAA
!-------------------------------------------------------------------------------
!
!*       0.    DECLARATIONS
!              ------------
!
USE modd_surfex_mpi, ONLY : nrank, npio
!
USE modd_data_cover_n, ONLY : data_cover_t
USE modd_diag_n, ONLY : diag_options_t
USE modd_surf_atm_n, ONLY : surf_atm_t
USE modd_ch_isba_n, ONLY : ch_isba_t, ch_isba_np_t
USE modd_surfex_n, ONLY : isba_diag_t
USE modd_dst_n, ONLY : dst_np_t
USE modd_gr_biog_n, ONLY : gr_biog_t
USE modd_isba_options_n, ONLY : isba_options_t
USE modd_isba_n, ONLY : isba_np_t, isba_p_t, isba_npe_t, isba_pe_t, isba_s_t
!
#ifdef SFX_ARO
USE modd_io_surf_aro,   ONLY : nblock
#endif
!
USE modd_xios, ONLY : lallow_add_dim, yground_layer_dim_name, &
                      ywground_layer_dim_name, ywiground_layer_dim_name, &
                      yswband_dim_name
!
USE modd_surf_par,   ONLY : xundef, nundef
!
USE modd_csts,       ONLY : xrholw, xtt, xlmtt
!
USE modd_dst_surf
!
USE modd_agri,     ONLY : lagrip
!
USE mode_diag
!
USE modi_init_io_surf_n
USE modi_write_surf
USE modi_end_io_surf_n
USE modi_write_field_1d_patch
!
#ifdef SFX_OL
USE modd_io_surf_ol, ONLY : ldef
USE modn_io_offline, ONLY : xtstep_output
#endif
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*       0.1   Declarations of arguments
!              -------------------------
!
!
TYPE(data_cover_t), INTENT(INOUT) :: DTCO
TYPE(diag_options_t), INTENT(INOUT) :: DUO
TYPE(surf_atm_t), INTENT(INOUT) :: U
TYPE(ch_isba_np_t), INTENT(INOUT) :: NCHI
TYPE(ch_isba_t), INTENT(INOUT) :: CHI
TYPE(isba_diag_t), INTENT(INOUT) :: ID
TYPE(dst_np_t), INTENT(INOUT) :: NDST
TYPE(gr_biog_t), INTENT(INOUT) :: GB
TYPE(isba_options_t), INTENT(INOUT) :: IO
TYPE(isba_s_t), INTENT(INOUT) :: S
TYPE(isba_np_t), INTENT(INOUT) :: NP
TYPE(isba_npe_t), INTENT(INOUT) :: NPE
!
 CHARACTER(LEN=6),  INTENT(IN)  :: HPROGRAM ! program calling
!
!*       0.2   Declarations of local variables
!              -------------------------------
!
TYPE(isba_p_t), POINTER :: PK
TYPE(isba_pe_t), POINTER :: PEK
!
INTEGER           :: IRESP          ! IRESP  : return-code if a problem appears
 CHARACTER(LEN=12) :: YRECFM         ! Name of the article to be write
 CHARACTER(LEN=100):: YCOMMENT       ! Comment string
 CHARACTER(LEN=2)  :: YNUM
!
LOGICAL           :: GRESET
INTEGER           :: JSV, JSW, JP, ISIZE
INTEGER           :: ISIZE_LMEB_PATCH   ! Number of patches where multi-energy balance should be applied
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
!-------------------------------------------------------------------------------
!
!         Initialisation for IO
!
IF (lhook) CALL dr_hook('WRITE_DIAG_SEB_ISBA_N',0,zhook_handle)
!
isize_lmeb_patch=count(io%LMEB_PATCH(:))
!
greset=.true.
#ifdef SFX_ARO
greset=(nblock>0)
#endif
#ifdef SFX_OL
IF (ldef) greset = .false.
#endif
!
#ifdef SFX_OL
IF (id%O%LSURF_BUDGET .AND. duo%LRESETCUMUL .AND. id%O%LSURF_BUDGETC .AND. .NOT.ldef) THEN
  !
  ! Output variables are not instantaneous but averaged over the output time step
  ! Fluxes by patch
  DO jp = 1,io%NPATCH
    CALL avg_diag_tstep_surf(xtstep_output, id%NDC%AL(jp), id%ND%AL(jp))
  ENDDO
  CALL avg_diag_tstep_surf(xtstep_output, id%DC, id%D)
  !  
  IF (id%DE%LSURF_EVAP_BUDGET) THEN
    DO jp = 1,io%NPATCH
      CALL avg_diag_tstep_evap(xtstep_output, id%NDEC%AL(jp), id%NDE%AL(jp))
    ENDDO
    CALL avg_diag_tstep_evap(xtstep_output, id%DEC, id%DE)
    !
    IF (id%DE%LWATER_BUDGET) THEN
      DO jp = 1,io%NPATCH
        CALL avg_diag_tstep_water(xtstep_output, id%NDEC%AL(jp), id%NDE%AL(jp))
      ENDDO
      CALL avg_diag_tstep_water(xtstep_output, id%DEC, id%DE)
    ENDIF
  ENDIF
  !
END IF
#endif
!
IF ( id%DM%LPROSNOW ) THEN
  CALL init_io_surf_n(dtco, u, hprogram,'NATURE','ISBA  ','WRITE','ISBA_PROGNOSTIC.OUT.nc')
ELSE
  CALL init_io_surf_n(dtco, u, hprogram,'NATURE','ISBA  ','WRITE','ISBA_DIAGNOSTICS.OUT.nc')
ENDIF
!
!-------------------------------------------------------------------------------
!
!*       2.     Richardson number :
!               -----------------
!
IF (id%O%N2M>=1) THEN
  !
  yrecfm='RI_ISBA'
  ycomment='Richardson number over tile nature'
  CALL write_surf(duo%CSELECT,hprogram,yrecfm,id%D%XRI(:),iresp,hcomment=ycomment)
  !
END IF
!
!*       3.     Energy fluxes :
!               -------------
!
IF (id%O%LSURF_BUDGET) THEN
  !
  yrecfm='TALB_ISBA'
  ycomment='total albedo over tile nature (-)'
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,id%D%XALBT(:),iresp,hcomment=ycomment)
  !  
  yrecfm='RN_ISBA'
  ycomment='Net radiation over tile nature'//' (W/m2)'
  CALL write_surf(duo%CSELECT,hprogram,yrecfm,id%D%XRN(:),iresp,hcomment=ycomment)
  !
  yrecfm='H_ISBA'
  ycomment='Sensible heat flux over tile nature'//' (W/m2)'
  CALL write_surf(duo%CSELECT,hprogram,yrecfm,id%D%XH(:),iresp,hcomment=ycomment)
  !
  yrecfm='LE_ISBA'
  ycomment='total latent heat flux over tile nature'//' (W/m2)'
  CALL write_surf(duo%CSELECT,hprogram,yrecfm,id%D%XLE(:),iresp,hcomment=ycomment)
  !
  yrecfm='LEI_ISBA'
  ycomment='sublimation latent heat flux over tile nature'//' (W/m2)'
  CALL write_surf(duo%CSELECT,hprogram,yrecfm,id%D%XLEI(:),iresp,hcomment=ycomment)
  !
  yrecfm='GFLUX_ISBA'
  ycomment='Ground flux over tile nature'//' (W/m2)'
  CALL write_surf(duo%CSELECT,hprogram,yrecfm,id%D%XGFLUX(:),iresp,hcomment=ycomment)
  !
  yrecfm='EVAP_ISBA'
  ycomment='total evaporative flux for tile nature'//' (Kg/m2/s)'
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,id%D%XEVAP(:),iresp,hcomment=ycomment)
  !
  yrecfm='SUBL_ISBA'
  ycomment='sublimation flux for tile nature'//' (Kg/m2/s)'
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,id%D%XSUBL(:),iresp,hcomment=ycomment)
  ! 
  IF (id%O%LRAD_BUDGET  .OR. (id%O%LSURF_BUDGETC .AND. .NOT.duo%LRESET_BUDGETC)) THEN
    !
    yrecfm='SWD_ISBA'
    ycomment='short wave downward radiation over tile nature'//' (W/m2)'
    CALL write_surf(duo%CSELECT,hprogram,yrecfm,id%D%XSWD(:),iresp,hcomment=ycomment)
    !
    yrecfm='SWU_ISBA'
    ycomment='short wave upward radiation over tile nature'//' (W/m2)'
    CALL write_surf(duo%CSELECT,hprogram,yrecfm,id%D%XSWU(:),iresp,hcomment=ycomment)
    !
    yrecfm='LWD_ISBA'
    ycomment='long wave downward radiation over tile nature'//' (W/m2)'
    CALL write_surf(duo%CSELECT,hprogram,yrecfm,id%D%XLWD(:),iresp,hcomment=ycomment)
    !
    yrecfm='LWU_ISBA'
    ycomment='long wave upward radiation over tile nature'//' (W/m2)'
    CALL write_surf(duo%CSELECT,hprogram,yrecfm,id%D%XLWU(:),iresp,hcomment=ycomment)
    !    
    IF (lallow_add_dim)  THEN
      !
      yrecfm='SWD_ISBA_'
      ycomment='X_Y_'//yrecfm//' (W/m2)'
      CALL write_surf(duo%CSELECT,&
           hprogram,yrecfm,id%D%XSWBD(:,:),iresp,hcomment=ycomment, hnam_dim=yswband_dim_name)
      !
      yrecfm='SWU_ISBA_'
      ycomment='X_Y_'//yrecfm//' (W/m2)'
      CALL write_surf(duo%CSELECT,&
           hprogram,yrecfm,id%D%XSWBU(:,:),iresp,hcomment=ycomment, hnam_dim=yswband_dim_name)
      !
    ELSE
      !
      DO jsw=1, SIZE(id%D%XSWBD,2)
        ynum=achar(48+jsw)
        !
        yrecfm='SWD_ISBA_'//ynum
        ycomment='short wave downward radiation over tile nature for spectral band'//ynum//' (W/m2)'
        CALL write_surf(duo%CSELECT,hprogram,yrecfm,id%D%XSWBD(:,jsw),iresp,hcomment=ycomment)
        !
        yrecfm='SWU_ISBA_'//ynum
        ycomment='short wave upward radiation over tile nature for spectral band'//ynum//' (W/m2)'
        CALL write_surf(duo%CSELECT,hprogram,yrecfm,id%D%XSWBU(:,jsw),iresp,hcomment=ycomment)
        !
      ENDDO
      !
    ENDIF
    !
  ENDIF
  !
  yrecfm='FMU_ISBA'
  ycomment='u component of wind stress'//' (Pa)'  
  CALL write_surf(duo%CSELECT,hprogram,yrecfm,id%D%XFMU(:),iresp,hcomment=ycomment)
  !
  yrecfm='FMV_ISBA'
  ycomment='v component of wind stress'//' (Pa)'  
  CALL write_surf(duo%CSELECT,hprogram,yrecfm,id%D%XFMV(:),iresp,hcomment=ycomment)
  !
END IF
!
!*       4.    Specific Energy fluxes :(for each patch)
!              ----------------------------------------
!
IF (id%DE%LSURF_EVAP_BUDGET) CALL write_evap_bud(id%DE,"_ISBA ",.false.)
!
!*       6.     parameters at 2 and 10 meters :
!               -------------------------------
!
IF (id%O%N2M>=1) THEN
  !
  yrecfm='T2M_ISBA'
  ycomment='X_Y_'//yrecfm//' (K)'
  CALL write_surf(duo%CSELECT,hprogram,yrecfm,id%D%XT2M(:),iresp,hcomment=ycomment)
  !
  yrecfm='T2MMIN_ISBA'
  ycomment='X_Y_'//yrecfm//' (K)'
  CALL write_surf(duo%CSELECT,hprogram,yrecfm,id%D%XT2M_MIN(:),iresp,hcomment=ycomment)
  IF(greset)id%D%XT2M_MIN(:)=xundef
  !
  yrecfm='T2MMAX_ISBA'
  ycomment='X_Y_'//yrecfm//' (K)'
  CALL write_surf(duo%CSELECT,hprogram,yrecfm,id%D%XT2M_MAX(:),iresp,hcomment=ycomment)
  IF(greset)id%D%XT2M_MAX(:)=-xundef
  !
  yrecfm='Q2M_ISBA'
  ycomment='X_Y_'//yrecfm//' (KG/KG)'
  CALL write_surf(duo%CSELECT,hprogram,yrecfm,id%D%XQ2M(:),iresp,hcomment=ycomment)
  !
  yrecfm='HU2M_ISBA'
  ycomment='X_Y_'//yrecfm//' (-)'
  CALL write_surf(duo%CSELECT,hprogram,yrecfm,id%D%XHU2M(:),iresp,hcomment=ycomment)
  !
  yrecfm='HU2MMIN_ISBA'
  ycomment='X_Y_'//yrecfm//' (-)'
  CALL write_surf(duo%CSELECT,hprogram,yrecfm,id%D%XHU2M_MIN(:),iresp,hcomment=ycomment)
  IF(greset)id%D%XHU2M_MIN(:)=xundef
  !
  yrecfm='HU2MMAX_ISBA'
  ycomment='X_Y_'//yrecfm//' (-)'
  CALL write_surf(duo%CSELECT,hprogram,yrecfm,id%D%XHU2M_MAX(:),iresp,hcomment=ycomment)
  IF(greset)id%D%XHU2M_MAX(:)=-xundef
  !
  yrecfm='ZON10M_ISBA'
  ycomment='X_Y_'//yrecfm//' (M/S)'
  CALL write_surf(duo%CSELECT,hprogram,yrecfm,id%D%XZON10M(:),iresp,hcomment=ycomment)
  !
  yrecfm='MER10M_ISBA'
  ycomment='X_Y_'//yrecfm//' (M/S)'
  CALL write_surf(duo%CSELECT,hprogram,yrecfm,id%D%XMER10M(:),iresp,hcomment=ycomment)
  !
  yrecfm='W10M_ISBA'
  ycomment='X_Y_'//yrecfm//' (M/S)'
  CALL write_surf(duo%CSELECT,hprogram,yrecfm,id%D%XWIND10M(:),iresp,hcomment=ycomment)
  !
  yrecfm='W10MMAX_ISBA'
  ycomment='X_Y_'//yrecfm//' (M/S)'
  CALL write_surf(duo%CSELECT,hprogram,yrecfm,id%D%XWIND10M_MAX(:),iresp,hcomment=ycomment)
  IF(greset)id%D%XWIND10M_MAX(:)=0.0
  !
  yrecfm='SFCO2_ISBA'
  ycomment='X_Y_'//yrecfm//' (M.kgCO2.S-1.kgAIR-1)'
  CALL write_surf(duo%CSELECT,hprogram,yrecfm,id%D%XSFCO2(:),iresp,hcomment=ycomment)
  !  
END IF
!----------------------------------------------------------------------------
!
!*       7.     Transfer coefficients
!               ---------------------
!
IF (id%O%LCOEF) THEN
  !
  yrecfm='CD_ISBA'
  ycomment='X_Y_'//yrecfm
  CALL write_surf(duo%CSELECT,hprogram,yrecfm,id%D%XCD(:),iresp,hcomment=ycomment)
  !
  yrecfm='CH_ISBA'
  ycomment='X_Y_'//yrecfm
  CALL write_surf(duo%CSELECT,hprogram,yrecfm,id%D%XCH(:),iresp,hcomment=ycomment)
  !
  yrecfm='CE_ISBA'
  ycomment='X_Y_'//yrecfm
  CALL write_surf(duo%CSELECT,hprogram,yrecfm,id%D%XCE(:),iresp,hcomment=ycomment)
  !
  yrecfm='Z0_ISBA'
  ycomment='X_Y_'//yrecfm//' (M)'
  CALL write_surf(duo%CSELECT,hprogram,yrecfm,id%D%XZ0(:),iresp,hcomment=ycomment)
  !
  yrecfm='Z0H_ISBA'
  ycomment='X_Y_'//yrecfm//' (M)'
  CALL write_surf(duo%CSELECT,hprogram,yrecfm,id%D%XZ0H(:),iresp,hcomment=ycomment)
  !
ENDIF
!
!----------------------------------------------------------------------------
!
!*       8.     Surface humidity
!               ----------------
IF (id%O%LSURF_VARS) THEN
  !
  yrecfm='QS_ISBA'
  ycomment='X_Y_'//yrecfm//' (KG/KG)'
  CALL write_surf(duo%CSELECT,hprogram,yrecfm,id%D%XQS(:),iresp,hcomment=ycomment)
  !
ENDIF
!
!----------------------------------------------------------------------------
!
isize = u%NSIZE_NATURE
!
!User want (or not) patch output
IF (id%O%LPATCH_BUDGET.AND.(io%NPATCH >1)) THEN
  !
  !*      10.     Richardson number (for each patch)
  !               -----------------
  !
  IF (id%O%N2M>=1) THEN
    !  
    yrecfm='RI_'
    ycomment='X_Y_'//yrecfm      
    DO jp=1,io%NPATCH
      CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
              np%AL(jp)%NR_P,id%ND%AL(jp)%XRI(:),isize,s%XWORK_WR)
    ENDDO
    !
  END IF
  !
  !*       11.     Energy fluxes :(for each patch)
  !                -------------
  !
  IF (id%O%LSURF_BUDGET) THEN
    !
    yrecfm='TALB_'
    ycomment='total albedo per patch'
    DO jp=1,io%NPATCH
      CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
              np%AL(jp)%NR_P,id%ND%AL(jp)%XALBT(:),isize,s%XWORK_WR)
    ENDDO      
    !
    yrecfm='RN_'
    ycomment='X_Y_'//yrecfm//' (W/m2)'
    DO jp=1,io%NPATCH
      CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
              np%AL(jp)%NR_P,id%ND%AL(jp)%XRN(:),isize,s%XWORK_WR)
    ENDDO      
    !
    yrecfm='H_'
    ycomment='X_Y_'//yrecfm//' (W/m2)'
    DO jp=1,io%NPATCH
        CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
            np%AL(jp)%NR_P,id%ND%AL(jp)%XH(:),isize,s%XWORK_WR)
    ENDDO      
    !
    yrecfm='LE_'
    ycomment='X_Y_'//yrecfm//' (W/m2)'
    DO jp = 1,io%NPATCH
      CALL  write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
              np%AL(jp)%NR_P,npe%AL(jp)%XLE(:),isize,s%XWORK_WR)  
    ENDDO      
    !
    yrecfm='LEI_'
    ycomment='X_Y_'//yrecfm//' (W/m2)'      
    DO jp=1,io%NPATCH
      CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
              np%AL(jp)%NR_P,id%ND%AL(jp)%XLEI(:),isize,s%XWORK_WR)
    ENDDO  
    !
    yrecfm='GFLUX_'
    ycomment='X_Y_'//yrecfm//' (W/m2)'      
    DO jp=1,io%NPATCH
      CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
              np%AL(jp)%NR_P,id%ND%AL(jp)%XGFLUX(:),isize,s%XWORK_WR)
    ENDDO   
    !
    yrecfm='EVAP_'
    ycomment='X_Y_'//yrecfm//' (W/m2)'      
    DO jp=1,io%NPATCH
      CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
              np%AL(jp)%NR_P,id%ND%AL(jp)%XEVAP(:),isize,s%XWORK_WR)
    ENDDO    
    !
    yrecfm='SUBL_'
    ycomment='X_Y_'//yrecfm//' (W/m2)'      
    DO jp=1,io%NPATCH
      CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
              np%AL(jp)%NR_P,id%ND%AL(jp)%XSUBL(:),isize,s%XWORK_WR)
    ENDDO   
    !     
    IF (id%O%LRAD_BUDGET .OR. (id%O%LSURF_BUDGETC .AND. .NOT.duo%LRESET_BUDGETC)) THEN
      !
      yrecfm='SWD_'
      ycomment='X_Y_'//yrecfm//' (W/m2)'        
      DO jp=1,io%NPATCH
      CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
                np%AL(jp)%NR_P,id%ND%AL(jp)%XSWD(:),isize,s%XWORK_WR)
      ENDDO       
      !
      yrecfm='SWU_'
      ycomment='X_Y_'//yrecfm//' (W/m2)'        
      DO jp=1,io%NPATCH
        CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
                np%AL(jp)%NR_P,id%ND%AL(jp)%XSWU(:),isize,s%XWORK_WR)
      ENDDO
      !
      yrecfm='LWD_'
      ycomment='X_Y_'//yrecfm//' (W/m2)'        
      DO jp=1,io%NPATCH
        CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
                np%AL(jp)%NR_P,id%ND%AL(jp)%XLWD(:),isize,s%XWORK_WR)
      ENDDO  
      !
      yrecfm='LWU_'
      ycomment='X_Y_'//yrecfm//' (W/m2)'        
      DO jp=1,io%NPATCH
        CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
                np%AL(jp)%NR_P,id%ND%AL(jp)%XLWU(:),isize,s%XWORK_WR)
      ENDDO
      !
      DO jsw=1, SIZE(id%D%XSWBD,2)
        ynum=achar(48+jsw)
        !
        yrecfm='SWD_'//ynum
        ycomment='X_Y_'//yrecfm//' (W/m2)'          
        DO jp=1,io%NPATCH
          CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
              np%AL(jp)%NR_P,id%ND%AL(jp)%XSWD(:),isize,s%XWORK_WR)
        ENDDO
        !
        yrecfm='SWU_'//ynum
        ycomment='X_Y_'//yrecfm//' (W/m2)'          
        DO jp=1,io%NPATCH
          CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
              np%AL(jp)%NR_P,id%ND%AL(jp)%XSWU(:),isize,s%XWORK_WR)
        ENDDO
        !
      ENDDO
      !
    ENDIF
    !
    yrecfm='FMU_'
    ycomment='X_Y_'//yrecfm//' (Pa)'      
    DO jp=1,io%NPATCH
      CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
              np%AL(jp)%NR_P,id%ND%AL(jp)%XFMU(:),isize,s%XWORK_WR)
    ENDDO
    !
    yrecfm='FMV_'
    ycomment='X_Y_'//yrecfm//' (Pa)'      
    DO jp=1,io%NPATCH
      CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
              np%AL(jp)%NR_P,id%ND%AL(jp)%XFMV(:),isize,s%XWORK_WR)
    ENDDO
    !
  END IF
  !
  !*       12.    Specific Energy fluxes :(for each patch)
  !               ----------------------------------------
  !
  IF (id%DE%LSURF_EVAP_BUDGET) CALL write_evap_bud_patch(id%NDE,'_ ',.false.)
  !
  !*       13.    surface temperature parameters at 2 and 10 meters (for each patch):
  !               -------------------------------------------------------------------
  !
  IF (id%O%N2M>=1.AND..NOT.io%LCANOPY) THEN
    !
    yrecfm='T2M_'
    ycomment='X_Y_'//yrecfm//' (K)'      
    DO jp=1,io%NPATCH
      CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
              np%AL(jp)%NR_P,id%ND%AL(jp)%XT2M(:),isize,s%XWORK_WR)
    ENDDO
    !
    yrecfm='T2MMIN_'
    ycomment='X_Y_'//yrecfm//' (K)' 
    DO jp=1,io%NPATCH
      CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
              np%AL(jp)%NR_P,id%ND%AL(jp)%XT2M_MIN(:),isize,s%XWORK_WR)
      IF (greset) id%ND%AL(jp)%XT2M_MIN(:)=xundef        
    ENDDO        
    !      
    yrecfm='T2MMAX_'
    ycomment='X_Y_'//yrecfm//' (K)'
    DO jp=1,io%NPATCH
      CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
              np%AL(jp)%NR_P,id%ND%AL(jp)%XT2M_MAX(:),isize,s%XWORK_WR)
      IF (greset) id%ND%AL(jp)%XT2M_MAX(:)=-xundef        
    ENDDO
    ! 
    yrecfm='Q2M_'
    ycomment='X_Y_'//yrecfm//' (KG/KG)'      
    DO jp=1,io%NPATCH
      CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
              np%AL(jp)%NR_P,id%ND%AL(jp)%XQ2M(:),isize,s%XWORK_WR)
    ENDDO
    !
    yrecfm='HU2M_'
    ycomment='X_Y_'//yrecfm//' (PERCENT)'      
    DO jp=1,io%NPATCH
      CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
              np%AL(jp)%NR_P,id%ND%AL(jp)%XHU2M(:),isize,s%XWORK_WR)
    ENDDO
    !
    yrecfm='ZON10M_'
    ycomment='X_Y_'//yrecfm//' (M/S)'      
    DO jp=1,io%NPATCH
      CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
              np%AL(jp)%NR_P,id%ND%AL(jp)%XZON10M(:),isize,s%XWORK_WR)
    ENDDO
    !
    yrecfm='MER10M_'
    ycomment='X_Y_'//yrecfm//' (M/S)'
    DO jp=1,io%NPATCH
      CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
              np%AL(jp)%NR_P,id%ND%AL(jp)%XMER10M(:),isize,s%XWORK_WR)
    ENDDO
    !    
    yrecfm='W10M_'
    ycomment='X_Y_'//yrecfm//' (M/S)'
    DO jp=1,io%NPATCH
      CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
              np%AL(jp)%NR_P,id%ND%AL(jp)%XWIND10M(:),isize,s%XWORK_WR)
    ENDDO
    !
  END IF
  !
ENDIF
!
!----------------------------------------------------------------------------
!
!*       9.     Diag of prognostic fields
!               -------------------------
!
IF (duo%LPROVAR_TO_DIAG) CALL provar_to_diag
!
!----------------------------------------------------------------------------
!
!*       15.     chemical diagnostics:
!               --------------------
!
IF (chi%SVI%NBEQ>0 .AND. chi%CCH_DRY_DEP=="WES89 ") THEN
  !
  DO jsv = 1,SIZE(chi%CCH_NAMES,1)
    yrecfm='DVNT'//trim(chi%CCH_NAMES(jsv))
    WRITE(ycomment,'(A13,I3.3)')'(m/s) DV_NAT_',jsv
    DO jp = 1,io%NPATCH
      CALL  write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
                                 np%AL(jp)%NR_P,nchi%AL(jp)%XDEP(:,jsv),isize,s%XWORK_WR)
    ENDDO
  END DO
  !
ENDIF
!
IF (chi%SVI%NBEQ>0 .AND. chi%LCH_BIO_FLUX) THEN
  !
  IF (ASSOCIATED(gb%XFISO)) THEN
    yrecfm='FISO'
    WRITE(ycomment,'(A21)')'FISO (molecules/m2/s)'
    CALL write_surf(duo%CSELECT,hprogram,yrecfm,gb%XFISO(:),iresp,hcomment=ycomment)
  END IF
  !
  IF (ASSOCIATED(gb%XFISO)) THEN
    yrecfm='FMONO'
    WRITE(ycomment,'(A22)')'FMONO (molecules/m2/s)'
    CALL write_surf(duo%CSELECT,hprogram,yrecfm,gb%XFMONO(:),iresp,hcomment=ycomment)
  END IF
  !
ENDIF
!
IF (chi%LCH_NO_FLUX) THEN
  IF (ASSOCIATED(gb%XNOFLUX)) THEN
    yrecfm='NOFLUX'
    WRITE(ycomment,'(A21)')'NOFLUX (molecules/m2/s)'
    CALL write_surf(duo%CSELECT, hprogram,yrecfm,gb%XNOFLUX(:),iresp,hcomment=ycomment)
  END IF
END IF
!
IF (chi%SVI%NDSTEQ > 0)THEN
  !
  DO jsv = 1,ndstmde ! for all dust modes
    WRITE(yrecfm,'(A5,I3.3)')'F_DST',jsv
    ycomment='X_Y_'//yrecfm//' (kg/m2/s)'
    DO jp = 1,io%NPATCH
      CALL  write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
                                 np%AL(jp)%NR_P,ndst%AL(jp)%XSFDST(:,jsv),isize,s%XWORK_WR)
    ENDDO    
  END DO
  !
ENDIF
!
!----------------------------------------------------------------------------
!
!*       5.    Cumulated Energy fluxes
!              -----------------------
!
 CALL end_io_surf_n(hprogram)
 CALL init_io_surf_n(dtco, u, hprogram,'NATURE','ISBA  ','WRITE','ISBA_DIAG_CUMUL.OUT.nc')
!
IF (id%O%LSURF_BUDGETC) THEN
  !
  CALL write_evap_bud(id%DEC,"C_ISBA",(id%O%LSURF_BUDGETC .AND. .NOT.duo%LRESET_BUDGETC))
  !
  IF(io%LGLACIER)THEN
    yrecfm='ICE_FC_ISBA'
    ycomment='X_Y_'//yrecfm//' (Kg/m2)'
    CALL write_surf(duo%CSELECT,hprogram,yrecfm,id%DEC%XICEFLUX(:),iresp,hcomment=ycomment)
  ENDIF
  !  
  yrecfm='RNC_ISBA'
  ycomment='X_Y_'//yrecfm//' (J/m2)'
  CALL write_surf(duo%CSELECT,hprogram,yrecfm,id%DC%XRN(:),iresp,hcomment=ycomment)
  !
  yrecfm='HC_ISBA'
  ycomment='X_Y_'//yrecfm//' (J/m2)'
  CALL write_surf(duo%CSELECT,hprogram,yrecfm,id%DC%XH(:),iresp,hcomment=ycomment)
  !
  yrecfm='LEC_ISBA'
  ycomment='X_Y_'//yrecfm//' (J/m2)'
  CALL write_surf(duo%CSELECT,hprogram,yrecfm,id%DC%XLE(:),iresp,hcomment=ycomment)
  !
  yrecfm='LEIC_ISBA'
  ycomment='X_Y_'//yrecfm//' (J/m2)'
  CALL write_surf(duo%CSELECT,hprogram,yrecfm,id%DC%XLEI(:),iresp,hcomment=ycomment)
  !
  yrecfm='GFLUXC_ISBA'
  ycomment='X_Y_'//yrecfm//' (J/m2)'
  CALL write_surf(duo%CSELECT,hprogram,yrecfm,id%DC%XGFLUX(:),iresp,hcomment=ycomment)
  !
  yrecfm='EVAPC_ISBA'
  ycomment='total evaporative flux for tile nature'//' (Kg/m2/s)'
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,id%DC%XEVAP(:),iresp,hcomment=ycomment)
  !
  yrecfm='SUBLC_ISBA'
  ycomment='sublimation flux for tile nature'//' (Kg/m2/s)'
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,id%DC%XSUBL(:),iresp,hcomment=ycomment)
  !  
  IF (id%O%LRAD_BUDGET .OR. (id%O%LSURF_BUDGETC .AND. .NOT.duo%LRESET_BUDGETC)) THEN
    !
    yrecfm='SWDC_ISBA'
    ycomment='X_Y_'//yrecfm//' (J/m2)'
    CALL write_surf(duo%CSELECT,hprogram,yrecfm,id%DC%XSWD(:),iresp,hcomment=ycomment)
    !
    yrecfm='SWUC_ISBA'
    ycomment='X_Y_'//yrecfm//' (J/m2)'
    CALL write_surf(duo%CSELECT,hprogram,yrecfm,id%DC%XSWU(:),iresp,hcomment=ycomment)
    !
    yrecfm='LWDC_ISBA'
    ycomment='X_Y_'//yrecfm//' (J/m2)'
    CALL write_surf(duo%CSELECT,hprogram,yrecfm,id%DC%XLWD(:),iresp,hcomment=ycomment)
    !
    yrecfm='LWUC_ISBA'
    ycomment='X_Y_'//yrecfm//' (J/m2)'
    CALL write_surf(duo%CSELECT,hprogram,yrecfm,id%DC%XLWU(:),iresp,hcomment=ycomment)
    !
  ENDIF
  !
  yrecfm='FMUC_ISBA'
  ycomment='X_Y_'//yrecfm//' (Pa.s)'  
  CALL write_surf(duo%CSELECT,hprogram,yrecfm,id%DC%XFMU(:),iresp,hcomment=ycomment)
  !
  yrecfm='FMVC_ISBA'
  ycomment='X_Y_'//yrecfm//' (Pa.s)'  
  CALL write_surf(duo%CSELECT,hprogram,yrecfm,id%DC%XFMV(:),iresp,hcomment=ycomment)
  !
ENDIF
!
!----------------------------------------------------------------------------
!
!User want (or not) patch output
IF (id%O%LPATCH_BUDGET.AND.(io%NPATCH >1)) THEN
    !
    !*       14.    Cumulated Energy fluxes :(for each patch)
    !               -----------------------------------------
    !
    IF (id%O%LSURF_BUDGETC) THEN
      !
      CALL write_evap_bud_patch(id%NDEC,'C_',(id%O%LSURF_BUDGETC .AND. .NOT.duo%LRESET_BUDGETC))
      !  
      IF(io%LGLACIER)THEN              
        yrecfm='ICE_FC_'
        ycomment='X_Y_'//yrecfm//' (Kg/m2)'
        DO jp=1,io%NPATCH
          CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
                np%AL(jp)%NR_P,id%NDEC%AL(jp)%XICEFLUX(:),isize,s%XWORK_WR)
        ENDDO
      ENDIF
      !      
      yrecfm='RNC_'
      ycomment='X_Y_'//yrecfm//' (J/m2)'
      DO jp=1,io%NPATCH
        CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
                np%AL(jp)%NR_P,id%NDC%AL(jp)%XRN(:),isize,s%XWORK_WR)
      ENDDO
      !     
      yrecfm='HC_'
      ycomment='X_Y_'//yrecfm//' (J/m2)'
      DO jp=1,io%NPATCH
        CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
                np%AL(jp)%NR_P,id%NDC%AL(jp)%XH(:),isize,s%XWORK_WR)
      ENDDO
      ! 
      yrecfm='LEC_'
      ycomment='X_Y_'//yrecfm//' (J/m2)'
      DO jp=1,io%NPATCH
        CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
                np%AL(jp)%NR_P,id%NDC%AL(jp)%XLE(:),isize,s%XWORK_WR)
      ENDDO
      !   
      yrecfm='LEIC_'
      ycomment='X_Y_'//yrecfm//' (J/m2)'
      DO jp=1,io%NPATCH
        CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
                np%AL(jp)%NR_P,id%NDC%AL(jp)%XLEI(:),isize,s%XWORK_WR)
      ENDDO
      !
      yrecfm='GFLUXC_'
      ycomment='X_Y_'//yrecfm//' (J/m2)'
      DO jp=1,io%NPATCH
        CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
                np%AL(jp)%NR_P,id%NDC%AL(jp)%XGFLUX(:),isize,s%XWORK_WR)
      ENDDO
      !     
      yrecfm='EVAPC_'
      ycomment='X_Y_'//yrecfm//' (W/m2)'
      DO jp=1,io%NPATCH
        CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
                np%AL(jp)%NR_P,id%NDC%AL(jp)%XEVAP(:),isize,s%XWORK_WR)
      ENDDO
      !    
      yrecfm='SUBLC_'
      ycomment='X_Y_'//yrecfm//' (W/m2)'
      DO jp=1,io%NPATCH
        CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
                np%AL(jp)%NR_P,id%NDC%AL(jp)%XSUBL(:),isize,s%XWORK_WR)
      ENDDO      
      !    
      IF (id%O%LRAD_BUDGET .OR. (id%O%LSURF_BUDGETC .AND. .NOT.duo%LRESET_BUDGETC)) THEN
        !
        yrecfm='SWDC_'
        ycomment='X_Y_'//yrecfm//' (J/m2)'
        DO jp=1,io%NPATCH
          CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
                np%AL(jp)%NR_P,id%NDC%AL(jp)%XSWD(:),isize,s%XWORK_WR)
        ENDDO
        !      
        yrecfm='SWUC_'
        ycomment='X_Y_'//yrecfm//' (J/m2)'
        DO jp=1,io%NPATCH
          CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
                np%AL(jp)%NR_P,id%NDC%AL(jp)%XSWU(:),isize,s%XWORK_WR)
        ENDDO
        !     
        yrecfm='LWDC_'
        ycomment='X_Y_'//yrecfm//' (J/m2)'
        DO jp=1,io%NPATCH
          CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
                np%AL(jp)%NR_P,id%NDC%AL(jp)%XLWD(:),isize,s%XWORK_WR)
        ENDDO
        !     
        yrecfm='LWUC_'
        ycomment='X_Y_'//yrecfm//' (J/m2)'
        DO jp=1,io%NPATCH
          CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
                np%AL(jp)%NR_P,id%NDC%AL(jp)%XLWU(:),isize,s%XWORK_WR)
        ENDDO
        !
      ENDIF
      ! 
      yrecfm='FMUC_'
      ycomment='X_Y_'//yrecfm//' (Pa.s)'   
      DO jp=1,io%NPATCH
          CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
                np%AL(jp)%NR_P,id%NDC%AL(jp)%XFMU(:),isize,s%XWORK_WR)
      ENDDO
      ! 
      yrecfm='FMVC_'
      ycomment='X_Y_'//yrecfm//' (Pa.s)'
       DO jp=1,io%NPATCH
          CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
                np%AL(jp)%NR_P,id%NDC%AL(jp)%XFMV(:),isize,s%XWORK_WR)
      ENDDO
      !
    ENDIF
    !
    !-------------------------------------------------------------------------------
ENDIF
!
!         End of IO
!
 CALL end_io_surf_n(hprogram)
!
!-------------------------------------------------------------------------------
!
!User want (or not) patch output
!-------------------------------------------------------------------------------
!
IF ( duo%LRESETCUMUL .AND. id%O%LSURF_BUDGETC ) THEN
  !
  DO jp = 1,io%NPATCH
    CALL init_surf_bud(id%NDC%AL(jp),0.)
  ENDDO
  !  
  IF (id%DE%LSURF_EVAP_BUDGET) THEN
    DO jp = 1,io%NPATCH
      CALL init_evap_bud(id%NDEC%AL(jp))
    ENDDO
    !
    IF (id%DE%LWATER_BUDGET) THEN
      DO jp = 1,io%NPATCH
        CALL init_water_bud(id%NDEC%AL(jp))
      ENDDO
    ENDIF
  ENDIF
  !
END IF
!
IF (lhook) CALL dr_hook('WRITE_DIAG_SEB_ISBA_N',1,zhook_handle)
!
CONTAINS
!
!-------------------------------------------------------------------------------
!
SUBROUTINE write_evap_bud(DEA,HTERM,OFLAG)
!
USE modd_diag_evap_isba_n, ONLY : diag_evap_isba_t
!
TYPE(diag_evap_isba_t) :: DEA
 CHARACTER(LEN=6), INTENT(IN) :: HTERM
LOGICAL, INTENT(IN) :: OFLAG
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
IF (lhook) CALL dr_hook('WRITE_DIAG_SEB_ISBA_N:WRITE_EVAP_BUD',0,zhook_handle)
!
yrecfm='LEG'//trim(hterm)
ycomment='bare ground evaporation for tile nature'//' (W/m2)'
 CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XLEG(:),iresp,hcomment=ycomment)
!
yrecfm='LEGI'//trim(hterm)
ycomment='bare ground sublimation for tile nature'//' (W/m2)'
 CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XLEGI(:),iresp,hcomment=ycomment)
!
yrecfm='LEV'//trim(hterm)
ycomment='total vegetation evaporation for tile nature'//' (W/m2)'
 CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XLEV(:),iresp,hcomment=ycomment)
  !
yrecfm='LES'//trim(hterm)
ycomment='snow sublimation for tile nature'//' (W/m2)'  
 CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XLES(:),iresp,hcomment=ycomment)
  !
IF(npe%AL(1)%TSNOW%SCHEME=='3-L' .OR. npe%AL(1)%TSNOW%SCHEME=='CRO')THEN  
  yrecfm='LESL'//trim(hterm)
  ycomment='liquid water evaporation over snow for tile nature'//' (W/m2)'
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XLESL(:),iresp,hcomment=ycomment)
  yrecfm='SNDRIF'//trim(hterm)
  ycomment='blowing snow sublimation for tile nature'//' (Kg/m2/s)'
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XSNDRIFT(:),iresp,hcomment=ycomment)    
ENDIF
  !  
yrecfm='LER'//trim(hterm)
ycomment='canopy direct evaporation for tile nature'//' (W/m2)'
 CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XLER(:),iresp,hcomment=ycomment)
  !
yrecfm='LETR'//trim(hterm)
ycomment='vegetation transpiration for tile nature'//' (W/m2)'
 CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XLETR(:),iresp,hcomment=ycomment)
  !
yrecfm='DRAIN'//trim(hterm)
ycomment='drainage for tile nature'//' (Kg/m2/s)'
 CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XDRAIN(:),iresp,hcomment=ycomment)
  !
IF(io%CRUNOFF=='SGH'.AND.io%CISBA=='DIF')THEN
  yrecfm='QSB'//trim(hterm)
  ycomment='lateral subsurface flow for tile nature'//' (Kg/m2/s)'
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XQSB(:),iresp,hcomment=ycomment)
ENDIF
  !
yrecfm='RUNOFF'//trim(hterm)
ycomment='runoff for tile nature'//' (Kg/m2/s)'
 CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XRUNOFF(:),iresp,hcomment=ycomment)
  !
IF(io%CHORT=='SGH'.OR.io%CISBA=='DIF')THEN
  yrecfm='HORTON'//trim(hterm)
  ycomment='horton runoff for tile nature'//' (Kg/m2/s)'
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XHORT(:),iresp,hcomment=ycomment)
ENDIF
  !
yrecfm='DRIVEG'//trim(hterm)
ycomment='X_Y_'//yrecfm//' (Kg/m2/s)'
 CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XDRIP(:),iresp,hcomment=ycomment)
  !
yrecfm='RRVEG'//trim(hterm)
ycomment='X_Y_'//yrecfm//' (Kg/m2/s)'
 CALL write_surf(duo%CSELECT, &
                  hprogram,yrecfm,dea%XRRVEG(:),iresp,hcomment=ycomment)
  !
yrecfm='SNOMLT'//trim(hterm)
ycomment='snow melting rate'//' (Kg/m2/s)'
 CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XMELT(:),iresp,hcomment=ycomment)
  !
IF(lagrip)THEN
  yrecfm='IRRIG'//trim(hterm)
  ycomment='irrigation rate'//' (Kg/m2/s)'
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XIRRIG_FLUX(:),iresp,hcomment=ycomment)
ENDIF  
! MEB STUFF
IF (isize_lmeb_patch>0) THEN
  yrecfm='LELIT'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XLELITTER(:),iresp,hcomment=ycomment)
    !
  yrecfm='LELITI'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XLELITTERI(:),iresp,hcomment=ycomment)
    !
  yrecfm='DRIPLIT'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XDRIPLIT(:),iresp,hcomment=ycomment)
    !
  yrecfm='RRLIT'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XRRLIT(:),iresp,hcomment=ycomment)
    !        
  yrecfm='LEV_CV'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XLEV_CV(:),iresp,hcomment=ycomment)  
    !
  yrecfm='LES_CV'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XLES_CV(:),iresp,hcomment=ycomment)
    !
  yrecfm='LETR_CV'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XLETR_CV(:),iresp,hcomment=ycomment)
    !
  yrecfm='LER_CV'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XLER_CV(:),iresp,hcomment=ycomment)
    !    
   yrecfm='LE_CV'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XLE_CV(:),iresp,hcomment=ycomment)
    !
  yrecfm='H_CV'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XH_CV(:),iresp,hcomment=ycomment)
    !
  yrecfm='MELT_CV'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (kg/m2/s)'
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XMELT_CV(:),iresp,hcomment=ycomment)
    !
  yrecfm='FRZ_CV'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (kg/m2/s)'
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XFRZ_CV(:),iresp,hcomment=ycomment)
    !  
  yrecfm='LETR_GV'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XLETR_GV(:),iresp,hcomment=ycomment)
    !
  yrecfm='LER_GV'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XLER_GV(:),iresp,hcomment=ycomment)
    !
  yrecfm='LE_GV'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XLE_GV(:),iresp,hcomment=ycomment)  
    !
  yrecfm='H_GV'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XH_GV(:),iresp,hcomment=ycomment)  
    !
  yrecfm='LE_GN'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XLE_GN(:),iresp,hcomment=ycomment)  
    !
  yrecfm='H_GN'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XH_GN(:),iresp,hcomment=ycomment)
    !
  yrecfm='SR_GN'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (kg/m2/s)'
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XSR_GN(:),iresp,hcomment=ycomment)  
    !
  yrecfm='SWDN_GN'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XSWDOWN_GN(:),iresp,hcomment=ycomment)
    !
  yrecfm='LWDN_GN'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XLWDOWN_GN(:),iresp,hcomment=ycomment)  
    !
  yrecfm='LE_CA'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XLE_CA(:),iresp,hcomment=ycomment)
    !
  yrecfm='H_CA'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XH_CA(:),iresp,hcomment=ycomment)  
    !
  yrecfm='SWNT_V'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XSWNET_V(:),iresp,hcomment=ycomment)
    !
  yrecfm='SWNT_G'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XSWNET_G(:),iresp,hcomment=ycomment)
    !
  yrecfm='SWNT_N'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XSWNET_N(:),iresp,hcomment=ycomment)
    !
  yrecfm='SWNT_NS'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XSWNET_NS(:),iresp,hcomment=ycomment)
    !
  yrecfm='LWNT_V'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XLWNET_V(:),iresp,hcomment=ycomment)
    !
  yrecfm='LWNT_G'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XLWNET_G(:),iresp,hcomment=ycomment)
    !
  yrecfm='LWNT_N'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XLWNET_N(:),iresp,hcomment=ycomment)
    !
ENDIF
  ! END MEB STUFF
  !
IF(io%LFLOOD)THEN
  !        
  yrecfm='IFLOOD'//trim(hterm)
  ycomment='flood soil infiltration (Kg/m2/s)'    
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XIFLOOD(:),iresp,hcomment=ycomment)
    !
  yrecfm='PFLOOD'//trim(hterm)
  ycomment='intercepted precipitation by floodplains (Kg/m2/s)'    
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XPFLOOD(:),iresp,hcomment=ycomment)
    !
  yrecfm='LEF'//trim(hterm)
  ycomment='total floodplains evaporation (W/m2)'   
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XLE_FLOOD(:),iresp,hcomment=ycomment)
    !
  yrecfm='LEIF'//trim(hterm)
  ycomment='solid floodplains evaporation (W/m2)'    
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XLEI_FLOOD(:),iresp,hcomment=ycomment)
    !
ENDIF
  !
IF(io%CPHOTO/='NON')THEN
    !
  yrecfm='GPP'//trim(hterm)
  ycomment='gross primary production over tile nature (kgCO2/m2/s)'
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XGPP(:),iresp,hcomment=ycomment)
    !
  IF (hterm(1:1)=="C") THEN
    yrecfm='RC_AUTO'//trim(hterm(2:))
  ELSE
    yrecfm='R_AUTO'//trim(hterm)
  ENDIF      
  ycomment='autotrophic respiration over tile nature (kgCO2/m2/s)'
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XRESP_AUTO(:),iresp,hcomment=ycomment)
    !
  IF (hterm(1:1)=="C") THEN
    yrecfm='RC_ECO'//trim(hterm(2:))
  ELSE
    yrecfm='R_ECO'//trim(hterm)
  ENDIF      
  ycomment='ecosystem respiration over tile nature (kgCO2/m2/s)'
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XRESP_ECO(:),iresp,hcomment=ycomment)
    !
ENDIF
!  
IF(id%DE%LWATER_BUDGET .OR. oflag)THEN 
    !
  yrecfm='RAINF'//trim(hterm)
  ycomment='input rainfall rate (Kg/m2/s)'
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XRAINFALL(:),iresp,hcomment=ycomment)
    !
  yrecfm='SNOWF'//trim(hterm)
  ycomment='input snowfall rate (Kg/m2/s)'
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XSNOWFALL(:),iresp,hcomment=ycomment)
    !
  yrecfm='DWG'//trim(hterm)
  ycomment='change in liquid soil moisture (Kg/m2/s)'
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XDWG(:),iresp,hcomment=ycomment)
    !
  yrecfm='DWGI'//trim(hterm)
  ycomment='change in solid soil moisture (Kg/m2/s)'
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XDWGI(:),iresp,hcomment=ycomment)
    !
  yrecfm='DWR'//trim(hterm)
  ycomment='change in water on canopy (Kg/m2/s)'
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XDWR(:),iresp,hcomment=ycomment)
    !
  yrecfm='DSWE'//trim(hterm)
  ycomment='change in snow water equivalent (Kg/m2/s)'
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XDSWE(:),iresp,hcomment=ycomment)
    !
  yrecfm='WATBUD'//trim(hterm)
  ycomment='isba water budget as residue (Kg/m2/s)'
  CALL write_surf(duo%CSELECT, hprogram,yrecfm,dea%XWATBUD(:),iresp,hcomment=ycomment)
    !
ENDIF
!
IF (lhook) CALL dr_hook('WRITE_DIAG_SEB_ISBA_N:WRITE_EVAP_BUD',1,zhook_handle)

!
END SUBROUTINE write_evap_bud
!
SUBROUTINE write_evap_bud_patch(NDEA,HTERM,OFLAG)
!
USE modd_diag_evap_isba_n, ONLY : diag_evap_isba_np_t
!
TYPE(diag_evap_isba_np_t) :: NDEA
 CHARACTER(LEN=2), INTENT(IN) :: HTERM
LOGICAL, INTENT(IN) :: OFLAG
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
IF (lhook) CALL dr_hook('WRITE_DIAG_SEB_ISBA_N:WRITE_EVAP_BUD_PATCH',0,zhook_handle)
!
yrecfm='LEG'//trim(hterm)
ycomment='X_Y_'//yrecfm//' (W/m2)'
DO jp=1,io%NPATCH
  CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XLEG(:),isize,s%XWORK_WR)
ENDDO
!
yrecfm='LEGI'//trim(hterm)
ycomment='X_Y_'//yrecfm//' (W/m2)'
DO jp=1,io%NPATCH
  CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XLEGI(:),isize,s%XWORK_WR)
ENDDO
!
yrecfm='LEV'//trim(hterm)
ycomment='X_Y_'//yrecfm//' (W/m2)'
DO jp=1,io%NPATCH
  CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XLEV(:),isize,s%XWORK_WR)
ENDDO
!
 
yrecfm='LES'//trim(hterm)
ycomment='X_Y_'//yrecfm//' (W/m2)'
DO jp=1,io%NPATCH
  CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XLES(:),isize,s%XWORK_WR)
ENDDO
!
IF(npe%AL(1)%TSNOW%SCHEME=='3-L' .OR. npe%AL(1)%TSNOW%SCHEME=='CRO')THEN 
  yrecfm='LESL'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  DO jp=1,io%NPATCH
  CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XLESL(:),isize,s%XWORK_WR)
ENDDO
  yrecfm='SNDRIF'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (kg/m2/s)'
  DO jp=1,io%NPATCH
  CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XSNDRIFT(:),isize,s%XWORK_WR)
ENDDO
ENDIF
! 
yrecfm='LER'//trim(hterm)
ycomment='X_Y_'//yrecfm//' (W/m2)'
DO jp=1,io%NPATCH
  CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XLER(:),isize,s%XWORK_WR)
ENDDO
!
yrecfm='LETR'//trim(hterm)
ycomment='X_Y_'//yrecfm//' (W/m2)'
DO jp=1,io%NPATCH
  CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XLETR(:),isize,s%XWORK_WR)
ENDDO
!
yrecfm='DRAIN'//trim(hterm)
ycomment='X_Y_'//yrecfm//' (Kg/m2/s)'
DO jp=1,io%NPATCH
  CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XDRAIN(:),isize,s%XWORK_WR)
ENDDO
!
IF(io%CRUNOFF=='SGH'.AND.io%CISBA=='DIF')THEN
  yrecfm='QSB'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (Kg/m2/s)'
  DO jp=1,io%NPATCH
    CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
          np%AL(jp)%NR_P,ndea%AL(jp)%XQSB(:),isize,s%XWORK_WR)
  ENDDO  
ENDIF
!
yrecfm='RUNOFF'//trim(hterm)
ycomment='X_Y_'//yrecfm//' (Kg/m2/s)'
DO jp=1,io%NPATCH
  CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XRUNOFF(:),isize,s%XWORK_WR)
ENDDO
!
IF(io%CHORT=='SGH'.OR.io%CISBA=='DIF')THEN
 
  yrecfm='HORTON'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (Kg/m2/s)'
  DO jp=1,io%NPATCH
  CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XHORT(:),isize,s%XWORK_WR)
  ENDDO

ENDIF
!
yrecfm='DRIVEG'//trim(hterm)
ycomment='X_Y_'//yrecfm//' (Kg/m2/s)'
DO jp=1,io%NPATCH
  CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XDRIP(:),isize,s%XWORK_WR)
ENDDO
!
yrecfm='RRVEG'//trim(hterm)
ycomment='X_Y_'//yrecfm//' (Kg/m2/s)'
DO jp=1,io%NPATCH
  CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XRRVEG(:),isize,s%XWORK_WR)
ENDDO
!
yrecfm='SNOMLT'//trim(hterm)
ycomment='X_Y_'//yrecfm//' (Kg/m2/s)'
DO jp=1,io%NPATCH
  CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XMELT(:),isize,s%XWORK_WR)
ENDDO
!
! MEB STUFF
IF (isize_lmeb_patch>0) THEN
  !
  yrecfm='LELIT'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  DO jp=1,io%NPATCH
    CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XLELITTER(:),isize,s%XWORK_WR)
  ENDDO
  !
  yrecfm='LELITI'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  DO jp=1,io%NPATCH
    CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XLELITTERI(:),isize,s%XWORK_WR)
  ENDDO
  !
  yrecfm='DRIPLIT'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  DO jp=1,io%NPATCH
    CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XDRIPLIT(:),isize,s%XWORK_WR)
  ENDDO
  !
  yrecfm='RRLIT'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  DO jp=1,io%NPATCH
    CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XRRLIT(:),isize,s%XWORK_WR)
  ENDDO
  !
  yrecfm='LEV_CV'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  DO jp=1,io%NPATCH
    CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XLEV_CV(:),isize,s%XWORK_WR)
  ENDDO
  !
  yrecfm='LES_CV'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  DO jp=1,io%NPATCH
    CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XLES_CV(:),isize,s%XWORK_WR)
  ENDDO
  !
  yrecfm='LETR_CV'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  DO jp=1,io%NPATCH
    CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XLETR_CV(:),isize,s%XWORK_WR)
  ENDDO
  !
  yrecfm='LER_CV'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  DO jp=1,io%NPATCH
    CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XLER_CV(:),isize,s%XWORK_WR)
  ENDDO
  !
  yrecfm='LE_CV'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  DO jp=1,io%NPATCH
    CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XLE_CV(:),isize,s%XWORK_WR)
  ENDDO
  !
  yrecfm='H_CV'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  DO jp=1,io%NPATCH
    CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XH_CV(:),isize,s%XWORK_WR)
  ENDDO  
  !
  yrecfm='MELT_CV'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (kg/m2/s)'
  DO jp=1,io%NPATCH
    CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XMELT_CV(:),isize,s%XWORK_WR)
  ENDDO
  !
  yrecfm='FRZ_CV'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (kg/m2/s)'
  DO jp=1,io%NPATCH
    CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XFRZ_CV(:),isize,s%XWORK_WR)
  ENDDO  
  !
  yrecfm='LE_GV'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  DO jp=1,io%NPATCH
    CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XLE_GV(:),isize,s%XWORK_WR)
  ENDDO  
  !
  yrecfm='H_GV'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  DO jp=1,io%NPATCH
    CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XH_GV(:),isize,s%XWORK_WR)
  ENDDO  
  !
  yrecfm='LE_GN'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  DO jp=1,io%NPATCH
    CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XLE_GN(:),isize,s%XWORK_WR)
  ENDDO 
  !
  yrecfm='H_GN'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  DO jp=1,io%NPATCH
    CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XH_GN(:),isize,s%XWORK_WR)
  ENDDO
  !
  yrecfm='SR_GN'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (kg/m2/s)'
  DO jp=1,io%NPATCH
    CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XSR_GN(:),isize,s%XWORK_WR)
  ENDDO  
  !
  yrecfm='SWDN_GN'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  DO jp=1,io%NPATCH
    CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XSWDOWN_GN(:),isize,s%XWORK_WR)
  ENDDO
  !
  yrecfm='LWDN_GN'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  DO jp=1,io%NPATCH
    CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XLWDOWN_GN(:),isize,s%XWORK_WR)
  ENDDO  
  !
  yrecfm='LE_CA'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  DO jp=1,io%NPATCH
    CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XLE_CA(:),isize,s%XWORK_WR)
  ENDDO
  !
  yrecfm='H_CA'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  DO jp=1,io%NPATCH
    CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XH_CA(:),isize,s%XWORK_WR)
  ENDDO  
  !
  yrecfm='SWNT_V'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  DO jp=1,io%NPATCH
    CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XSWNET_V(:),isize,s%XWORK_WR)
  ENDDO
  !
  yrecfm='SWNT_G'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  DO jp=1,io%NPATCH
    CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XSWNET_G(:),isize,s%XWORK_WR)
  ENDDO
  !
  yrecfm='SWNT_N'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  DO jp=1,io%NPATCH
    CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XSWNET_N(:),isize,s%XWORK_WR)
  ENDDO
  !
  yrecfm='SWNT_NS'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  DO jp=1,io%NPATCH
    CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XSWNET_NS(:),isize,s%XWORK_WR)
  ENDDO
  !
  yrecfm='LWNT_V'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  DO jp=1,io%NPATCH
    CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XLWNET_V(:),isize,s%XWORK_WR)
  ENDDO
  !
  yrecfm='LWNT_G'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  DO jp=1,io%NPATCH
    CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XLWNET_G(:),isize,s%XWORK_WR)
  ENDDO
  !
  yrecfm='LWNT_N'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  DO jp=1,io%NPATCH
    CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XLWNET_N(:),isize,s%XWORK_WR)
  ENDDO

ENDIF
! END MEB STUFF
!
IF(lagrip)THEN
  yrecfm='IRRIG'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (Kg/m2/s)'
  DO jp=1,io%NPATCH
    CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XIRRIG_FLUX(:),isize,s%XWORK_WR)
  ENDDO
ENDIF
!
IF(io%LFLOOD)THEN
  !
  yrecfm='IFLOOD'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (Kg/m2/s)'
  DO jp=1,io%NPATCH
    CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XIFLOOD(:),isize,s%XWORK_WR)
  ENDDO
  !
  yrecfm='PFLOOD'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (Kg/m2/s)'
  DO jp=1,io%NPATCH
    CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XPFLOOD(:),isize,s%XWORK_WR)
  ENDDO
  !

  yrecfm='LEF'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  DO jp=1,io%NPATCH
  CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XLE_FLOOD(:),isize,s%XWORK_WR)
ENDDO
  !

  yrecfm='LEIF'//trim(hterm)
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  DO jp=1,io%NPATCH
    CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XLEI_FLOOD(:),isize,s%XWORK_WR)
  ENDDO
  !
ENDIF
!
IF(io%CPHOTO/='NON')THEN
  !
 
  yrecfm='GPP'//trim(hterm)
  ycomment='gross primary production per patch (kgCO2/m2/s)'
  DO jp=1,io%NPATCH
    CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XGPP(:),isize,s%XWORK_WR)
  ENDDO
  !
  IF (hterm(1:1)=="C") THEN
    yrecfm='RC_AUTO'//trim(hterm(2:))
  ELSE
    yrecfm='R_AUTO'//trim(hterm)
  ENDIF
  ycomment='autotrophic respiration per patch (kgCO2/m2/s)'
    DO jp=1,io%NPATCH
    CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XRESP_AUTO(:),isize,s%XWORK_WR)
  ENDDO
  !
  IF (hterm(1:1)=="C") THEN
    yrecfm='RC_ECO'//trim(hterm(2:))
  ELSE
    yrecfm='R_ECO'//trim(hterm)
  ENDIF  
  ycomment='ecosystem respiration per patch (kgCO2/m2/s)'
  DO jp=1,io%NPATCH
    CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XRESP_ECO(:),isize,s%XWORK_WR)
  ENDDO
  !
ENDIF
!
IF(id%DE%LWATER_BUDGET .OR. oflag)THEN 
  !
  yrecfm='DWG'//trim(hterm)
  ycomment='change in liquid soil moisture per patch (Kg/m2/s)'
    DO jp=1,io%NPATCH
    CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XDWG(:),isize,s%XWORK_WR)
  ENDDO
  !
  yrecfm='DWGI'//trim(hterm)
  ycomment='change in solid soil moisture per patch (Kg/m2/s)'
    DO jp=1,io%NPATCH
    CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XDWGI(:),isize,s%XWORK_WR)
  ENDDO
  !
  yrecfm='DWR'//trim(hterm)
  ycomment='change in water on canopy per patch (Kg/m2/s)'
  DO jp=1,io%NPATCH
    CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XDWR(:),isize,s%XWORK_WR)
  ENDDO
  !
  yrecfm='DSWE'//trim(hterm)
  ycomment='change in snow water equivalent per patch (Kg/m2/s)'
  DO jp=1,io%NPATCH
    CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XDSWE(:),isize,s%XWORK_WR)
  ENDDO
  !
  yrecfm='WATBUD'//trim(hterm)
  ycomment='isba water budget as residue per patch (Kg/m2/s)'
  DO jp=1,io%NPATCH
    CALL write_field_1d_patch(duo%CSELECT,hprogram,yrecfm,ycomment,jp,&
        np%AL(jp)%NR_P,ndea%AL(jp)%XWATBUD(:),isize,s%XWORK_WR)
  ENDDO
  !
ENDIF
!
IF (lhook) CALL dr_hook('WRITE_DIAG_SEB_ISBA_N:WRITE_EVAP_BUD_PATCH',1,zhook_handle)
!
END SUBROUTINE write_evap_bud_patch
!
!
SUBROUTINE provar_to_diag
!
REAL, DIMENSION(U%NSIZE_NATURE)             :: ZPATCH, ZWORK
REAL, DIMENSION(U%NSIZE_NATURE,SIZE(NPE%AL(1)%XWG,2)) :: ZWG
REAL, DIMENSION(U%NSIZE_NATURE,SIZE(NPE%AL(1)%XWG,2)) :: ZWGI
REAL, DIMENSION(U%NSIZE_NATURE,SIZE(NPE%AL(1)%XTG,2)) :: ZTG
REAL, DIMENSION(U%NSIZE_NATURE,SIZE(NP%AL(1)%XDG,2)) :: ZDG_TOT
REAL, DIMENSION(U%NSIZE_NATURE,SIZE(NP%AL(1)%XDG,2)) :: ZDG
!
REAL, DIMENSION(U%NSIZE_NATURE,IO%NNBIOMASS)   :: ZBIOMASS
REAL, DIMENSION(U%NSIZE_NATURE,IO%NNSOILCARB)  :: ZSOILCARB
REAL, DIMENSION(U%NSIZE_NATURE,IO%NNLITTLEVS)  :: ZLIGNIN_STRUC
REAL, DIMENSION(U%NSIZE_NATURE,IO%NNLITTER,IO%NNLITTLEVS)  :: ZLITTER
!
 CHARACTER(LEN=4 ) :: YLVL
REAL              :: ZMISS
INTEGER           :: JL, JP, JJ, INI, IWORK, IDEPTH, IMASK
!
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
IF (lhook) CALL dr_hook('WRITE_DIAG_SEB_ISBA_N:PROVAR_TO_DIAG',0,zhook_handle)
!
ini=u%NSIZE_NATURE
!
! * soil temperatures (K)
!
IF(io%LTEMP_ARP)THEN
  iwork=io%NTEMPLAYER_ARP
ELSEIF(io%CISBA/='DIF')THEN
  iwork=2
ELSE
  iwork=io%NGROUND_LAYER
ENDIF
!
ztg(:,:)=0.0
DO jp=1,io%NPATCH
  pk => np%AL(jp)
  pek => npe%AL(jp)
  DO jl=1,iwork
    DO jj=1,pk%NSIZE_P
      imask = pk%NR_P(jj)
      ztg(imask,jl) = ztg(imask,jl) + pk%XPATCH(jj) * pek%XTG(jj,jl)
      ENDDO
   ENDDO
ENDDO
!
IF (lallow_add_dim)  THEN 
  yrecfm='TG_ISBA' ; 
  ycomment='Soil temperature (K)'
  CALL write_surf(duo%CSELECT, &
                  hprogram,yrecfm,ztg(:,:),iresp,ycomment,hnam_dim=yground_layer_dim_name)
ELSE
  DO jl=1,iwork
    WRITE(ylvl,'(I4)') jl
    yrecfm='TG'//adjustl(ylvl(:len_trim(ylvl)))
    yrecfm=yrecfm(:len_trim(yrecfm))//'_ISBA'
    ycomment='X_Y_'//yrecfm//' (K)'
    CALL write_surf(duo%CSELECT,hprogram,yrecfm,ztg(:,jl),iresp,hcomment=ycomment)
  END DO
ENDIF
!
! * Compute soil liquid and ice water content (kg/m2 and m3/m3) 
!
zwg(:,:)=0.0
zwgi(:,:)=0.0
zdg_tot(:,:)=0.0
!
IF(io%CISBA=='DIF')THEN
  !
  DO jp=1,io%NPATCH
    pk => np%AL(jp)
    pek => npe%AL(jp)

    DO jj=1,pk%NSIZE_P
      imask = pk%NR_P(jj)
      DO jl=1,io%NGROUND_LAYER
!
!       liquid and ice water content
        idepth= pk%NWG_LAYER(jj)
        IF(jl<=idepth)THEN    
          zwg(imask,jl) = zwg(imask,jl)+ pk%XPATCH(jj) * pek%XWG (jj,jl) * pk%XDZG(jj,jl)
          zwgi(imask,jl)= zwgi(imask,jl)+ pk%XPATCH(jj) * pek%XWGI(jj,jl) * pk%XDZG(jj,jl)
          zdg_tot(imask,jl)= zdg_tot(imask,jl)+ pk%XPATCH(jj)* pk%XDZG(jj,jl)
        ENDIF
!                      
      ENDDO
    ENDDO
  ENDDO
!  
ELSE
  !
  DO jp=1,io%NPATCH
    pk => np%AL(jp)
    pek => npe%AL(jp)

    zdg(1:pk%NSIZE_P,1) = pk%XDG(:,1)
    zdg(1:pk%NSIZE_P,2) = pk%XDG(:,2)
    IF(io%CISBA=='3-L')THEN
      zdg(1:pk%NSIZE_P,3) = pk%XDG(:,3)-pk%XDG(:,2)
    ENDIF

    DO jj=1,pk%NSIZE_P
      imask = pk%NR_P(jj)
      DO jl=1,io%NGROUND_LAYER
        
        zwg(imask,jl) = zwg(imask,jl)+ pk%XPATCH(jj) *pek%XWG(jj,jl)*  zdg(jj,jl)
        zwgi(imask,jl)= zwgi(imask,jl)+ pk%XPATCH(jj) *pek%XWGI(jj,jl)* zdg(jj,jl)
        zdg_tot(imask,jl)=zdg_tot(imask,jl)+pk%XPATCH(jj)*zdg(jj,jl)
      ENDDO
    ENDDO
  ENDDO
!  
ENDIF
!
IF(hprogram=='AROME '.OR.hprogram=='FA    ')THEN
  zmiss=0.0
ELSE
  zmiss=xundef
ENDIF
!
WHERE(zdg_tot(:,:)>0.0)
      zwg(:,:)=zwg(:,:)/zdg_tot(:,:)
      zwgi(:,:)=zwgi(:,:)/zdg_tot(:,:)
ELSEWHERE
      zwg(:,:)=zmiss
      zwgi(:,:)=zmiss    
ENDWHERE
!
! * soil liquid water content (m3/m3) and soil moisture (kg/m2)
!
IF (lallow_add_dim)  THEN 
  yrecfm='WG_ISBA' ; 
  ycomment='Soil liquid water content (m3/m3)' 
  CALL write_surf(duo%CSELECT, &
                  hprogram,yrecfm,zwg(:,:),iresp,ycomment,hnam_dim=ywground_layer_dim_name)
ELSE
  DO jl=1,io%NGROUND_LAYER
    WRITE(ylvl,'(I4)') jl
    yrecfm='WG'//adjustl(ylvl(:len_trim(ylvl)))
    yrecfm=yrecfm(:len_trim(yrecfm))//'_ISBA'
    ycomment='Soil liquid water content (m3/m3)' 
    CALL write_surf(duo%CSELECT,hprogram,yrecfm,zwg(:,jl),iresp,hcomment=ycomment)
  END DO
ENDIF
!
! * soil ice water content (m3/m3) and soil ice mass (kg/m2)
!
iwork=io%NGROUND_LAYER
IF(io%CISBA/='DIF')THEN
  iwork=2 ! No ice in the FR 3-layers
ENDIF
!
IF (lallow_add_dim)  THEN 
  yrecfm='WGI_ISBA' ; 
  ycomment='Soil solid water content (m3/m3)' 
  CALL write_surf(duo%CSELECT, &
                  hprogram,yrecfm,zwgi(:,:),iresp,ycomment,hnam_dim=ywiground_layer_dim_name)
ELSE
  DO jl=1,iwork
    WRITE(ylvl,'(I4)') jl
    yrecfm='WGI'//adjustl(ylvl(:len_trim(ylvl)))
    yrecfm=yrecfm(:len_trim(yrecfm))//'_ISBA'
    ycomment='Soil solid water content (m3/m3)' 
    CALL write_surf(duo%CSELECT,hprogram,yrecfm,zwgi(:,jl),iresp,hcomment=ycomment) 
  END DO
ENDIF
!
! * water intercepted on leaves (kg/m2)
!
zwork(:)=0.0
DO jp=1,io%NPATCH
  DO jj=1,np%AL(jp)%NSIZE_P
    imask = np%AL(jp)%NR_P(jj)
    zwork(imask) = zwork(imask) + np%AL(jp)%XPATCH(jj) * npe%AL(jp)%XWR(jj)
  ENDDO
ENDDO
!
yrecfm='WR_ISBA'
ycomment='X_Y_'//yrecfm//' (kg/m2)'
 CALL write_surf(duo%CSELECT,hprogram,yrecfm,zwork(:),iresp,hcomment=ycomment)
!
! * Glacier ice storage (semi-prognostic) (kg/m2)
!
IF(io%LGLACIER)THEN
  !
  zwork(:)=0.0
  DO jp=1,io%NPATCH
    DO jj=1,np%AL(jp)%NSIZE_P
      imask = np%AL(jp)%NR_P(jj)
      zwork(imask) = zwork(imask) + np%AL(jp)%XPATCH(jj) * npe%AL(jp)%XICE_STO(jj)
    ENDDO    
  ENDDO    
  !
  yrecfm='ICE_STO_ISBA'
  ycomment='X_Y_'//yrecfm//' (kg/m2)'
  CALL write_surf(duo%CSELECT,hprogram,yrecfm,zwork(:),iresp,hcomment=ycomment)
  !
ENDIF
!
! * Snow albedo (-) 
!
zpatch(:) = 0.0
zwork(:) = 0.0
DO jp=1,io%NPATCH
  pk => np%AL(jp)
  pek => npe%AL(jp)

  DO jj=1,pk%NSIZE_P
    imask = pk%NR_P(jj)

    IF(pek%TSNOW%ALB(jj)/=xundef)THEN
      zwork(imask) = zwork(imask) + pk%XPATCH(jj) * pek%TSNOW%ALB(jj)
      zpatch(imask) = zpatch(imask) + pk%XPATCH(jj)
    ENDIF
  ENDDO
ENDDO
!
WHERE(zpatch(:)>0.0)
  zwork(:) = zwork(:) / zpatch(:)
ELSEWHERE
  zwork(:) = xundef
ENDWHERE
!
yrecfm='ASN_ISBA'
ycomment='X_Y_'//yrecfm//' (-)'
 CALL write_surf(duo%CSELECT,hprogram,yrecfm,zwork(:),iresp,hcomment=ycomment)
!  
IF(npe%AL(1)%TSNOW%SCHEME=='3-L' .OR. npe%AL(1)%TSNOW%SCHEME=='CRO')THEN
  !
  ! * Snow reservoir (kg/m2) by layer
  !
  DO jl = 1,npe%AL(1)%TSNOW%NLAYER
    !
    zwork(:)=0.0
    DO jp=1,io%NPATCH
      pk => np%AL(jp)
      pek => npe%AL(jp)    
      DO jj=1,pk%NSIZE_P
        imask = pk%NR_P(jj)
        zwork(imask) = zwork(imask) + pk%XPATCH(jj) * pek%TSNOW%WSNOW(jj,jl)
      ENDDO
    ENDDO
    !
    WRITE(ylvl,'(I4)') jl
    yrecfm='WSN_'//adjustl(ylvl(:len_trim(ylvl)))
    yrecfm=yrecfm(:len_trim(yrecfm))//'_ISBA'
    ycomment='X_Y_'//yrecfm//' (kg/m2)'
    CALL write_surf(duo%CSELECT,hprogram,yrecfm,zwork(:),iresp,hcomment=ycomment)
    !
  ENDDO
  !
  ! * Snow depth (m)
  !
  DO jl = 1,npe%AL(1)%TSNOW%NLAYER
    !
    zwork(:)=0.0
    DO jp=1,io%NPATCH
      pk => np%AL(jp)
      pek => npe%AL(jp)
      DO jj=1,pk%NSIZE_P
        imask = pk%NR_P(jj)
        zwork(imask) = zwork(imask) + pk%XPATCH(jj) * &
                pek%TSNOW%WSNOW(jj,jl)/pek%TSNOW%RHO(jj,jl)
      ENDDO
    ENDDO
    !
    WRITE(ylvl,'(I4)') jl
    yrecfm='DSN_'//adjustl(ylvl(:len_trim(ylvl)))
    yrecfm=yrecfm(:len_trim(yrecfm))//'_ISBA'
    ycomment='X_Y_'//yrecfm//' (kg/m2)'
    CALL write_surf(duo%CSELECT,hprogram,yrecfm,zwork(:),iresp,hcomment=ycomment)
    !
  ENDDO
  !
  ! * Snow temperature (k)
  !  
  IF(hprogram=='AROME '.OR.hprogram=='FA    ')THEN
    zmiss=xtt
  ELSE
    zmiss=xundef
  ENDIF
  !  
  DO jl = 1,npe%AL(1)%TSNOW%NLAYER
    !
    zwork(:) = 0.0
    zpatch(:) = 0.0
    DO jp=1,io%NPATCH
      pk => np%AL(jp)
      pek => npe%AL(jp)    
      DO jj=1,pk%NSIZE_P
        imask = pk%NR_P(jj)
        IF(pek%TSNOW%WSNOW(jj,jl)>0.)THEN
          zwork(imask) = zwork(imask) + pk%XPATCH(jj) * pek%TSNOW%TEMP(jj,jl) 
          zpatch(imask) = zpatch(imask) + pk%XPATCH(jj)
        ENDIF
      ENDDO
    ENDDO
    !
    WHERE(zpatch(:)>0.0)
      zwork(:) = zwork(:) / zpatch(:)
    ELSEWHERE
      zwork(:) = zmiss
    ENDWHERE
    !
    WRITE(ylvl,'(I4)') jl
    yrecfm='TSN_'//adjustl(ylvl(:len_trim(ylvl)))
    yrecfm=yrecfm(:len_trim(yrecfm))//'_ISBA'
    ycomment='X_Y_'//yrecfm//' (K)'
    CALL write_surf(duo%CSELECT,hprogram,yrecfm,zwork(:),iresp,hcomment=ycomment)
    !
  ENDDO
  !
  ! * Snow age (day)
  !    
  DO jl = 1,npe%AL(1)%TSNOW%NLAYER
    !
    zwork(:) = 0.0
    zpatch(:) = 0.0
    DO jp=1,io%NPATCH
      pk => np%AL(jp)
      pek => npe%AL(jp)    
      DO jj=1,pk%NSIZE_P
        imask = pk%NR_P(jj)    
        IF(pek%TSNOW%WSNOW(jj,jl)>0.)THEN    
          zwork(imask) = zwork(imask) + pk%XPATCH(jj) * pek%TSNOW%AGE(jj,jl) 
          zpatch(imask) = zpatch(imask) + pk%XPATCH(jj)
          ENDIF
       ENDDO
    ENDDO
    !
    WHERE(zpatch(:)>0.0)
      zwork(:) = zwork(:) / zpatch(:)
    ENDWHERE
    !
    WRITE(ylvl,'(I4)') jl
    yrecfm='AGSN_'//adjustl(ylvl(:len_trim(ylvl)))
    yrecfm=yrecfm(:len_trim(yrecfm))//'_ISBA'
    ycomment='X_Y_'//yrecfm//' (day_since_snowfall)'
    CALL write_surf(duo%CSELECT,hprogram,yrecfm,zwork(:),iresp,hcomment=ycomment)
    !
  ENDDO
  !
ENDIF
!
! * Isba-Ags biomass reservoir
!
IF(io%CPHOTO=='NIT'.OR.io%CPHOTO=='NCB')THEN
!
  zbiomass(:,:)=0.0
  DO jp=1,io%NPATCH
    pk => np%AL(jp)
    pek => npe%AL(jp)    
    DO jj=1,pk%NSIZE_P
      imask = pk%NR_P(jj)      
      DO jl=1,io%NNBIOMASS
        zbiomass(imask,jl) = zbiomass(imask,jl) + pk%XPATCH(jj) * pek%XBIOMASS(jj,jl)
      ENDDO
    ENDDO
  ENDDO
!
  DO jl = 1,io%NNBIOMASS
    WRITE(ylvl,'(I4)') jl
    yrecfm='BIOM'//adjustl(ylvl(:len_trim(ylvl)))
    yrecfm=yrecfm(:len_trim(yrecfm))//'_ISBA'
    ycomment='X_Y_'//yrecfm//' (kgDM/m2)'
    CALL write_surf(duo%CSELECT,hprogram,yrecfm,zbiomass(:,jl),iresp,hcomment=ycomment)  
  ENDDO
!
ENDIF
!
! * Isba-CC carbon reservoir
!
IF(io%CRESPSL=='CNT')THEN
!
  zlitter(:,:,:)=0.0
  zlignin_struc(:,:)=0.0
  DO jp=1,io%NPATCH
    pk => np%AL(jp)
    pek => npe%AL(jp)    
    DO jj=1,pk%NSIZE_P
      imask = pk%NR_P(jj) 
      DO jl=1,io%NNLITTLEVS
        zlitter(imask,1,jl) = zlitter(imask,1,jl) + pk%XPATCH(jj) * pek%XLITTER(jj,1,jl)
        zlitter(imask,2,jl) = zlitter(imask,2,jl) + pk%XPATCH(jj) * pek%XLITTER(jj,2,jl)
        zlignin_struc(imask,jl) = zlignin_struc(imask,jl) + pk%XPATCH(jj) * pek%XLIGNIN_STRUC(jj,jl)
      ENDDO
    ENDDO
  ENDDO
!       
  DO jl=1,io%NNLITTLEVS
     WRITE(ylvl,'(I4)') jl
     yrecfm='LIT1_'//adjustl(ylvl(:len_trim(ylvl)))
     yrecfm=yrecfm(:len_trim(yrecfm))//'_ISBA'
     ycomment='X_Y_'//yrecfm//' (gC/m2)'
     CALL write_surf(duo%CSELECT,hprogram,yrecfm,zlitter(:,1,jl),iresp,hcomment=ycomment)  
     WRITE(ylvl,'(I4)') jl
     yrecfm='LIT2_'//adjustl(ylvl(:len_trim(ylvl)))
     yrecfm=yrecfm(:len_trim(yrecfm))//'_ISBA'
     ycomment='X_Y_'//yrecfm//' (gC/m2)'
     CALL write_surf(duo%CSELECT, hprogram,yrecfm,zlitter(:,2,jl),iresp,hcomment=ycomment)
     WRITE(ylvl,'(I4)') jl
     yrecfm='LIGSTR'//adjustl(ylvl(:len_trim(ylvl)))
     yrecfm=yrecfm(:len_trim(yrecfm))//'_ISBA'
     ycomment='X_Y_'//yrecfm//' (-)'
     CALL write_surf(duo%CSELECT,hprogram,yrecfm,zlignin_struc(:,jl),iresp,hcomment=ycomment)      
  END DO
!
  zsoilcarb(:,:)=0.0
  DO jp=1,io%NPATCH
    pk => np%AL(jp)
    pek => npe%AL(jp)    
    DO jj=1,pk%NSIZE_P
      imask = pk%NR_P(jj)   
      DO jl=1,io%NNSOILCARB
        zsoilcarb(imask,jl) = zsoilcarb(imask,jl) + pk%XPATCH(jj) * pek%XSOILCARB(jj,jl)
      ENDDO
    ENDDO
  ENDDO
!
  DO jl = 1,io%NNSOILCARB
    WRITE(ylvl,'(I4)') jl
    yrecfm='SCARB'//adjustl(ylvl(:len_trim(ylvl)))
    yrecfm=yrecfm(:len_trim(yrecfm))//'_ISBA'
    ycomment='X_Y_'//yrecfm//' (gC/m2)'
    CALL write_surf(duo%CSELECT, hprogram,yrecfm,zsoilcarb(:,jl),iresp,hcomment=ycomment)  
  ENDDO
!
ENDIF
!
IF (lhook) CALL dr_hook('WRITE_DIAG_SEB_ISBA_N:PROVAR_TO_DIAG',1,zhook_handle)
!
END SUBROUTINE provar_to_diag
!
END SUBROUTINE write_diag_seb_isba_n