!OPTIONS XOPT(NOEVAL)

 & PTRBGA,PVOBGA,PSTBGA,PAEOPS,PAEOPL,PAEOPU,&

 & PAEOPF,&

 & PAEOPD,PTRPT,PAEADK,PAEADM,  PAEROS      )

!**** *SUAERV15* - PARAM. FOR THE VERTICAL DISTRIBUTIONS OF AEROSOLS.

!****              FROZEN VERSION (CYCLE 15) OF SUAERV

!     PURPOSE.

!     --------

!          THIS ROUTINE COMPUTES THE VALUES *PVDAEN* (*N=*S,*L,*U OR *D

!     FOR SEA,LAND,URBAN OR DESERT) OF A SURFACE-NORMALISED VERTICAL

!     DISTRIBUTION OF AEROSOLS' OPTICAL DEPHTS FROM THE ARGUMENT *PETAH*

!     (VERTICAL COORDINATE) AT *KLEVP1* LEVELS. IT ALSO SETS VALUES FOR

!     NON-GEOGRAPHICALLY WEIGHTED TOTAL OPTICAL DEPTHS (AT 0.55 E-06

!     WAVE-LENGTH) *PAEOPN* FOR THE SAME FOUR TYPES AND SIMILEAR OPTICAL

!     DEPHTS DIVIDED BY PRESSURE FOR BACKGROUND WELL-MIXED AEROSOLS

!     OF THREE TYPES *PMNBGA* (*MN*=*TR*,*VO* OR *ST* FOR TROPOSPHERIC,

!     VOLCANIC (STRATOSPHERIC ASHES) OR STRATOSPHERIC (SULFURIC TYPE)).

!     IT FINALLY SET VALUES FOR THE POWER TO BE APPLIED TO A TEMPERATURE

!     RATIO SMALLER THAN ONE IN ORDER TO OBTAIN AN IDEX ONE IN THE

!     STRATOSPHERE AND ZERO IN THE TROPOSPHERE WITH A RELATIVELY SMOOTH

!     TRANSITION (*PTRPT*), AS WELL AS FOR ADSORPTION COEFFICIENTS FOR

!     WATER TO THE THREE TYPE OF TROPOSPHERIC AEROSOLS (*PAEADK*) WITH

!     A MINIMUM VALUE (IN THE WHOLE ATMOSPHERE) FOR THE SUM OF THE

!     PRODUCTS OF *PAEADK* BY THE OPTICAL DEPTHS DIVIDED BY PRESURE

!     THICKNESS: *PAEADM*.

!**   INTERFACE.

!     ----------

!          *SUAERV15* IS CALLED FROM *SUECRAD15*.

!          THERE ARE SIXTEEN DUMMY ARGUMENTS: *PETAH* IS THE VERTICAL

!     COORDINATE.

!                                             *PVDAEN* (*N=*S,*L,*U OR

!     *D) ARE THE NORMALISED VERTICAL DISTRIBUTIONS.

!                                             *KLEVP1* IS THE NUMBER OF

!     LEVELS.

!                                             *PMNBGA* (*MN*=*TR*,*VO*

!     OR *ST*) ARE THE BACKGROUND OPTICAL DEPTHS DIVIDED BY PRESSURE.

!                                             *PAEOPN* (*N=*S,*L,*U OR

!     *D) ARE THE TOTAL OPTICAL DEPHTS FOR THE VERTICALLY VARYING

!     AEROSOLS.

!                                             *PTRPT* IS THE TEMPERATURE

!     EXPONENT FOR THE STRATOSPHERIC DEFINITION.

!                                             *PAEADK* (1,2,3) AND

!     AND *PAEADM* ARE THE CONSTANTS FOR THE DEFINITION OF THE QUANTITY

!     OF WATER VAPOUR THAT WILL BE ADSORBED TO THE DRY AEROSOLS TO FORM

!     MOIST AEROSOLS.

!     METHOD.

!     -------

!          STRAIGHTFORWARD, EQUIVALENT HEIGTHS ARE GIVEN IN METERS (8434

!     FOR THE ATMOSPHERE) AND TROPOSPHERIC AND STRATOSPHERIC PRESSURE

!     BOUNDARY VALUES ARE SET AT 101325 AND 19330 *PASCAL.

!     EXTERNALS.

!     ----------

!          NONE.

!     REFERENCE.

!     ----------

!          NONE.

!     AUTHOR

!     ------

!        96-11: Ph. Dandin. Meteo-France

!     J.-J. MORCRETTE  E.C.M.W.F.     91/07/14   ADAPTATION TO I.F.S.

!     J.-F. GELEYN     E.C.M.W.F.     04/11/82.

!     MODIFICATIONS

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

!        M.Hamrud      01-Oct-2003 CY28 Cleaning

!        A.Alias       29-09-2005 Sulfate aerosols (Hu Rong Ming)

!                                 PVDAEF and PAEOPF added

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

USE PARKIND1  ,ONLY : JPIM     ,JPRB

USE YOMHOOK   ,ONLY : LHOOK,   DR_HOOK

IMPLICIT NONE

INTEGER(KIND=JPIM),INTENT(IN)    :: KLEV

REAL(KIND=JPRB)   ,INTENT(IN)    :: PETAH(KLEV+1)

REAL(KIND=JPRB)   ,INTENT(OUT)   :: PVDAES(KLEV+1)

REAL(KIND=JPRB)   ,INTENT(OUT)   :: PVDAEL(KLEV+1)

REAL(KIND=JPRB)   ,INTENT(OUT)   :: PVDAEU(KLEV+1)

REAL(KIND=JPRB)   ,INTENT(OUT)   :: PVDAED(KLEV+1)

REAL(KIND=JPRB)   ,INTENT(OUT)   :: PVDAEF(KLEV+1)

REAL(KIND=JPRB)   ,INTENT(OUT)   :: PTRBGA

REAL(KIND=JPRB)   ,INTENT(OUT)   :: PVOBGA

REAL(KIND=JPRB)   ,INTENT(OUT)   :: PSTBGA

REAL(KIND=JPRB)   ,INTENT(OUT)   :: PAEOPS

REAL(KIND=JPRB)   ,INTENT(OUT)   :: PAEOPL

REAL(KIND=JPRB)   ,INTENT(OUT)   :: PAEOPU

REAL(KIND=JPRB)   ,INTENT(OUT)   :: PAEOPF

REAL(KIND=JPRB)   ,INTENT(OUT)   :: PAEOPD

REAL(KIND=JPRB)   ,INTENT(OUT)   :: PTRPT

REAL(KIND=JPRB)   ,INTENT(OUT)   :: PAEADK(3)

REAL(KIND=JPRB)   ,INTENT(OUT)   :: PAEADM

REAL(KIND=JPRB)   ,INTENT(OUT)   :: PAEROS

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

INTEGER(KIND=JPIM) :: JK

REAL(KIND=JPRB) :: ZHSD, ZHSF, ZHSL, ZHSS, ZHSU

REAL(KIND=JPRB) :: ZHOOK_HANDLE

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

!*         1.     COMPUTATIONS.

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

ZHSS=MAX(1.0_JPRB,8434._JPRB/1000._JPRB)

ZHSL=MAX(1.0_JPRB,8434._JPRB/1000._JPRB)

ZHSU=MAX(1.0_JPRB,8434._JPRB/1000._JPRB)

ZHSD=MAX(1.0_JPRB,8434._JPRB/3000._JPRB)

ZHSF=MAX(1.0_JPRB,8434._JPRB/1000._JPRB)

ENDIF

ENDDO

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