flag_diag_update.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 flag_diag_update (CHF, CHI, CHS, CHT, CHW, DGEI, DGF, DGI, DGMF, DGMI, DGMTO, &
                             dgo, dgs, dgsi, dgu, dgt, dgw, i, u, &
                             ofrac, odiag_grid, k2m, osurf_budget, orad_budget, ocoef,  &
                            osurf_vars, kbeq, kdsteq, odiag_ocean, odiag_seaice,       &
                            owater_profile,                                            &
                            osurf_evap_budget, oflood,  opgd_isba, och_no_flux_isba,   &
                            osurf_misc_budget_isba, opgd_teb, osurf_misc_budget_teb    )
!     ############################################################
!
!!****  *FLAG_DIAG_UPDATE* - routine to modify selection of output fields
!!
!!    PURPOSE
!!    -------
!!
!!**  METHOD
!!    ------
!!
!!    EXTERNAL
!!    --------
!!
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!
!!    AUTHOR
!!    ------
!!      P. Le Moigne   *Meteo France*   
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    02/2008 
!
!       B.Decharme  10/2009 flag to desactivate writing of pgd 
!!      Modified    04/2013, P. Le Moigne: FLake chemistry
!!      Modified    01/2014, S. Senesi   : introduce sea-ice model
!-------------------------------------------------------------------------------
!
!*       0.    DECLARATIONS
!              ------------
!
!
!
!
USE modd_ch_flake_n, ONLY : ch_flake_t
USE modd_ch_isba_n, ONLY : ch_isba_t
USE modd_ch_seaflux_n, ONLY : ch_seaflux_t
USE modd_ch_teb_n, ONLY : ch_teb_t
USE modd_ch_watflux_n, ONLY : ch_watflux_t
USE modd_diag_evap_isba_n, ONLY : diag_evap_isba_t
USE modd_diag_flake_n, ONLY : diag_flake_t
USE modd_diag_isba_n, ONLY : diag_isba_t
USE modd_diag_misc_flake_n, ONLY : diag_misc_flake_t
USE modd_diag_misc_isba_n, ONLY : diag_misc_isba_t
USE modd_diag_misc_teb_option_n, ONLY : diag_misc_teb_options_t
USE modd_diag_ocean_n, ONLY : diag_ocean_t
USE modd_diag_seaflux_n, ONLY : diag_seaflux_t
USE modd_diag_seaice_n, ONLY : diag_seaice_t
USE modd_diag_surf_atm_n, ONLY : diag_surf_atm_t
USE modd_diag_teb_n, ONLY : diag_teb_t
USE modd_diag_watflux_n, ONLY : diag_watflux_t
USE modd_isba_n, ONLY : isba_t
USE modd_surf_atm_n, ONLY : surf_atm_t
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*       0.1   Declarations of arguments
!              -------------------------
!
!
TYPE(ch_flake_t), INTENT(INOUT) :: chf
TYPE(ch_isba_t), INTENT(INOUT) :: chi
TYPE(ch_seaflux_t), INTENT(INOUT) :: chs
TYPE(ch_teb_t), INTENT(INOUT) :: cht
TYPE(ch_watflux_t), INTENT(INOUT) :: chw
TYPE(diag_evap_isba_t), INTENT(INOUT) :: dgei
TYPE(diag_flake_t), INTENT(INOUT) :: dgf
TYPE(diag_isba_t), INTENT(INOUT) :: dgi
TYPE(diag_misc_flake_t), INTENT(INOUT) :: dgmf
TYPE(diag_misc_isba_t), INTENT(INOUT) :: dgmi
TYPE(diag_misc_teb_options_t), INTENT(INOUT) :: dgmto
TYPE(diag_ocean_t), INTENT(INOUT) :: dgo
TYPE(diag_seaflux_t), INTENT(INOUT) :: dgs
TYPE(diag_seaice_t), INTENT(INOUT) :: dgsi
TYPE(diag_surf_atm_t), INTENT(INOUT) :: dgu
TYPE(diag_teb_t), INTENT(INOUT) :: dgt
TYPE(diag_watflux_t), INTENT(INOUT) :: dgw
TYPE(isba_t), INTENT(INOUT) :: i
TYPE(surf_atm_t), INTENT(INOUT) :: u
!
LOGICAL, INTENT(IN) :: ofrac
LOGICAL, INTENT(IN) :: odiag_grid
INTEGER, INTENT(IN) :: k2m
LOGICAL, INTENT(IN) :: osurf_budget
LOGICAL, INTENT(IN) :: orad_budget
LOGICAL, INTENT(IN) :: ocoef
LOGICAL, INTENT(IN) :: osurf_vars
!
INTEGER, INTENT(IN) :: kbeq
INTEGER, INTENT(IN) :: kdsteq
!
LOGICAL, INTENT(IN) :: odiag_ocean
LOGICAL, INTENT(IN) :: odiag_seaice
LOGICAL, INTENT(IN) :: owater_profile
LOGICAL, INTENT(IN) :: osurf_evap_budget
LOGICAL, INTENT(IN) :: oflood
LOGICAL, INTENT(IN) :: osurf_misc_budget_isba
LOGICAL, INTENT(IN) :: och_no_flux_isba
LOGICAL, INTENT(IN) :: opgd_isba
LOGICAL, INTENT(IN) :: osurf_misc_budget_teb
LOGICAL, INTENT(IN) :: opgd_teb
!
!*       0.2   Declarations of local variables
!              -------------------------------
!
REAL(KIND=JPRB) :: zhook_handle
!-------------------------------------------------------------------------------
IF (lhook) CALL dr_hook('FLAG_DIAG_UPDATE',0,zhook_handle)
!
dgu%LFRAC = ofrac
dgu%LDIAG_GRID = odiag_grid
dgu%N2M = k2m
dgu%LSURF_BUDGET = osurf_budget
dgu%LRAD_BUDGET = orad_budget
dgu%LCOEF = ocoef
dgu%LSURF_VARS = osurf_vars
!
dgu%LSURF_BUDGETC = .false.
IF (.NOT.dgu%LRESET_BUDGETC) dgu%LSURF_BUDGETC = .true.
!
IF (u%CSEA   =='SEAFLX') CALL flag_sea_update(kbeq, odiag_ocean, odiag_seaice, &
                                            dgu%LSURF_BUDGETC, dgu%LRAD_BUDGET, k2m)
IF (u%CWATER =='FLAKE ') CALL flag_flake_update(kbeq, owater_profile, dgu%LSURF_BUDGETC, dgu%LRAD_BUDGET, k2m)
IF (u%CWATER =='WATFLX') CALL flag_water_update(kbeq, dgu%LSURF_BUDGETC, dgu%LRAD_BUDGET, k2m)
IF (u%CNATURE=='ISBA  ') CALL flag_isba_update(kbeq, kdsteq, osurf_evap_budget, oflood, &
                                opgd_isba, och_no_flux_isba, osurf_misc_budget_isba,&
                                dgu%LSURF_BUDGETC, dgu%LRAD_BUDGET, k2m)
IF (u%CTOWN  =='TEB   ') CALL flag_teb_update(kbeq, opgd_teb, osurf_misc_budget_teb, dgu%LRAD_BUDGET, k2m)
!
IF (lhook) CALL dr_hook('FLAG_DIAG_UPDATE',1,zhook_handle)
!-------------------------------------------------------------------------------
 CONTAINS
!
SUBROUTINE flag_sea_update(KBEQ, ODIAG_OCEAN, ODIAG_SEAICE, OSURF_BUDGETC, ORAD_BUDGET, K2M)
!
!
IMPLICIT NONE
!
INTEGER, INTENT(IN) :: kbeq
LOGICAL, INTENT(IN) :: odiag_ocean
LOGICAL, INTENT(IN) :: odiag_seaice
LOGICAL, INTENT(IN) :: osurf_budgetc
LOGICAL, INTENT(IN) :: orad_budget
INTEGER, INTENT(IN) :: k2m
!
REAL(KIND=JPRB) :: zhook_handle
!
IF (lhook) CALL dr_hook('FLAG_DIAG_UPDATE:FLAG_SEA_UPDATE',0,zhook_handle)
 chs%SVS%NBEQ = kbeq
dgo%LDIAG_OCEAN = odiag_ocean
dgsi%LDIAG_SEAICE = odiag_seaice
dgs%LSURF_BUDGETC = osurf_budgetc
dgs%LRAD_BUDGET = orad_budget
dgs%N2M = k2m
IF (lhook) CALL dr_hook('FLAG_DIAG_UPDATE:FLAG_SEA_UPDATE',1,zhook_handle)
!
END SUBROUTINE flag_sea_update
!
SUBROUTINE flag_water_update(KBEQ, OSURF_BUDGETC, ORAD_BUDGET, K2M)
!
!
IMPLICIT NONE
!
INTEGER, INTENT(IN) :: kbeq
LOGICAL, INTENT(IN) :: osurf_budgetc
LOGICAL, INTENT(IN) :: orad_budget
INTEGER, INTENT(IN) :: k2m
!
REAL(KIND=JPRB) :: zhook_handle
!
IF (lhook) CALL dr_hook('FLAG_DIAG_UPDATE:FLAG_WATER_UPDATE',0,zhook_handle)
 chw%SVW%NBEQ = kbeq
dgw%LSURF_BUDGETC = osurf_budgetc
dgw%LRAD_BUDGET = orad_budget
dgw%N2M = k2m
IF (lhook) CALL dr_hook('FLAG_DIAG_UPDATE:FLAG_WATER_UPDATE',1,zhook_handle)
!
END SUBROUTINE flag_water_update
!
SUBROUTINE flag_flake_update(KBEQ, OWATER_PROFILE, OSURF_BUDGETC, ORAD_BUDGET, K2M)
!
!
IMPLICIT NONE
!
INTEGER, INTENT(IN) :: kbeq
LOGICAL, INTENT(IN) :: owater_profile
LOGICAL, INTENT(IN) :: osurf_budgetc
LOGICAL, INTENT(IN) :: orad_budget
INTEGER, INTENT(IN) :: k2m
!
REAL(KIND=JPRB) :: zhook_handle
!
IF (lhook) CALL dr_hook('FLAG_DIAG_UPDATE:FLAG_FLAKE_UPDATE',0,zhook_handle)
 chf%SVF%NBEQ = kbeq
dgmf%LWATER_PROFILE = owater_profile
dgf%LSURF_BUDGETC = osurf_budgetc
dgf%LRAD_BUDGET = orad_budget
dgf%N2M = k2m
IF (lhook) CALL dr_hook('FLAG_DIAG_UPDATE:FLAG_FLAKE_UPDATE',1,zhook_handle)
!
END SUBROUTINE flag_flake_update
!
!
SUBROUTINE flag_isba_update(KBEQ, KDSTEQ, OSURF_EVAP_BUDGET, OFLOOD, &
                            opgd, och_no_flux, osurf_misc_budget, &
                            osurf_budgetc, orad_budget, k2m)
!
!
IMPLICIT NONE
!
INTEGER, INTENT(IN) :: kbeq
INTEGER, INTENT(IN) :: kdsteq
LOGICAL, INTENT(IN) :: osurf_evap_budget
LOGICAL, INTENT(IN) :: oflood
LOGICAL, INTENT(IN) :: opgd
LOGICAL, INTENT(IN) :: och_no_flux
LOGICAL, INTENT(IN) :: osurf_misc_budget
LOGICAL, INTENT(IN) :: osurf_budgetc
LOGICAL, INTENT(IN) :: orad_budget
INTEGER, INTENT(IN) :: k2m
!
REAL(KIND=JPRB) :: zhook_handle
!
IF (lhook) CALL dr_hook('FLAG_DIAG_UPDATE:FLAG_ISBA_UPDATE',0,zhook_handle)
 chi%SVI%NBEQ = kbeq
 chi%SVI%NDSTEQ = kdsteq
dgei%LSURF_EVAP_BUDGET = osurf_evap_budget
i%LFLOOD = oflood
dgi%LPGD = opgd
 chi%LCH_NO_FLUX = och_no_flux
dgmi%LSURF_MISC_BUDGET = osurf_misc_budget
dgei%LSURF_BUDGETC = osurf_budgetc
dgi%LPATCH_BUDGET = osurf_budgetc
dgi%LRAD_BUDGET = orad_budget
dgi%N2M = k2m
IF (lhook) CALL dr_hook('FLAG_DIAG_UPDATE:FLAG_ISBA_UPDATE',1,zhook_handle)
!
END SUBROUTINE flag_isba_update
!
SUBROUTINE flag_teb_update(KBEQ, OPGD, OSURF_MISC_BUDGET, ORAD_BUDGET, K2M)
!
!
IMPLICIT NONE
!
INTEGER, INTENT(IN) :: kbeq
LOGICAL, INTENT(IN) :: opgd
LOGICAL, INTENT(IN) :: osurf_misc_budget
LOGICAL, INTENT(IN) :: orad_budget
INTEGER, INTENT(IN) :: k2m
!
REAL(KIND=JPRB) :: zhook_handle
!
IF (lhook) CALL dr_hook('FLAG_DIAG_UPDATE:FLAG_TEB_UPDATE',0,zhook_handle)
 cht%SVT%NBEQ = kbeq
dgt%LPGD = opgd
dgmto%LSURF_MISC_BUDGET = osurf_misc_budget
dgt%LRAD_BUDGET = orad_budget
dgt%N2M = k2m
IF (lhook) CALL dr_hook('FLAG_DIAG_UPDATE:FLAG_TEB_UPDATE',1,zhook_handle)
!
END SUBROUTINE flag_teb_update
!
END SUBROUTINE flag_diag_update