suphy2.F90

Go to the documentation of this file.

!OPTIONS XOPT(NOEVAL)
SUBROUTINE suphy2(KULOUT)

!**** *SUPHY2*   - Initialize common YOMPHY2 physics controlling
!                  constants

!     Purpose.
!     --------
!           Initialize YOMPHY2, the common that contains the parameters
!           for the control part of the physics of the model.

!**   Interface.
!     ----------
!        *CALL* *SUPHY2(KULOUT)

!        Explicit arguments :
!        --------------------
!        KULOUT : Logical unit for the output

!        Implicit arguments :
!        --------------------
!        COMMON YOMPHY2

!     Method.
!     -------
!        See documentation

!     Externals.
!     ----------

!     Reference.
!     ----------
!        Documentation ARPEGE

!     Author.
!     -------
!        J.-F. Geleyn .
!        Original : 90-9-1

!     Modifications.
!     --------------
!        R. EL Khatib : 93-04-02 Set-up defaults controled by LECMWF
!        J.-F. Geleyn : 93-08-19 New cloudiness diagnostics.
!        J.-F. Geleyn : 95-04-10 Anti-fibril. Girard-Delage.
!        P. Marquet   : 97-02-18 Value of VETAF=VAH/VP00+VBH.
!        J.M. Piriou  : 97-04-17 XMULAF default value.
!        E. Bazile    : 98-03-10 Introduce XMUCVPP.
!        W. Owcarz    : 2000-03-27 Set a default value for TSPHY
!        R. EL Khatib : 2000-06-13 RIPBLC
!        R. EL Khatib : 2000-08-21 Turbulent gusts setup
!        J.M. Piriou  : 2002-01-10 set default values to operational ones.
!        Modified by R. EL Khatib : 02-03-29 Control XMULAF<0 ; add LMULAF
!        Modified by D. Banciu    : 02-12-09 Introduction of XDAMP
!        M.Hamrud      01-Oct-2003 CY28 Cleaning
!     ------------------------------------------------------------------

USE parkind1  ,ONLY : jpim     ,jprb
USE yomhook   ,ONLY : lhook,   dr_hook

! Ce qui concerne NULNAM commente par MPL le 15.04.09
!USE YOMLUN   , ONLY : NULNAM
USE yomct0b  , ONLY : lecmwf
! commente par MPL 25.11.08
!USE YOMGEM   , ONLY : VALH     ,VBH
USE yomdim   , ONLY : nflevg
USE yomphy2  , ONLY : ntshm    ,ntsml    ,xmucvpp  ,lmulaf   ,&
 & xmulaf   ,xdamp    ,hclp     ,htcls    ,&
 & ripblc   ,&
 & lraftur  ,gz0raf   ,facraf   ,&
 & hvcls    ,htshm    ,htsml    ,&
 & tsphy  

IMPLICIT NONE

INTEGER(KIND=JPIM),INTENT(IN)    :: KULOUT 
INTEGER(KIND=JPIM) :: JLEV

REAL(KIND=JPRB) :: ZVETAF
REAL(KIND=JPRB) :: ZHOOK_HANDLE

#include "abor1.intfb.h"
#include "posnam.intfb.h"

#include "namphy2.h"
!     ------------------------------------------------------------------

!*       1.    Set default values.
!              -------------------

!        1.1 Set implicit default values

IF (lhook) CALL dr_hook('SUPHY2',0,zhook_handle)
xmulaf=-1.75_jprb
xmucvpp=0._jprb
xdamp=0._jprb
hclp=1500._jprb
htcls=2._jprb
hvcls=10._jprb
htshm=0.450_jprb
htsml=0.785_jprb
tsphy=1._jprb
ripblc=0.5_jprb
lraftur=.false.
gz0raf=10.0_jprb
facraf=15.0_jprb
lmulaf=.false.

!        1.2 Modify default values according to LECMWF

IF (lecmwf) THEN
ELSE
  lraftur=.true.
ENDIF

!     Remark : values for TSPHY, NTSHM/ML are calculated and not set up.

!     ------------------------------------------------------------------

!*       2.    Modify default values.
!              ----------------------

! Ce qui concerne NAMPHY2 commente par MPL le 15.04.09
!CALL POSNAM(NULNAM,'NAMPHY2')
!READ(NULNAM,NAMPHY2)
!     ------------------------------------------------------------------

!*       3.    Compute cloud transition indexes.
!              ---------------------------------

ntshm=0
ntsml=0
! commente par MPL 25.11.08
!DO JLEV=1,NFLEVG
!  ZVETAF=(VALH(JLEV)+VBH(JLEV)+VALH(JLEV-1)+VBH(JLEV-1))*0.5_JPRB
!  IF (ZVETAF <= HTSHM) THEN
!    NTSHM=JLEV
!  ENDIF
!  IF (ZVETAF <= HTSML) THEN
!    NTSML=JLEV
!  ENDIF
!ENDDO

!     ------------------------------------------------------------------

!*       4.    Print final values.
!              -------------------

WRITE(unit=kulout,fmt='('' COMMON YOMPHY2 '')')
WRITE(unit=kulout,fmt='('' XMUCVPP = '',E10.4,'' XMULAF = '',E10.4 &
 & ,'' XDAMP = '',E10.4 &
 & ,'' LMULAF = '',L2,/,'' HTCLS = '',E10.4 &
 & ,'' HVCLS = '',E10.4,'' HCLP = '',E10.4,/&
 & ,'' RIPBLC  = '',F8.4 &
 & ,'' LRAFTUR = '',L2,'' GZ0RAF = '',E10.4,'' FACRAF = '',E10.4 &
 & ,'' HTSHM = '',F8.4,'' NTSHM = '',I3,'' HTSML = '',F8.4 &
 & ,'' NTSML = '',I3 &
 & )')&
 & xmucvpp,xmulaf,xdamp,lmulaf,&
 & htcls,hvcls,hclp,&
 & ripblc,&
 & lraftur,gz0raf,facraf,&
 & htshm,ntshm,htsml,ntsml  

!*       5.    Control
!              -------

IF (xmulaf > 0.0_jprb) THEN
  WRITE(kulout,*) 'XMULAF SHOULD BE NEGATIVE'
  CALL abor1('SUPHY2 : ABOR1 CALLED')
ENDIF

IF ((xdamp /= 0.0_jprb).AND.(xmucvpp /= 0.0_jprb)) THEN
  WRITE(unit=kulout,fmt='(A)') 'INCONSISTENCY BETWEEN XDAMP AND XMUCVPP !'
  CALL abor1('XDAMP/=0. IMPLIES XMUCVPP=0.!...')
ENDIF

!     ------------------------------------------------------------------

IF (lhook) CALL dr_hook('SUPHY2',1,zhook_handle)
END SUBROUTINE suphy2