diag_seaice_initn.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 diag_seaice_init_n (DGO, DI, DIC, DGMSI, OREAD_BUDGETC, S, &
                                      HPROGRAM,KLU,KSW)
!     #####################
!
!!****  *DIAG_SEAFLUX_INIT_n* - routine to initialize SEAFLUX diagnostic variables
!!
!!    PURPOSE
!!    -------
!!
!!**  METHOD
!!    ------
!!
!!    EXTERNAL
!!    --------
!!
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!
!!    AUTHOR
!!    ------
!!      V. Masson   *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    01/2004 
!!      Modified    01/2006 : sea flux parameterization.
!!      Modified    08/2009 : cumulative sea flux 
!!      B. decharme 04/2013 : Add EVAP and SUBL diag
!!      S.Senesi    01/2014 : introduce fractional seaice 
!-------------------------------------------------------------------------------
!
!*       0.    DECLARATIONS
!              ------------
!
USE mode_diag
!
USE modd_diag_n, ONLY : diag_t, diag_options_t
USE modd_diag_misc_seaice_n, ONLY : diag_misc_seaice_t
USE modd_seaflux_n, ONLY : seaflux_t
!
#ifdef SFX_OL
USE modn_io_offline,     ONLY : lrestart
#endif
USE modd_surf_par,       ONLY : xundef
USE modd_sfx_oasis,      ONLY : lcpl_sea,lcpl_seaice
!
USE modi_read_surf
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*       0.1   Declarations of arguments
!              -------------------------
!
TYPE(diag_options_t), INTENT(INOUT) :: DGO
TYPE(diag_t), INTENT(INOUT) :: DI
TYPE(diag_t), INTENT(INOUT) :: DIC
TYPE(diag_misc_seaice_t), INTENT(INOUT) :: DGMSI
LOGICAL, INTENT(IN) :: OREAD_BUDGETC
TYPE(seaflux_t), INTENT(INOUT) :: S
!
INTEGER, INTENT(IN) :: KLU   ! size of arrays
INTEGER, INTENT(IN) :: KSW   ! number of SW spectral bands
 CHARACTER(LEN=6), INTENT(IN):: HPROGRAM  ! program calling
!
!*       0.2   Declarations of local variables
!              -------------------------------
!
INTEGER           :: IVERSION
INTEGER           :: IRESP          ! IRESP  : return-code if a problem appears
 CHARACTER(LEN=12) :: YREC           ! Name of the article to be read
!
REAL(KIND=JPRB)   :: ZHOOK_HANDLE
!
!-------------------------------------------------------------------------------
!
!* surface energy budget
!
IF (lhook) CALL dr_hook('DIAG_SEAICE_INIT_N',0,zhook_handle)
!
 CALL alloc_bud(dgo,di,klu,ksw)
!
!* cumulative surface energy budget
!
#ifdef SFX_OL
IF (dgo%LSURF_BUDGETC .OR. (lrestart .AND. .NOT.dgo%LRESET_BUDGETC)) THEN
#else
IF (dgo%LSURF_BUDGETC .OR. .NOT.dgo%LRESET_BUDGETC) THEN
#endif
  !
  CALL alloc_surf_bud(dic,klu,klu,ksw)
  CALL init_surf_bud(dic,0.)
  !
ELSE
  CALL alloc_surf_bud(dic,0,0,0)
ENDIF
!
!* Seaice model diagnostics init 
!
IF (dgmsi%LDIAG_MISC_SEAICE) THEN
  ALLOCATE(dgmsi%XSIT(klu))
  ALLOCATE(dgmsi%XSND(klu)) 
  ALLOCATE(dgmsi%XMLT(klu))  
  dgmsi%XSIT=xundef
  dgmsi%XSND=xundef
  dgmsi%XMLT=xundef
ELSE
  ALLOCATE(dgmsi%XSIT  (0))
  ALLOCATE(dgmsi%XSND  (0))
  ALLOCATE(dgmsi%XMLT  (0))
ENDIF
!
IF(lcpl_seaice.OR.s%LHANDLE_SIC)THEN
  ALLOCATE(s%XCPL_SEAICE_SNET(klu))
  ALLOCATE(s%XCPL_SEAICE_HEAT(klu))
  ALLOCATE(s%XCPL_SEAICE_EVAP(klu))
  s%XCPL_SEAICE_SNET(:) = 0.0
  s%XCPL_SEAICE_HEAT(:) = 0.0
  s%XCPL_SEAICE_EVAP(:) = 0.0
ELSE
  ALLOCATE(s%XCPL_SEAICE_SNET(0))
  ALLOCATE(s%XCPL_SEAICE_HEAT(0))
  ALLOCATE(s%XCPL_SEAICE_EVAP(0))
ENDIF
!
IF (lhook) CALL dr_hook('DIAG_SEAICE_INIT_N',1,zhook_handle)
!
!-------------------------------------------------------------------------------
!
END SUBROUTINE diag_seaice_init_n