suaerv15.F90

Go to the documentation of this file.

!OPTIONS XOPT(NOEVAL)
SUBROUTINE suaerv15 ( KLEV  , PETAH,&
 & pvdaes,pvdael,pvdaeu,pvdaed,pvdaef,&
 & 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.
!                 -------------

IF (lhook) CALL dr_hook('SUAERV15',0,zhook_handle)
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)
pvdaes(1)=0._jprb
pvdael(1)=0._jprb
pvdaeu(1)=0._jprb
pvdaed(1)=0._jprb
pvdaef(1)=0._jprb
IF(petah(1) /= 0.0_jprb) THEN
  pvdaes(1)=petah(1)**zhss
  pvdael(1)=petah(1)**zhsl
  pvdaeu(1)=petah(1)**zhsu
  pvdaed(1)=petah(1)**zhsd
  pvdaef(1)=petah(1)**zhsf
ENDIF
DO jk=2,klev+1
  pvdaes(jk)=petah(jk)**zhss
  pvdael(jk)=petah(jk)**zhsl
  pvdaeu(jk)=petah(jk)**zhsu
  pvdaed(jk)=petah(jk)**zhsd
  pvdaef(jk)=petah(jk)**zhsf
ENDDO
ptrbga=0.03_jprb/(101325._jprb-19330._jprb)
pvobga=0.007_jprb/19330._jprb
pstbga=0.045_jprb/19330._jprb
paeops=0.05_jprb
paeopl=0.2_jprb
paeopu=0.1_jprb
paeopd=1.9_jprb
paeopf=0.22_jprb
ptrpt=30._jprb
paeadk(1)=+.3876e-03_jprb
paeadk(2)=+.6693e-02_jprb
paeadk(3)=+.8563e-03_jprb
paeadm=2.6e-10_jprb

paeros=0.1462e-16_jprb

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

IF (lhook) CALL dr_hook('SUAERV15',1,zhook_handle)
END SUBROUTINE suaerv15