source: LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/suinit.F90 @ 5441

SUBROUTINE SUINIT(klon,klev)
#ifdef DOC

!     **** *SUINIT* - SCM initialization.

!     Purpose.
!     --------

!     **   Interface.
!     ----------

!     Explicit arguments :    None.
!     --------------------

!     Implicit arguments :    None.
!     --------------------

!     Method.
!     -------

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

!     Reference.
!     ----------

!     Author.
!     -------
!     Eric Bazile, Francois Bouyssel et Jean-Marcel Piriou

!     Modifications.
!     --------------
!     Original :97-02-01
!               Jozef Vivoda, SHMI: calling sequence as in 3D model
!                                   and ECMWF setup
!     2001-11-27 P. Marquet : several printout on listing (NULOUT=15)

!     ------------------------------------------------------------------
#endif

USE PARKIND1  ,ONLY : JPIM     ,JPRB
!#include "tsmbkind.h"

USE PARDIM, ONLY : JPMXLE
USE YOMCT0B  , ONLY : LECMWF
USE YOMRIP   , ONLY : NINDAT   ,NSSSSS
USE YOMDIM  
USE YOMDPHY  
! MPL 29042010: NDLNPR,RHYDR0 non initialises et pour ne pas mettre tout sudyn.F90
USE YOMDYN  , ONLY : TSTEP , NDLNPR , RHYDR0        ! MPL 29042010
!USE YOMEVOL  , ONLY : TECH     ,FREQFS   ,FREQFE   , FREQDDH
!USE YOMCT0   , ONLY : LFROG
! quelques ajouts qui viennent de suallo
USE YOMGEM   , ONLY : VDELA    , VDELB   ,VC       ,NLOEN    ,NLOENG  ,NGPTOT
USE YOMSTA   , ONLY : STZ      ,STPREH   ,STPRE    ,STPHI    ,STTEM   ,STDEN
USE YOEAERD  , ONLY : CVDAES   ,CVDAEL   ,CVDAEU   ,CVDAED
USE YOEOVLP  , ONLY : RA1OVLP
USE YOECLD   , ONLY : CETA
USE YOECND   , ONLY : CEVAPCU
USE YOMTOPH  , ONLY : RMESOU   ,RMESOT   ,RMESOQ
USE YOMGC    , ONLY : GEMU     ,GELAM    ,GELAT    ,GECLO    ,GESLO    ,GM       ,GAW


IMPLICIT NONE
LOGICAL LLTRACE, LLDEBUG
integer klon,klev
CHARACTER*200 CFICP
CHARACTER*200 CFLUX
CHARACTER*200 CLIST
CHARACTER*200 CFDDH
CHARACTER*80 CNMEXP


LLTRACE=.TRUE.
LLDEBUG=.TRUE.

!     ------------------------
!     *    READ NAMELISTS.
!     ------------------------

!----------------------------------------------------------------
! Elements indispensables de SUNAM pour faire tourner RRTM dans LMDZ
!-------------------------------------------------------------------
CFICP='Profile'
CFLUX='Output'
CLIST='Listing'
CFDDH='DHFDL'
CNMEXP='SCM'
TSTEP=450
! MPL 29042010 - RHYDR0 - upper boundary contition for hydrostatic
RHYDR0=1._JPRB
! MPL 29042010
! NDLNPR : NDLNPR=0: conventional formulation of delta, i.e. ln(P(l)/P(l-1)).
!          NDLNPR=1: formulation of delta used in non hydrostatic model,
NDLNPR=0
PRINT *,'SUINIT: RHYDR0 NDLNPR',RHYDR0,NDLNPR

!----------------------------------------------------------------
! Elements indispensables de SUDIM pour faire tourner RRTM dans LMDZ
!-------------------------------------------------------------------
NDLON=klon
NFLEVG=klev
NPROMA=klon

!-------------------------------------------------------------------
!JV    Initialize constants
!     ---------------------
!JV
IF (LLTRACE)  WRITE(*,*) " coucou SUINIT : avant SUCST"
WRITE(*,FMT='('' ---------------- '')')
WRITE(*,FMT='(''     SUCST : '')')
WRITE(*,FMT='('' ---------------- '')')
NINDAT=20090408      !!!!! A REVOIR (MPL)
NSSSSS=0  ! LMDZ demarre tjrs a 00h -- MPL 15.04.09
CALL SUCST(6,NINDAT,NSSSSS,1)
PRINT *,'SUINIT: NINDAT, NSSSSS',NINDAT, NSSSSS

IF (LLDEBUG) THEN
WRITE(*,FMT='(''  SUINIT / apres : SUCST '')')
ENDIF


!     ------------------------
!     *    ALLOCATES RECUPERES DE SUALLO
!     ------------------------
ALLOCATE(VDELA  (MAX(JPMXLE,NFLEVG)))
ALLOCATE(VDELB  (MAX(JPMXLE,NFLEVG)))
VDELB = 0  !ym missing init
ALLOCATE( VC      (NFLEVG) )
VC = 0    !ym missing init
ALLOCATE( NLOEN   (NPROMA) )
ALLOCATE( NLOENG   (NPROMA) )
ALLOCATE( STZ     (NFLEVG) )
ALLOCATE( CVDAES  (NFLEVG+1))
ALLOCATE( CVDAEL  (NFLEVG+1))
ALLOCATE( CVDAEU  (NFLEVG+1))
ALLOCATE( CVDAED  (NFLEVG+1))
ALLOCATE(RA1OVLP(NFLEVG))

ALLOCATE(STPREH(0:NFLEVG)) ! Nouvel ajout MPL 22062010
ALLOCATE(STPRE(NFLEVG))
ALLOCATE(STPHI(NFLEVG))
ALLOCATE(STTEM(NFLEVG))
ALLOCATE(STDEN(NFLEVG))

ALLOCATE(CETA(NFLEVG))    ! Nouvel ajout MPL 28062010
ALLOCATE(CEVAPCU(NFLEVG))
ALLOCATE(RMESOU(NFLEVG))
ALLOCATE(RMESOT(NFLEVG))
ALLOCATE(RMESOQ(NFLEVG))

!     ------------------------
!     *    ALLOCATES RECUPERES DE SUGEM2
!     ------------------------

ALLOCATE(GEMU   (NGPTOT)) ! Nouvel ajout MPL 28062010
ALLOCATE(GELAM  (NGPTOT))
ALLOCATE(GELAT  (NGPTOT))
ALLOCATE(GECLO  (NGPTOT))
ALLOCATE(GESLO  (NGPTOT))
ALLOCATE(GM     (NGPTOT))
ALLOCATE(GAW    (NGPTOT))
!  
!     ------------------------------------------------------------------

END SUBROUTINE SUINIT