write_diag_seb_surf_atmn.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_surf_atm_n (DTCO, DGO, D, DC, U, HGRID, HPROGRAM)
!     #################################
!
!!****  *WRITE_DIAG_SEB_SURF_ATM_n* - writes surface diagnostics
!!
!!    PURPOSE
!!    -------
!!
!!
!!**  METHOD
!!    ------
!!
!!
!!    REFERENCE
!!    ---------
!!
!!
!!    AUTHOR
!!    ------
!!      V. Masson   *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    01/2004
!!      Modified    01/2006 : sea flux parameterization.
!!      Modified    08/2009 : cumulated diag
!!      Juan        6/12/2011: parallel bug , remove local ANY(XAVG_ZON10M) test
!!      B. Decharme  06/13   Add QS, evap and sublimation diags
!-------------------------------------------------------------------------------
!
!*       0.    DECLARATIONS
!              ------------
!
!
USE modd_data_cover_n, ONLY : data_cover_t
USE modd_diag_n, ONLY : diag_t, diag_options_t
USE modd_surf_atm_n, ONLY : surf_atm_t
!
USE modd_surf_par, ONLY : xundef
!
USE modd_xios, ONLY : lallow_add_dim, yswband_dim_name
!
USE modi_init_io_surf_n
USE modi_write_surf
USE modi_end_io_surf_n
USE modi_sum_on_all_procs
!
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) :: DGO
TYPE(diag_t), INTENT(INOUT) :: D
TYPE(diag_t), INTENT(INOUT) :: DC
TYPE(surf_atm_t), INTENT(INOUT) :: U
 CHARACTER(LEN=*), INTENT(IN) :: HGRID
!
 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=12) :: YRECFM         ! Name of the article to be read
 CHARACTER(LEN=100):: YCOMMENT       ! Comment string
 CHARACTER(LEN=2)  :: YNUM
!
INTEGER           :: JSW
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
!-------------------------------------------------------------------------------
!
!         Initialisation for IO
!
IF (lhook) CALL dr_hook('WRITE_DIAG_SEB_SURF_ATM_N',0,zhook_handle)
 CALL init_io_surf_n(dtco, u, hprogram,'FULL  ','SURF  ','WRITE','SURF_ATM_DIAGNOSTICS.OUT.nc')
!
!*       1.     Richardson number :
!               -----------------
!
IF (dgo%N2M>=1) THEN
  !        
  yrecfm='RI'
  ycomment='X_Y_'//yrecfm
  !
  CALL write_surf(dgo%CSELECT,hprogram,yrecfm,d%XRI(:),iresp,hcomment=ycomment)
  !
ENDIF
!
!*       2.     parameters at surface, 2 and 10 meters :
!               ----------------------------------------
!
IF (dgo%N2M>=1.OR.dgo%LSURF_BUDGET.OR.dgo%LSURF_BUDGETC) THEN
  !
  yrecfm='TS'
  ycomment='X_Y_'//yrecfm//' (K)'
  CALL write_surf(dgo%CSELECT,hprogram,yrecfm,d%XTS(:),iresp,hcomment=ycomment)
  !
  yrecfm='TSRAD'
  ycomment='X_Y_'//yrecfm//' (K)'
  CALL write_surf(dgo%CSELECT,hprogram,yrecfm,d%XTRAD(:),iresp,hcomment=ycomment)
  !
  yrecfm='EMIS'
  ycomment='X_Y_'//yrecfm//' (-)'
  CALL write_surf(dgo%CSELECT,hprogram,yrecfm,d%XEMIS(:),iresp,hcomment=ycomment)
  !
  yrecfm='SFCO2'
  ycomment='X_Y_'//yrecfm//' (M.kgCO2.S-1.kgAIR-1)'
  CALL write_surf(dgo%CSELECT,hprogram,yrecfm,d%XSFCO2(:),iresp,hcomment=ycomment)
  !
ENDIF
!
IF (dgo%N2M>=1) THEN
  !
  yrecfm='T2M'
  ycomment='X_Y_'//yrecfm//' (K)'
  CALL write_surf(dgo%CSELECT,hprogram,yrecfm,d%XT2M(:),iresp,hcomment=ycomment)
  !
  yrecfm='T2MMIN'
  ycomment='X_Y_'//yrecfm//' (K)'
  CALL write_surf(dgo%CSELECT,hprogram,yrecfm,d%XT2M_MIN(:),iresp,hcomment=ycomment)
  !
  yrecfm='T2MMAX'
  ycomment='X_Y_'//yrecfm//' (K)'
  CALL write_surf(dgo%CSELECT,hprogram,yrecfm,d%XT2M_MAX(:),iresp,hcomment=ycomment)
  !
  yrecfm='Q2M'
  ycomment='X_Y_'//yrecfm//' (KG/KG)'
  CALL write_surf(dgo%CSELECT,hprogram,yrecfm,d%XQ2M(:),iresp,hcomment=ycomment)
  !
  yrecfm='HU2M'
  ycomment='X_Y_'//yrecfm//' (-)'
  CALL write_surf(dgo%CSELECT,hprogram,yrecfm,d%XHU2M(:),iresp,hcomment=ycomment)
  !
  yrecfm='HU2MMIN'
  ycomment='X_Y_'//yrecfm//' (-)'
  CALL write_surf(dgo%CSELECT,hprogram,yrecfm,d%XHU2M_MIN(:),iresp,hcomment=ycomment)
  !
  yrecfm='HU2MMAX'
  ycomment='X_Y_'//yrecfm//' (-)'
  CALL write_surf(dgo%CSELECT,hprogram,yrecfm,d%XHU2M_MAX(:),iresp,hcomment=ycomment)
  !
  IF ( sum_on_all_procs(hprogram,hgrid,d%XZON10M(:)/= xundef) > 0. ) THEN
    !
    yrecfm='ZON10M'
    ycomment='X_Y_'//yrecfm//' (M/S)'
    CALL write_surf(dgo%CSELECT,hprogram,yrecfm,d%XZON10M(:),iresp,hcomment=ycomment)
    !
    yrecfm='MER10M'
    ycomment='X_Y_'//yrecfm//' (M/S)'
    CALL write_surf(dgo%CSELECT,hprogram,yrecfm,d%XMER10M(:),iresp,hcomment=ycomment)
    !
    yrecfm='W10M'
    ycomment='X_Y_'//yrecfm//' (M/S)'
    CALL write_surf(dgo%CSELECT,hprogram,yrecfm,d%XWIND10M(:),iresp,hcomment=ycomment)
    !
    yrecfm='W10MMAX'
    ycomment='X_Y_'//yrecfm//' (M/S)'
    CALL write_surf(dgo%CSELECT,hprogram,yrecfm,d%XWIND10M_MAX(:),iresp,hcomment=ycomment)
    !
  ENDIF
  !
  IF (dgo%L2M_MIN_ZS) THEN
    !
    yrecfm='T2M_MIN_ZS'
    ycomment='X_Y_'//yrecfm//' (K)'
    CALL write_surf(dgo%CSELECT,hprogram,yrecfm,d%XT2M_MIN_ZS(:),iresp,hcomment=ycomment)
    !
    yrecfm='Q2M_MIN_ZS'
    ycomment='X_Y_'//yrecfm//' (KG/KG)'
    CALL write_surf(dgo%CSELECT,hprogram,yrecfm,d%XQ2M_MIN_ZS(:),iresp,hcomment=ycomment)
    !
    yrecfm='HU2M_MIN_ZS'
    ycomment='X_Y_'//yrecfm//' (KG/KG)'
    CALL write_surf(dgo%CSELECT,hprogram,yrecfm,d%XHU2M_MIN_ZS(:),iresp,hcomment=ycomment)
    !
  END IF
  !
END IF
!
!*       3.     Energy fluxes :
!               -------------
!
IF (dgo%LSURF_BUDGET) THEN
  !
  yrecfm='RN'
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  CALL write_surf(dgo%CSELECT,hprogram,yrecfm,d%XRN(:),iresp,hcomment=ycomment)
  !
  yrecfm='H'
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  CALL write_surf(dgo%CSELECT,hprogram,yrecfm,d%XH(:),iresp,hcomment=ycomment)
  !
  yrecfm='LE'
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  CALL write_surf(dgo%CSELECT,hprogram,yrecfm,d%XLE(:),iresp,hcomment=ycomment)
  !
  yrecfm='LEI'
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  CALL write_surf(dgo%CSELECT,hprogram,yrecfm,d%XLEI(:),iresp,hcomment=ycomment)
  !
  yrecfm='GFLUX'
  ycomment='X_Y_'//yrecfm//' (W/m2)'
  CALL write_surf(dgo%CSELECT,hprogram,yrecfm,d%XGFLUX(:),iresp,hcomment=ycomment)
  !
  yrecfm='EVAP'
  ycomment='X_Y_'//yrecfm//' (kg/m2/s)'
  !
  CALL write_surf(dgo%CSELECT,hprogram,yrecfm,d%XEVAP(:),iresp,hcomment=ycomment)
  !
  yrecfm='SUBL'
  ycomment='X_Y_'//yrecfm//' (kg/m2/s)'
  !
  CALL write_surf(dgo%CSELECT,hprogram,yrecfm,d%XSUBL(:),iresp,hcomment=ycomment)
  !
  IF (dgo%LRAD_BUDGET) THEN
    !         
    yrecfm='SWD'
    ycomment='X_Y_'//yrecfm//' (W/m2)'
    CALL write_surf(dgo%CSELECT,hprogram,yrecfm,d%XSWD(:),iresp,hcomment=ycomment)
    !
    yrecfm='SWU'
    ycomment='X_Y_'//yrecfm//' (W/m2)'
    CALL write_surf(dgo%CSELECT,hprogram,yrecfm,d%XSWU(:),iresp,hcomment=ycomment)
    !
    yrecfm='LWD'
    ycomment='X_Y_'//yrecfm//' (W/m2)'
    CALL write_surf(dgo%CSELECT,hprogram,yrecfm,d%XLWD(:),iresp,hcomment=ycomment)
    !
    yrecfm='LWU'
    ycomment='X_Y_'//yrecfm//' (W/m2)'
    CALL write_surf(dgo%CSELECT,hprogram,yrecfm,d%XLWU(:),iresp,hcomment=ycomment)
    !
    IF (lallow_add_dim)  THEN 
      !
      yrecfm='SWD_'
      ycomment='X_Y_'//yrecfm//' (W/m2)'
      CALL write_surf(dgo%CSELECT,&
           hprogram,yrecfm,d%XSWBD(:,:),iresp,hcomment=ycomment, hnam_dim=yswband_dim_name)
      !
      yrecfm='SWU_'
      ycomment='X_Y_'//yrecfm//' (W/m2)'
      CALL write_surf(dgo%CSELECT,&
           hprogram,yrecfm,d%XSWBD(:,:),iresp,hcomment=ycomment, hnam_dim=yswband_dim_name)
       !
    ELSE
      !    
      DO jsw=1, SIZE(d%XSWBD,2)
        ynum=achar(48+jsw)
        !
        yrecfm='SWD_'//ynum
        ycomment='X_Y_'//yrecfm//' (W/m2)'
        CALL write_surf(dgo%CSELECT,hprogram,yrecfm,d%XSWBD(:,jsw),iresp,hcomment=ycomment)
        !
        yrecfm='SWU_'//ynum
        ycomment='X_Y_'//yrecfm//' (W/m2)'
        CALL write_surf(dgo%CSELECT,hprogram,yrecfm,d%XSWBU(:,jsw),iresp,hcomment=ycomment)
        !
      ENDDO
      !
    ENDIF
    !
  ENDIF
  !
  yrecfm='FMUNOSSO'
  ycomment='X_Y_'//yrecfm//' (kg/ms2)'
  CALL write_surf(dgo%CSELECT,hprogram,yrecfm,d%XFMU(:),iresp,hcomment=ycomment)
  !
  yrecfm='FMVNOSSO'
  ycomment='X_Y_'//yrecfm//' (kg/ms2)'
  CALL write_surf(dgo%CSELECT,hprogram,yrecfm,d%XFMV(:),iresp,hcomment=ycomment)
  !
  yrecfm='FMU'
  ycomment='X_Y_'//yrecfm//' (kg/ms2)'
  CALL write_surf(dgo%CSELECT,hprogram,yrecfm,d%XSSO_FMU(:),iresp,hcomment=ycomment)
  !
  yrecfm='FMV'
  ycomment='X_Y_'//yrecfm//' (kg/ms2)'
  CALL write_surf(dgo%CSELECT,hprogram,yrecfm,d%XSSO_FMV(:),iresp,hcomment=ycomment)
  !
END IF
!
! * Cumulated diag
!
IF (dgo%LSURF_BUDGETC) THEN
  !
  yrecfm='RNC'
  ycomment='X_Y_'//yrecfm//' (J/m2)'
  CALL write_surf(dgo%CSELECT,hprogram,yrecfm,dc%XRN(:),iresp,hcomment=ycomment)
  !
  yrecfm='HC'
  ycomment='X_Y_'//yrecfm//' (J/m2)'
  CALL write_surf(dgo%CSELECT,hprogram,yrecfm,dc%XH(:),iresp,hcomment=ycomment)
  !
  yrecfm='LEC'
  ycomment='X_Y_'//yrecfm//' (J/m2)'
  CALL write_surf(dgo%CSELECT,hprogram,yrecfm,dc%XLE(:),iresp,hcomment=ycomment)
  !
  yrecfm='LEIC'
  ycomment='X_Y_'//yrecfm//' (J/m2)'
  CALL write_surf(dgo%CSELECT,hprogram,yrecfm,dc%XLEI(:),iresp,hcomment=ycomment)
  !
  yrecfm='GFLUXC'
  ycomment='X_Y_'//yrecfm//' (J/m2)'
  CALL write_surf(dgo%CSELECT,hprogram,yrecfm,dc%XGFLUX(:),iresp,hcomment=ycomment)
  !
  yrecfm='EVAPC'
  ycomment='X_Y_'//yrecfm//' (kg/m2)'
  !
  CALL write_surf(dgo%CSELECT,hprogram,yrecfm,dc%XEVAP(:),iresp,hcomment=ycomment)
  !
  yrecfm='SUBLC'
  ycomment='X_Y_'//yrecfm//' (kg/m2)'
  !
  CALL write_surf(dgo%CSELECT,hprogram,yrecfm,dc%XSUBL(:),iresp,hcomment=ycomment)
  !
  IF (dgo%LRAD_BUDGET .OR. (dgo%LSURF_BUDGETC .AND. .NOT.dgo%LRESET_BUDGETC)) THEN
    !        
    yrecfm='SWDC'
    ycomment='X_Y_'//yrecfm//' (J/m2)'
    CALL write_surf(dgo%CSELECT,hprogram,yrecfm,dc%XSWD(:),iresp,hcomment=ycomment)
    !
    yrecfm='SWUC'
    ycomment='X_Y_'//yrecfm//' (J/m2)'
    CALL write_surf(dgo%CSELECT,hprogram,yrecfm,dc%XSWU(:),iresp,hcomment=ycomment)
    !
    yrecfm='LWDC'
    ycomment='X_Y_'//yrecfm//' (J/m2)'
    CALL write_surf(dgo%CSELECT,hprogram,yrecfm,dc%XLWD(:),iresp,hcomment=ycomment)
    !
    yrecfm='LWUC'
    ycomment='X_Y_'//yrecfm//' (J/m2)'
    CALL write_surf(dgo%CSELECT,hprogram,yrecfm,dc%XLWU(:),iresp,hcomment=ycomment)
    !
  ENDIF
  !
  yrecfm='FMUC'
  ycomment='X_Y_'//yrecfm//' (kg/ms)'
  CALL write_surf(dgo%CSELECT,hprogram,yrecfm,dc%XFMU(:),iresp,hcomment=ycomment)
  !
  yrecfm='FMVC'
  ycomment='X_Y_'//yrecfm//' (kg/ms)'
  CALL write_surf(dgo%CSELECT,hprogram,yrecfm,dc%XFMV(:),iresp,hcomment=ycomment)
  !
END IF
!
!
!*       4.     Transfer coefficients
!               ---------------------
!
IF (dgo%LCOEF) THEN
  !
  yrecfm='CD'
  ycomment='X_Y_'//yrecfm
  CALL write_surf(dgo%CSELECT,hprogram,yrecfm,d%XCD(:),iresp,hcomment=ycomment)
  !
  yrecfm='CH'
  ycomment='X_Y_'//yrecfm
  CALL write_surf(dgo%CSELECT,hprogram,yrecfm,d%XCH(:),iresp,hcomment=ycomment)
  !
  yrecfm='CE'
  ycomment='X_Y_'//yrecfm
  CALL write_surf(dgo%CSELECT,hprogram,yrecfm,d%XCE(:),iresp,hcomment=ycomment)
  !
  yrecfm='Z0'
  ycomment='X_Y_'//yrecfm
  CALL write_surf(dgo%CSELECT,hprogram,yrecfm,d%XZ0(:),iresp,hcomment=ycomment)
  !
  yrecfm='Z0H'
  ycomment='X_Y_'//yrecfm
  CALL write_surf(dgo%CSELECT,hprogram,yrecfm,d%XZ0H(:),iresp,hcomment=ycomment)
  !
  yrecfm='UREF'
  ycomment='X_Y_'//yrecfm
  CALL write_surf(dgo%CSELECT,hprogram,yrecfm,d%XUREF(:),iresp,hcomment=ycomment)
  !
  yrecfm='ZREF'
  ycomment='X_Y_'//yrecfm
  CALL write_surf(dgo%CSELECT,hprogram,yrecfm,d%XZREF(:),iresp,hcomment=ycomment)
  !
END IF
!
!
!*       5.     Surface humidity
!               ----------------
!
IF (dgo%LSURF_VARS) THEN
!
yrecfm='QS'
ycomment='X_Y_'//yrecfm//' (kg/kg)'
!
 CALL write_surf(dgo%CSELECT,hprogram,yrecfm,d%XQS(:),iresp,hcomment=ycomment)
!
ENDIF
!
!-------------------------------------------------------------------------------
!
!         End of IO
!
 CALL end_io_surf_n(hprogram)
IF (lhook) CALL dr_hook('WRITE_DIAG_SEB_SURF_ATM_N',1,zhook_handle)
!
!
END SUBROUTINE write_diag_seb_surf_atm_n