write_diag_misc_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_misc_isba_n (DTCO, DGU, U, DGI, DGMI, I, &
                                         hprogram)
!     #################################
!
!!****  *WRITE_DIAG_MISC_ISBA* - writes the ISBA diagnostic fields
!!
!!    PURPOSE
!!    -------
!!
!!
!!**  METHOD
!!    ------
!!
!!    REFERENCE
!!    ---------
!!
!!
!!    AUTHOR
!!    ------
!!      P. Le Moigne   *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    10/2004
!!      B. Decharme    2008  Total Albedo, Total SWI and Floodplains
!!      B. Decharme 06/2009  key to write (or not) patch result
!!      A.L. Gibelin 04/09 : Add respiration diagnostics
!!      A.L. Gibelin 05/09 : Add carbon spinup
!!      A.L. Gibelin 07/09 : Suppress RDK and transform GPP as a diagnostic
!!      D. Carrer    04/11 : Add FAPAR and effective LAI
!!      B. Decharme  09/2012 : suppress NWG_LAYER (parallelization problems)
!!      B. Decharme  09/12 : Carbon fluxes in diag_evap
!!      B. Decharme  09/12   New diag for DIF:
!!                           F2 stress
!!                           Root zone swi, wg and wgi
!!                           swi, wg and wgi comparable to ISBA-FR-DG2 and DG3 layers
!!                           active layer thickness over permafrost
!!                           frozen layer thickness over non-permafrost
!!      B. Decharme  06/13   All snow outputs noted SN
!!                           XTSRAD_NAT instead of XAVG_TSRAD
!!                           delete NWG_SIZE
!!                           water table depth
!!
!-------------------------------------------------------------------------------
!
!*       0.    DECLARATIONS
!              ------------
!
USE modd_data_cover_n, ONLY : data_cover_t
USE modd_diag_surf_atm_n, ONLY : diag_surf_atm_t
USE modd_surf_atm_n, ONLY : surf_atm_t
USE modd_diag_isba_n, ONLY : diag_isba_t
USE modd_diag_misc_isba_n, ONLY : diag_misc_isba_t
USE modd_isba_n, ONLY : isba_t
!
USE modd_surf_par,        ONLY :   nundef, xundef
!
USE modd_assim, ONLY : lassim, cassim_isba, nvar
!                                 
USE modd_agri,            ONLY :   lagrip
!
USE modi_init_io_surf_n
USE modi_write_surf
USE modi_end_io_surf_n
!
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_surf_atm_t), INTENT(INOUT) :: dgu
TYPE(surf_atm_t), INTENT(INOUT) :: u
TYPE(diag_isba_t), INTENT(INOUT) :: dgi
TYPE(diag_misc_isba_t), INTENT(INOUT) :: dgmi
TYPE(isba_t), INTENT(INOUT) :: i
!
 CHARACTER(LEN=6),  INTENT(IN)  :: hprogram ! program calling
!
!*       0.2   Declarations of local variables
!              -------------------------------
!
INTEGER           :: iresp          ! IRESP  : return-code if a problem appears
 CHARACTER(LEN=1) :: yvar
 CHARACTER(LEN=12) :: yrecfm         ! Name of the article to be read
 CHARACTER(LEN=100):: ycomment       ! Comment string
 CHARACTER(LEN=2)  :: ylvl
 CHARACTER(LEN=20) :: yform
!
INTEGER           :: jlayer, jj, idepth, jvar
!
REAL(KIND=JPRB) :: zhook_handle
!
!-------------------------------------------------------------------------------
!
!         Initialisation for IO
!
IF (lhook) CALL dr_hook('WRITE_DIAG_MISC_ISBA_N',0,zhook_handle)
 CALL init_io_surf_n(dtco, dgu, u, &
                     hprogram,'NATURE','ISBA  ','WRITE')
!
!-------------------------------------------------------------------------------
!
IF (dgmi%LSURF_MISC_BUDGET) THEN
  !
  !*       2.     Miscellaneous fields :
  !
  !-------------------------------------------------------------------------------
  !
  !        2.1    Halstead coefficient
  !               --------------------
  !
  yrecfm='HV_ISBA'
  ycomment='Halstead coefficient averaged over tile nature (-)'
  CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XAVG_HV(:),iresp,hcomment=ycomment)
  !
  !        2.2    Snow fractions
  !               --------------
  !
  yrecfm='PSNG_ISBA'
  ycomment='snow fraction over ground averaged over tile nature (-)'
  CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XAVG_PSNG(:),iresp,hcomment=ycomment)
  !
  yrecfm='PSNV_ISBA'
  ycomment='snow fraction over vegetation averaged over tile nature (-)'
  CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XAVG_PSNV(:),iresp,hcomment=ycomment)
  !
  yrecfm='PSN_ISBA'
  ycomment='total snow fraction averaged over tile nature (-)'
  CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XAVG_PSN(:),iresp,hcomment=ycomment)
  !
  !        2.3    Total Albedo and surface temperature
  !               ------------------------------------
  !
  yrecfm='TALB_ISBA'
  ycomment='total albedo over tile nature (-)'
  CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XAVG_ALBT(:),iresp,hcomment=ycomment)
  !
  IF (i%TSNOW%SCHEME=='3-L' .OR. i%TSNOW%SCHEME=='CRO') THEN
    !        
    yrecfm='TS_ISBA'
    ycomment='total surface temperature (isba+snow) over tile nature'
    CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgi%XAVG_TS(:),iresp,hcomment=ycomment)
    !
    yrecfm='TSRAD_ISBA'
    ycomment='total radiative surface temperature (isba+snow) over tile nature'
    CALL write_surf(dgu, u, &
                  hprogram,yrecfm,i%XTSRAD_NAT(:),iresp,hcomment=ycomment)
    !
  END IF
  !
  !        2.4    Soil Wetness Index, Water content and active layer depth
  !               --------------------------------------------------------
  !  
  IF(i%CISBA=='DIF')THEN
    DO jlayer = 1,i%NGROUND_LAYER
     DO jj=1,SIZE(i%NWG_LAYER,1)
        idepth=maxval(i%NWG_LAYER(jj,:),i%NWG_LAYER(jj,:)/=nundef)
        IF(jlayer>idepth)THEN  
          dgmi%XAVG_SWI (jj,jlayer) = xundef
          dgmi%XAVG_TSWI(jj,jlayer) = xundef
        ENDIF
      ENDDO 
    ENDDO
  ENDIF         
  !
  DO jlayer=1,i%NGROUND_LAYER
    !
    WRITE(ylvl,'(I2)') jlayer
    !
    yrecfm='SWI'//adjustl(ylvl(:len_trim(ylvl)))
    yrecfm=yrecfm(:len_trim(yrecfm))//'_ISBA'
    yform='(A29,I1.1,A4)'
    IF (jlayer >= 10)  yform='(A29,I2.2,A4)'
    WRITE(ycomment,fmt=yform) 'soil wetness index for layer ',jlayer,' (-)'
    CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XAVG_SWI(:,jlayer),iresp,hcomment=ycomment)
    !
    yrecfm='TSWI'//adjustl(ylvl(:len_trim(ylvl)))
    yrecfm=yrecfm(:len_trim(yrecfm))//'_ISBA'
    yform='(A29,I1.1,A4)'
    IF (jlayer >= 10)  yform='(A29,I2.2,A4)'
    WRITE(ycomment,fmt=yform) 'total swi (liquid+solid) for layer ',jlayer,' (-)'
    CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XAVG_TSWI(:,jlayer),iresp,hcomment=ycomment)
    !
  END DO
  !
  yrecfm='SWI_T_ISBA'
  ycomment='soil wetness index over the soil column (-)'
  CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XSOIL_SWI(:),iresp,hcomment=ycomment)
  !
  yrecfm='TSWI_T_ISBA'
  ycomment='total soil wetness index over the soil column (-)'
  CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XSOIL_TSWI(:),iresp,hcomment=ycomment)
  !
  yrecfm='WGTOT_T_ISBA'
  ycomment='total water content (liquid+solid) over the soil column (kg/m2)'
  CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XSOIL_TWG(:),iresp,hcomment=ycomment)
  !
  yrecfm='WGI_T_ISBA'
  ycomment='total ice content (solid) over the soil column (kg/m2)'
  CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XSOIL_TWGI(:),iresp,hcomment=ycomment)
  !
  yrecfm='WGTOT_ISBA'
  ycomment='total volumetric water content (liquid+solid) over the soil column (m3/m3)'
  CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XSOIL_WG(:),iresp,hcomment=ycomment)
  !
  yrecfm='WGI_ISBA'
  ycomment='total volumetric ice content (solid) over the soil column (m3/m3)'
  CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XSOIL_WGI(:),iresp,hcomment=ycomment)
  !
  IF(i%CISBA=='DIF') THEN
    !
    IF (dgmi%LSURF_MISC_DIF)THEN
      !
      yrecfm='TSWI_D2_ISBA'
      ycomment='total soil wetness index over comparable FR-DG2 reservoir (-)'
      CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XFRD2_TSWI(:),iresp,hcomment=ycomment)
      !
      yrecfm='WG_D2_ISBA'
      ycomment='liquid water content over comparable FR-DG2 reservoir (m3/m3)'
      CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XFRD2_TWG(:),iresp,hcomment=ycomment)
      !
      yrecfm='WGI_D2_ISBA'
      ycomment='ice content over comparable FR-DG2 reservoir (m3/m3)'
      CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XFRD2_TWGI(:),iresp,hcomment=ycomment)  
      !
      yrecfm='TSWI_D3_ISBA'
      ycomment='total soil wetness index over comparable FR-DG3 reservoir (-)'
      CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XFRD3_TSWI(:),iresp,hcomment=ycomment)
      !
      yrecfm='WG_D3_ISBA'
      ycomment='liquid water content over comparable FR-DG3 reservoir (m3/m3)'
      CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XFRD3_TWG(:),iresp,hcomment=ycomment)
      !
      yrecfm='WGI_D3_ISBA'
      ycomment='ice content over comparable FR-DG3 reservoir (m3/m3)'
      CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XFRD3_TWGI(:),iresp,hcomment=ycomment)  
      !
    ENDIF
    !
    yrecfm='ALT_ISBA'
    ycomment='active layer thickness over permafrost (m)'
    CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XAVG_ALT(:),iresp,hcomment=ycomment)
    !
    yrecfm='FLT_ISBA'
    ycomment='frozen layer thickness over non-permafrost (m)'
    CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XAVG_FLT(:),iresp,hcomment=ycomment)
    !
  ENDIF
  !
  !        2.5    Snow outputs
  !               -------------
  !
  yrecfm='WSN_T_ISBA'
  ycomment='Total_snow_reservoir (kg/m2)'
  CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XAVG_TWSNOW(:),iresp,hcomment=ycomment)
  !
  yrecfm='DSN_T_ISBA'
  ycomment='Total_snow_depth (m)'
  CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XAVG_TDSNOW(:),iresp,hcomment=ycomment)
  !
  yrecfm='TSN_T_ISBA'
  ycomment='Total_snow_temperature (K)'
  CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XAVG_TTSNOW(:),iresp,hcomment=ycomment)
  !
  !        2.6    SGH scheme
  !               ----------
  !
  IF(i%CRUNOFF=='SGH '.OR.i%CRUNOFF=='DT92')THEN     
    yrecfm='FSAT_ISBA'
    ycomment='Soil saturated fraction (-)'
    CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XAVG_FSAT(:),iresp,hcomment=ycomment)
  ENDIF
  !
  IF(i%CRAIN=='SGH ')THEN
    yrecfm='MUF_ISBA'
    ycomment='fraction of the grid cell reached by the rainfall (-)'
    CALL write_surf(dgu, u, &
                  hprogram,yrecfm,i%XMUF(:),iresp,hcomment=ycomment)
  ENDIF
  !
  !        2.7    Flooding scheme
  !               ---------------
  !
  IF(i%LFLOOD)THEN
    !
    yrecfm='FFG_ISBA'
    ycomment='flood fraction over ground averaged over tile nature (-)'
    CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XAVG_FFG(:),iresp,hcomment=ycomment)
    !
    yrecfm='FFV_ISBA'
    ycomment='flood fraction over vegetation averaged over tile nature (-)'
    CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XAVG_FFV(:),iresp,hcomment=ycomment)
    !
    yrecfm='FF_ISBA'
    ycomment='total flood fraction averaged over tile nature (-)'
    CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XAVG_FF(:),iresp,hcomment=ycomment)
    !
    yrecfm='FFLOOD_ISBA'
    ycomment='Grdi-cell potential flood fraction (-)'
    CALL write_surf(dgu, u, &
                  hprogram,yrecfm,i%XFFLOOD(:),iresp,hcomment=ycomment)
    !
    yrecfm='PIFLOOD_ISBA'
    ycomment='Grdi-cell Potential_floodplain_infiltration (kg/m2/s)'
    CALL write_surf(dgu, u, &
                  hprogram,yrecfm,i%XPIFLOOD(:),iresp,hcomment=ycomment)
    !
  ENDIF
  !
  !        2.8    Total LAI
  !               ---------
  !
  IF(i%CPHOTO/='NON'.OR.i%NPATCH>1)THEN        
    yrecfm='LAI_ISBA'
    ycomment='leaf area index (m2/m2)'
    CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XAVG_LAI(:),iresp,hcomment=ycomment)
  ENDIF
  !
  !        2.9    Water table depth
  !               -----------------
  !
  IF(i%LWTD)THEN
    !
    yrecfm='FWTD_ISBA'
    ycomment='grid-cell fraction of water table to rise'
    CALL write_surf(dgu, u, &
                  hprogram,yrecfm,i%XFWTD(:),iresp,hcomment=ycomment)
    !
    yrecfm='WTD_ISBA'
    ycomment='water table depth from RRM model or observation (m)'
    CALL write_surf(dgu, u, &
                  hprogram,yrecfm,i%XWTD(:),iresp,hcomment=ycomment)
    !
  ENDIF
  !*       3.     Miscellaneous fields for each patch :
  !               -------------------------------------
  !
  !----------------------------------------------------------------------------
  !User wants (or not) patch output
  IF(dgi%LPATCH_BUDGET)THEN
    !----------------------------------------------------------------------------
    !
    !        3.1    Soil Wetness Index and active layer depth
    !               -----------------------------------------   
    !
    DO jlayer=1,i%NGROUND_LAYER
      !
      WRITE(ylvl,'(I2)') jlayer
      !
      yrecfm='SWI'//adjustl(ylvl(:len_trim(ylvl)))
      yform='(A39,I1.1,A4)'
      IF (jlayer >= 10)  yform='(A39,I2.2,A4)'
      WRITE(ycomment,fmt=yform) 'soil wetness index per patch for layer ',jlayer,' (-)'
      CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XSWI(:,jlayer,:),iresp,hcomment=ycomment)
      !
      yrecfm='TSWI'//adjustl(ylvl(:len_trim(ylvl)))
      yform='(A39,I1.1,A4)'
      IF (jlayer >= 10)  yform='(A39,I2.2,A4)'
      WRITE(ycomment,fmt=yform) 'total swi (liquid+solid) per patch for layer ',jlayer,' (-)'
      CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XTSWI(:,jlayer,:),iresp,hcomment=ycomment)
      !
    END DO
    !
    IF(i%CISBA=='DIF')THEN
      !
      yrecfm='ALT_P'
      ycomment='active layer thickness over permafrost per patch (m)'
      CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XALT(:,:),iresp,hcomment=ycomment)
      !
      yrecfm='FLT_P'
      ycomment='frozen layer thickness over non-permafrost per patch (m)'
      CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XFLT(:,:),iresp,hcomment=ycomment) 
      !
    ENDIF
    !    
    !        3.2    Snow fractions
    !               --------------
    !
    yrecfm='PSNG_P'
    ycomment='snow fraction per patch over ground '
    CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XDPSNG(:,:),iresp,hcomment=ycomment)
    !
    yrecfm='PSNV_P'
    ycomment='snow fraction per patch over vegetation'
    CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XDPSNV(:,:),iresp,hcomment=ycomment)
    !
    yrecfm='PSN_P'
    ycomment='total snow fraction per patch'
    CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XDPSN(:,:),iresp,hcomment=ycomment)
    !
    !        3.3    SGH scheme
    !               ----------
    !
    IF(i%CRUNOFF=='DT92')THEN     
      yrecfm='FSAT_P'
      ycomment='Soil saturated fraction per patch (-)'
      CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XDFSAT(:,:),iresp,hcomment=ycomment)
    ENDIF
    !
    !        3.3    Flood fractions
    !               --------------
    !
    IF(i%LFLOOD)THEN
      !        
      yrecfm='FFG_P'
      ycomment='flood fraction per patch over ground '
      CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XDFFG(:,:),iresp,hcomment=ycomment)
      !
      yrecfm='FFV_P'
      ycomment='flood fraction per patch over vegetation'
      CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XDFFV(:,:),iresp,hcomment=ycomment)
      !
      yrecfm='FF_P'
      ycomment='total flood fraction per patch'
      CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XDFF(:,:),iresp,hcomment=ycomment)
      !
    ENDIF
    !
    !        3.4    Total Albedo
    !               ------------
    !
    yrecfm='TALB'
    ycomment='total albedo per patch'
    !
    CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XALBT(:,:),iresp,hcomment=ycomment)
    !
    IF (i%TSNOW%SCHEME=='3-L' .OR. i%TSNOW%SCHEME=='CRO') THEN
      yrecfm='TS_P'
      ycomment='total surface temperature (isba+snow) per patch'
      CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgi%XTS(:,:),iresp,hcomment=ycomment)
      yrecfm='TSRAD_P'
      ycomment='total radiative surface temperature (isba+snow) per patch'
      CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgi%XTSRAD(:,:),iresp,hcomment=ycomment)
    ENDIF
    !
    !        3.5    Halstead coefficient
    !               --------------------
    !
    yrecfm='HV'
    ycomment='Halstead coefficient per patch'
    CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XHV(:,:),iresp,hcomment=ycomment)
    !
    !        3.6  Snow outputs 
    !        -----------------
    !
    yrecfm='WSN_T_P'
    ycomment='X_Y_WSNOW_TOT (kg/m2) per patch'
    CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XTWSNOW(:,:),iresp,hcomment=ycomment)
    !
    yrecfm='DSN_T_P'
    ycomment='X_Y_DSNOW_TOT (m) per patch'
    CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XTDSNOW(:,:),iresp,hcomment=ycomment)
    !
    yrecfm='TSN_T_P'
    ycomment='X_Y_TSNOW_TOT (k) per patch'
    CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XTTSNOW(:,:),iresp,hcomment=ycomment)
    !
    IF (i%TSNOW%SCHEME=='3-L' .OR. i%TSNOW%SCHEME=='CRO') THEN
      !
      DO jlayer=1,i%TSNOW%NLAYER
        !
        WRITE(ylvl,'(I2)') jlayer
        !
        yrecfm='SNOWLIQ'//adjustl(ylvl(:len_trim(ylvl)))
        yform='(A17,I1.1,A4)'
        IF (jlayer >= 10)  yform='(A17,I2.2,A4)'
        WRITE(ycomment,fmt=yform) 'snow liquid water',jlayer,' (m)'
        CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XSNOWLIQ(:,jlayer,:),iresp,hcomment=ycomment)
        !
        yrecfm='SNOWTEMP'//adjustl(ylvl(:len_trim(ylvl)))
        yform='(A16,I1.1,A4)'
        IF (jlayer >= 10)  yform='(A16,I2.2,A4)'
        WRITE(ycomment,fmt=yform) 'snow temperature',jlayer,' (K)'
        CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XSNOWTEMP(:,jlayer,:),iresp,hcomment=ycomment)
        !
      END DO
      !        
    ENDIF
    !
  END IF
  !
  IF (lagrip) THEN
    !
    !        2.8    Irrigation threshold
    !               --------------------
    !
    yrecfm='IRRISEUIL'
    ycomment='irrigation threshold per patch'
    CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XSEUIL(:,:),iresp,hcomment=ycomment)
    !
  ENDIF
  !
  IF (i%LTR_ML) THEN
    !
    yrecfm='FAPAR'
    ycomment='FAPAR (-)'
    CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XFAPAR(:,:),iresp,hcomment=ycomment)
    !
    yrecfm='FAPIR'
    ycomment='FAPIR (-)'
    CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XFAPIR(:,:),iresp,hcomment=ycomment)
    !
    yrecfm='FAPAR_BS'
    ycomment='FAPAR_BS (-)'
    CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XFAPAR_BS(:,:),iresp,hcomment=ycomment)
    !
    yrecfm='FAPIR_BS'
    ycomment='FAPIR_BS (-)'
    CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XFAPIR_BS(:,:),iresp,hcomment=ycomment)
    !
    yrecfm='DFAPARC'
    ycomment='DFAPARC (-)'
    CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XDFAPARC(:,:),iresp,hcomment=ycomment)
    !
    yrecfm='DFAPIRC'
    ycomment='DFAPIRC (-)'
    CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XDFAPIRC(:,:),iresp,hcomment=ycomment)
    !
    yrecfm='DLAI_EFFC'
    ycomment='DLAI_EFFC (m2/m2)'
    CALL write_surf(dgu, u, &
                  hprogram,yrecfm,dgmi%XDLAI_EFFC(:,:),iresp,hcomment=ycomment)
    !
  ENDIF
  ! 
  IF (lassim .AND. cassim_isba=="EKF  ") THEN
    !
    DO jvar = 1,nvar
      WRITE(yvar,fmt='(I1.1)') jvar
      yrecfm="ANAL_INCR"//yvar
      ycomment="by patch"
      CALL write_surf(dgu, u, &
                      hprogram,yrecfm,i%XINCR(:,i%NPATCH*(jvar-1)+1:i%NPATCH*jvar),iresp,hcomment=ycomment)
    ENDDO
    !
  ENDIF
  !
ENDIF
!         End of IO
!
 CALL end_io_surf_n(hprogram)
IF (lhook) CALL dr_hook('WRITE_DIAG_MISC_ISBA_N',1,zhook_handle)
!
END SUBROUTINE write_diag_misc_isba_n