source: LMDZ5/trunk/libf/phymar/suaerv.F90 @ 4038

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

Inclusion de la physique de MAR


Integration of MAR physics

File size: 4.9 KB
Line 
1SUBROUTINE SUAERV &
2 &(  KLEV  , PETAH &
3 &,  PVDAES, PVDAEL, PVDAEU, PVDAED &
4 &,  PTRBGA, PVOBGA, PSTBGA, PAEOPS, PAEOPL, PAEOPU &
5 &,  PAEOPD, PTRPT , PAEADK, PAEADM, PAEROS &
6 &)
7
8
9!**** *SUAERV* - PARAMETERS FOR THE VERTICAL DISTRIBUTIONS OF AEROSOLS.
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!          *SUAERV* IS CALLED FROM *SUECRAD*.
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!     J.-F. GELEYN     E.C.M.W.F.     04/11/82.
74
75!     MODIFICATIONS
76!     -------------
77!     J.-J. MORCRETTE  E.C.M.W.F.     91/07/14   ADAPTATION TO I.F.S.
78!     JJ Morcrette 960903  block data SUAERH moved after end of SUAERV
79
80!     ------------------------------------------------------------------
81
82
83
84#include "tsmbkind.h"
85
86IMPLICIT NONE
87
88
89!     DUMMY INTEGER SCALARS
90INTEGER_M :: KLEV
91
92!     DUMMY REAL SCALARS
93REAL_B :: PAEADM
94REAL_B :: PAEOPD
95REAL_B :: PAEOPL
96REAL_B :: PAEOPS
97REAL_B :: PAEOPU
98REAL_B :: PAEROS
99REAL_B :: PSTBGA
100REAL_B :: PTRBGA
101REAL_B :: PTRPT
102REAL_B :: PVOBGA
103
104
105!     ------------------------------------------------------------------
106REAL_B :: PETAH(KLEV+1)
107
108REAL_B :: PVDAES(KLEV+1),PVDAEL(KLEV+1),PVDAEU(KLEV+1),PVDAED(KLEV+1)
109REAL_B :: PAEADK(3)
110
111!     LOCAL INTEGER SCALARS
112INTEGER_M :: JK
113
114!     LOCAL REAL SCALARS
115REAL_B :: ZHSD, ZHSL, ZHSS, ZHSU
116
117
118!     ------------------------------------------------------------------
119
120
121!*         1.     COMPUTATIONS.
122!                 -------------
123
124
125ZHSS=MAX(_ONE_,8434._JPRB/1000._JPRB)
126ZHSL=MAX(_ONE_,8434._JPRB/1000._JPRB)
127ZHSU=MAX(_ONE_,8434._JPRB/1000._JPRB)
128ZHSD=MAX(_ONE_,8434._JPRB/3000._JPRB)
129PVDAES(1)=0._JPRB
130PVDAEL(1)=0._JPRB
131PVDAEU(1)=0._JPRB
132PVDAED(1)=0._JPRB
133
134! Martin CONTROL
135!PRINT*,'Control SUAERV'
136!PRINT*,'=PVDAES(1)=',PVDAES(1)
137!PRINT*,'PETAH(1)=',PETAH(1)
138!PRINT*,'PETAH=',PETAH
139!PRINT*,'ZHSS=',ZHSS
140!PRINT*,'End Control SUAERV'
141! Martin END CONTROL
142
143IF(PETAH(1) /= _ZERO_) THEN
144  PVDAES(1)=PETAH(1)**ZHSS
145  PVDAEL(1)=PETAH(1)**ZHSL
146  PVDAEU(1)=PETAH(1)**ZHSU
147  PVDAED(1)=PETAH(1)**ZHSD
148ENDIF
149DO JK=2,KLEV+1
150  PVDAES(JK)=PETAH(JK)**ZHSS
151  PVDAEL(JK)=PETAH(JK)**ZHSL
152  PVDAEU(JK)=PETAH(JK)**ZHSU
153  PVDAED(JK)=PETAH(JK)**ZHSD
154ENDDO
155PTRBGA= 0.03_JPRB/(101325._JPRB-19330._JPRB)
156PVOBGA= 0.007_JPRB/19330._JPRB
157PSTBGA= 0.045_JPRB/19330._JPRB
158PAEOPS= 0.05_JPRB
159PAEOPL= 0.2_JPRB
160PAEOPU= 0.1_JPRB
161PAEOPD= 1.9_JPRB
162PTRPT = 30._JPRB
163PAEADK(1)=+.3876E-03_JPRB
164PAEADK(2)=+.6693E-02_JPRB
165PAEADK(3)=+.8563E-03_JPRB
166PAEADM   = 2.6E-10_JPRB
167
168PAEROS   = 0.1462E-16_JPRB
169
170!     ------------------------------------------------------------------
171
172RETURN
173END SUBROUTINE SUAERV
Note: See TracBrowser for help on using the repository browser.