suinit.F90

Go to the documentation of this file.

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)))
ALLOCATE( vc(nflevg) )
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