write_diag_seb_seafluxn.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_seaflux_n (DTCO, DGU, U, CHS, DGS, S, &
                                           hprogram)
!     #################################
!
!!****  *WRITE_DIAG_SEB_SEAFLUX_n* - write the SEAFLUX diagnostic fields
!!
!!    PURPOSE
!!    -------
!!
!!
!!**  METHOD
!!    ------
!!
!!    REFERENCE
!!    ---------
!!
!!
!!    AUTHOR
!!    ------
!!      V. Masson   *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    01/2004
!!      Modified    01/2006 : sea flux parameterization.
!!      Modified    08/2009 : cumulated diag
!!      B. Decharme 06/2013 : Add evap and sublimation diag
!!                            Delete LPROVAR_TO_DIAG here
!!      S.Senesi    01/2014 : add diags on seaice 
!!      S. Belamari 06/2014 : Introduce GRESET 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_seaflux_n, ONLY : ch_seaflux_t
USE modd_diag_seaflux_n, ONLY : diag_seaflux_t
USE modd_seaflux_n, ONLY : seaflux_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_seaflux_t), INTENT(INOUT) :: chs
TYPE(diag_seaflux_t), INTENT(INOUT) :: dgs
TYPE(seaflux_t), INTENT(INOUT) :: s
!
 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
LOGICAL           :: gmisc
!
REAL(KIND=JPRB) :: zhook_handle
!
!-------------------------------------------------------------------------------
!
IF (lhook) CALL dr_hook('WRITE_DIAG_SEB_SEAFLUX_N',0,zhook_handle)
!
!         Initialisation for IO
!
greset=.true.
#ifdef SFX_ARO
greset=(nblock>0)
#endif
!
 CALL init_io_surf_n(dtco, dgu, u, &
                    hprogram,'SEA   ','SEAFLX','WRITE')
!
!
!*       1.     Surface temperature :
!               ---------------------
!
gmisc=(dgs%N2M>=1.OR.dgs%LSURF_BUDGET.OR.dgs%LSURF_BUDGETC)
!
IF (gmisc.AND.s%LHANDLE_SIC) THEN
    !
    yrecfm='TS_SEA'
    ycomment='X_Y_'//yrecfm//' (K)'
    !
    CALL write_surf(dgu, u, &
                    hprogram,yrecfm,dgs%XTS(:),iresp,hcomment=ycomment)
    !
    yrecfm='TSRAD_SEA'
    ycomment='X_Y_'//yrecfm//' (K)'
    !
    CALL write_surf(dgu, u, &
                    hprogram,yrecfm,dgs%XTSRAD(:),iresp,hcomment=ycomment)
    !
ENDIF
!
!*       2.     Richardson number :
!               -----------------
!
IF (dgs%N2M>=1) THEN
   !
   yrecfm='RI_SEA'
   ycomment='X_Y_'//yrecfm
   !
   CALL write_surf(dgu, u, &
                    hprogram,yrecfm,dgs%XRI(:),iresp,hcomment=ycomment)
   !
ENDIF
 !
 !*       3.     Energy fluxes :
 !               -------------
 !
IF (dgs%LSURF_BUDGET) THEN
   !
   yrecfm='RN_SEA'
   ycomment='X_Y_'//yrecfm//' (W/m2)'
   !
   CALL write_surf(dgu, u, &
                    hprogram,yrecfm,dgs%XRN(:),iresp,hcomment=ycomment)
   !
   yrecfm='H_SEA'
   ycomment='X_Y_'//yrecfm//' (W/m2)'
   !
   CALL write_surf(dgu, u, &
                    hprogram,yrecfm,dgs%XH(:),iresp,hcomment=ycomment)
   !
   yrecfm='LE_SEA'
   ycomment='X_Y_'//yrecfm//' (W/m2)'
   !
   CALL write_surf(dgu, u, &
                    hprogram,yrecfm,dgs%XLE(:),iresp,hcomment=ycomment)
   !
   yrecfm='LEI_SEA'
   ycomment='X_Y_'//yrecfm//' (W/m2)'
   CALL write_surf(dgu, u, &
                    hprogram,yrecfm,dgs%XLE_ICE(:),iresp,hcomment=ycomment) 
   !
   yrecfm='GFLUX_SEA'
   ycomment='X_Y_'//yrecfm//' (W/m2)'
   !
   CALL write_surf(dgu, u, &
                    hprogram,yrecfm,dgs%XGFLUX(:),iresp,hcomment=ycomment)
   !
   yrecfm='EVAP_SEA'
   ycomment='X_Y_'//yrecfm//' (kg/m2/s)'
   !
   CALL write_surf(dgu, u, &
                    hprogram,yrecfm,dgs%XEVAP(:),iresp,hcomment=ycomment)
   !
   yrecfm='SUBL_SEA'
   ycomment='X_Y_'//yrecfm//' (kg/m2/s)'
   !
   CALL write_surf(dgu, u, &
                    hprogram,yrecfm,dgs%XSUBL(:),iresp,hcomment=ycomment)
   !
   IF (dgs%LRAD_BUDGET) THEN
      !
      yrecfm='SWD_SEA'
      ycomment='X_Y_'//yrecfm//' (W/m2)'
      !
      CALL write_surf(dgu, u, &
                    hprogram,yrecfm,dgs%XSWD(:),iresp,hcomment=ycomment)
      !
      yrecfm='SWU_SEA'
      ycomment='X_Y_'//yrecfm//' (W/m2)'
      !
      CALL write_surf(dgu, u, &
                    hprogram,yrecfm,dgs%XSWU(:),iresp,hcomment=ycomment)
      !
      yrecfm='LWD_SEA'
      ycomment='X_Y_'//yrecfm//' (W/m2)'
      !
      CALL write_surf(dgu, u, &
                    hprogram,yrecfm,dgs%XLWD(:),iresp,hcomment=ycomment)
      !
      yrecfm='LWU_SEA'
      ycomment='X_Y_'//yrecfm//' (W/m2)'
      !
      CALL write_surf(dgu, u, &
                    hprogram,yrecfm,dgs%XLWU(:),iresp,hcomment=ycomment)
      !
      DO jsw=1, SIZE(dgs%XSWBD,2)
         ynum=achar(48+jsw)
         !
         yrecfm='SWD_SEA_'//ynum
         ycomment='X_Y_'//yrecfm//' (W/m2)'
         !
         CALL write_surf(dgu, u, &
                    hprogram,yrecfm,dgs%XSWBD(:,jsw),iresp,hcomment=ycomment)
         !
         yrecfm='SWU_SEA_'//ynum
         ycomment='X_Y_'//yrecfm//' (W/m2)'
         !
         CALL write_surf(dgu, u, &
                    hprogram,yrecfm,dgs%XSWBU(:,jsw),iresp,hcomment=ycomment)
         !
      ENDDO
      !
   ENDIF
   !
   yrecfm='FMU_SEA'
   ycomment='X_Y_'//yrecfm//' (kg/ms2)'
   !
   CALL write_surf(dgu, u, &
                    hprogram,yrecfm,dgs%XFMU(:),iresp,hcomment=ycomment)
   !
   yrecfm='FMV_SEA'
   ycomment='X_Y_'//yrecfm//' (kg/ms2)'
   !
   CALL write_surf(dgu, u, &
                    hprogram,yrecfm,dgs%XFMV(:),iresp,hcomment=ycomment)
   !
END IF
!
IF (dgs%LSURF_BUDGETC) THEN
   !
   yrecfm='RNC_SEA'
   ycomment='X_Y_'//yrecfm//' (J/m2)'
   !
   CALL write_surf(dgu, u, &
                    hprogram,yrecfm,dgs%XRNC(:),iresp,hcomment=ycomment)
   !
   yrecfm='HC_SEA'
   ycomment='X_Y_'//yrecfm//' (J/m2)'
   !
   CALL write_surf(dgu, u, &
                    hprogram,yrecfm,dgs%XHC(:),iresp,hcomment=ycomment)
   !
   yrecfm='LEC_SEA'
   ycomment='X_Y_'//yrecfm//' (J/m2)'
   !
   CALL write_surf(dgu, u, &
                    hprogram,yrecfm,dgs%XLEC(:),iresp,hcomment=ycomment)
   !
   yrecfm='LEIC_SEA'
   ycomment='X_Y_'//yrecfm//' (W/m2)'
   CALL write_surf(dgu, u, &
                    hprogram,yrecfm,dgs%XLEC_ICE(:),iresp,hcomment=ycomment) 
   !
   yrecfm='GFLUXC_SEA'
   ycomment='X_Y_'//yrecfm//' (J/m2)'
   !
   CALL write_surf(dgu, u, &
                    hprogram,yrecfm,dgs%XGFLUXC(:),iresp,hcomment=ycomment)
   !
   yrecfm='EVAPC_SEA'
   ycomment='X_Y_'//yrecfm//' (kg/m2)'
   !
   CALL write_surf(dgu, u, &
                    hprogram,yrecfm,dgs%XEVAPC(:),iresp,hcomment=ycomment)
   !
   yrecfm='SUBLC_SEA'
   ycomment='X_Y_'//yrecfm//' (kg/m2)'
   !
   CALL write_surf(dgu, u, &
                    hprogram,yrecfm,dgs%XSUBLC(:),iresp,hcomment=ycomment)
   !
   IF (dgs%LRAD_BUDGET .OR. (dgs%LSURF_BUDGETC .AND. .NOT.dgu%LRESET_BUDGETC)) THEN
      !
      yrecfm='SWDC_SEA'
      ycomment='X_Y_'//yrecfm//' (J/m2)'
      !
      CALL write_surf(dgu, u, &
                    hprogram,yrecfm,dgs%XSWDC(:),iresp,hcomment=ycomment)
      !
      yrecfm='SWUC_SEA'
      ycomment='X_Y_'//yrecfm//' (J/m2)'
      !
      CALL write_surf(dgu, u, &
                    hprogram,yrecfm,dgs%XSWUC(:),iresp,hcomment=ycomment)
      !
      yrecfm='LWDC_SEA'
      ycomment='X_Y_'//yrecfm//' (J/m2)'
      !
      CALL write_surf(dgu, u, &
                    hprogram,yrecfm,dgs%XLWDC(:),iresp,hcomment=ycomment)
      !
      yrecfm='LWUC_SEA'
      ycomment='X_Y_'//yrecfm//' (J/m2)'
      !
      CALL write_surf(dgu, u, &
                    hprogram,yrecfm,dgs%XLWUC(:),iresp,hcomment=ycomment)
      !
   ENDIF
   !
   yrecfm='FMUC_SEA'
   ycomment='X_Y_'//yrecfm//' (kg/ms)'
   !
   CALL write_surf(dgu, u, &
                    hprogram,yrecfm,dgs%XFMUC(:),iresp,hcomment=ycomment)
   !
   yrecfm='FMVC_SEA'
   ycomment='X_Y_'//yrecfm//' (kg/ms)'
   !
   CALL write_surf(dgu, u, &
                    hprogram,yrecfm,dgs%XFMVC(:),iresp,hcomment=ycomment)
   !
END IF
!
IF (dgs%LSURF_BUDGET.OR.dgs%LSURF_BUDGETC) THEN
!
  yrecfm='TALB_SEA'
  ycomment='total albedo over tile sea (-)'
  CALL write_surf(dgu, u, &
                    hprogram,yrecfm,dgs%XALBT(:),iresp,hcomment=ycomment)
!        
ENDIF
!
!*       4.     transfer coefficients
!               ---------------------
!
IF (dgs%LCOEF) THEN
   !
   yrecfm='CD_SEA'
   ycomment='X_Y_'//yrecfm//' (W/s2)'
   !
   CALL write_surf(dgu, u, &
                    hprogram,yrecfm,dgs%XCD(:),iresp,hcomment=ycomment)
   !
   yrecfm='CH_SEA'
   ycomment='X_Y_'//yrecfm//' (W/s)'
   !
   CALL write_surf(dgu, u, &
                    hprogram,yrecfm,dgs%XCH(:),iresp,hcomment=ycomment)
   !
   yrecfm='CE_SEA'
   ycomment='X_Y_'//yrecfm//' (W/s/K)'
   !
   CALL write_surf(dgu, u, &
                    hprogram,yrecfm,dgs%XCE(:),iresp,hcomment=ycomment)
   !
   yrecfm='Z0_SEA'
   ycomment='X_Y_'//yrecfm//' (M)'
   !
   CALL write_surf(dgu, u, &
                    hprogram,yrecfm,dgs%XZ0(:),iresp,hcomment=ycomment)
   !
   yrecfm='Z0H_SEA'
   ycomment='X_Y_'//yrecfm//' (M)'
   !
   CALL write_surf(dgu, u, &
                    hprogram,yrecfm,dgs%XZ0H(:),iresp,hcomment=ycomment)
   !
END IF
!
!
!*       5.     Surface humidity
!               ----------------
!
IF (dgs%LSURF_VARS) THEN
   !
   yrecfm='QS_SEA'
   ycomment='X_Y_'//yrecfm//' (KG/KG)'
   !
   CALL write_surf(dgu, u, &
                    hprogram,yrecfm,dgs%XQS(:),iresp,hcomment=ycomment)
   !
ENDIF
!
!
!*       6.     parameters at 2 and 10 meters :
!               -----------------------------
!
IF (dgs%N2M>=1) THEN
   !
   yrecfm='T2M_SEA'
   ycomment='X_Y_'//yrecfm//' (K)'
   !
   CALL write_surf(dgu, u, &
                    hprogram,yrecfm,dgs%XT2M(:),iresp,hcomment=ycomment)
   !
   yrecfm='T2MMIN_SEA'
   ycomment='X_Y_'//yrecfm//' (K)'
   !
   CALL write_surf(dgu, u, &
                    hprogram,yrecfm,dgs%XT2M_MIN(:),iresp,hcomment=ycomment)
   IF(greset)dgs%XT2M_MIN(:)=xundef
   !
   yrecfm='T2MMAX_SEA'
   ycomment='X_Y_'//yrecfm//' (K)'
   !
   CALL write_surf(dgu, u, &
                    hprogram,yrecfm,dgs%XT2M_MAX(:),iresp,hcomment=ycomment)
   IF(greset)dgs%XT2M_MAX(:)=0.0
   !
   yrecfm='Q2M_SEA'
   ycomment='X_Y_'//yrecfm//' (KG/KG)'
   !
   CALL write_surf(dgu, u, &
                    hprogram,yrecfm,dgs%XQ2M(:),iresp,hcomment=ycomment)
   !
   yrecfm='HU2M_SEA'
   ycomment='X_Y_'//yrecfm//' (-)'
   !
   CALL write_surf(dgu, u, &
                    hprogram,yrecfm,dgs%XHU2M(:),iresp,hcomment=ycomment)
   !
   yrecfm='HU2MMIN_SEA'
   ycomment='X_Y_'//yrecfm//' (-)'
   !
   CALL write_surf(dgu, u, &
                    hprogram,yrecfm,dgs%XHU2M_MIN(:),iresp,hcomment=ycomment)
   IF(greset)dgs%XHU2M_MIN(:)=xundef
   !
   yrecfm='HU2MMAX_SEA'
   ycomment='X_Y_'//yrecfm//' (-)'
   !
   CALL write_surf(dgu, u, &
                    hprogram,yrecfm,dgs%XHU2M_MAX(:),iresp,hcomment=ycomment)
   IF(greset)dgs%XHU2M_MAX(:)=-xundef
   !
   yrecfm='ZON10M_SEA'
   ycomment='X_Y_'//yrecfm//' (M/S)'
   !
   CALL write_surf(dgu, u, &
                    hprogram,yrecfm,dgs%XZON10M(:),iresp,hcomment=ycomment)
   !
   yrecfm='MER10M_SEA'
   ycomment='X_Y_'//yrecfm//' (M/S)'
   !
   CALL write_surf(dgu, u, &
                    hprogram,yrecfm,dgs%XMER10M(:),iresp,hcomment=ycomment)
   !
   yrecfm='W10M_SEA'
   ycomment='X_Y_'//yrecfm//' (M/S)'
   !
   CALL write_surf(dgu, u, &
                    hprogram,yrecfm,dgs%XWIND10M(:),iresp,hcomment=ycomment)
   !
   yrecfm='W10MMAX_SEA'
   ycomment='X_Y_'//yrecfm//' (M/S)'
   !
   CALL write_surf(dgu, u, &
                    hprogram,yrecfm,dgs%XWIND10M_MAX(:),iresp,hcomment=ycomment)
   IF(greset)dgs%XWIND10M_MAX(:)=0.0
   !
END IF
!
!
!*       7.     chemical diagnostics:
!               --------------------
!
IF (chs%SVS%NBEQ>0 .AND. chs%CCH_DRY_DEP=="WES89 ") THEN
   DO jsv = 1,SIZE(chs%CCH_NAMES,1)
      yrecfm='DV_SEA_'//trim(chs%CCH_NAMES(jsv))
      WRITE(ycomment,'(A13,I3.3)')'(m/s) DV_SEA_',jsv
      CALL write_surf(dgu, u, &
                    hprogram,yrecfm,chs%XDEP(:,jsv),iresp,hcomment=ycomment)
   END DO
ENDIF
!
!------------------------------------------------------------------------------
!
!         End of IO
!
 CALL end_io_surf_n(hprogram)
!
IF (lhook) CALL dr_hook('WRITE_DIAG_SEB_SEAFLUX_N',1,zhook_handle)
!
!
END SUBROUTINE write_diag_seb_seaflux_n