write_diag_seb_flaken.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_flake_n (DTCO, DGU, U, CHF, DGF, &
                                         hprogram)
!     #################################
!
!!****  *WRITE_DIAG_SEB_FLAKE_n* - writes FLAKE diagnostics
!!
!!    PURPOSE
!!    -------
!!
!!
!!**  METHOD
!!    ------
!!          
!!    REFERENCE
!!    ---------
!!
!!
!!    AUTHOR
!!    ------
!!      V. Masson   *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    01/2004
!!      Modified    01/2006 : sea flux parameterization.
!!      P.LeMoigne    04/2013 : Add accumulated diagnostics
!!      Modified    04/2013, P. Le Moigne: FLake chemistry
!!      S. Belamari 06/2014 : Introduce NBLOCK to avoid errors due to NBLOCK=0
!!                            when coupled with ARPEGE/ALADIN/AROME
!!      B. Decharme 02/2016 : NBLOCK instead of LCOUNTW for compilation in AAA
!-------------------------------------------------------------------------------
!
!*       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_ch_flake_n, ONLY : ch_flake_t
USE modd_diag_flake_n, ONLY : diag_flake_t
!
#ifdef SFX_ARO
USE modd_io_surf_aro,   ONLY : nblock
#endif
!
USE modd_surf_par,      ONLY : xundef
!
!
!
!
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(ch_flake_t), INTENT(INOUT) :: chf
TYPE(diag_flake_t), INTENT(INOUT) :: dgf
!
 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
!
LOGICAL           :: greset
INTEGER           :: jsv, jsw
REAL(KIND=JPRB)   :: zhook_handle
!-------------------------------------------------------------------------------
!
!         Initialisation for IO
!
IF (lhook) CALL dr_hook('WRITE_DIAG_SEB_FLAKE_N',0,zhook_handle)
!
greset=.true.
#ifdef SFX_ARO
greset=(nblock>0)
#endif
!
 CALL init_io_surf_n(dtco, dgu, u, &
                     hprogram,'WATER ','FLAKE ','WRITE')
!
!
!*       2.     Richardson number :
!               -----------------
!
IF (dgf%N2M>=1) THEN

yrecfm='RI_WAT'
ycomment='Bulk-Richardson number for water'
!
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,dgf%XRI(:),iresp,hcomment=ycomment)
!
END IF
!
!*       3.     Energy fluxes :
!               -------------
!
IF (dgf%LSURF_BUDGET) THEN

yrecfm='RN_WAT'
ycomment='net radiation for water'//' (W/m2)'
!
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,dgf%XRN(:),iresp,hcomment=ycomment)
!
yrecfm='H_WAT'
ycomment='sensible heat flux for water'//' (W/m2)'
!
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,dgf%XH(:),iresp,hcomment=ycomment)
!
yrecfm='LE_WAT'
ycomment='total latent heat flux for water'//' (W/m2)'
!
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,dgf%XLE(:),iresp,hcomment=ycomment)
!
yrecfm='LEI_WAT'
ycomment='sublimation latent heat flux for water-ice'//' (W/m2)'
!
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,dgf%XLEI(:),iresp,hcomment=ycomment)
!
yrecfm='GFLUX_WAT'
ycomment='conduction flux for water'//' (W/m2)'
!
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,dgf%XGFLUX(:),iresp,hcomment=ycomment)
!
yrecfm='EVAP_WAT'
ycomment='X_Y_'//yrecfm//' (kg/m2/s)'
!
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,dgf%XEVAP(:),iresp,hcomment=ycomment)
!
yrecfm='SUBL_WAT'
ycomment='X_Y_'//yrecfm//' (kg/m2/s)'
!
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,dgf%XSUBL(:),iresp,hcomment=ycomment)
!
IF (dgf%LRAD_BUDGET) THEN
!
   yrecfm='SWD_WAT'
   ycomment='short wave downward radiation for water'//' (W/m2)'
   !
   CALL write_surf(dgu, u, &
                 hprogram,yrecfm,dgf%XSWD(:),iresp,hcomment=ycomment)
   !
   yrecfm='SWU_WAT'
   ycomment='short wave upward radiation for water'//' (W/m2)'
   !
   CALL write_surf(dgu, u, &
                 hprogram,yrecfm,dgf%XSWU(:),iresp,hcomment=ycomment)
   !
   yrecfm='LWD_WAT'
   ycomment='downward long wave radiation'//' (W/m2)'
   !
   CALL write_surf(dgu, u, &
                 hprogram,yrecfm,dgf%XLWD(:),iresp,hcomment=ycomment)
   !
   yrecfm='LWU_WAT'
   ycomment='upward long wave radiation'//' (W/m2)'
   !
   CALL write_surf(dgu, u, &
                 hprogram,yrecfm,dgf%XLWU(:),iresp,hcomment=ycomment)
   !       
   DO jsw=1, SIZE(dgf%XSWBD,2)
      ynum=achar(48+jsw)
      !
      yrecfm='SWD_WAT_'//ynum
      ycomment='downward short wave radiation by spectral band '//' (W/m2)'
      !
      CALL write_surf(dgu, u, &
                 hprogram,yrecfm,dgf%XSWBD(:,jsw),iresp,hcomment=ycomment)
      !
      yrecfm='SWU_WAT_'//ynum
      ycomment='upward short wave radiation by spectral band'//' (W/m2)'
      !
      CALL write_surf(dgu, u, &
                 hprogram,yrecfm,dgf%XSWBU(:,jsw),iresp,hcomment=ycomment)
      !
   ENDDO
!
ENDIF
!
yrecfm='FMU_WAT'
ycomment='u-component of momentum flux for water'//' (kg/ms2)'
!
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,dgf%XFMU(:),iresp,hcomment=ycomment)
yrecfm='FMV_WAT'
ycomment='v-component of momentum flux for water'//' (kg/ms2)'
!
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,dgf%XFMV(:),iresp,hcomment=ycomment)
!
END IF
!
IF (dgf%LSURF_BUDGET.OR.dgf%LSURF_BUDGETC) THEN
!
  yrecfm='TALB_WAT'
  ycomment='total albedo over tile water (-)'
  CALL write_surf(dgu, u, &
                 hprogram,yrecfm,dgf%XALBT(:),iresp,hcomment=ycomment)
!
  yrecfm='WSN_WAT'
  ycomment='snow water equivalent over tile water (-)'
  CALL write_surf(dgu, u, &
                 hprogram,yrecfm,dgf%XSWE(:),iresp,hcomment=ycomment)
!        
ENDIF
!
IF (dgf%LSURF_BUDGETC) THEN
!
yrecfm='RNC_WAT'
ycomment='X_Y_'//yrecfm//' (J/m2)'
!
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,dgf%XRNC(:),iresp,hcomment=ycomment)
!
yrecfm='HC_WAT'
ycomment='X_Y_'//yrecfm//' (J/m2)'
!
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,dgf%XHC(:),iresp,hcomment=ycomment)
!
yrecfm='LEC_WAT'
ycomment='X_Y_'//yrecfm//' (J/m2)'
!
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,dgf%XLEC(:),iresp,hcomment=ycomment)
!
yrecfm='LEIC_WAT'
ycomment='X_Y_'//yrecfm//' (J/m2)'
!
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,dgf%XLEIC(:),iresp,hcomment=ycomment)
!
yrecfm='GFLUXC_WAT'
ycomment='X_Y_'//yrecfm//' (J/m2)'
!
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,dgf%XGFLUXC(:),iresp,hcomment=ycomment)
!
yrecfm='EVAPC_WAT'
ycomment='X_Y_'//yrecfm//' (kg/m2)'
!
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,dgf%XEVAPC(:),iresp,hcomment=ycomment)
!
yrecfm='SUBLC_WAT'
ycomment='X_Y_'//yrecfm//' (kg/m2)'
!
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,dgf%XSUBLC(:),iresp,hcomment=ycomment)
!
IF (dgf%LRAD_BUDGET .OR. (dgf%LSURF_BUDGETC .AND. .NOT.dgu%LRESET_BUDGETC)) THEN
!
   yrecfm='SWDC_WAT'
   ycomment='X_Y_'//yrecfm//' (J/m2)'
   !
   CALL write_surf(dgu, u, &
                 hprogram,yrecfm,dgf%XSWDC(:),iresp,hcomment=ycomment)
   !
   yrecfm='SWUC_WAT'
   ycomment='X_Y_'//yrecfm//' (J/m2)'
   !
   CALL write_surf(dgu, u, &
                 hprogram,yrecfm,dgf%XSWUC(:),iresp,hcomment=ycomment)
   !
   yrecfm='LWDC_WAT'
   ycomment='X_Y_'//yrecfm//' (J/m2)'
   !
   CALL write_surf(dgu, u, &
                 hprogram,yrecfm,dgf%XLWDC(:),iresp,hcomment=ycomment)
   !
   yrecfm='LWUC_WAT'
   ycomment='X_Y_'//yrecfm//' (J/m2)'
   !
   CALL write_surf(dgu, u, &
                 hprogram,yrecfm,dgf%XLWUC(:),iresp,hcomment=ycomment)
!
ENDIF
!
yrecfm='FMUC_WAT'
ycomment='X_Y_'//yrecfm//' (kg/ms)'
!
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,dgf%XFMUC(:),iresp,hcomment=ycomment)
!
yrecfm='FMVC_WAT'
ycomment='X_Y_'//yrecfm//' (kg/ms)'
!
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,dgf%XFMVC(:),iresp,hcomment=ycomment)
!
END IF
!
!
!*       4.     Transfer coefficients
!               ---------------------
!
IF (dgf%LCOEF) THEN

yrecfm='CD_WAT'
ycomment='drag coefficient for wind over water (W/s2)'
!
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,dgf%XCD(:),iresp,hcomment=ycomment)
!
yrecfm='CH_WAT'
ycomment='drag coefficient for heat (W/s)'
!
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,dgf%XCH(:),iresp,hcomment=ycomment)
!
yrecfm='CE_WAT'
ycomment='drag coefficient for vapor (W/s/K)'
!
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,dgf%XCE(:),iresp,hcomment=ycomment)
!
yrecfm='Z0_WAT'
ycomment='roughness length over water (m)'

 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,dgf%XZ0(:),iresp,hcomment=ycomment)
!
yrecfm='Z0H_WAT'
ycomment='thermal roughness length over water (m)'
!
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,dgf%XZ0H(:),iresp,hcomment=ycomment)
!
END IF
!
!
!*       5.     Surface humidity
!               ----------------
!
IF (dgf%LSURF_VARS) THEN

yrecfm='QS_WAT'
ycomment='specific humidity over water'//' (KG/KG)'
!
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,dgf%XQS(:),iresp,hcomment=ycomment)
!
ENDIF
!

!
!*       6.     parameters at 2 and 10 meters :
!               -----------------------------
!
IF (dgf%N2M>=1) THEN
!
yrecfm='T2M_WAT'
ycomment='2 meters temperature'//' (K)'
!
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,dgf%XT2M(:),iresp,hcomment=ycomment)
!
yrecfm='T2MMIN_WAT'
ycomment='X_Y_'//yrecfm//' (K)'
!
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,dgf%XT2M_MIN(:),iresp,hcomment=ycomment)
IF(greset)dgf%XT2M_MIN(:)=xundef
!
yrecfm='T2MMAX_WAT'
ycomment='X_Y_'//yrecfm//' (K)'
!
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,dgf%XT2M_MAX(:),iresp,hcomment=ycomment)
IF(greset)dgf%XT2M_MAX(:)=0.0
!
yrecfm='Q2M_WAT'
ycomment='2 meters specific humidity'//' (KG/KG)'
!
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,dgf%XQ2M(:),iresp,hcomment=ycomment)
!
yrecfm='HU2M_WAT'
ycomment='2 meters relative humidity'//' (KG/KG)'
!
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,dgf%XHU2M(:),iresp,hcomment=ycomment)
!
yrecfm='HU2MMIN_WAT'
ycomment='X_Y_'//yrecfm//' (-)'
!
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,dgf%XHU2M_MIN(:),iresp,hcomment=ycomment)
IF(greset)dgf%XHU2M_MIN(:)=xundef
!
yrecfm='HU2MMAX_WAT'
ycomment='X_Y_'//yrecfm//' (-)'
!
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,dgf%XHU2M_MAX(:),iresp,hcomment=ycomment)
IF(greset)dgf%XHU2M_MAX(:)=-xundef
!
yrecfm='ZON10M_WAT'
ycomment='10 meters zonal wind'//' (M/S)'
!
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,dgf%XZON10M(:),iresp,hcomment=ycomment)
!
yrecfm='MER10M_WAT'
ycomment='10 meters meridian wind'//' (M/S)'
!
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,dgf%XMER10M(:),iresp,hcomment=ycomment)
!
yrecfm='W10M_WAT'
ycomment='X_Y_'//yrecfm//' (M/S)'
!
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,dgf%XWIND10M(:),iresp,hcomment=ycomment)
!
yrecfm='W10MMAX_WAT'
ycomment='X_Y_'//yrecfm//' (M/S)'
!
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,dgf%XWIND10M_MAX(:),iresp,hcomment=ycomment)
IF(greset)dgf%XWIND10M_MAX(:)=0.0
!
END IF
!
!
!*       7.     chemical diagnostics:
!               --------------------
!
IF (chf%SVF%NBEQ>0 .AND. chf%CCH_DRY_DEP=="WES89 ") THEN
  DO jsv = 1,SIZE(chf%CCH_NAMES,1)
    yrecfm='DV_WAT_'//trim(chf%CCH_NAMES(jsv))
    WRITE(ycomment,'(A13,I3.3)')'(m/s) DV_WAT_',jsv
    CALL write_surf(dgu, u, &
                 hprogram,yrecfm,chf%XDEP(:,jsv),iresp,hcomment=ycomment)
  END DO
ENDIF
!
!
!*       8.     prognostic variable diagnostics:
!               --------------------------------
!
IF(dgu%LPROVAR_TO_DIAG)THEN
!
  yrecfm='TS_WATER'
  ycomment='TS_WATER (K)'
  CALL write_surf(dgu, u, &
                 hprogram,yrecfm,dgf%XDIAG_TS(:),iresp,hcomment=ycomment)
!
ENDIF
!
!-------------------------------------------------------------------------------
!
!         End of IO
!
 CALL end_io_surf_n(hprogram)
IF (lhook) CALL dr_hook('WRITE_DIAG_SEB_FLAKE_N',1,zhook_handle)
!
!
END SUBROUTINE write_diag_seb_flake_n