source: LMDZ5/branches/LMDZ5_SPLA/libf/phylmd/rrtm/suaerv15.F90

Last change on this file was 1990, checked in by Laurent Fairhead, 11 years ago

Corrections à la version r1989 pour permettre la compilation avec RRTM
Inclusion de la licence CeCILL_V2 pour RRTM


Changes to revision r1989 to enable RRTM code compilation
RRTM part put under CeCILL_V2 licence

  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
File size: 5.7 KB
Line 
1!OPTIONS XOPT(NOEVAL)
2SUBROUTINE SUAERV15 ( KLEV  , PETAH,&
3 & PVDAES,PVDAEL,PVDAEU,PVDAED,PVDAEF,&
4 & PTRBGA,PVOBGA,PSTBGA,PAEOPS,PAEOPL,PAEOPU,&
5 & PAEOPF,&
6 & PAEOPD,PTRPT,PAEADK,PAEADM,  PAEROS      ) 
7
8!**** *SUAERV15* - PARAM. FOR THE VERTICAL DISTRIBUTIONS OF AEROSOLS.
9!****              FROZEN VERSION (CYCLE 15) OF SUAERV
10
11!     PURPOSE.
12!     --------
13
14!          THIS ROUTINE COMPUTES THE VALUES *PVDAEN* (*N=*S,*L,*U OR *D
15!     FOR SEA,LAND,URBAN OR DESERT) OF A SURFACE-NORMALISED VERTICAL
16!     DISTRIBUTION OF AEROSOLS' OPTICAL DEPHTS FROM THE ARGUMENT *PETAH*
17!     (VERTICAL COORDINATE) AT *KLEVP1* LEVELS. IT ALSO SETS VALUES FOR
18!     NON-GEOGRAPHICALLY WEIGHTED TOTAL OPTICAL DEPTHS (AT 0.55 E-06
19!     WAVE-LENGTH) *PAEOPN* FOR THE SAME FOUR TYPES AND SIMILEAR OPTICAL
20!     DEPHTS DIVIDED BY PRESSURE FOR BACKGROUND WELL-MIXED AEROSOLS
21!     OF THREE TYPES *PMNBGA* (*MN*=*TR*,*VO* OR *ST* FOR TROPOSPHERIC,
22!     VOLCANIC (STRATOSPHERIC ASHES) OR STRATOSPHERIC (SULFURIC TYPE)).
23!     IT FINALLY SET VALUES FOR THE POWER TO BE APPLIED TO A TEMPERATURE
24!     RATIO SMALLER THAN ONE IN ORDER TO OBTAIN AN IDEX ONE IN THE
25!     STRATOSPHERE AND ZERO IN THE TROPOSPHERE WITH A RELATIVELY SMOOTH
26!     TRANSITION (*PTRPT*), AS WELL AS FOR ADSORPTION COEFFICIENTS FOR
27!     WATER TO THE THREE TYPE OF TROPOSPHERIC AEROSOLS (*PAEADK*) WITH
28!     A MINIMUM VALUE (IN THE WHOLE ATMOSPHERE) FOR THE SUM OF THE
29!     PRODUCTS OF *PAEADK* BY THE OPTICAL DEPTHS DIVIDED BY PRESURE
30!     THICKNESS: *PAEADM*.
31
32!**   INTERFACE.
33!     ----------
34
35!          *SUAERV15* IS CALLED FROM *SUECRAD15*.
36!          THERE ARE SIXTEEN DUMMY ARGUMENTS: *PETAH* IS THE VERTICAL
37!     COORDINATE.
38!                                             *PVDAEN* (*N=*S,*L,*U OR
39!     *D) ARE THE NORMALISED VERTICAL DISTRIBUTIONS.
40!                                             *KLEVP1* IS THE NUMBER OF
41!     LEVELS.
42!                                             *PMNBGA* (*MN*=*TR*,*VO*
43!     OR *ST*) ARE THE BACKGROUND OPTICAL DEPTHS DIVIDED BY PRESSURE.
44!                                             *PAEOPN* (*N=*S,*L,*U OR
45!     *D) ARE THE TOTAL OPTICAL DEPHTS FOR THE VERTICALLY VARYING
46!     AEROSOLS.
47!                                             *PTRPT* IS THE TEMPERATURE
48!     EXPONENT FOR THE STRATOSPHERIC DEFINITION.
49!                                             *PAEADK* (1,2,3) AND
50!     AND *PAEADM* ARE THE CONSTANTS FOR THE DEFINITION OF THE QUANTITY
51!     OF WATER VAPOUR THAT WILL BE ADSORBED TO THE DRY AEROSOLS TO FORM
52!     MOIST AEROSOLS.
53
54!     METHOD.
55!     -------
56
57!          STRAIGHTFORWARD, EQUIVALENT HEIGTHS ARE GIVEN IN METERS (8434
58!     FOR THE ATMOSPHERE) AND TROPOSPHERIC AND STRATOSPHERIC PRESSURE
59!     BOUNDARY VALUES ARE SET AT 101325 AND 19330 *PASCAL.
60
61!     EXTERNALS.
62!     ----------
63
64!          NONE.
65
66!     REFERENCE.
67!     ----------
68
69!          NONE.
70
71!     AUTHOR
72!     ------
73!        96-11: Ph. Dandin. Meteo-France
74!     J.-J. MORCRETTE  E.C.M.W.F.     91/07/14   ADAPTATION TO I.F.S.
75!     J.-F. GELEYN     E.C.M.W.F.     04/11/82.
76
77!     MODIFICATIONS
78!     -------------
79!        M.Hamrud      01-Oct-2003 CY28 Cleaning
80!        A.Alias       29-09-2005 Sulfate aerosols (Hu Rong Ming)
81!                                 PVDAEF and PAEOPF added
82!     ------------------------------------------------------------------
83
84USE PARKIND1  ,ONLY : JPIM     ,JPRB
85USE YOMHOOK   ,ONLY : LHOOK,   DR_HOOK
86
87IMPLICIT NONE
88
89INTEGER(KIND=JPIM),INTENT(IN)    :: KLEV
90REAL(KIND=JPRB)   ,INTENT(IN)    :: PETAH(KLEV+1)
91REAL(KIND=JPRB)   ,INTENT(OUT)   :: PVDAES(KLEV+1)
92REAL(KIND=JPRB)   ,INTENT(OUT)   :: PVDAEL(KLEV+1)
93REAL(KIND=JPRB)   ,INTENT(OUT)   :: PVDAEU(KLEV+1)
94REAL(KIND=JPRB)   ,INTENT(OUT)   :: PVDAED(KLEV+1)
95REAL(KIND=JPRB)   ,INTENT(OUT)   :: PVDAEF(KLEV+1)
96REAL(KIND=JPRB)   ,INTENT(OUT)   :: PTRBGA
97REAL(KIND=JPRB)   ,INTENT(OUT)   :: PVOBGA
98REAL(KIND=JPRB)   ,INTENT(OUT)   :: PSTBGA
99REAL(KIND=JPRB)   ,INTENT(OUT)   :: PAEOPS
100REAL(KIND=JPRB)   ,INTENT(OUT)   :: PAEOPL
101REAL(KIND=JPRB)   ,INTENT(OUT)   :: PAEOPU
102REAL(KIND=JPRB)   ,INTENT(OUT)   :: PAEOPF
103REAL(KIND=JPRB)   ,INTENT(OUT)   :: PAEOPD
104REAL(KIND=JPRB)   ,INTENT(OUT)   :: PTRPT
105REAL(KIND=JPRB)   ,INTENT(OUT)   :: PAEADK(3)
106REAL(KIND=JPRB)   ,INTENT(OUT)   :: PAEADM
107REAL(KIND=JPRB)   ,INTENT(OUT)   :: PAEROS
108!     ------------------------------------------------------------------
109INTEGER(KIND=JPIM) :: JK
110
111REAL(KIND=JPRB) :: ZHSD, ZHSF, ZHSL, ZHSS, ZHSU
112REAL(KIND=JPRB) :: ZHOOK_HANDLE
113
114!     ------------------------------------------------------------------
115
116!*         1.     COMPUTATIONS.
117!                 -------------
118
119IF (LHOOK) CALL DR_HOOK('SUAERV15',0,ZHOOK_HANDLE)
120ZHSS=MAX(1.0_JPRB,8434._JPRB/1000._JPRB)
121ZHSL=MAX(1.0_JPRB,8434._JPRB/1000._JPRB)
122ZHSU=MAX(1.0_JPRB,8434._JPRB/1000._JPRB)
123ZHSD=MAX(1.0_JPRB,8434._JPRB/3000._JPRB)
124ZHSF=MAX(1.0_JPRB,8434._JPRB/1000._JPRB)
125PVDAES(1)=0._JPRB
126PVDAEL(1)=0._JPRB
127PVDAEU(1)=0._JPRB
128PVDAED(1)=0._JPRB
129PVDAEF(1)=0._JPRB
130IF(PETAH(1) /= 0.0_JPRB) THEN
131  PVDAES(1)=PETAH(1)**ZHSS
132  PVDAEL(1)=PETAH(1)**ZHSL
133  PVDAEU(1)=PETAH(1)**ZHSU
134  PVDAED(1)=PETAH(1)**ZHSD
135  PVDAEF(1)=PETAH(1)**ZHSF
136ENDIF
137DO JK=2,KLEV+1
138  PVDAES(JK)=PETAH(JK)**ZHSS
139  PVDAEL(JK)=PETAH(JK)**ZHSL
140  PVDAEU(JK)=PETAH(JK)**ZHSU
141  PVDAED(JK)=PETAH(JK)**ZHSD
142  PVDAEF(JK)=PETAH(JK)**ZHSF
143ENDDO
144PTRBGA=0.03_JPRB/(101325._JPRB-19330._JPRB)
145PVOBGA=0.007_JPRB/19330._JPRB
146PSTBGA=0.045_JPRB/19330._JPRB
147PAEOPS=0.05_JPRB
148PAEOPL=0.2_JPRB
149PAEOPU=0.1_JPRB
150PAEOPD=1.9_JPRB
151PAEOPF=0.22_JPRB
152PTRPT=30._JPRB
153PAEADK(1)=+.3876E-03_JPRB
154PAEADK(2)=+.6693E-02_JPRB
155PAEADK(3)=+.8563E-03_JPRB
156PAEADM=2.6E-10_JPRB
157
158PAEROS=0.1462E-16_JPRB
159
160!     ------------------------------------------------------------------
161
162IF (LHOOK) CALL DR_HOOK('SUAERV15',1,ZHOOK_HANDLE)
163END SUBROUTINE SUAERV15
Note: See TracBrowser for help on using the repository browser.