ini_rain_ice.F90

!MNH_LIC Copyright 1995-2019 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
!-----------------------------------------------------------------
!     ######spl
      SUBROUTINE INI_RAIN_ICE ( KLUOUT, HCLOUD )
!     ###########################################################
!
!!****  *INI_RAIN_ICE * - initialize the constants necessary for the warm and
!!                        cold microphysical schemes.
!!
!!    PURPOSE
!!    -------
!!      The purpose of this routine is to initialize the constants used to
!!    resolve the mixed phase microphysical scheme. The collection kernels of
!!    the precipitating particles are recomputed if necessary if some parameters
!!    defining the ice categories have been modified. The number of small
!!    time steps leading to stable scheme for the rain, ice, snow and ggraupeln
!!    sedimentation is also computed (time-splitting technique).
!!
!!**  METHOD
!!    ------
!!      The constants are initialized to their numerical values and the number
!!    of small time step is computed by dividing the 2* Deltat time interval of
!!    the Leap-frog scheme so that the stability criterion for the rain
!!    sedimentation is fulfilled for a Raindrop maximal fall velocity equal
!!    VTRMAX. The parameters defining the collection kernels are read and are
!!    checked against the new ones. If any change occurs, these kernels are
!!    recomputed and their numerical values are written in the output listing.
!!
!!    EXTERNAL
!!    --------
!!      GAMMA    :  gamma function
!!
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!      Module MODD_CST
!!        XPI                  !
!!        XP00                 ! Reference pressure
!!        XRD                  ! Gaz constant for dry air
!!        XRHOLW               ! Liquid water density
!!      Module MODD_REF
!!        XTHVREFZ             ! Reference virtual pot.temp. without orography
!!      Module MODD_PARAMETERS
!!        JPVEXT               !
!!      Module MODD_RAIN_ICE_DESCR
!!      Module MODD_RAIN_ICE_PARAM
!!
!!    REFERENCE
!!    ---------
!!      Book2 of documentation ( routine INI_RAIN_ICE )
!!
!!    AUTHOR
!!    ------
!!      J.-P. Pinty      * Laboratoire d'Aerologie*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    04/12/95
!!      J.-P. Pinty 05/04/96 Add automatic control and regeneration of the
!!                           collection kernels
!!      J.-P. Pinty 10/05/96 Correction of ZRATE and computations of RIM
!!      J.-P. Pinty 24/11/97 Sedimentation of ice made for Columns and bug for XAG
!!      J.-P. Lafore 23/11/98 Back to Lin et al. 83 formulation for RIAUTS
!!                            with a Critical ice content set to .5 g/Kg
!!      N. Asencio  13/08/98 parallel code: PDZMIN is computed outside in ini_modeln
!!      J.-P. Lafore 12/8/98 In case of nesting microphysics constants of
!!                           MODD_RAIN_ICE_PARAM are computed only once.
!!                           Only KSPLTR is computed for each model.
!!      J. Stein    20/04/99 remove 2 unused local variables
!!      G Molinie   21/05/99 Bug in XEXRCFRI and XRCFRI
!!      J.-P. Pinty 24/06/00 Bug in RCRIMS
!!      J.-P. Pinty 24/12/00 Update hail case
!!      J.-P. Chaboureau & J.-P. Pinty
!!                  24/03/01 Update XCRIAUTI for cirrus cases
!!      J.-P. Pinty 24/11/01 Update ICE3/ICE4 options
!!      S. Riette 2016-11: new ICE3/ICE4 options
!!      P. Wautelet 22/01/2019 bug correction: incorrect write
!  P. Wautelet 26/04/2019: replace non-standard FLOAT function by REAL function
!
!-------------------------------------------------------------------------------
!
!*       0.    DECLARATIONS
!              ------------
!
USE MODD_CST
USE MODD_LUNIT
USE MODD_PARAMETERS
USE MODD_PARAM_ICE
USE MODD_RAIN_ICE_DESCR
USE MODD_RAIN_ICE_PARAM
USE MODD_REF
!
USE MODI_GAMMA
USE MODI_GAMMA_INC
USE MODE_RRCOLSS, ONLY: RRCOLSS
USE MODE_RZCOLX, ONLY: RZCOLX
USE MODE_RSCOLRG, ONLY: RSCOLRG
USE MODE_READ_XKER_RACCS, ONLY: READ_XKER_RACCS
USE MODE_READ_XKER_SDRYG, ONLY: READ_XKER_SDRYG
USE MODE_READ_XKER_RDRYG, ONLY: READ_XKER_RDRYG
USE MODE_READ_XKER_SWETH, ONLY: READ_XKER_SWETH
USE MODE_READ_XKER_GWETH, ONLY: READ_XKER_GWETH
USE MODE_READ_XKER_RWETH, ONLY: READ_XKER_RWETH
!
USE PARKIND1, ONLY : JPRB
USE YOMHOOK , ONLY : LHOOK, DR_HOOK
!
IMPLICIT NONE
!
!*       0.1   Declarations of dummy arguments :
!
!
INTEGER,                 INTENT(IN) :: KLUOUT   ! Logical unit number for prints
!
CHARACTER (LEN=4), INTENT(IN)       :: HCLOUD    ! Indicator of the cloud scheme
!
!
!
!*       0.2   Declarations of local variables :
!
INTEGER :: IKB                ! Coordinates of the first physical
                              ! points along z
INTEGER :: J1,J2              ! Internal loop indexes
REAL :: ZT                    ! Work variable
REAL :: ZVTRMAX               ! Raindrop maximal fall velocity
REAL :: ZRHO00                ! Surface reference air density
REAL :: ZE, ZRV               ! Work array for ZRHO00 computation
REAL :: ZRATE                 ! Geometrical growth of Lbda in the tabulated
                              ! functions and kernels
REAL :: ZBOUND                ! XDCSLIM*Lbda_s: upper bound for the partial
                              ! integration of the riming rate of the aggregates
REAL :: ZEGS, ZEGR, ZEHS, &   ! Bulk collection efficiencies
      & ZEHG, ZEHR
!
INTEGER :: IND                ! Number of interval to integrate the kernels
REAL :: ZESR                  ! Mean efficiency of rain-aggregate collection
REAL :: ZFDINFTY              ! Factor used to define the "infinite" diameter
!
!
!
LOGICAL  :: GFLAG   ! Logical flag for printing the constatnts on the output
                    ! listing
REAL     :: ZCONC_MAX ! Maximal concentration for snow
REAL     :: ZGAMC,ZGAMC2 ! parameters
                    ! involving various moments of the generalized gamma law
REAL     :: ZFACT_NUCL! Amplification factor for the minimal ice concentration
REAL     :: ZXR     ! Value of x_r in N_r = C_r lambda_r ** x_r
!
INTEGER  :: KND
INTEGER  :: KACCLBDAS,KACCLBDAR,KDRYLBDAG,KDRYLBDAS,KDRYLBDAR
INTEGER  :: KWETLBDAS,KWETLBDAG,KWETLBDAR,KWETLBDAH
REAL     :: PALPHAR,PALPHAS,PALPHAG,PALPHAH
REAL     :: PNUR,PNUS,PNUG,PNUH
REAL     :: PBR,PBS,PBG
REAL     :: PCR,PCS,PCG,PCH
REAL     :: PDR,PDS,PDG,PDH
REAL     :: PESR,PEGS,PEGR,PEHS,PEHG,PEHR
REAL     :: PFDINFTY
REAL     :: PACCLBDAS_MAX,PACCLBDAR_MAX,PACCLBDAS_MIN,PACCLBDAR_MIN
REAL     :: PDRYLBDAG_MAX,PDRYLBDAS_MAX,PDRYLBDAG_MIN,PDRYLBDAS_MIN
REAL     :: PDRYLBDAR_MAX,PDRYLBDAR_MIN
REAL     :: PWETLBDAS_MAX,PWETLBDAG_MAX,PWETLBDAS_MIN,PWETLBDAG_MIN
REAL     :: PWETLBDAR_MAX,PWETLBDAH_MAX,PWETLBDAR_MIN,PWETLBDAH_MIN
!
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!-------------------------------------------------------------------------------
!
IF (LHOOK) CALL DR_HOOK('INI_RAIN_ICE',0,ZHOOK_HANDLE)
!
!*       0.     FUNCTION STATEMENTS
!               -------------------
!
!
!*       0.1    p_moment of the Generalized GAMMA function
!
!
!
!
IF (ALLOCATED(XRTMIN)) THEN       ! In case of nesting microphysics constants of
                                  ! MODD_RAIN_ICE_PARAM are computed only once,
                                  ! but if INI_RAIN_ICE has been called already
                                  ! one must change the XRTMIN size.
  DEALLOCATE(XRTMIN)
END IF
!
!-------------------------------------------------------------------------------
!
!*       2.     CHARACTERISTICS OF THE SPECIES
!               ------------------------------
!
!
!*       2.1    Cloud droplet and Raindrop characteristics
!
XAC = (XPI/6.0)*XRHOLW
XBC = 3.0
XCC = XRHOLW*XG/(18.0*1.7E-5) ! Stokes flow (Pruppacher p 322 for T=273K)
XDC = 2.0