source: LMDZ6/branches/DYNAMICO-conv/libf/phymar/suaersn.F90 @ 5452

Last change on this file since 5452 was 2089, checked in by Laurent Fairhead, 10 years ago

Inclusion de la physique de MAR

Integration of MAR physics
File size: 7.5 KB
Line 
1SUBROUTINE SUAERSN (KTSW,KSW)
2
3!**** *SUAERS*   - INITIALIZE COMMON YOEAER
4
5!     PURPOSE.
6!     --------
7!           INITIALIZE YOEAER, 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 YOEAER
21
22!     METHOD.
23!     -------
24!        SEE DOCUMENTATION
25
26!     EXTERNALS.
27!     ----------
28
29!     REFERENCE.
30!     ----------
31!        ECMWF RESEARCH DEPARTMENT DOCUMENTATION OF THE "IFS MODEL"
32
33!     AUTHOR.
34!     -------
35!        JEAN-JACQUES MORCRETTE *ECMWF*
36
37!     MODIFICATIONS.
38!     --------------
39!        ORIGINAL : 88-02-15
40!        96-01-27  JJ Morcrette  Various spectral resolutions
41!        99-05-25  JJMorcrette   Revised aerosol optical properties
42!        00-10-25  JJMorcrette   6 spectral intervals
43
44!     ------------------------------------------------------------------
45
46#include "tsmbkind.h"
47
48USE YOESW    , ONLY : RTAUA     ,RPIZA    ,RCGA
49
50!      ----------------------------------------------------------------
51
52IMPLICIT NONE
53
54
55!     DUMMY INTEGER SCALARS
56INTEGER_M :: KSW
57INTEGER_M :: KTSW
58
59REAL_B :: ZTAUA2(2,6)  ,ZPIZA2(2,6)  ,ZCGA2(2,6)
60REAL_B :: ZTAUA4(4,6)  ,ZPIZA4(4,6)  ,ZCGA4(4,6)
61REAL_B :: ZTAUA6(6,6)  ,ZPIZA6(6,6)  ,ZCGA6(6,6)
62
63!     LOCAL INTEGER SCALARS
64INTEGER_M :: JAER, JNU
65
66!      ----------------------------------------------------------------
67
68!*       1.    SHORTWAVE COEFFICIENTS
69!              ----------------------
70!=======================================================================
71!-- The (old) five aerosol types were respectively:
72
73!  1/ continental average (+desert)       2/ maritime
74!  3/ urban                               4/ volcanic active
75!  5/ stratospheric background
76
77!-- old values were not spectrally defined:
78! ZTAU2  = .730719, .912819, .725059, .745405, .682188
79! ZPIZA2 = .872212, .982545, .623143, .944887, .997975
80! ZCGA2  = .647596, .739002, .580845, .662657, .624246
81!=======================================================================
82
83!-- The six aerosol types are respectively:
84
85!  1/ continental average                 2/ maritime
86!  3/ desert                              4/ urban
87!  5/ volcanic active                     6/ stratospheric background
88
89! The quantities given are:
90! TAU : ratio of average optical thickness in interval to that at 0.55
91!       micron
92! PIZA: average single scattering albedo
93! CGA : average asymmetry factor
94
95! computed from Hess and Koepke (con, mar, des, urb)
96!          from Bonnel et al.   (vol, str)
97
98
99!        1.1   TWO SPECTRAL INTERVALS (0.25-0.69-4.00microns)
100
101ZTAUA2(1, :)= (/&
102 &1.69446_JPRB , 1.11855_JPRB , 1.09212_JPRB , 1.72145_JPRB , 1.03858_JPRB , 1.12044_JPRB /)
103ZTAUA2(2, :)= (/&
104 &0.40174_JPRB , 0.89383_JPRB , 0.89546_JPRB , 0.40741_JPRB , 0.51143_JPRB , 0.32646_JPRB /)
105 
106ZPIZA2(1, :)= (/&
107 &.9148907_JPRB, .9956173_JPRB, .7504584_JPRB, .8131335_JPRB, .9401905_JPRB, .9999999_JPRB/)
108ZPIZA2(2, :)= (/&
109 &.8814597_JPRB, .9920407_JPRB, .9239428_JPRB, .7546879_JPRB, .9515548_JPRB, .9938563_JPRB/)
110 
111ZCGA2(1, :)= (/&
112 &0.729019_JPRB, 0.803129_JPRB, 0.784592_JPRB, 0.712208_JPRB, .7008249_JPRB, .7270548_JPRB/)
113ZCGA2(2, :)= (/&
114 &0.663224_JPRB, 0.793746_JPRB, 0.696315_JPRB, 0.652612_JPRB, .6608509_JPRB, .6318786_JPRB/)
115
116
117!        1.2   FOUR SPECTRAL INTERVALS (0.25-0.69-1.19-2.38-4.00microns)
118
119ZTAUA4(1, :)= (/&
120 &1.69446_JPRB , 1.11855_JPRB , 1.09212_JPRB , 1.72145_JPRB , 1.03858_JPRB , 1.12044_JPRB /)
121ZTAUA4(2, :)= (/&
122 &0.52838_JPRB , 0.93285_JPRB , 0.93449_JPRB , 0.53078_JPRB , 0.67148_JPRB , 0.46608_JPRB /)
123ZTAUA4(3, :)= (/&
124 &0.20543_JPRB , 0.84642_JPRB , 0.84958_JPRB , 0.21673_JPRB , 0.28270_JPRB , 0.10915_JPRB /)
125ZTAUA4(4, :)= (/&
126 &0.10849_JPRB , 0.66699_JPRB , 0.65255_JPRB , 0.11600_JPRB , 0.06529_JPRB , 0.04468_JPRB /)
127 
128ZPIZA4(1, :)= (/&
129 &.9148907_JPRB, .9956173_JPRB, .7504584_JPRB, .8131335_JPRB, .9401905_JPRB, .9999999_JPRB/)
130ZPIZA4(2, :)= (/&
131 &.8970131_JPRB, .9984940_JPRB, .9245594_JPRB, .7768385_JPRB, .9532763_JPRB, .9999999_JPRB/)
132ZPIZA4(3, :)= (/&
133 &.8287144_JPRB, .9949396_JPRB, .9279543_JPRB, .6765051_JPRB, .9467578_JPRB, .9955938_JPRB/)
134ZPIZA4(4, :)= (/&
135 &.5230504_JPRB, .7868518_JPRB, .8531531_JPRB, .4048149_JPRB, .8748231_JPRB, .2355667_JPRB/)
136 
137ZCGA4(1, :)= (/&
138 &0.729019_JPRB, 0.803129_JPRB, 0.784592_JPRB, 0.712208_JPRB, .7008249_JPRB, .7270548_JPRB/)
139ZCGA4(2, :)= (/&
140 &0.668431_JPRB, 0.788530_JPRB, 0.698682_JPRB, 0.657422_JPRB, .6735182_JPRB, .6519706_JPRB/)
141ZCGA4(3, :)= (/&
142 &0.636342_JPRB, 0.802467_JPRB, 0.691305_JPRB, 0.627497_JPRB, .6105750_JPRB, .4760794_JPRB/)
143ZCGA4(4, :)= (/&
144 &0.700610_JPRB, 0.818871_JPRB, 0.702399_JPRB, 0.689886_JPRB, .4629866_JPRB, .1907639_JPRB/)
145
146
147!        1.3   SIX SPECTRAL INTERVALS (0.185-0.25-0.44-0.69-1.19-2.38-4.00microns)
148
149ZTAUA6(1, :)= (/&
150 &1.69446_JPRB , 1.11855_JPRB , 1.09212_JPRB , 1.72145_JPRB , 1.03858_JPRB , 1.12044_JPRB /)
151ZTAUA6(2, :)= (/&
152 &1.69446_JPRB , 1.11855_JPRB , 1.09212_JPRB , 1.72145_JPRB , 1.03858_JPRB , 1.12044_JPRB /)
153ZTAUA6(3, :)= (/&
154 &1.69446_JPRB , 1.11855_JPRB , 1.09212_JPRB , 1.72145_JPRB , 1.03858_JPRB , 1.12044_JPRB /)
155ZTAUA6(4, :)= (/&
156 &0.52838_JPRB , 0.93285_JPRB , 0.93449_JPRB , 0.53078_JPRB , 0.67148_JPRB , 0.46608_JPRB /)
157ZTAUA6(5, :)= (/&
158 &0.20543_JPRB , 0.84642_JPRB , 0.84958_JPRB , 0.21673_JPRB , 0.28270_JPRB , 0.10915_JPRB /)
159ZTAUA6(6, :)= (/&
160 &0.10849_JPRB , 0.66699_JPRB , 0.65255_JPRB , 0.11600_JPRB , 0.06529_JPRB , 0.04468_JPRB /)
161 
162ZPIZA6(1, :)= (/&
163 &.9148907_JPRB, .9956173_JPRB, .7504584_JPRB, .8131335_JPRB, .9401905_JPRB, .9999999_JPRB/)
164ZPIZA6(2, :)= (/&
165 &.9148907_JPRB, .9956173_JPRB, .7504584_JPRB, .8131335_JPRB, .9401905_JPRB, .9999999_JPRB/)
166ZPIZA6(3, :)= (/&
167 &.9148907_JPRB, .9956173_JPRB, .7504584_JPRB, .8131335_JPRB, .9401905_JPRB, .9999999_JPRB/)
168ZPIZA6(4, :)= (/&
169 &.8970131_JPRB, .9984940_JPRB, .9245594_JPRB, .7768385_JPRB, .9532763_JPRB, .9999999_JPRB/)
170ZPIZA6(5, :)= (/&
171 &.8287144_JPRB, .9949396_JPRB, .9279543_JPRB, .6765051_JPRB, .9467578_JPRB, .9955938_JPRB/)
172ZPIZA6(6, :)= (/&
173 &.5230504_JPRB, .7868518_JPRB, .8531531_JPRB, .4048149_JPRB, .8748231_JPRB, .2355667_JPRB/)
174 
175ZCGA6(1, :)= (/&
176 &0.729019_JPRB, 0.803129_JPRB, 0.784592_JPRB, 0.712208_JPRB, .7008249_JPRB, .7270548_JPRB/)
177ZCGA6(2, :)= (/&
178 &0.729019_JPRB, 0.803129_JPRB, 0.784592_JPRB, 0.712208_JPRB, .7008249_JPRB, .7270548_JPRB/)
179ZCGA6(3, :)= (/&
180 &0.729019_JPRB, 0.803129_JPRB, 0.784592_JPRB, 0.712208_JPRB, .7008249_JPRB, .7270548_JPRB/)
181ZCGA6(4, :)= (/&
182 &0.668431_JPRB, 0.788530_JPRB, 0.698682_JPRB, 0.657422_JPRB, .6735182_JPRB, .6519706_JPRB/)
183ZCGA6(5, :)= (/&
184 &0.636342_JPRB, 0.802467_JPRB, 0.691305_JPRB, 0.627497_JPRB, .6105750_JPRB, .4760794_JPRB/)
185ZCGA6(6, :)= (/&
186 &0.700610_JPRB, 0.818871_JPRB, 0.702399_JPRB, 0.689886_JPRB, .4629866_JPRB, .1907639_JPRB/)
187
188
189!      ----------------------------------------------------------------
190
191IF (KSW == 2) THEN
192  DO JNU=1,KSW
193    DO JAER=1,6
194      RTAUA(JNU,JAER)=ZTAUA2(JNU,JAER)
195      RPIZA(JNU,JAER)=ZPIZA2(JNU,JAER)
196      RCGA(JNU,JAER) =ZCGA2 (JNU,JAER)
197    ENDDO
198  ENDDO
199ELSEIF (KSW == 4) THEN
200  DO JNU=1,KSW
201    DO JAER=1,6
202      RTAUA(JNU,JAER)=ZTAUA4(JNU,JAER)
203      RPIZA(JNU,JAER)=ZPIZA4(JNU,JAER)
204      RCGA(JNU,JAER) =ZCGA4 (JNU,JAER)
205    ENDDO
206  ENDDO
207ELSEIF (KSW == 6) THEN
208  DO JNU=1,KSW
209    DO JAER=1,6
210      RTAUA(JNU,JAER)=ZTAUA6(JNU,JAER)
211      RPIZA(JNU,JAER)=ZPIZA6(JNU,JAER)
212      RCGA(JNU,JAER) =ZCGA6 (JNU,JAER)
213    ENDDO
214  ENDDO
215ELSE
216  STOP 'SUAERSN: WRONG NUMBER OF SPECTRAL INTERVALS'
217ENDIF
218
219!      ----------------------------------------------------------------
220
221RETURN
222END SUBROUTINE SUAERSN
Note: See TracBrowser for help on using the repository browser.