dslt_init_modes.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 dslt_init_modes (KEQ, KSV_BEG, KSV_END, OVARSIG, ORGFIX, &
                            kmdebeg, kmde)
!!    ###########################################
!!
!!*** *DSLT_INIT_MODES*
!!
!!    PURPOSE
!!    -------
!!    Find the number of dust modes to be transported
!!    Each mode needs 3 moments to be described, so logically, the number of modes is
!!    The number of dust tracers divided by 3
!!     
!!
!!    REFERENCE
!!    ---------
!!    Modified dst_init_names (march 2005)    
!!
!!    AUTHOR
!!    ------
!!    Alf Grini <alf.grini@cnrm.meteo.fr>
!!
!!    MODIFICATIONS
!!    -------------
!!
!!    EXTERNAL
!!    --------
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!-------------------------------------------------------------------------------
!
!*       0.     DECLARATIONS
!               ------------
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
USE modi_abor1_sfx
!
IMPLICIT NONE
!
!*      0.1    declarations of arguments
!

INTEGER,                         INTENT(IN) :: keq       ! number of dust variables
INTEGER,                         INTENT(IN) :: ksv_beg   ! First number of dust tracer
INTEGER,                         INTENT(IN) :: ksv_end   ! Last number of dust tracer
LOGICAL,                         INTENT(IN) :: ovarsig   ! type of standard deviation (fixed or variable)
LOGICAL,                         INTENT(IN) :: orgfix    ! type of mean radius
INTEGER,                         INTENT(OUT) :: kmdebeg  ! Place in scalar list of dustmass in first mode
INTEGER,                         INTENT(OUT) :: kmde     ! Number of dust modes
REAL(KIND=JPRB) :: zhook_handle


!Check if you have a multiple of 3 dust related variables, and 
!Set the number of modes to the number of dust related variables
!divided by 3
IF (lhook) CALL dr_hook('DSLT_INIT_MODES',0,zhook_handle)
!
kmdebeg = ksv_beg
kmde    = ksv_end - ksv_beg + 1
!
IF (ovarsig) THEN !case three moments by modes
  IF(mod(kmde,3).NE.0.) THEN
    CALL abor1_sfx('DST_INIT_MODES: (1) WRONG NUMBER OF DUST VARIABLES')
  ELSE
    kmde = kmde / 3
  ENDIF
ELSE IF (.NOT.orgfix) THEN ! case two moment by modes
  IF(mod(kmde,2).ne.0.)THEN
   CALL abor1_sfx('DST_INIT_MODES: (1) WRONG NUMBER OF DUST VARIABLES')
  ELSE
    kmde = kmde / 2
  END IF
END IF
!
IF (lhook) CALL dr_hook('DSLT_INIT_MODES',1,zhook_handle)
!
END SUBROUTINE dslt_init_modes