source: LMDZ5/branches/testing/libf/phylmd/rrtm/suaer15.F90 @ 5431

Last change on this file since 5431 was 1999, checked in by Laurent Fairhead, 11 years ago

Merged trunk changes r1920:1997 into testing branch

  • 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: 4.4 KB
Line 
1SUBROUTINE SUAER15
2
3!**** *SUAER15*   - INITIALIZE COMMON YOMAER15
4
5!     PURPOSE.
6!     --------
7!           INITIALIZE YOMAER15, THE COMMON THAT CONTAINS THE
8!           RADIATIVE CHARACTERISTICS OF THE AEROSOLS
9
10!**   INTERFACE.
11!     ----------
12!              -----        -----
13
14!        EXPLICIT ARGUMENTS :
15!        --------------------
16!        NONE
17
18!        IMPLICIT ARGUMENTS :
19!        --------------------
20!        COMMON YOMAER15
21
22!     METHOD.
23!     -------
24!        SEE DOCUMENTATION
25
26!     EXTERNALS.
27!     ----------
28
29!     REFERENCE.
30!     ----------
31
32!     04-06: F. Bouyssel. Meteo-France
33!     05-09: A. Alias - PAESOD (black_carbon) is added to PAELAN (oraganic)
34!                     - the sulfate is put in a separate type
35!                     (see Hu Ron Ming work) - P.Marquet
36
37!=======================================================================
38!-- The (old) five aerosol types were respectively:
39
40!  1/ continental average (+desert)       2/ maritime
41!  3/ urban                               4/ volcanic active
42!  5/ stratospheric background
43
44!=======================================================================
45
46!-- The (new) six aerosol types are respectively:
47
48!  1/ continental average                 2/ maritime
49!  3/ desert                              4/ volcanic active
50!  5/ stratospheric background            6/ sulfate
51
52!     ------------------------------------------------------------------
53
54USE PARKIND1  ,ONLY : JPIM     ,JPRB
55USE YOMHOOK   ,ONLY : LHOOK,   DR_HOOK
56
57USE YOMRAD15 , ONLY : LNEWAER15
58USE YOMAER15 , ONLY : RAER15    ,TAUA15   ,RPIZA15  ,RCGA15
59
60IMPLICIT NONE
61
62!     ------------------------------------------------------------------
63
64REAL(KIND=JPRB) :: ZHOOK_HANDLE
65IF (LHOOK) CALL DR_HOOK('SUAER15',0,ZHOOK_HANDLE)
66IF (.NOT.LNEWAER15) THEN
67 
68  RAER15=RESHAPE((/&
69   & .038520_JPRB, .037196_JPRB, .040532_JPRB, .054934_JPRB, .038520_JPRB ,&
70   & .12613_JPRB , .18313_JPRB , .10357_JPRB , .064106_JPRB, .126130_JPRB ,&
71   & .012579_JPRB, .013649_JPRB, .018652_JPRB, .025181_JPRB, .012579_JPRB ,&
72   & .011890_JPRB, .016142_JPRB, .021105_JPRB, .028908_JPRB, .011890_JPRB ,&
73   & .013792_JPRB, .026810_JPRB, .052203_JPRB, .066338_JPRB, .013792_JPRB ,&
74   & .012579_JPRB, .013649_JPRB, .018652_JPRB, .025181_JPRB, .012579_JPRB /)&
75   & ,SHAPE=(/5,6/)) 
76
77  TAUA15=RESHAPE((/&
78   & .730719_JPRB,.730719_JPRB,&
79   & .912819_JPRB,.912819_JPRB,&
80   & .725059_JPRB,.725059_JPRB,&
81   & .745405_JPRB,.745405_JPRB,&
82   & .682188_JPRB,.682188_JPRB,&
83   & 1.35059_JPRB,.725059_JPRB /)&
84   & ,SHAPE=(/2,6/)) 
85
86  RPIZA15=RESHAPE((/&
87   & .872212_JPRB,.872212_JPRB,&
88   & .982545_JPRB,.982545_JPRB,&
89   & .623143_JPRB,.623143_JPRB,&
90   & .944887_JPRB,.944887_JPRB,&
91   & .997975_JPRB,.997975_JPRB,&
92   & .999999_JPRB,.999999_JPRB /)&
93   & ,SHAPE=(/2,6/)) 
94
95  RCGA15=RESHAPE((/&
96   & .647596_JPRB,.647596_JPRB,&
97   & .739002_JPRB,.739002_JPRB,&
98   & .580845_JPRB,.580845_JPRB,&
99   & .662657_JPRB,.662657_JPRB,&
100   & .624246_JPRB,.624246_JPRB,&
101   & .680845_JPRB,.680845_JPRB /)&
102   & ,SHAPE=(/2,6/)) 
103
104ELSE
105
106  RAER15( :, 1)= (/&
107   & .036271_JPRB, .030153_JPRB, .017343_JPRB, .015002_JPRB, .008806_JPRB /) 
108  RAER15( :, 2)= (/&
109   & .026561_JPRB, .032657_JPRB, .017977_JPRB, .014210_JPRB, .016775_JPRB /) 
110  RAER15( :, 3)= (/&
111   & .014897_JPRB, .016359_JPRB, .019789_JPRB, .030777_JPRB, .013341_JPRB /) 
112  RAER15( :, 4)= (/&
113   & .011890_JPRB, .016142_JPRB, .021105_JPRB, .028908_JPRB, .011890_JPRB /) 
114  RAER15( :, 5)= (/&
115   & .013792_JPRB, .026810_JPRB, .052203_JPRB, .066338_JPRB, .013792_JPRB /) 
116  RAER15( :, 6)= (/&
117   & .012579_JPRB, .013649_JPRB, .018652_JPRB, .025181_JPRB, .012579_JPRB /)
118
119  TAUA15(1, :)= (/&
120   & 1.69446_JPRB, 1.11855_JPRB, 1.09212_JPRB, 1.03858_JPRB, 1.12044_JPRB, 1.35059_JPRB  /) 
121  TAUA15(2, :)= (/&
122   & 0.40174_JPRB, 0.89383_JPRB, 0.89546_JPRB, 0.51143_JPRB, 0.32646_JPRB, .725059_JPRB /) 
123 
124  RPIZA15(1, :)= (/&
125   & .9148907_JPRB, .9956173_JPRB, .7504584_JPRB, .9401905_JPRB, .9999999_JPRB, .999999_JPRB/) 
126  RPIZA15(2, :)= (/&
127   & .8814597_JPRB, .9920407_JPRB, .9239428_JPRB, .9515548_JPRB, .9938563_JPRB, .999999_JPRB/) 
128 
129  RCGA15(1, :)= (/&
130   & 0.729019_JPRB, 0.803129_JPRB, 0.784592_JPRB, .7008249_JPRB, .7270548_JPRB, .680845_JPRB/) 
131  RCGA15(2, :)= (/&
132   & 0.663224_JPRB, 0.793746_JPRB, 0.696315_JPRB, .6608509_JPRB, .6318786_JPRB, .680845_JPRB/) 
133
134ENDIF
135
136!      ----------------------------------------------------------------
137
138IF (LHOOK) CALL DR_HOOK('SUAER15',1,ZHOOK_HANDLE)
139END SUBROUTINE SUAER15
Note: See TracBrowser for help on using the repository browser.