source: LMDZ5/trunk/libf/phylmd/rrtm/suecrad15.F90 @ 4521

!OPTIONS XOPT(NOEVAL)
SUBROUTINE SUECRAD15 (KULOUT, KLEV, PETAH )

!**** *SUECRAD15*   - INITIALIZE COMMONS YOMRxx15 CONTROLLING RADIATION
!****                 FROZEN VERSION (CYCLE 15) OF SUECRAD

!     PURPOSE.
!     --------
!           INITIALIZE YOMRAD15, THE COMMON THAT CONTROLS THE
!           RADIATION OF THE MODEL, AND YOMRDU15 THAT INCLUDES
!           ADJUSTABLE PARAMETERS FOR RADIATION COMPUTATIONS

!**   INTERFACE.
!     ----------
!        CALL *SUECRAD15* FROM *SUPHEC*
!              ---------        ------

!        EXPLICIT ARGUMENTS :
!        --------------------
!        NONE

!        IMPLICIT ARGUMENTS :
!        --------------------
!        COMMONS YOMRAD15, YOMRDU15

!     METHOD.
!     -------
!        SEE DOCUMENTATION

!     EXTERNALS.
!     ----------
!        SUAER, SUAERH, SUAERV, SULW, SUSW, SUCLD, SUOCST, SUSAT

!     REFERENCE.
!     ----------
!        ECMWF Research Department documentation of the IFS

!     AUTHOR.
!     -------
!        96-11: Ph. Dandin. Meteo-France
!        ORIGINAL : 88-12-15 BY JEAN-JACQUES MORCRETTE  *ECMWF*

!     MODIFICATIONS.
!     --------------
!        R. El Khatib 01-02-02 proper initialization of NFRRC moved in SUCFU
!        M.Hamrud      01-Oct-2003 CY28 Cleaning
!        F. Bouyssel 27-09-04 initialisation of NSW
!        A. Alias    29-09-05 Sulfate aerosols (Hu Rong Ming)
!     ------------------------------------------------------------------

USE PARKIND1  ,ONLY : JPIM     ,JPRB
USE YOMHOOK   ,ONLY : LHOOK,   DR_HOOK

USE YOMCT0   , ONLY : NPRINTLEV
USE YOMDIM   , ONLY : NDLON    ,NSMAX   ,NGPBLKS  ,NFLEVG   ,NPROMA
USE YOMDYN   , ONLY : TSTEP
! Ce qui concerne NULNAM commente par MPL le 15.04.09
!USE YOMLUN   , ONLY : NULNAM
USE YOMCST   , ONLY : RDAY     ,RG       ,RCPD
USE YOMPHY   , ONLY : LRAYFM15
USE YOEPHY   , ONLY : LEPHYS   ,LERADI
USE YOMRAD15 , ONLY : NAER15   ,NFLUX15  ,NMODE15  ,NRAD15   ,&
 & NRADFR15 ,NRADPFR15,NRADPLA15,NRINT15  ,NRADNFR15,&
 & NRADSFR15,NOVLP15  ,NRPROMA15,NRADF2C15,NRADC2F15,&
 & LERAD6H15,LERADHS15,LRADAER15,LNEWAER15
USE YOERAD   , ONLY : NAER     ,NTSW
!USE YOERAD   , ONLY : NAER     ,NSW      ,NTSW
! NSW mis dans .def MPL 20140211
USE YOMRDU15 , ONLY : NUAER15  ,NTRAER15 ,RCDAY15  ,R10E15   ,&
 & REELOG15 ,REPSC15  ,REPSCO15 ,REPSCQ15 ,REPSCT15 ,&
 & REPSCW15 ,DIFF15  
USE YOMAERD15, ONLY : CVDAES15 ,CVDAEL15 ,CVDAEU15 ,CVDAED15 ,&
 & CVDAEF15 ,&
 & RCAEOPS15,RCAEOPL15,RCAEOPU15,RCAEOPD15,RCTRBGA15,&
 & RCAEOPF15,&
 & RCVOBGA15,RCSTBGA15,RCTRPT15 ,RCAEADM15,RCAEROS15,&
 & RCAEADK15  
USE YOMPRAD  , ONLY : LODBGRADI,LODBGRADL
USE YOMRADF  , ONLY : EMTD     ,EMTU      ,TRSW    ,RMOON

IMPLICIT NONE

include "clesphys.h"

INTEGER(KIND=JPIM),INTENT(IN)    :: KLEV
INTEGER(KIND=JPIM),INTENT(IN)    :: KULOUT 
REAL(KIND=JPRB)   ,INTENT(IN)    :: PETAH(KLEV+1)
LOGICAL :: LLP

#include "namrad15.h"
!      ----------------------------------------------------------------

LOGICAL :: LLMESS

INTEGER(KIND=JPIM) :: IRADFR, IST1HR, IST6HR


REAL(KIND=JPRB) :: ZSTPHR, ZTSTEP
REAL(KIND=JPRB) :: ZHOOK_HANDLE

#include "posnam.intfb.h"
#include "suaer15.intfb.h"
#include "suaerv15.intfb.h"
#include "suecradi15.intfb.h"
#include "suecradl.intfb.h"
#include "sulw15.intfb.h"
#include "surdi15.intfb.h"
#include "susat.intfb.h"
#include "susw15.intfb.h"

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

!*         1.       SET DEFAULT VALUES.
!                   -------------------

!*         1.1      PRESET INDICES IN *YOMRAD15*
!                   --------------------------

IF (LHOOK) CALL DR_HOOK('SUECRAD15',0,ZHOOK_HANDLE)
LLMESS=.FALSE.
LERAD6H15=.TRUE.
LERADHS15=.TRUE.
LRADAER15=.TRUE.
LNEWAER15=.FALSE.
NAER15   =1
NAER=0
NFLUX15  =6
NMODE15  =0
NRAD15   =1
NRADFR15 =-3
NRADPFR15=36
NRADPLA15=15
NRINT15  =4
NRADF2C15=1
NRADC2F15=1
NUAER15  = 24
NTRAER15 = 15
NOVLP15 = 1
NSW=2
NTSW=2
IF(NSMAX >= 106) THEN
  NRPROMA15 = 80
ELSEIF(NSMAX == 63) THEN
  NRPROMA15=48
ELSE
  NRPROMA15=20
ENDIF

!*         1.3      SET SECURITY PARAMETERS
!                   -----------------------

REPSC15  = 1.E-12_JPRB
REPSCO15 = 1.E-12_JPRB
REPSCQ15 = 1.E-12_JPRB
REPSCT15 = 1.E-12_JPRB
REPSCW15 = 1.E-12_JPRB
REELOG15 = 1.E-12_JPRB

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

!*         2.       READ VALUES OF RADIATION CONFIGURATION
!                   --------------------------------------

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

!     INITIALISE DATA STRUCTURES REQUIRED FOR RADIATION INTERPOLATION

LODBGRADI=.FALSE.
CALL SUECRADI15

IF( LLMESS )THEN

!     INITIALISE DATA STRUCTURES REQUIRED FOR RADIATION COURSE GRID
!     LOAD BALANCING

  LODBGRADL=.FALSE.
! CALL SUECRADL   ! MPL 1.12.08
  CALL ABOR1('JUSTE APRES CALL SUECRADL COMMENTE')
ENDIF

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

!*       3.    INITIALIZE RADIATION COEFFICIENTS.
!              ----------------------------------

RCDAY15  = RDAY * RG / RCPD
DIFF15   = 1.66_JPRB
R10E15   = 0.4342945_JPRB

CALL SURDI15

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

!*       4.    INITIALIZE RADIATION ABSORPTION COEFFICIENTS
!              --------------------------------------------

CALL SULW15
CALL SUSW15

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

!*       5.    INITIALIZE AEROSOL OPTICAL PARAMETERS AND DISTRIBUTION
!              ------------------------------------------------------

!     INITIALIZATION DONE IN BLOCK DATA SUAERHBD

!- optical properties
CALL SUAER15

!     CALL SUAERH

CALL SUAERV15 ( KLEV  , PETAH,&
 & CVDAES15, CVDAEL15, CVDAEU15, CVDAED15, CVDAEF15,&
 & RCTRBGA15, RCVOBGA15, RCSTBGA15, RCAEOPS15, RCAEOPL15, RCAEOPU15,&
 & RCAEOPF15,&
 & RCAEOPD15, RCTRPT15 , RCAEADK15, RCAEADM15, RCAEROS15        )  

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

!*       6.    INITIALIZE SATELLITE GEOMETRICAL/RADIOMETRIC PARAMETERS
!              -------------------------------------------------------

IF (LEPHYS) THEN
  IF (NMODE15 > 1) THEN
    CALL SUSAT
  ENDIF
ENDIF

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

!*       7.    INITIALIZE CLIMATOLOGICAL OZONE DISTRIBUTION
!              --------------------------------------------
!                  (not done here!!!  called from APLPAR as it depends
!                     on model pressure levels!)

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

!*       8.    SET UP MODEL CONFIGURATION FOR TIME-SPACE INTERPOLATION
!              -------------------------------------------------------

ZTSTEP=MAX(TSTEP,1.0_JPRB)
ZSTPHR=3600._JPRB/ZTSTEP
IRADFR=NRADFR15
IF(NRADFR15 < 0) THEN
  NRADFR15=-NRADFR15*ZSTPHR+0.5_JPRB
ENDIF
NRADPFR15=NRADPFR15*NRADFR15
IF (MOD(NRADPLA15,2) == 0.AND. NRADPLA15 /= 0) THEN
  NRADPLA15=NRADPLA15+1
ENDIF

IST1HR=ZSTPHR+0.05_JPRB
IST6HR=6._JPRB*ZSTPHR+0.05_JPRB
IF (MOD(3600._JPRB,ZTSTEP) > 0.1_JPRB) THEN
  IST1HR=IST1HR+1
  DO WHILE (MOD(IST6HR,IST1HR) /= 0)
    IST1HR=IST1HR+1
  ENDDO
ENDIF
NRADSFR15=IST1HR
NRADNFR15=NRADFR15

IF(LRAYFM15) THEN
  NRPROMA15=NDLON+6+(1-MOD(NDLON,2))
ENDIF



!*       9.    ALLOCATE WORK ARRAYS
!               --------------------

LLP = NPRINTLEV >= 1

  ALLOCATE(EMTD(NPROMA,NFLEVG+1,NGPBLKS))
  IF(LLP)WRITE(UNIT=KULOUT,FMT=9) 'EMTD     ',SIZE(EMTD     ),SHAPE(EMTD     )
  ALLOCATE(TRSW(NPROMA,NFLEVG+1,NGPBLKS))
  IF(LLP)WRITE(UNIT=KULOUT,FMT=9) 'TRSW     ',SIZE(TRSW     ),SHAPE(TRSW     )
  ALLOCATE(EMTU(NPROMA,NFLEVG+1,NGPBLKS))
  IF(LLP)WRITE(UNIT=KULOUT,FMT=9) 'EMTU     ',SIZE(EMTU     ),SHAPE(EMTU     )
  ALLOCATE(RMOON(NPROMA,NGPBLKS))
  IF(LLP)WRITE(UNIT=KULOUT,FMT=9) 'RMOON    ',SIZE(RMOON    ),SHAPE(RMOON    )

9 FORMAT(1X,'ARRAY ',A10,' ALLOCATED ',8I8)

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

!*       9.    PRINT FINAL VALUES.
!              -------------------




IF (LLP) THEN
  WRITE(UNIT=KULOUT,FMT='('' COMMON YOMRAD15 '')')
  WRITE(UNIT=KULOUT,FMT='('' LERADI  = '',L5 &
   & ,'' LERAD6H15 = '',L5)')&
   & LERADI,LERAD6H15  
  WRITE(UNIT=KULOUT,FMT='('' NRADFR15  = '',I2 &
   & ,'' NRADPFR15 = '',I3 &
   & ,'' NRADPLA15 = '',I2 &
   & ,'' NRINT15   = '',I1 &
   & ,'' NRPROMA15 = '',I5 &
   & ,'' NRADF2C15 = '',I1 &
   & ,'' NRADC2F15 = '',I1 &
   & )')&
   & NRADFR15,NRADPFR15,NRADPLA15,NRINT15,&
   & NRPROMA15,NRADF2C15,NRADC2F15  
   
  WRITE(UNIT=KULOUT,FMT='('' LERADHS15= '',L5,'' LRADAER15= '',L5 &
   & ,'' LNEWAER15= '',L5 &
   & ,'' NMODE15 = '',I1 &
   & ,'' NAER15  = '',I1 &
   & ,'' NFLUX15 = '',I2 &
   & ,'' NRAD15  = '',I2 &
   & )')&
   & LERADHS15,LRADAER15,LNEWAER15,NMODE15,NAER15,NFLUX15,NRAD15  
  WRITE(KULOUT,FMT='('' WARNING! CLOUD OVERLAP ASSUMPTION IS''&
   & ,'' NOVLP15   = '',I2 &
   & )')&
   & NOVLP15  

  WRITE(UNIT=KULOUT,FMT='('' MODULE YOERAD '')')
  WRITE(UNIT=KULOUT,FMT='('' NSW = '',I2, '' NTSW = '',I2)') NSW,NTSW
ENDIF

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

IF (LHOOK) CALL DR_HOOK('SUECRAD15',1,ZHOOK_HANDLE)
END SUBROUTINE SUECRAD15